r - R 栅格函数将不接受 WKT 格式的 crs
问题描述
我正在尝试生成一个栅格并为其分配一个 CRS 投影。但是,我的 CRS 采用新的 WKT 格式,并且 raster() 函数要求我提供 proj4string。到目前为止,这是我的代码:
library(sf)
library(raster)
crs_dem <- st_crs(
'PROJCS["NAD_1983_2011_StatePlane_California_II_FIPS_0402",
GEOGCS["GCS_NAD_1983_2011",
DATUM["D_NAD_1983_2011",
SPHEROID["GRS_1980",6378137.0,298.257222101]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Lambert_Conformal_Conic"],
PARAMETER["False_Easting",2000000.0],
PARAMETER["False_Northing",500000.0],
PARAMETER["Central_Meridian",-122.0],
PARAMETER["Standard_Parallel_1",38.33333333333334],
PARAMETER["Standard_Parallel_2",39.83333333333334],
PARAMETER["Latitude_Of_Origin",37.66666666666666],
UNIT["Meter",1.0]]')
ext <- extent(1895000, 1935000, 579500, 616500)
grid <- raster(ext, resolution = c(40,40), crs = crs(dem))
上面的代码生成一个 crs = NA 的栅格。我也试过用 crs(grid)<- 和 projection(grid)<- 分配它,但没有运气。如何让我的特定 CRS 文件与此栅格相关联?
解决方案
@slamballais 的回答成功了!昨晚我还通过反复试验找到了另一种(不太干净)的方法,所以这里有两种解决方案。
选项1:
test <- sp::CRS(crs_dem$input)
grid <- raster(ext, resolution = c(40,40), crs = test)
选项 2:
library(dplyr)
aoi <- ext %>%
as('SpatialPolygons') %>%
st_as_sf %>%
st_set_crs(crs_dem)
grid <- raster(ext, resolution = c(40,40), crs = projection(aoi))
推荐阅读
- c# - 将连接服务添加到 nuget 包
- java - Firebase -- 文档引用必须有偶数个段,但 Rezepte 有 1
- react-native - React Navigation 5 如何映射抽屉堆栈屏幕?
- r - 合并和左连接不断在 R 中重复行项目
- scala - 覆盖 toString 时如何使用条件语句
- python - 将 PCollection 作为 Apache Beam 的侧输入传递时出现 KeyError
- node.js - 在 nodejs 中运行 rust 程序
- oracle - ORA-20000: 找不到 AOP 服务器
- c# - 比较表示日期的整数
- css - 无法更改我的 Wordpress 标题上的 CSS