Hi John,<br><br>Glad to hear you got output.  Unfortunately, no, I don&#39;t have a dissolve function written that can be plugged into this.  I had started to write one last year, as dissolving raster extents was something I needed to do occasionally, but I always ended up doing that step in ArcGIS instead.<br>
<br>Roger<br>--<br><br><div class="gmail_quote">On Thu, Apr 9, 2009 at 4:43 AM, John Zastrow <span dir="ltr">&lt;<a href="mailto:jcz@northredoubt.com">jcz@northredoubt.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="im">Roger André wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi John,<br>
<br>
I wrote a very un-optimized tool to do this a little while ago.  It&#39;s a Python script that you can probably run in FWTools.  Just be sure that you put it in your FWTools/bin directory, and then that you create an appropriate environment setting script with the same name in FWTools/bin_safe.  Just do a search for gdal_merge.py and you&#39;ll see what I mean.<br>

<br>
You&#39;ll need to do some editing of the script to specify the pixel values you want to vectorize, but that should be pretty straight-forward.<br>
<br>
Good luck, and let me know if you need any help getting it to run.<br>
<br>
Roger<br>
--<br>
<br>
</blockquote>
<br>
<br></div>
Hi Roger,<br>
<br>
Thank you very much. Your script did produce a shape output. Do you know if there is a method to do the &quot;dissolve&quot; on the polygons either inline with the vectorizing or as secondary step? Googling turns up people searching for such a tool in the GDAL/OGR world but, I did not see any final conclusions.<br>

<br>
Thank you,<br><font color="#888888">
John<br>
<br>
</font></blockquote></div><br>