微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

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

如何解决在 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 举报,一经查实,本站将立刻删除。