<div>hi</div><div>I've a method for producing map overlay with Grads in mercator projection, but is not so straightforward</div><div><br></div><div>i'm starting from rotated lat-lon or lambert conformal preprojected data in grib files</div>
<div>1) remove pdef directives in ctl files and sobstitute the "linear" progression in xdef and ydef with "levels" directive with the list of latitudes (ydef) and longitudes (xdef) of any single grid nodes</div>
<div>This steps in my case is producing a huge .ctl files (the grib ctl files describing a 1600x2000 grid); a perl text-manipulation script is doing this work</div><div>2) plot a tiff image with "set gxout kml", </div>
<div>3) converting and rescaling with imagemagick the *tiff image in PNG</div><div><br></div>
<div><br></div><div><br><br><div class="gmail_quote">2012/6/5 Jonas Kaufmann <span dir="ltr"><<a href="mailto:me@j0nes.de" target="_blank">me@j0nes.de</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello,<br>
<br>
do you have an idea how to use 'set gxout geotiff' to generate<br>
displayable PNG images for use on a map software? Or do you have any<br>
other hints how Mercator projection in GrADS works?<br>
<div><br>
Thank you for your help!<br>
<br>
Jonas<br>
<br>
</div><div><div>On Fri, Jun 1, 2012 at 10:55 AM, Jonas Kaufmann <<a href="mailto:me@j0nes.de" target="_blank">me@j0nes.de</a>> wrote:<br>
> Hello,<br>
><br>
> thank you for your answer! Too bad there is no support for Mercator<br>
> projection available (yet).<br>
><br>
> I just had a look at the GeoTIFF output, but I did not succeed in<br>
> creating a PNG (or whatever) file that I could overlay on OSM or<br>
> Google Maps. The tif file gets created, but it is looking almost<br>
> completely black when viewing with an image viewer (but I think this<br>
> is expected behaviour?). I have tried to convert it to PNG using the<br>
> GDAL tools with gdal_translate, but I only get grayscale noise images.<br>
><br>
> Can someone explain the correct workflow to go from 'set gxout<br>
> geotiff' to a displayable PNG image?<br>
><br>
> Thank you for your help!<br>
><br>
><br>
> Jonas<br>
><br>
><br>
> On Thu, May 31, 2012 at 5:05 PM, Chris Slocum<br>
> <<a href="mailto:cslocum@atmos.colostate.edu" target="_blank">cslocum@atmos.colostate.edu</a>> wrote:<br>
>> Hi Jonas,<br>
>><br>
>> I too have been looking for a Mercator projection to<br>
>> use with something like Open Street Map. To date,<br>
>> I do not think there is a public extension for a Mercator<br>
>> projection but the idea has been bounced around.<br>
>><br>
>> I know that Weather Underground had someone<br>
>> develop the projection for their WunderMap so it<br>
>> is possible to do.<br>
>><br>
>> In the interim, you could look at the GeoTIFF option.<br>
>><br>
>> -Chris<br>
>><br>
>><br>
>><br>
>> On Thu, May 31, 2012 at 6:46 AM, Jonas Kaufmann <<a href="mailto:me@j0nes.de" target="_blank">me@j0nes.de</a>> wrote:<br>
>>><br>
>>> Hello,<br>
>>><br>
>>> I am trying to generate an image with GrADS that should cover the<br>
>>> whole world. The image should be based on an forecast from the GFS<br>
>>> model. I want to use this image as a layer on top of Open Street Map.<br>
>>><br>
>>> I can already create an image without borders, labels etc that I can<br>
>>> put on my map. However it seems that the projection still is wrong,<br>
>>> e.g. the continent borders drawn by GrADS do not match the borders on<br>
>>> the map. I think this might be related to the map projection used, in<br>
>>> OSM this should be Mercator projection which is not available in<br>
>>> GrADS. The complete script can be found at the bottom of the mail.<br>
>>><br>
>>> Is there a way to make the forecast map fit on OSM?<br>
>>><br>
>>><br>
>>> Thank you for your help!<br>
>>><br>
>>><br>
>>> Jonas<br>
>>><br>
>>><br>
>>> -------<br>
>>><br>
>>> 'open gfs'<br>
>>><br>
>>> 'set gxout shaded'<br>
>>> 'set parea 0 11 0 8.5'<br>
>>> 'set mpdraw off'<br>
>>> 'set frame off'<br>
>>> 'set xlab off'<br>
>>> 'set ylab off'<br>
>>> 'set mproj latlon'<br>
>>> 'set grads off'<br>
>>><br>
>>> 'set mpdset hires'<br>
>>> 'set map 1 1 1'<br>
>>><br>
>>> 'set lon -180 180'<br>
>>> 'set lat -90 90'<br>
>>><br>
>>> 'd mag(ugrd10m, vgrd10m)'<br>
>>><br>
>>> 'printim world.png x512 y512'<br>
>>> 'quit'<br>
>>> _______________________________________________<br>
>>> gradsusr mailing list<br>
>>> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
>>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>><br>
>><br>
>><br>
>> _______________________________________________<br>
>> gradsusr mailing list<br>
>> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>><br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>------------------------------------------------------------------------<br><font size="1"><b><span style="font-family:tahoma,sans-serif">Bertozzi Angelo</span></b><br style="font-family:tahoma,sans-serif">
<span style="font-family:tahoma,sans-serif">Meteorologist - HPCC admin & research</span><br style="font-family:tahoma,sans-serif"><span style="font-family:tahoma,sans-serif">Atmospheric and Environmental numerical modeling division - head</span><br style="font-family:tahoma,sans-serif">
<span style="font-family:tahoma,sans-serif">Mobile <a href="tel:345-71.335.07" value="+393457133507" target="_blank">345-71.335.07</a> - Office <a href="tel:0546-046207" value="+390546046207" target="_blank">0546-046207</a></span><br style="font-family:tahoma,sans-serif">
<b style="font-family:tahoma,sans-serif">METEOCENTER (Sms Srl)</b></font><br style="font-family:tahoma,sans-serif">
<font size="1"><span style="font-family:tahoma,sans-serif">Parco delle Arti e delle Scienze "TORRICELLI" - Via Albert Einstein, 8 - 48018 FAENZA (RA) - ITALY</span></font><br><br>-------------------------------------------------------------------------<br>
<font size="1">Le informazioni contenute in questa e-mail e negli eventuali allegati sono riservate e destinate esclusivamente alla persona sopraindicata. Qualora non foste il destinatario, siete pregati di distruggere questo messaggio e notificarci il problema immediatamente. In ogni caso, non dovrete spedire a terzi, copiare, usare o diffondere il contenuto di questa e-mail e degli eventuali allegati. Si ricorda che la diffusione l'utilizzo e/o la conservazione dei dati ricevuti per errore costituiscono violazione alle disposizioni del D.lgs. n. 196/2003 (Codice in materia di protezione dei dati personali) oltre a costituire violazione di carattere penale ai sensi dell'art. 616 Codice penale (Violazione, sottrazione e soppressione di corrispondenza).</font><br>
</div>