如何将投影的netcdf文件转换为光栅以在R中打印

如何将投影的netcdf文件转换为光栅以在R中打印,r,netcdf,r-raster,map-projections,R,Netcdf,R Raster,Map Projections,我有一个netcdf雷击文件,其中包含x、y和Lat Lon的变量,xy数据的坐标参考系是albers_Continary_equal_area,还有许多附加变量假东距、假北距等,与投影一起。我想将文件转换为光栅,我需要填写一个crs参数。我希望光栅以投影格式显示,但不知道如何正确设置crs条目的格式 这是否可能(光栅是否仅适用于lat-lon数据?) 如果是,我如何转换netcdf文件中的信息以使用raster命令中的crs参数 我试图从netcdf文件中获取属性,但不清楚如何使用它们 我已经

我有一个netcdf雷击文件,其中包含x、y和Lat Lon的变量,xy数据的坐标参考系是albers_Continary_equal_area,还有许多附加变量假东距、假北距等,与投影一起。我想将文件转换为光栅,我需要填写一个crs参数。我希望光栅以投影格式显示,但不知道如何正确设置crs条目的格式

这是否可能(光栅是否仅适用于lat-lon数据?)

如果是,我如何转换netcdf文件中的信息以使用raster命令中的crs参数

我试图从netcdf文件中获取属性,但不清楚如何使用它们

我已经使用下面的代码从crs变量获取属性,但是

grid_mapping_name <- ncatt_get(lightning, "crs", "grid_mapping_name")
standard_parallel <- ncatt_get(lightning, "crs", "standard_parallel")
longitude_of_central_meridian <- ncatt_get(lightning, "crs", 
"longitude_of_central_meridian")
latitude_of_projection_origin <- ncatt_get(lightning, "crs", 
"latitude_of_projection_origin")
false_easting <- ncatt_get(lightning, "crs", "false_easting")
false_northing <- ncatt_get(lightning, "crs", "false_northing")
grid\u mapping\u名称使用
proj4string
尝试
CRS(“+proj=aea+lat_1=29.5+lat_2=45.5+lat_0=37.5+lon_0=-96+x_0=0+y_0=0+ellps=GRS80+datum=NAD83+units=m+no_defs”)

您要查找的是名为
proj4string
,它是对点的坐标参考系进行编码的字符串

例如,
CRS(“+proj=longlat+datum=WGS84”)
sp
包中的CRS对象,其
proj4string
为“+proj=longlat+datum=WGS84”。许多数据都来自于熟悉的长/宽CRS

你的数据是在阿尔伯斯圆锥面积相等的地方,我不得不用谷歌搜索项目。以下是我所做的,以供将来参考。我在谷歌上搜索了“albers_Continary_equal_area proj4string”,这让我大吃一惊。然后我抓起proj4string。就这样

使用
proj4string
尝试
CRS(“+proj=aea+lat_1=29.5+lat_2=45.5+lat_0=37.5+lon_0=-96+x_0=0+y_0=0+ellps=GRS80+datum=NAD83+units=m+no_defs”)

您要查找的是名为
proj4string
,它是对点的坐标参考系进行编码的字符串

例如,
CRS(“+proj=longlat+datum=WGS84”)
sp
包中的CRS对象,其
proj4string
为“+proj=longlat+datum=WGS84”。许多数据都来自于熟悉的长/宽CRS

你的数据是在阿尔伯斯圆锥面积相等的地方,我不得不用谷歌搜索项目。以下是我所做的,以供将来参考。我在谷歌上搜索了“albers_Continary_equal_area proj4string”,这让我大吃一惊。然后我抓起proj4string。就这样

我愿意

library(raster)
b <- brick("file.nc")
库(光栅)
b我愿意

library(raster)
b <- brick("file.nc")
库(光栅)
B