Fatal编程技术网
  • r/
  • R中光栅对象坐标系的转换

R中光栅对象坐标系的转换

r

R中光栅对象坐标系的转换,r,R,我正在使用Sentinel 2图像数据进行NDVI监测 我有一个区域的村级边界，输入为R中的SpatialPolygonsDataFrame。当我试图裁剪出感兴趣的区域时，到达NDVI光栅层后，它表示范围不重叠 我发现我的区域变量的CRS与我的NDVI图层（utm）不同（longlat） >crs（ndvi） CRS参数： +项目=utm+区域=42+基准=WGS84+单位=m+无定义+ellps=WGS84+牵引力=0,0,0 >crs（区域） CRS参数： +项目=longlat+datum

我正在使用Sentinel 2图像数据进行NDVI监测

我有一个区域的村级边界，输入为R中的
SpatialPolygonsDataFrame
。当我试图裁剪出感兴趣的区域时，到达NDVI光栅层后，它表示范围不重叠

我发现我的区域变量的CRS与我的NDVI
图层（utm）
不同（longlat）

>crs（ndvi）
CRS参数：
+项目=utm+区域=42+基准=WGS84+单位=m+无定义+ellps=WGS84+牵引力=0,0,0
>crs（区域）
CRS参数：
+项目=longlat+datum=WGS84+no_defs+ellps=WGS84+towgs84=0,0,0

>r您可以使用“区域”上的
rgdal:：spTranform
创建与“ndvi”具有相同crs的对象


库（rgdal）
规则

> crs(ndvi)
CRS arguments:
 +proj=utm +zone=42 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0

> crs(region)
CRS arguments:
 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0

> r <- crop(ndvi,extent(region))
Error in .local(x, y, ...) : extents do not overlap`



library(rgdal)
reg <- spTransform(region, crs(ndvi))



r <- crop(ndvi, reg)



m <- mask(r, reg)