Thanks Doug, this is great feedback. It is folks like you, and the others who responded, that making working w/ open sources tools a pretty awesome alternative.<br><br>Cheers, <br>Derek<br><br><div class="gmail_quote">On Fri, Mar 23, 2012 at 7:50 AM, <span dir="ltr"><<a href="mailto:Doug_Newcomb@fws.gov">Doug_Newcomb@fws.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br><font face="sans-serif">Derek,</font>
<br><font face="sans-serif"> You
can do this in GRASS ( using gdal in the background) by doing the
following:</font>
<br><font face="sans-serif">1) create a workspace in the same projection
as your existing raster ( when you first start GRASS it will ask
you to set up a workspace with projection information)</font>
<br><font face="sans-serif">2) import the existing raster
layer into GRASS, ( r.in.gdal,<a href="http://grass.fbk.eu/grass64/manuals/html64_user/r.in.gdal.html" target="_blank">http://grass.fbk.eu/grass64/manuals/html64_user/r.in.gdal.html</a>)</font>
<br><font face="sans-serif">3)set the region to match the imported
raster layer in extent and resolution (g.region, </font><a href="http://grass.fbk.eu/grass64/manuals/html64_user/g.region.html" target="_blank"><font face="sans-serif">http://grass.fbk.eu/grass64/manuals/html64_user/g.region.html)
</font></a>
<br><font face="sans-serif">4) import the county vector layer (v.in.ogr,
</font><a href="http://grass.fbk.eu/grass64/manuals/html64_user/v.in.ogr.html" target="_blank"><font face="sans-serif">http://grass.fbk.eu/grass64/manuals/html64_user/v.in.ogr.html
)</font></a>
<br><font face="sans-serif">5) convert the vector layer to a raster
layer(v.to.rast, </font><a href="http://grass.fbk.eu/grass64/manuals/html64_user/v.to.rast.html" target="_blank"><font face="sans-serif">http://grass.fbk.eu/grass64/manuals/html64_user/v.to.rast.html)
</font></a>
<br>
<br><font face="sans-serif">When you create the new raster it will
use the existing alignment and resolution of the region ( based off of
the existing raster) in the placement of the cells. Even though
I reference the tabular commands in the above process, it can be easily
done through both GRASS gui's</font>
<br>
<br><font face="sans-serif">Doug</font>
<br><font face="sans-serif"> </font>
<br>
<br>
<br><font face="sans-serif">Doug Newcomb
<br>
USFWS<br>
Raleigh, NC<br>
<a href="tel:919-856-4520%20ext.%2014" value="+19198564520" target="_blank">919-856-4520 ext. 14</a> <a href="mailto:doug_newcomb@fws.gov" target="_blank">doug_newcomb@fws.gov</a><br>
---------------------------------------------------------------------------------------------------------<br>
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.</font>
<br>
<br>
<br>
<p></p><table width="100%">
<tbody><tr valign="top">
<td width="40%"><font face="sans-serif" size="1"><b>Derek Morgan <<a href="mailto:jdmorgan@unca.edu" target="_blank">jdmorgan@unca.edu</a>></b>
</font>
<br><font face="sans-serif" size="1">Sent by: <a href="mailto:gdal-dev-bounces@lists.osgeo.org" target="_blank">gdal-dev-bounces@lists.osgeo.org</a></font>
<p><font face="sans-serif" size="1">03/22/2012 10:12 PM</font>
</p></td><td width="59%">
<table width="100%">
<tbody><tr valign="top">
<td>
<div align="right"><font face="sans-serif" size="1">To</font></div>
</td><td><font face="sans-serif" size="1"><a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a></font>
</td></tr><tr valign="top">
<td>
<div align="right"><font face="sans-serif" size="1">cc</font></div>
</td><td>
</td></tr><tr valign="top">
<td>
<div align="right"><font face="sans-serif" size="1">Subject</font></div>
</td><td><font face="sans-serif" size="1">[gdal-dev] shifting (or snapping) one
raster to another</font></td></tr></tbody></table>
<br>
<table>
<tbody><tr valign="top">
<td>
</td><td></td></tr></tbody></table>
<br></td></tr></tbody></table>
<br>
<br>
<br><tt><font><div><div class="h5">Hello,<br>
<br>
I am in a jam with ArcGIS, and am hoping there is a easier way to do<br>
this with GDAL. Here is the situation: I am converting a vector<br>
dataset of county polygons to raster at a pixel resolution of 231 m.<br>
However, I want this to line up with a different raster dataset I have<br>
at a much larger extent, but also the same resolution 231 m. As you<br>
probably guessed the two datasets are slightly askew as they have<br>
different extents. I would like to align the newly created county<br>
raster with the larger raster data set by shifting (or snapping) it to<br>
the closest pixels. Hopefully, this makes enough sense for someone
to<br>
let me know if this can be done in GDAL.<br>
<br>
Thanks ahead of time,<br>
Derek<br></div></div><div class="im">
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
</div></font></tt><a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank"><tt><font>http://lists.osgeo.org/mailman/listinfo/gdal-dev<br>
</font></tt></a>
<br></blockquote></div><br>