<div>Tom,</div><div>Would the GRASS 3D Raster (Voxel) functions meet your needs?</div><div><br></div><div><a href="http://grass.osgeo.org/grass64/manuals/html64_user/raster3D.html">http://grass.osgeo.org/grass64/manuals/html64_user/raster3D.html</a></div>

<div><br></div><div>Doug</div><br><div class="gmail_quote">On Wed, Nov 14, 2012 at 12:11 AM, Tom Roche <span dir="ltr"><<a href="mailto:Tom_Roche@pobox.com" target="_blank">Tom_Roche@pobox.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
summary: I'd appreciate advice regarding tools and methods for<br>
transforming values from an unprojected global 3D spatial grid to a<br>
projected 3D grid with different horizontal and vertical resolution, or<br>
pointers to other resources to consult.<br>
<br>
details:<br>
<br>
I have output from a global atmospheric model that I'd like to use as<br>
initial and boundary conditions for a regional model. The global input<br>
netCDF has dimensions=2.5° lon x 1.875° lat x 56 vertical levels. The<br>
regional model runs over North America using a 12-km grid projected<br>
Lambert Comformal Conic (LCC), with 34 vertical levels. Since its top<br>
height is less than that of the global input, the extents of the output<br>
domain are fully contained within the input domain.<br>
<br>
Each box or voxel defined by the global input grid contains an estimate<br>
for the N2O concentration for that volume. From those I want to compute<br>
the concentrations for each output gridbox volume. I'd appreciate your<br>
recommendations for tools that can handle this usecase. Particularly, is<br>
this something that gdalwarp can do?<br>
<br>
FWIW, the best tool I've seen so far is R package=gstat, but (IIUC)<br>
<br>
- gstat expects projected input. I'm not sure if I can work around that<br>
  for this usecase.<br>
<br>
- as the name implies, 'gstat' is doing geostatistics, e.g., variogram-<br>
  and covariance-based modeling. I'm not sure either how to setup the<br>
  distance weighting for my scenario, and, frankly, I remain unconvinced<br>
  that a statistical approach is necessary for this application (though<br>
  it may be a sufficient or the best-available approach). This may be<br>
  due to my statistical ignorance, however.<br>
<br>
your assistance is appreciated, Tom Roche <<a href="mailto:Tom_Roche@pobox.com">Tom_Roche@pobox.com</a>><br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a></blockquote></div><br><br clear="all"><div><br></div>-- <br><div>Doug Newcomb</div><div>USFWS</div>

<div>Raleigh, NC</div><div>919-856-4520 ext. 14 <a href="mailto:doug_newcomb@fws.gov" target="_blank">doug_newcomb@fws.gov</a></div><div>---------------------------------------------------------------------------------------------------------</div>

<div>The opinions I express are my own and are not representative of the official policy of the U.S.Fish and Wildlife Service or Dept. of the Interior.   Life is too short for undocumented, proprietary data formats.</div>

<br>