这应该有效:
> ll = projectRaster(wrf_temporary,crs="+init=epsg:4326")
> plot(ll[[1]])
并产生这个:
乍一看是正确的,但格林威治子午线(经度=0)并不穿过伦敦。
如果分辨率设置为其他值(4 或 12),则不会出现此问题。当使用分辨率=36 调用时,您会得到一个比其他值更大的栅格,但具有相同的投影字符串。这不可能是正确的。例如,以下绘图上的 (0,0) 坐标应该是地球上的同一点,因为它们声称是相同的投影,但事实并非如此。
我认为分辨率=12是正确的。这是转换为经纬度的栅格,并覆盖了欧盟矢量覆盖范围:
完美排列。
所以我会说它是一个错误 - 或者getRaster
猜测投影并弄错,或者在裁剪后没有应用投影更改,或者它使用的任何服务对投影撒谎。
getRaster
正在猜测。投影信息位于下载的 NetCDF 文件中。以另一种格式。对于分辨率 = 36 文件,投影信息部分如下:
int Lambert_Conformal ;
Lambert_Conformal:standard_parallel = 43., 43. ;
Lambert_Conformal:latitude_of_projection_origin = 24.2280006408691 ;
Lambert_Conformal:false_easting = 2182.62935 ;
Lambert_Conformal:false_northing = -269.65597 ;
Lambert_Conformal:grid_mapping_name = "lambert_conformal_conic" ;
Lambert_Conformal:longitude_of_central_meridian = -14.1000003814697 ;
Lambert_Conformal:_CoordinateTransformType = "Projection" ;
Lambert_Conformal:_CoordinateAxisTypes = "GeoX GeoY" ;
这与 R 包使用的分辨率 = 12 投影有点不同。因此,制作一个合规的:
p = "+proj=lcc +lat_1=43 +lat_2=43 +lat_0=24.2280006408691 +lon_0=-14.1000003814697 +x_0=2182629.35 +y_0=-269655.97 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=km +no_defs"
wrf <- getRaster('temp', day = testDay, resolution=36)
projection(wrf) = p
然后我用 EU 覆盖进行测试......
请注意,这次我重新预测了欧盟。正在做projectRaster
到 latlong 也可以工作:
> wrfLL = projectRaster(wrf, crs="+init=epsg:4326")
> plot(wrfLL[[1]])
> abline(v=0)
格林威治子午线处于其应有的位置。