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 [[email protected]] Sendt: 28. oktober 2013 16:55 Til: Jerl Simpson Cc: Simon Richard Proud; [email protected] 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 <[email protected]<mailto:[email protected]>> 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: [email protected]<mailto:[email protected]> [[email protected]<mailto:[email protected]>] på vegne af Jerl Simpson [[email protected]<mailto:[email protected]>] Sendt: 28. oktober 2013 16:17 Til: [email protected]<mailto:[email protected]> 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 [email protected]<mailto:[email protected]> 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 [email protected]<mailto:[email protected]> 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
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
