Geek asks and answers

How to get latitude and longitude for an eumetsat image pixel using python?

I have an image of 3712x3712 pixels in size of the geostationary satellite eumetsat. There is black on the ground, so the image looks like this:

Example image that is not an eumetsat image

For each pixel of the earth I would like to get its latitude and longitude. I know that pyproj and I were able to create a projection like this:

sat = pyproj.Proj('+proj=geos +lon_0 +h=035785831.0 +x_0=0 +y_0=0')

but getting a latlon pixel (using sat(x,y,inverse=True), where xand yare the pixel coordinates on the image) is obviously impossible, since the projection does not know the dimension (3712x3712) of my image.

What am I missing?

+5
source share
1 answer

I think you are using the correct projection library and settings.

( ) - eumetsat 3 .

, lon/lat x, y (-81 , 81 - , . eumetsat , http://www.eumetsat.int/):

import pyproj
sat = pyproj.Proj('+proj=geos +lon_0 +h=035785831.0 +x_0=0 +y_0=0')
x,y =  sat( 81.299, 0, radians = False, errcheck = True)
print (x * 2.0 / 3712.0 ) / 1000.0

2.927, , eumetsat.

, (, ), x/y lat/lon . , .

, , . ). .

+2

More articles:


All Articles