Hi Jerl,<br><br>  If you have python with matplotlib and the gdal bindings installed, you can use gdal to read the DEM (which I assume you are getting gdal_contour to draw contours from), and then matplotlibs contour-function to draw the labelled contours. You can draw  on top of an imported map/image-file. I know that matplotlib will read at least png&#39;s, and I think also other image formats (with the right packages installed).<br>
<br> Cheers,<br> Simon<br><br><br><div class="gmail_quote">On Mon, May 16, 2011 at 10:08 PM, Jerl Simpson <span dir="ltr">&lt;<a href="mailto:jsimpson@wxtrends.com">jsimpson@wxtrends.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Ok, thank you.<br>
<br>
Does anyone else have any suggestions?<br>
<br>
<br>
<br>
On Mon, May 16, 2011 at 3:06 PM, Chaitanya kumar CH<br>
<div><div></div><div class="h5">&lt;<a href="http://chaitanya.ch" target="_blank">chaitanya.ch</a>@<a href="http://gmail.com" target="_blank">gmail.com</a>&gt; wrote:<br>
&gt; Jerl,<br>
&gt;<br>
&gt; I&#39;m afraid gdal_rasterize can&#39;t do the labelling. You will need a<br>
&gt; specialized renderer for that. I&#39;m not familiar with any of the desktop<br>
&gt; applications that do this. Perhaps, qgis has a plugin.<br>
&gt;<br>
&gt; On Tue, May 17, 2011 at 1:30 AM, Jerl Simpson &lt;<a href="mailto:jsimpson@wxtrends.com">jsimpson@wxtrends.com</a>&gt; wrote:<br>
&gt;&gt;<br>
&gt;&gt; I already have the resultant image, it looks great.   Just drawing<br>
&gt;&gt; lines on a map.  I&#39;m using a blank image so I can more clearly see<br>
&gt;&gt; what&#39;s going on.  It&#39;s just a black and white image, and that&#39;s all I<br>
&gt;&gt; want.<br>
&gt;&gt; I&#39;m trying to get the contour labeled, so I know which one I&#39;m<br>
&gt;&gt; viewing.  So on the 1000 contour, I should see 1000 on the line<br>
&gt;&gt; somewhere, and if it&#39;s a 547 contour, I should see 547 on the line<br>
&gt;&gt; somewhere.<br>
&gt;&gt;<br>
&gt;&gt; Is that possible?<br>
&gt;&gt;<br>
&gt;&gt; Thanks,<br>
&gt;&gt;<br>
&gt;&gt; Jerl<br>
&gt;&gt;<br>
&gt;&gt; On Mon, May 16, 2011 at 2:41 PM, Chaitanya kumar CH<br>
&gt;&gt; &lt;<a href="http://chaitanya.ch" target="_blank">chaitanya.ch</a>@<a href="http://gmail.com" target="_blank">gmail.com</a>&gt; wrote:<br>
&gt;&gt; &gt; Jerl,<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; The contour attribute and the Z values are same here. Just one of them<br>
&gt;&gt; &gt; is<br>
&gt;&gt; &gt; enough.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; With the Z values near the magnitude of 1000, the burn value of 255 is<br>
&gt;&gt; &gt; not<br>
&gt;&gt; &gt; going to make much of a difference. Use it if you want to offset any<br>
&gt;&gt; &gt; negative Z values.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; You should mention contour.shp instead of just &#39;contour&#39; for<br>
&gt;&gt; &gt; &lt;src_datasource&gt;.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; Make sure BLANK.tif doesn&#39;t have different extents/projection.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; Your final image will have a data type of Float64. You may have to run<br>
&gt;&gt; &gt; some<br>
&gt;&gt; &gt; contrast stretching operations on it.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; On Tue, May 17, 2011 at 12:04 AM, Jerl Simpson &lt;<a href="mailto:jsimpson@wxtrends.com">jsimpson@wxtrends.com</a>&gt;<br>
&gt;&gt; &gt; wrote:<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; Hi guys and gals:<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; I have something that&#39;s evading me.  I have a basic ESRI Shape file<br>
&gt;&gt; &gt;&gt; (created using gdal_contour) that I run through gdal_rasterize to<br>
&gt;&gt; &gt;&gt; generate very nice looking contour map.  But I want the value of the<br>
&gt;&gt; &gt;&gt; contour in there.<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; I&#39;ve been trying different things with the -a option, but I can&#39;t get<br>
&gt;&gt; &gt;&gt; it to work out.<br>
&gt;&gt; &gt;&gt; Here&#39;s the latest thing I&#39;ve tried for gdal_rasterize:<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; gdal_rasterize -a contour -3d -b 1 -b 2 -b 3 -b 4 -burn 255 -burn 255<br>
&gt;&gt; &gt;&gt; -burn 255 -burn 255 -burn 255 -l contour contour BLANK.tif<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; I put an extra burn value in there (have also tried with out) to see<br>
&gt;&gt; &gt;&gt; if maybe it was just going in as 0.  BLANK.tif is just a black image.<br>
&gt;&gt; &gt;&gt; I&#39;m sure I&#39;m missing something simple, just can&#39;t see it.  I thought<br>
&gt;&gt; &gt;&gt; the -a would burn the value of &quot;contour&quot; into the raster.<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; Thanks for any guidance!!<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; Here&#39;s the data from my shapefile:<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; INFO: Open of `contour/&#39;<br>
&gt;&gt; &gt;&gt;      using driver `ESRI Shapefile&#39; successful.<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; Layer name: contour<br>
&gt;&gt; &gt;&gt; Geometry: 3D Line String<br>
&gt;&gt; &gt;&gt; Feature Count: 157<br>
&gt;&gt; &gt;&gt; Extent: (-0.500000, -89.830132) - (359.500000, 89.144137)<br>
&gt;&gt; &gt;&gt; Layer SRS WKT:<br>
&gt;&gt; &gt;&gt; GEOGCS[&quot;GCS_WGS_1984&quot;,<br>
&gt;&gt; &gt;&gt;    DATUM[&quot;WGS_1984&quot;,<br>
&gt;&gt; &gt;&gt;        SPHEROID[&quot;WGS_1984&quot;,6378137,298.257223563]],<br>
&gt;&gt; &gt;&gt;    PRIMEM[&quot;Greenwich&quot;,0],<br>
&gt;&gt; &gt;&gt;    UNIT[&quot;Degree&quot;,0.017453292519943295]]<br>
&gt;&gt; &gt;&gt; ID: Integer (8.0)<br>
&gt;&gt; &gt;&gt; contour: Real (12.3)<br>
&gt;&gt; &gt;&gt; OGRFeature(contour):0<br>
&gt;&gt; &gt;&gt;  ID (Integer) = 0<br>
&gt;&gt; &gt;&gt;  contour (Real) =     1000.000<br>
&gt;&gt; &gt;&gt;  LINESTRING (87.765060240963862 80.0 1000,.....<br>
&gt;&gt; &gt;&gt; _______________________________________________<br>
&gt;&gt; &gt;&gt; gdal-dev mailing list<br>
&gt;&gt; &gt;&gt; <a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
&gt;&gt; &gt;&gt; <a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; --<br>
&gt;&gt; &gt; Best regards,<br>
&gt;&gt; &gt; Chaitanya kumar CH.<br>
&gt;&gt; &gt; /tʃaɪθənjə/ /kʊmɑr/<br>
&gt;&gt; &gt; <a href="tel:%2B91-9494447584" value="+919494447584">+91-9494447584</a><br>
&gt;&gt; &gt; 17.2416N 80.1426E<br>
&gt;&gt; &gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; --<br>
&gt;&gt; This communication is privileged and may contain confidential<br>
&gt;&gt; information.  It&#39;s intended only for the use of the person or entity<br>
&gt;&gt; named above. If you are not the intended recipient, do not distribute<br>
&gt;&gt; or copy this communication. If you have received this communication in<br>
&gt;&gt; error, please notify the sender immediately and return the original to<br>
&gt;&gt; the email address above. © Copyright 2011 Weather Trends<br>
&gt;&gt; International, Inc.<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; Best regards,<br>
&gt; Chaitanya kumar CH.<br>
&gt; /tʃaɪθənjə/ /kʊmɑr/<br>
&gt; <a href="tel:%2B91-9494447584" value="+919494447584">+91-9494447584</a><br>
&gt; 17.2416N 80.1426E<br>
&gt;<br>
<br>
<br>
<br>
</div></div>--<br>
<div class="im">This communication is privileged and may contain confidential<br>
information.  It&#39;s intended only for the use of the person or entity<br>
named above. If you are not the intended recipient, do not distribute<br>
or copy this communication. If you have received this communication in<br>
error, please notify the sender immediately and return the original to<br>
the email address above. © Copyright 2011 Weather Trends<br>
International, Inc.<br>
</div><div><div></div><div class="h5">_______________________________________________<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><br>
</div></div></blockquote></div><br>