如何解决为R
我尝试将此df (c.data)
转换为RasterBrick并将crs
设置为epsg:3035
;带有代码的LAmbert Equal Area Europe。该代码创建了栅格,但未设置CRS
c.data
# A tibble: 361,478 x 6
x y pop women mean_age hh_size
<int> <int> <int> <int> <int> <int>
1 4334500 2684500 NA NA NA NA
2 4335500 2684500 NA NA NA NA
3 4336500 2684500 NA NA NA NA
4 4337500 2684500 NA NA NA NA
5 4338500 2684500 NA NA NA NA
6 4339500 2684500 NA NA NA NA
7 4334500 2685500 NA NA NA NA
8 4335500 2685500 NA NA NA NA
9 4336500 2685500 NA NA NA NA
10 4337500 2685500 NA NA NA NA
# ... with 361,468 more rows
c.data_ras <- rasterFromXYZ(c.data,crs = st_crs(3035)$proj4string)
c.data_ras
class : RasterBrick
dimensions : 868,642,557256,4 (nrow,ncol,ncell,nlayers)
resolution : 1000,1000 (x,y)
extent : 4031000,4673000,2684000,3552000 (xmin,xmax,ymin,ymax)
crs : NA
source : memory
names : pop,women,mean_age,hh_size
min values : 1,1,1
max values : 6,5,5
相反,我收到此警告
警告消息:
In .newCRS(value) :
+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs is not a valid PROJ.4 CRS string
请,我该怎么办?
谢谢
感谢@lime和@sequoia,我尝试了两个建议,但我仍然遇到相同的错误
c.data_ras <- rasterFromXYZ(c.data,crs = st_crs("+init=epsg:3035")$proj4string)
Warning message:
In .newCRS(value) :
+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs is not a valid PROJ.4 CRS string
和
c.data_ras <- rasterFromXYZ(c.data,crs = st_crs("epsg:3035")$proj4string)
Warning message:
In .newCRS(value) :
+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs is not a valid PROJ.4 CRS string
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。