<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 fpstyle="1" ocsi="0" style="">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Oops, sorry! I forgot to include the -a_ullr variables. I have:
<div>
<div><span style="font-size: 10pt;">up_left_x</span>=-5570248.832537</div>
<div><span style="font-size: 10pt;">up_left_y</span>=5570244.832537</div>
<div><span style="font-size: 10pt;">low_rig_x</span>=5567248.429179</div>
<div><span style="font-size: 10pt;">low_rig_y</span>=-5567252.429179</div>
<div><br>
</div>
<div>That should cover the full disk and seems to work with all satellites - you just need to vary the longitude to change between GOES/MSG/MTSAT etc.</div>
<div>Cheers,</div>
<div>Simon</div>
<div><br>
</div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF419720" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>Fra:</b> Jerl Simpson [jsimpson@wxtrends.com]<br>
<b>Sendt:</b> 28. oktober 2013 16:41<br>
<b>Til:</b> Simon Richard Proud<br>
<b>Cc:</b> gdal-dev@lists.osgeo.org<br>
<b>Emne:</b> Re: [gdal-dev] Satellite projections<br>
</font><br>
</div>
<div></div>
<div>
<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 face="Tahoma" color="#000000"><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>
<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><u></u><u></u></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><u></u><u></u></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><u></u><u></u></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)"><u></u> </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>
</div>
</div>
</div>
</div>
<br clear=all>
</body>
</html>