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