[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[python #IRW-840489]: Fulldisk lon/lat -> x/y data point conversion



Hi,

The conversion is described in this document:

https://www.goes-r.gov/users/docs/PUG-L1b-vol3.pdf

You'll want to look at section 5.1.2.8, "Navigation of Image Data". If you're 
doing this in Python, you may just want to rely on the calculations in PyProj 
(using the "geos" transform https://proj.org/operations/projections/geos.html) 
or CartoPy (using the Geostationary projection class 
https://scitools.org.uk/cartopy/docs/latest/crs/projections.html#geostationary).
 The one caveat for either of these two approaches is that you'll need to 
multiply your x,y values by the satellite height to get the proper coordinates 
to pass into those projection classes.

Hope that helps,

Ryan

> Hello-
> 
> I am working on a python script to display satellite data based on a central 
> longitude and latitude. However, the data grid for the full disk satellite is 
> not uniform, and my best guess still leaves me very inaccurate.
> 
> Can you please provide the conversation between lon/lat and x/y? To be 
> specific, I am working with imagery that can be found here.
> 
> https://thredds.ucar.edu/thredds/catalog/satellite/goes/east/products/CloudAndMoistureImagery/FullDisk/catalog.xml
> 
> I hope this makes sense. Please let me know if you need more information from 
> me.


Ticket Details
===================
Ticket ID: IRW-840489
Department: Support Python
Priority: Low
Status: Closed
===================
NOTE: All email exchanges with Unidata User Support are recorded in the Unidata 
inquiry tracking system and then made publicly available through the web.  If 
you do not want to have your interactions made available in this way, you must 
let us know in each email you send to us.