ext<-范围(b)
lon<-seq(ext@xmin,ext@xmax,res(b)[1])
所以基本上你要创建一个序列向量,从x/y最小到最大,间隔为砖块的分辨率。
仅作说明:
X<-光栅(分辨率=C(40,40))
绘图(X)
col='red',pch='*',cex=5,add=t)
所以上面的方法给出了红色星号表示的纬度。如果您需要蓝色的,可以使用xyfromcellwhich返回光栅单元的坐标。
# create testraster
x <- raster(resolution=c(40,40))
x[]<- 1:ncell(x)
# plot
plot(x)
# add corner coordinates
plot(SpatialPoints(cbind(rep(extent(x)@xmin,10),seq(extent(x)@ymin,extent(x)@ymax,res(x)[2])),proj4string = crs(x)),
col='red',pch='*',cex=5,add=T)
# add cell centers
plot(SpatialPoints(xyFromCell(x,cellFromRowCol(x,1:nrow(x),1)),proj4string = crs(x)),
col='blue',pch='*',cex=5,add=T)
xyFromCell
它返回光栅单元的坐标。
