Преобразование положения пикселя в географические координаты (широта, долгота) для изображения SAR Sentinel-1

Question

Как я могу получить географические координаты из положения x,y на снимке со спутника Sentinel-1 Synthetic Aperture Radar (SAR)?

Например, я могу получить доступ к информации загруженного изображения sg как

from snappy import ProductIO
from snappy import PixelPos 
path='path_name'
product = ProductIO.readProduct(path)
sg = product.getSceneGeoCoding()

Но как я могу получить широту и долготу для желаемой точки (х, у) в sg используя механизм привязки ESA в Python?

Answer 1

Используя пользовательскую функцию ниже, мы можем легко конвертировать любую точку внутри нашего изображения sg к его координатам (широта, долгота):

def LatLon_from_XY(ProductSceneGeoCoding, x, y):
    #From x,y position in satellite image (SAR), get the Latitude and Longitude
    geopos = ProductSceneGeoCoding.getGeoPos(PixelPos(x, y), None)
    latitude = geopos.getLat()
    longitude = geopos.getLon()
    return latitude, longitude

Например, для данного sg произведение, мы можем получить координаты пикселя (х =12000, у =2000) как

latitude, longitude = LatLon_from_XY(sg, 12000, 2000)