<font class="Apple-style-span" face="arial, helvetica, sans-serif">thanks for that... </font><div><font class="Apple-style-span" face="arial, helvetica, sans-serif"><br></font></div><div><font class="Apple-style-span" face="arial, helvetica, sans-serif">What I&#39;m looking for is something like <span class="Apple-style-span" style="white-space: pre-wrap; ">ReadRaster but using the coordinates system used by the raster... (it could be lat,long etc...)</span></font><div>

<font class="Apple-style-span" face="arial, helvetica, sans-serif"><span class="Apple-style-span" style="white-space: pre-wrap; "><br></span></font></div><div><font class="Apple-style-span" face="arial, helvetica, sans-serif"><span class="Apple-style-span" style="white-space: pre-wrap; ">otherwise I would need to do something like</span></font></div>

<div><font class="Apple-style-span" face="arial, helvetica, sans-serif"><span class="Apple-style-span" style="white-space: pre-wrap; "><br></span></font></div><div><font class="Apple-style-span" face="arial, helvetica, sans-serif"><span class="Apple-style-span" style="white-space: pre-wrap; ">cs=read_raster_coordinate_system()</span></font></div>

<div><font class="Apple-style-span" face="arial, helvetica, sans-serif"><span class="Apple-style-span" style="white-space: pre-wrap; ">(x,y)= cs.convert_from_cs_to_pixel_xy(x,y)</span></font></div><div><font class="Apple-style-span" face="arial, helvetica, sans-serif"><span class="Apple-style-span" style="white-space: pre-wrap;">pixel=ReadRaster </span></font></div>

<div><font class="Apple-style-span" face="arial, helvetica, sans-serif"><span class="Apple-style-span" style="white-space: pre-wrap;"><br></span></font></div><div><font class="Apple-style-span" face="arial, helvetica, sans-serif"><span class="Apple-style-span" style="white-space: pre-wrap;">is there any single line command to do that ?<br>

</span></font><br><div class="gmail_quote">2011/1/12 Jorge Arévalo <span dir="ltr">&lt;<a href="mailto:jorge.arevalo@deimos-space.com">jorge.arevalo@deimos-space.com</a>&gt;</span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">

<div><div></div><div class="h5">On Wed, Jan 12, 2011 at 4:57 PM, Sebastian E. Ovide<br>
&lt;<a href="mailto:sebastian.ovide@gmail.com">sebastian.ovide@gmail.com</a>&gt; wrote:<br>
&gt; Hi All,<br>
&gt; what would be the simpler way to read a pixel value from a raster using<br>
&gt; python ?<br>
&gt; something like  pixelvalue=read(x,y) where x,y and in the raster coordinate<br>
&gt; system<br>
&gt; thanks<br>
&gt;<br>
&gt; --<br>
&gt; Sebastian E. Ovide<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
</div></div>&gt; _______________________________________________<br>
&gt; gdal-dev mailing list<br>
&gt; <a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
&gt; <a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
&gt;<br>
<br>
Hi Sebastian,<br>
<br>
This may help you:<br>
<a href="https://svn.osgeo.org/postgis/trunk/raster/scripts/python/pixval.py" target="_blank">https://svn.osgeo.org/postgis/trunk/raster/scripts/python/pixval.py</a><br>
<br>
Best regards,<br>
<font color="#888888"><br>
--<br>
Jorge Arévalo<br>
Internet &amp; Mobilty Division, DEIMOS<br>
<a href="mailto:jorge.arevalo@deimos-space.com">jorge.arevalo@deimos-space.com</a><br>
<a href="http://mobility.grupodeimos.com/" target="_blank">http://mobility.grupodeimos.com/</a><br>
<a href="http://gis4free.wordpress.com" target="_blank">http://gis4free.wordpress.com</a><br>
<a href="http://geohash.org/ezjqgrgzz0g" target="_blank">http://geohash.org/ezjqgrgzz0g</a><br>
</font></blockquote></div><br><br clear="all"><br>-- <br>Sebastian E. Ovide<br><br><br><br><br>
</div></div>