On Wed, Dec 8, 2010 at 2:31 PM, Tatiana.Gonzalez <span dir="ltr">&lt;<a href="mailto:Tatiana.Gonzalez@noaa.gov">Tatiana.Gonzalez@noaa.gov</a>&gt;</span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">






<div bgcolor="#ffffff" text="#000000">







<p class="MsoNormal"><small>Hi,<br>
Is there a way to project grads in native satellite image? I
have this image (see <span lang="EN-GB">10Z12DEC2008_wind.png)
but
what I really need is something like this (see sat_image.png)</span></small></p>
<p class="MsoNormal"><small><span lang="EN-GB"> </span></small></p>
<p class="MsoNormal"><small>Can anybody give me some ideas on how to
solve this problem?<br></small></p></div></blockquote><div><br></div><div>I&#39;d done this sort of thing in pygrads, and there are some examples plotting fields on a geostationary image. More information here:</div><div>
<br></div><div>                  <a href="http://opengrads.org/wiki/index.php?title=Python_Interface_to_GrADS">http://opengrads.org/wiki/index.php?title=Python_Interface_to_GrADS</a></div><div><br></div><div>This requires some python. For example, to plot contours on a conus image:  </div>
<div><br></div><div><span class="Apple-style-span" style="font-family: sans-serif; font-size: 13px; line-height: 19px; "><pre style="padding-top: 1em; padding-right: 1em; padding-bottom: 1em; padding-left: 1em; border-top-width: 1px; border-right-width: 1px; border-bottom-width: 1px; border-left-width: 1px; border-top-style: dashed; border-right-style: dashed; border-bottom-style: dashed; border-left-style: dashed; border-top-color: rgb(47, 111, 171); border-right-color: rgb(47, 111, 171); border-bottom-color: rgb(47, 111, 171); border-left-color: rgb(47, 111, 171); color: black; background-color: rgb(249, 249, 249); line-height: 1.1em; ">
       sat = Image.open(&quot;goes_conus.jpg&quot;)
       opts=(-109.65,14.23,-44.65,62.23,-75.) # image info
       ga.basemap(&#39;geos&#39;, opts=opts, resolution=&#39;l&#39;) 
       ga.contour(&#39;h500&#39;,bgim=sat,N=12,Map=False)
       title(&quot;500 hPa Heights&quot;)</pre></span></div><div><br></div><div>The resulting plot is here: <a href="http://sourceforge.net/dbimage.php?id=153948">dbimage.php</a></div><div><br></div><div>    Arlindo</div>
</div><br>-- <br>Arlindo da Silva<br><a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a><br>