<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style type="text/css" id="owaParaStyle"></style><style type="text/css"></style><style type="text/css"></style>
</head>
<body bgcolor="#FFFFFF" fpstyle="1" ocsi="0" style="">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi Reinoud,
<div>To be honest I've never had success in using GDAL for GOES in LRIT format (only for reprojecting when it is already in Geotiff). In that case I would recommend using PyTroll instead. It works well for all geostationary data: <a href="http://pytroll.org/" target="_blank" style="font-size: 10pt;">http://pytroll.org/</a></div>
<div><br>
</div>
<div>Cheers,</div>
<div>Simon</div>
<div><br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF300424" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>Fra:</b> Reinoud Bokhorst [reinoud.bokhorst@bmtargoss.com]<br>
<b>Sendt:</b> 28. oktober 2013 16:55<br>
<b>Til:</b> Jerl Simpson<br>
<b>Cc:</b> Simon Richard Proud; gdal-dev@lists.osgeo.org<br>
<b>Emne:</b> Re: [gdal-dev] Satellite projections<br>
</font><br>
</div>
<div></div>
<div>
<div class="moz-cite-prefix">Dear Jerl, Simon,<br>
Apologize for mixing into this thread but I am currently investigating how to process EUMETCAST data to more usuable formats like NetCDF and GeoTIFF.  The data stream contains MSG HRIT data and GOES15/17 + MTSAT2 LRIT data. I'll also need to reproject like
 you do so this post already proves useful!<br>
<br>
I can read the MSG data trough gdal/meteosatlib but I can't read the GOES and MTSAT2 LRIT data yet. Do you perhaps have some useful pointers to software to use? Online documentation is rather scarce...<br>
<br>
Thanks Reinoud<br>
<br>
<br>
<br>
On 10/28/2013 04:41 PM, Jerl Simpson wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">Thanks Simon!
<div><br>
</div>
<div>I'll give it a shot.  Meteosat is definitely a pain, and it's one of the satellites I'm working with.</div>
<div><br>
</div>
<div>If you don't mind me asking, what do you typically have in your -a_ullr designation for gdal_translate?  Or are you pulling out sections and it varies?  I'd like to get the full image, if possible.</div>
<div><br>
</div>
<div>Thanks again!</div>
<div><br>
</div>
<div>Jerl</div>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Mon, Oct 28, 2013 at 10:33 AM, Simon Richard Proud <span dir="ltr">
<<a href="mailto:srp@geo.ku.dk" target="_blank">srp@geo.ku.dk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
            .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div>
<div style="direction:ltr; font-size:10pt; font-family:Tahoma">Hi Jerl,
<div>We work with Meteosat data in full disk projection and it's a bit of a headache! Try these commands in GDAL and see if they work for you:</div>
<div>
<div><br>
</div>
<div>gdal_translate -a_srs  "+proj=geos +a=6378169 +b=6356583.8 +lon_0=9.49 +h=35785831" -a_ullr $up_left_x $up_left_y $low_rig_x $low_rig_y $INFILE $TEMPFILE</div>
<div>gdalwarp -s_srs '+proj=geos +lon_0=9.49 +h=35785831 +x_0=0.0' -t_srs '+proj=latlong +datum=WGS84' -tr 0.01 0.01 -te 0 20  70 70 -order 3 <span style="font-size:10pt">$TEMPFILE $OUTFILE</span></div>
<div><span style="font-size:10pt"><br>
</span></div>
<div><span style="font-size:10pt">Just remember to change the +lon_0 to whatever you need for GOES. Also change the -tr and -te for gdal :-) You can also add +lat_0 but to be honest I never found it made a difference (perhaps it is ignored?)</span></div>
<div><span style="font-size:10pt"><br>
</span></div>
<div><span style="font-size:10pt">Cheers,</span></div>
<div><span style="font-size:10pt">Simon</span></div>
<div><br>
</div>
<div style="font-size:16px; font-family:Times New
                    Roman">
<hr>
<div style="direction:ltr"><font color="#000000" face="Tahoma"><b>Fra:</b> <a href="mailto:gdal-dev-bounces@lists.osgeo.org" target="_blank">
gdal-dev-bounces@lists.osgeo.org</a> [<a href="mailto:gdal-dev-bounces@lists.osgeo.org" target="_blank">gdal-dev-bounces@lists.osgeo.org</a>] på vegne af Jerl Simpson [<a href="mailto:jsimpson@wxtrends.com" target="_blank">jsimpson@wxtrends.com</a>]<br>
<b>Sendt:</b> 28. oktober 2013 16:17<br>
<b>Til:</b> <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<b>Emne:</b> [gdal-dev] Satellite projections<br>
</font><br>
</div>
<div>
<div class="h5">
<div>
<div dir="ltr">Hi guys:<br>
<div><br>
</div>
<div>I have a few data sets that come from the GOES satellites (Full Disk) and a couple others.  I can use one of the SSEC's programs to convert from their format to a netCDF.  No useful SRS info comes along in the conversion.</div>
<div><br>
</div>
<div>I've applied an SRS to them using gdal_translate and PROJ.4 that looks like:<br>
<br>
"+proj=ortho +lon_0=-79.8 +lat_0=.208"  (lon/lat is the position of one of the satellites)<br>
</div>
<div><br>
</div>
<div>(I've also tried):</div>
<div>+proj=ortho +datum=WGS84 +ellps=WGS84 +lon_0=-79.8 +lat_0=.208<br>
</div>
<div><br>
</div>
<div>I'm wanting to get them into EPSG:4326.  When I do this, I end up with corner coordinates of a very small area just below the satellites position.  </div>
<div><br>
</div>
<div>Does anyone have an suggestions that might help me warp this into position?  I'm sure there's something obvious that I'm missing.</div>
<div>I found this page, and the copied/pasted bit about the projection, but my attempts have not been successful.</div>
<div><br>
</div>
<div><a href="http://www.gdal.org/frmt_msg.html" target="_blank">http://www.gdal.org/frmt_msg.html</a><br>
</div>
<div>
<div><i>Georeference and Projection</i></div>
<div><i><br>
</i></div>
<div><i>The images are using the Geostationary Satellite View projection. Most GIS packages don't recognize this projection (we only know of ILWIS that does have this projection), but gdalwarp.exe can be used to re-project the images.</i></div>
</div>
<div><i><br>
</i></div>
<div>Thanks!</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br clear="all">
</div>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div>
<p style="margin:0px; font-family:arial,sans-serif; font-size:13px; background-color:rgb(255,255,255)">
<font color="#1f497d">Jerl Simpson</font></p>
<p style="margin:0px; font-family:arial,sans-serif; font-size:13px; background-color:rgb(255,255,255)">
<font color="#1f497d">Sr. Systems Engineer</font></p>
<p style="margin:0px; color:rgb(34,34,34); font-family:arial,sans-serif; font-size:13px; background-color:rgb(255,255,255)">
<span style="color:rgb(31,73,125)">Weather Trends International</span></p>
<p style="margin:0px; color:rgb(34,34,34); font-family:arial,sans-serif; font-size:13px; background-color:rgb(255,255,255)">
<span style="color:rgb(31,73,125)"><a href="mailto:jsimpson@wxtrends.com" target="_blank">jsimpson@wxtrends.com</a></span></p>
<p style="margin:0px; color:rgb(34,34,34); font-family:arial,sans-serif; font-size:13px; background-color:rgb(255,255,255)">
<span style="color:rgb(31,73,125)">Corporate – <a href="http://www.wxtrends.com/" style="color:rgb(17,85,204)" target="_blank">http://www.wxtrends.com</a></span></p>
<p style="margin:0px; color:rgb(34,34,34); font-family:arial,sans-serif; font-size:13px; background-color:rgb(255,255,255)">
<span style="color:rgb(31,73,125)">B2B – <a href="http://www.wt360business.com/" style="color:rgb(17,85,204)" target="_blank">http://www.wt360business.com</a></span></p>
<p style="margin:0px; color:rgb(34,34,34); font-family:arial,sans-serif; font-size:13px; background-color:rgb(255,255,255)">
<span style="color:rgb(31,73,125)">B2C – <a href="http://www.wt360.com/" style="color:rgb(17,85,204)" target="_blank">http://www.wt360.com</a></span></p>
<p style="margin:0px; color:rgb(34,34,34); font-family:arial,sans-serif; font-size:13px; background-color:rgb(255,255,255)">
<span style="color:rgb(31,73,125)"> </span></p>
</div>
<div><br>
</div>
<div><font color="#ff0000">This communication is privileged and may contain confidential information. It's intended only for the use of the person or entity named above. If you are not the intended recipient, do not distribute or copy this communication. If
 you have received this communication in error, please notify the sender immediately and return the original to the email address above. © Copyright 2013 Weather Trends International, Inc.</font></div>
</div>
<br clear="all">
______________________________________________________________________<br>
This email has been scanned by the Symantec Email Security.cloud service.<br>
For more information please visit <a class="moz-txt-link-freetext" href="http://www.symanteccloud.com" target="_blank">
http://www.symanteccloud.com</a><br>
______________________________________________________________________<br>
<br>
<fieldset class="mimeAttachmentHeader" target="_blank"></fieldset> <br>
<pre>_______________________________________________
gdal-dev mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a></pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">-- 
Reinoud Bokhorst
Advisor Oceanography & ICT
BMT ARGOSS

Tel.:         +31 527-242299
Direct:     +31 527-656421
Mobile:    +31 6-54714427

<a class="moz-txt-link-abbreviated" href="http://www.bmt.org/linkedin" target="_blank">www.bmt.org/linkedin</a> | <a class="moz-txt-link-abbreviated" href="http://www.twitter.com/bmtgroup" target="_blank">www.twitter.com/bmtgroup</a> | <a class="moz-txt-link-abbreviated" href="http://www.facebook.com/bmtgroup" target="_blank">www.facebook.com/bmtgroup</a> | <a class="moz-txt-link-abbreviated" href="http://www.youtube.com/bmtgroupltd" target="_blank">www.youtube.com/bmtgroupltd</a>

BMT ARGOSS b.v., Voorsterweg 28, 8316 PT, Marknesse, The Netherlands

(Postal address : PO Box 61, 8325 ZH, Vollenhove, The Netherlands)

Registered in The Netherlands, Registered no. 39060160.

Unless otherwise agreed by BMT ARGOSS in writing, all  work, services, goods or products supplied by BMT ARGOSS shall be subject to and governed by BMT ARGOSS’s own terms and conditions which are available for inspection from BMT ARGOSS on request.

E-mail confidentiality notice and disclaimer:

The contents of this e-mail and any attachments are intended for the use of the mail addressee(s) shown. If you are not that person, you are not allowed to read it, to take any action based upon it or to copy it, forward, distribute or disclose the contents of it and you should please delete it from your system. BMT ARGOSS does not accept liability for any errors or omissions in the context of this e-mail or its attachments which arise as a result of internet transmission, nor accept liability for statements which are those of the author and clearly not made on behalf of BMT ARGOSS.

Please consider the environmental impacts of printing this email, and only do so if really necessary</pre>
</div>
</div>
</div>
</div>
<br clear=all> 
</body>
</html>