Hi,
I'm trying to import multidimensional netCDF files slicing the bands
according to time. I want to use the new temopral module in GRASS70
for further processing.
My process are as follows:
1. Download the netCDF usinf ncks, i.e.:
$ ncks -O -v r -d time,${MIN},${MAX} -d lon,115.00,155.00 -d lat,4.00,22.00 \
http://disc2.nascom.nasa.gov:80/dods/3B42RT_V7_rainrate \
3hourly_trmm_2.nc
This is a 3hourly rainfal data.
2. Running ncdump, the data are as follows:
$ ncdump -ct 3hourly_trmm_2.nc
netcdf 3hourly_trmm_2 {
dimensions:
lat = 72 ;
lon = 160 ;
time = 89 ;
variables:
double lat(lat) ;
lat:grads_dim = "y" ;
lat:grads_mapping = "linear" ;
lat:grads_size = "480" ;
lat:units = "degrees_north" ;
lat:long_name = "latitude" ;
lat:minimum = -59.875 ;
lat:maximum = 59.875 ;
lat:resolution = 0.25f ;
double lon(lon) ;
lon:grads_dim = "x" ;
lon:grads_mapping = "linear" ;
lon:grads_size = "1440" ;
lon:units = "degrees_east" ;
lon:long_name = "longitude" ;
lon:minimum = 0.125 ;
lon:maximum = 359.875 ;
lon:resolution = 0.25f ;
float r(time, lat, lon) ;
r:_FillValue = -99999.f ;
r:missing_value = -99999.f ;
r:long_name = "hourly rain rate (mm/hr) " ;
double time(time) ;
time:grads_dim = "t" ;
time:grads_mapping = "linear" ;
time:grads_size = "131078" ;
time:grads_min = "00z01mar2000" ;
time:grads_step = "3hr" ;
time:units = "days since 1-1-1 00:00:0.0" ;
time:long_name = "time" ;
time:minimum = "00z01mar2000" ;
time:maximum = "15z08jan2045" ;
time:resolution = 0.125f ;
// global attributes:
:title = "Three Hourly TRMM and Other Satellite Rainfall (3B42RT)" ;
:Conventions = "COARDS\n",
"GrADS" ;
:dataType = "Grid" ;
:history = "Tue Nov 11 15:30:44 2014: ncks -O -v r -d
time,41945,42033 -d lon,115.00,155.00 -d lat,4.00,22.00
http://disc2.nascom.nasa.gov:80/dods/3B42RT_V7_rainrate
3hourly_trmm_2.nc\n",
"Mon Nov 10 23:05:02 GMT-00:00 2014 : imported by GrADS
Data Server 2.0" ;
:NCO = "4.4.2" ;
data:
lat = 4.125, 4.375, 4.625, 4.875, 5.125, 5.375, 5.625, 5.875, 6.125, 6.375,
6.625, 6.875, 7.125, 7.375, 7.625, 7.875, 8.125, 8.375, 8.625, 8.875,
... ;
lon = 115.125, 115.375, 115.625, 115.875, 116.125, 116.375, 116.625,
116.875, 117.125, 117.375, 117.625, 117.875, 118.125, 118.375, 118.625,
... ;
time = "2014-07-09 03", "2014-07-09 06", "2014-07-09 09", "2014-07-09 12",
"2014-07-09 15", "2014-07-09 18", "2014-07-09 21", "2014-07-10",
... ;
}
3. I want to import it in a GRASS database where the raster filename
correspond to the date/time.
script.run_command('r.in.gdal', input='3hourly_trmm_2.nc', output='g',
flags='o')
script.list_pairs(type='rast')
[('g.1', 'PERMANENT'),
...
('g.85', 'PERMANENT')]
As you can see, the band/raster names were just an autoincrement of
the netCDF layers.
Any advice on doing this in Python? Thanks!
--
cheers,
maning
------------------------------------------------------
"Freedom is still the most radical idea of all" -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
------------------------------------------------------
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
