不匹配 R 中的栅格和多边形数据答案

【问题标题】：Unmatch raster and polygon data in R不匹配 R 中的栅格和多边形数据
【发布时间】：2018-09-10 01:50:14
【问题描述】：

多边形（从 shapefile 读取）是：

class       : SpatialPolygonsDataFrame 
features    : 653 
extent      : -7115213, 4895580, 1368240, 7805331  (xmin, xmax, ymin, 
ymax)
coord. ref. : +proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 
+x_0=0 +y_0=0 +datum=NAD27 +units=us-ft +no_defs +ellps=clrk66 
+nadgrids=@conus,@alaska,@ntv2_0.gsb,@ntv1_can.dat 
variables   : 6
names       : cat,   NAME,      AREA_MI,      lccount,        lcsum,      
lcmean 
min values  :   1, Alaska, 1.006402e+00, 1.000000e+01, 0.000000e+00,  
0.00000000 
max values  :  99, Alaska, 9.945810e-01, 9.900000e+01, 9.960000e+02, 
12.00000000

栅格（高程）数据是

class       : RasterLayer 
dimensions  : 6800, 9200, 62560000  (nrow, ncol, ncell)
resolution  : 1305.521, 946.6311  (x, y)
extent      : -7115213, 4895580, 1368240, 7805331  (xmin, xmax, ymin, 
ymax)
coord. ref. : +proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 
+x_0=0 +y_0=0 +datum=NAD27 +units=us-ft +no_defs +ellps=clrk66 
+nadgrids=@conus,@alaska,@ntv2_0.gsb,@ntv1_can.dat 
data source : /Users/hong/Documents/GitHub/Alaska/akshd300m.tif 
names       : akshd300m 
values      : 0, 255  (min, max)

有没有人可以解决这个问题？我用 tmap 来绘制这个图

【问题讨论】：

标签： r raster

【解决方案1】：

虽然两个数据集的坐标参考系统都报告为+proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 +x_0=0 +y_0=0 +datum=NAD27 +units=us-ft +no_defs +ellps=clrk66，但显然这不可能。通过将已知多边形转换为 crs，应该很容易找出哪个是正确的。例如

library(raster)
library(rgdal)
usa <- getData('GADM', country='USA', level=1)
ak <- usa[usa$NAME_1 == 'Alaska', ]
ak_asa <- spTransform(ak, '+proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 +x_0=0 +y_0=0 +datum=NAD27 +units=us-ft +no_defs +ellps=clrk66')

并且也在绘制那个。确定哪一个是正确的（希望两者中的一个是正确的），然后您可以尝试追溯另一个错误的地方。根据您提供的信息，我们无法为您解决。

【讨论】：