[gdal-dev] Satellite projections

Tue Oct 29 06:50:00 PDT 2013

Thanks for both your tips, I'll give it a try.

Reinoud

________________________________
From: Simon Richard Proud [srp at geo.ku.dk]
Sent: Monday, October 28, 2013 5:05 PM
To: Reinoud Bokhorst; Jerl Simpson
Cc: gdal-dev at lists.osgeo.org
Subject: SV: [gdal-dev] Satellite projections

Hi Reinoud,
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: http://pytroll.org/

Cheers,
Simon

________________________________
Fra: Reinoud Bokhorst [reinoud.bokhorst at bmtargoss.com]
Sendt: 28. oktober 2013 16:55
Til: Jerl Simpson
Cc: Simon Richard Proud; gdal-dev at lists.osgeo.org
Emne: Re: [gdal-dev] Satellite projections

Dear Jerl, Simon,
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!

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...

Thanks Reinoud

On 10/28/2013 04:41 PM, Jerl Simpson wrote:
Thanks Simon!

I'll give it a shot.  Meteosat is definitely a pain, and it's one of the satellites I'm working with.

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.

Thanks again!

Jerl

On Mon, Oct 28, 2013 at 10:33 AM, Simon Richard Proud <srp at geo.ku.dk<mailto:srp at geo.ku.dk>> wrote:
Hi Jerl,
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:

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
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 $TEMPFILE $OUTFILE

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?)

Cheers,
Simon

________________________________
Fra: gdal-dev-bounces at lists.osgeo.org<mailto:gdal-dev-bounces at lists.osgeo.org> [gdal-dev-bounces at lists.osgeo.org<mailto:gdal-dev-bounces at lists.osgeo.org>] på vegne af Jerl Simpson [jsimpson at wxtrends.com<mailto:jsimpson at wxtrends.com>]
Sendt: 28. oktober 2013 16:17
Til: gdal-dev at lists.osgeo.org<mailto:gdal-dev at lists.osgeo.org>
Emne: [gdal-dev] Satellite projections

Hi guys:

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.

I've applied an SRS to them using gdal_translate and PROJ.4 that looks like:

"+proj=ortho +lon_0=-79.8 +lat_0=.208"  (lon/lat is the position of one of the satellites)

(I've also tried):
+proj=ortho +datum=WGS84 +ellps=WGS84 +lon_0=-79.8 +lat_0=.208

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.

Does anyone have an suggestions that might help me warp this into position?  I'm sure there's something obvious that I'm missing.
I found this page, and the copied/pasted bit about the projection, but my attempts have not been successful.

http://www.gdal.org/frmt_msg.html
Georeference and Projection

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.

Thanks!

--

Jerl Simpson

Sr. Systems Engineer

Weather Trends International

jsimpson at wxtrends.com<mailto:jsimpson at wxtrends.com>

Corporate – http://www.wxtrends.com<http://www.wxtrends.com/>

B2B – http://www.wt360business.com<http://www.wt360business.com/>

B2C – http://www.wt360.com<http://www.wt360.com/>

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.

______________________________________________________________________
This email has been scanned by the Symantec Email Security.cloud service.
For more information please visit http://www.symanteccloud.com
______________________________________________________________________

_______________________________________________
gdal-dev mailing list
gdal-dev at lists.osgeo.org<mailto:gdal-dev at lists.osgeo.org>
http://lists.osgeo.org/mailman/listinfo/gdal-dev

--
Reinoud Bokhorst
Advisor Oceanography & ICT
BMT ARGOSS

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

www.bmt.org/linkedin<http://www.bmt.org/linkedin> | www.twitter.com/bmtgroup<http://www.twitter.com/bmtgroup> | www.facebook.com/bmtgroup<http://www.facebook.com/bmtgroup> | www.youtube.com/bmtgroupltd<http://www.youtube.com/bmtgroupltd>

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

______________________________________________________________________
This email has been scanned by the Symantec Email Security.cloud service.
For more information please visit http://www.symanteccloud.com
______________________________________________________________________

______________________________________________________________________
This email has been scanned by the Symantec Email Security.cloud service.
For more information please visit http://www.symanteccloud.com
______________________________________________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20131029/8ef6df04/attachment-0001.html>