Блог

Когда я создаю растр земной поверхности на основе wrld_simpl (или любого другого слоя окружающей среды, исходящего из worldclim), всегда появляются какие-то «невозможные» NAS на суше. Почему это должно произойти? Мне нужна идеальная маска мировой суши, чтобы извлекать записи, которые не упали в океан. Тем не менее, существует много записей на суше, и они все еще считаются NA.

Мой сценарий выглядит следующим образом:

 require(raster)
require(maptools)
data(wrld_simpl)
x=read.csv("https://www.dropbox.com/s/ncvu64r2fxgfd4e/NAlocations.csv?dl=0")
r=raster(ncols=360,nrows=(180))
extent(r)=extent(wrld_simpl)
r=rasterize(wrld_simpl,r,wrld_simpl$AREA)
plot(r)
x=x[-which(is.na(extract(r,x$lon,x$lat))),]# This should eliminate all locations on land.
points(x$lon,x$lat, col="red", cex=.3)            
  

Как это возможно? И было бы ли это способом создать чистый растр для world land?

Прямой read.csv доступ из dropbox у меня не работает. Если я сделаю

 z <- extract(r, x)
# NOT z <- extract(r, x[,1], x[,2]) !!! 

i <- which(is.na(z))
points(x[i,])
  

Я вижу кучу точек в воде у побережья Мозамбика.