如何將投影的netcdf文件轉換為柵格以在R中進行繪圖

Question

我有一個雷擊的netcdf文件，其中包含x，y和緯度Lon的變量。xy數據的坐標參考系統為albers_conical_equal_area，並且還有許多其他變量，例如虛假東移，虛假北移等。投影。 我想將文件轉換為柵格，並且需要填寫crs參數。 我想以投影格式顯示柵格，但不知道如何正確格式化crs條目。

這是否有可能（柵格僅適用於經緯度數據嗎？）

如果是這樣，我如何轉換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")

我不確定如何將此信息應用於raster命令

strikes <- raster(t(lightning.array)
        , xmn=min(x)
        , xmx=max(x)
        , ymn=min(y)
        , ymx=max(y)
        , crs=CRS("missing code here ")
)

我使用的文件在美國各地是一天的雷擊，並且我預計會導入很多這樣的文件，以便在整個感興趣的地區匯總全年的雷擊，因此我想弄清楚這個項目的問題首先在此文件上。

Answer 1

使用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圓錐形相等的區域中，我不得不用Google進行搜索。 這是我所做的工作，以供將來參考。 我用谷歌搜索“ albers_conical_equal_area proj4string”，這把我帶到了這里 。 然后我抓住了proj4string。 而已！

Answer 2

我會做

library(raster)
b <- brick("file.nc")

那應該滿足您的所有要求，甚至更多。