Hi Guys,
+1, I recommend gdalwarp as well. We've found it to be an effective
regridding tool, especially for uniform grids.
Cheers,
Chris
On Jul 9, 2012, at 7:03 AM, Kokron, Daniel S. (GSFC-610.1)[Computer Sciences
Corporation] wrote:
> I recommend gdalwarp from the GDAL package.
>
> Daniel Kokron
> NASA Ames (ARC-TN)
> SciCon group
> 301-286-3959
>
> ________________________________________
> From: netcdfgroup-bounces@xxxxxxxxxxxxxxxx
> [netcdfgroup-bounces@xxxxxxxxxxxxxxxx] On Behalf Of Tom Roche
> [Tom_Roche@xxxxxxxxx]
> Sent: Friday, July 06, 2012 5:43 PM
> To: netcdfgroup@xxxxxxxxxxxxxxxx
> Subject: [netcdfgroup] spatial "regridding"
>
> I'd appreciate advice regarding tools and methods for transforming
> values "contained" in one 2-d netCDF spatial grid in one {projection,
> resolution} to another, or pointers to other resources to consult.
> Apologies if the following is too verbose (I'm guessing this is a
> common usecase) or uses the wrong vocabulary (please correct where
> needed), but here's what I mean:
>
> I want to take N2O emissions 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 gridded lat-lon at
> differing resolutions (e.g., 0.1x0.1°). So I need to "regrid" the
> latter emissions. How to do? on linux, preferably programmatically.
> Details:
>
> EDGAR 4.2 gives the mass of N2O emitted from each 0.1x0.1° lat-lon
> gridcell over the globe. From that I can subset the portion that
> intersects the AQMEII North American grid
>
> https://github.com/TomRoche/cornbeltN2O/wiki/AQMEII-North-American-domain
>
> (Pretty sure I can do that with NCO.) I can then use R package=M3
>
> http://cran.r-project.org/web/packages/M3/
>
> to project the nodes of the EDGAR subgrid to LCC. Note that, @ 40° N
> (middle latitude of this subgrid), each EDGAR gridcell's extents are
> ~11 km N-S x ~8.5 km E-W.
>
> Over that "LCC EDGAR" subgrid I would (conceptually) overlay a 12x12
> km LCC CMAQ grid. (Which of course will deviate slightly from those
> exact dimensions--we have the exact grid centers, nodes, and areas.)
> Given
>
> 1 the LCC EDGAR grid's extents are larger than the CMAQ grid's extents
> (which I can make true)
>
> 2 each CMAQ gridcell is larger than each LCC EDGAR gridcell, over the
> CMAQ grid's extents (which I believe is true, but should check)
>
> 3 each gridcell in either has no internal structure, so emissions
> allocated to each gridcell are presumed to be "smeared evenly" by
> area over the gridcell
>
> each CMAQ gridcell will overlay between 2 and n LCC EDGAR gridcells
> (guessing that n=4 in this case), and I can reallocate emissions
> area-proportionally from each overlaid LCC EDGAR gridcell to the CMAQ
> gridcell. This presumes that I can, programmatically,
>
> 1 identify each underlay (LCC EDGAR) gridcell that an overlay (CMAQ)
> gridcell intersects
>
> and, for each of the underlay gridcells U_i (with i in 1..n)
> intersected by a given overlay gridcell O,
>
> 2 identify the emissions attributed to U_i e_{U_i} (easy, just lookup)
>
> 3 identify the total area of U_i a_{U_i} (should be an easy
> computation, given lookup of its corners/gridnodes)
>
> 4 identify the area of intersection I_i between U_i and O (I don't
> know how to do that)
>
> whereupon I can attribute to O emissions e_O =
>
> \sum_{i=1}^{n} \frac{I_i U_i}{a_{U_i}}
>
> This sounds doable, provided task 4 can be done programmatically, but
> it also sounds error-prone, and I've got a lot of other work to do on
> this project, so I'm more than willing to reuse existing code (esp R,
> since that would fit with my other scripts).
>
> Your coding recommendations are appreciated, as would be pointers to
> helpful resources. Thanks in advance, and feel free to forward,
> Tom Roche <Tom_Roche@xxxxxxxxx>
>
> _______________________________________________
> netcdfgroup mailing list
> netcdfgroup@xxxxxxxxxxxxxxxx
> For list information or to unsubscribe, visit:
> http://www.unidata.ucar.edu/mailing_lists/
> --
> This message has been scanned for viruses and dangerous
> content by GMAO's MailScanner, and is believed to be clean.
>
> _______________________________________________
> netcdfgroup mailing list
> netcdfgroup@xxxxxxxxxxxxxxxx
> For list information or to unsubscribe, visit:
> http://www.unidata.ucar.edu/mailing_lists/
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Chris Mattmann, Ph.D.
Senior Computer Scientist
NASA Jet Propulsion Laboratory Pasadena, CA 91109 USA
Office: 171-266B, Mailstop: 171-246
Email: chris.a.mattmann@xxxxxxxx
WWW: http://sunset.usc.edu/~mattmann/
Phone: +1 (818) 354-8810
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Adjunct Assistant Professor, Computer Science Department
University of Southern California, Los Angeles, CA 90089 USA
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
