[gdal-dev] [newbie] regridding netCDF

Tom Roche Tom_Roche at pobox.com
Mon Jul 9 16:16:29 PDT 2012


summary: new user with limited geoinformatic background (including
jargon familiarity) wants to know if any of the GDAL tools (e.g.,
gdaltransform) can "regrid" netCDF data, such that data attributed to
each gridcell in a source dataset at a given {projection, resolution}
tuple is appropriately allocated to each gridcell in a target dataset
at a different {projection, resolution}. If so, please give or point
to an example.

details:

I'd appreciate advice suitability of GDAL tools for transforming values
contained in one 2-d netCDF spatial grid in one {projection, resolution}
to another, or pointers to other resources to consult. I see GDAL
supports netCDF as a "raster format"

http://gdal.org/formats_list.html

but don't know if the following counts as "image processing" (which is
what I tend to associate GDAL with), or if that matters. Further
apologies in advance if the following is too verbose (I'm guessing this
is a common usecase--or hoping :-) or uses the wrong vocabulary (please
correct where needed), but here's what I mean:

I want to take emissions of gaseous N2O from various sources in the
contiguous US (CONUS) and input them to a particular atmospheric model
(CMAQ). My group has other inputs (notably, meteorology) for our
usecase on a 12 km x 12 km Lambert conformal conic (LCC) grid, and the
model wants everything in one format, so I need to get our N2O
emissions on that 12x12 LCC. However, most of my N2O emissions will be
coming from global inventories (EDGAR, GEIA, GFED), which are all
gridded lat-lon, but at differing resolutions (e.g., 0.1x0.1°). So
I need to "regrid" the latter emissions, such that the mass emitted in
each lat-lon gridcell in the source dataset is appropriately
attributed (probably by spatial averaging) to one or more gridcells in
the 12-km LCC target dataset (the extents of which are a subset of the
global source dataset).

How to do programmatically (driven by bash, python, R) on linux?

Your coding recommendations are appreciated, as would be pointers to
helpful resources, esp example code. TIA, and feel free to forward,
Tom Roche <Tom_Roche at pobox.com>


More information about the gdal-dev mailing list