2016-03-25 7 views
0

私が取り組んでいる課題については、発生点を使って分布モデルを実行しています。表示されるプロットは大きすぎるエクステントを示しており、xlimとylimを変更して、プロットがオカレンスポイント(オレンジ色の領域)のみを表示するようにしたいと考えています。私はそれをいくつかの方法で操作しようとしましたが、私が望む程度には焦点を当てません。私が本当に理解していないのは、xlimとylimが指定しているものだと思います。xlimとylimを操作してマップをプロットする方法

sparrow= gbif("ammodramus", "maritimus*", geo=FALSE) 
sparrow <- subset(sparrow, !is.na(lon) & !is.na(lat)) 
library(maptools) 
data(wrld_simpl) 
plot(wrld_simpl, xlim=c(-80,70), ylim=c(-60,10), axes=TRUE, col="lightyellow") 
points(sparrow$lon, sparrow$lat, col="orange", pch=20, cex=0.75) 
points(acgeo$lon, acgeo$lat, col="red", cex=0.75) 
+0

私は最近maptoolsを使っていませんが、xlimとylimはプロットの緯度と経度を指定する必要があります。私はあなたのコードを実行していないが、私はあなたが南アメリカからアフリカにすべてをマッピングしていると思われる。 – Dave2e

+0

あなたの例では 'acgeo'と書いてあります。おそらく' sparrow'を意味しています。また、 'library(dismo)'で始めると、再現性があります。 –

答えて

0

あなたは

library(sp) 
coordinates(sparrow) = ~lon + lat 
proj4string(sparrow) = proj4string(wrld_simpl) 
plot(sparrow, col = 'red', cex = .75, axes = TRUE) 
plot(wrld_simpl, add = TRUE) 

を行うと、プロットがsparrowからxlimylimを取ります。

+0

これはうまくいきました、ありがとう! – user5899223

関連する問題