[R-sig-Geo] Selecting variable within nc file in raster() or brick()
Agustin Lobo
alobolistas at gmail.com
Fri Jan 9 12:36:55 CET 2015
Hi!
For nc (v3) files having only 1 field variable along time such as
open.ncdf(file.path(ncdir,nom1))
[1] "file /media/alobo/FREECOM
HDD/BIGFILES/IMASCOARSESTOCK/SMOS/_DELSMH001DIBEA.nc has 4
dimensions:"
[1] "time Size: 257"
[1] "lat Size: 1100"
[1] "lon Size: 1600"
[1] "time_strlen Size: 20"
[1] "------------------------"
[1] "file /media/alobo/FREECOM
HDD/BIGFILES/IMASCOARSESTOCK/SMOS/_DELSMH001DIBEA.nc has 1 variables:"
[1] "float HighRes_Soil_Moisture[lon,lat,time] Longname:1 km Soil
Moisture Missval:-999"
I can simply use
nc1 <- brick(file.path(ncdir,nom1))
But, what in the case of having several variables in the nc file, i.e.:
open.ncdf(file.path(ncdir,nom5))
[1] "file /media/alobo/FREECOM
HDD/BIGFILES/IMASCOARSESTOCK/SMOS/_NRTSM003D025A_ALL_2013.nc has 4
dimensions:"
[1] "time Size: 365"
[1] "lat Size: 44"
[1] "lon Size: 62"
[1] "time_strlen Size: 20"
[1] "------------------------"
[1] "file /media/alobo/FREECOM
HDD/BIGFILES/IMASCOARSESTOCK/SMOS/_NRTSM003D025A_ALL_2013.nc has 4
variables:"
[1] "int L2_Points[lon,lat,time] Longname:Number of Measures Missval:0"
[1] "float SM[lon,lat,time] Longname:Surface soil moisture Missval:-999"
[1] "float SM_DQX[lon,lat,time] Longname:Surface soil moisture_dqx
Missval:-999"
[1] "float VARIANCE_SM[lon,lat,time] Longname:Variance of Surface
soil moisture Missval:-999"
Is there a way to select the desired variable within brick() ?
Or should I use get.var.ncdf() and then transform to raster?
(the only problem in this case is that I have to rotate etc, while
brick() does it all in the first case)
Thanks
Agus
More information about the R-sig-Geo
mailing list