Hi Kyle:
I used the FedData to get nlcd for several States. Pretty cool.
Now I am trying to download DAYMET data and, to be honest, I do not know why this problem is occurring but I thought it would be useful to check it is not the package. The script downloads the data, but at the time it starts to extract I get the error. I would really really appreciate if you could take a look. The error does not occur as a result of an obvious problem in one single tile.
Two things. I am using the package hosted in Github and secondly, there is a message in the nc's. I leave the script
library(rgdal)
library(gdalUtils)
library(raster)
library(maptools)
library(sp)
library(shapefiles)
library(devtools)
library(scales)
library(ncdf4)
library(FedData)
#Downloaded online from us census https://www.census.gov/geographies/mapping-files/time-series/geo/carto-boundary-file.html
states=readOGR("/your_path/", layer = "cb_2018_us_state_5m") #s
PA=states[states@data$STUSPS=="PA",]
plot(PA)
# The projection of the DAYMET data is presented here https://daac.ornl.gov/DAYMET/guides/Daymet_V3_CFMosaics.html#datacharact
lambert.conformal.conic <- CRS("+proj=lcc +lat_1=25 +lat_2=60 +lat_0=42.5 +lon_0=-100 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs")
PA.proj=spTransform(PA, lambert.conformal.conic)
test=get_daymet(template=PA.proj,
label = paste0("day_met", "PA"),
element="prcp",
years = c(2010),
force.redo=T) #FedData package
test=nc_open("your_path_to_Files/RAW/DAYMET/11931/prcp.nc")
test
Check the test dimension characteristics. It says: "[1] 'vobjtovarid4: **** WARNING **** I was asked to get a varid for dimension named nv BUT this dimension HAS NO DIMVAR! Code will probably fail at this point'"
Thanks in advance
Diego
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4