在 R 中正确地在美国投影 NetCDF 文件？

问题描述

我正在尝试使用来自此网站的每小时降水 (MPE) 数据：http://convection.meas.ncsu.edu:8080/thredds/catalog/ncep/mpe/netcdf/stage4/201809/catalog.html

但是，我真的很纠结如何让投影在美国上空正确对齐。你能提供的任何帮助都会很棒。根据我的理解，投影是极地立体投影，但它不容易转化为更常见的投影。 NWS 和 NOAA 在他们的网站上没有明确说明如何解决这个问题。找了我校园的 GIS 专家，他也无法帮助解决这个问题。

library(tidyverse)
library(raster)
library(rgdal)
library(sp)
library(mapview)
library(here)

file_list <- list.files("201809")


r<-raster::brick("201809//ST4_2018091613_01h.nc",varname="APCP_SFC")[[1]]

# create variable for projection of imported raster
prj <- projection(r)

# create a variable for the extent as imported
ext <- extent(r)

wgs<-"+proj=longlat +datum=wgs84 +ellps=wgs84 +no_defs"
wgsraster <- projectRaster(r,crs=wgs)

I came across another post on here asking a similar question but it didn't quite solve the problem. You can see how the data,especially in the northern half of the country,just isn't aligning correctly

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）