Snappy - Неправильное преобразование из пикселя в лат длинный - PullRequest
Новая
       8

Snappy - Неправильное преобразование из пикселя в лат длинный

0 голосов
/ 12 сентября 2018

Я использую snappy, чтобы попытаться выяснить, каковы координаты x, y, чтобы обрезать изображение.

Я провел несколько тестов с использованием функций snappy, но заметил, что тамс ними что-то не так.Я преобразовал X, Y из изображения в широту и долготу, а затем, используя эти координаты, попытался снова преобразовать их в X, Y, но я не получил тот же результат.

Идеяили конечная цель - получить координаты LatLong у геоджона, прочитать их и, используя эти координаты, получить X, Y на изображении.

Обратите внимание, что путь относится к файлу TIFF.

from snappy import ProductIO
from snappy import PixelPos, GeoPos
import numpy as np

path='/home/.../x.tiff'
###############################################################################
product = ProductIO.readProduct(path)
sg = product.getSceneGeoCoding()

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

latitude, longitude = LatLon_from_XY(sg, 11048, 1365)

print('LatLong from PixelPosition')
print(latitude)
print(longitude)
### 38.3976151718
### -5.47978868123

###############################################################################

def getPixelPosFromLatLong(source, lat,lon):
    if sg.canGetPixelPos() is not True:
        raise Exception('Cant''t get Pixel Position from this source')
    else:
        pos = GeoPos(lat,lon)
        pixpos = sg.getPixelPos(pos,None)
        X = np.round(pixpos.getX())
        Y = np.round(pixpos.getY())
    return [X,Y]

[X,Y] = getPixelPosFromLatLong(path, 38.3976151718, -5.47978868123)

print('Pixel Position from LatLong')
print(X)
print(Y)
### 10715.0
### 1143.0

Есть ли другой способ получить пиксель X, Y из изображения, используя lat long?

1 Ответ

0 голосов
/ 13 сентября 2018

Для тех из вас, кто интересуется той же проблемой, я нашел решение:

Изменить путь к ссылке не на файл .tiff, а на папку .SAFE

Затем я снова пишу свои функции, и теперь это выглядит так: 

from snappy import ProductIO
from snappy import PixelPos, GeoPos
import numpy as np

def LatLon_from_XY(ProductSceneGeoCoding, x, y):
    geoPos = ProductSceneGeoCoding.getGeoPos(PixelPos(x,y),None)
    lat = geoPos.getLat()
    lon = geoPos.getLon()
    return lat,lon

def XY_from_LatLon(ProductSceneGeoCoding, latitude, longitude):
    pixelPos = ProductSceneGeoCoding.getPixelPos(GeoPos(latitude, longitude),None)
    x = np.round(pixelPos.getX())
    y = np.round(pixelPos.getY())
    return x,y

###############################################################################

#Notice that the path does not refer to any file but to the folder .SAFE
path = '/.../X.SAFE'
product = ProductIO.readProduct(path)
sg = product.getSceneGeoCoding()
originalX = 13000
originalY = 13000

print('Original X,Y: ', originalX,originalY)

lat,lon = LatLon_from_XY(sg, originalX, originalY)
print(lat,lon)

x,y = XY_from_LatLon(sg,lat,lon)
print(x,y)

originalLat = 37.36475504265766
originalLon = -5.972416873450527

print('Original Lat,Lon: ', originalLat, originalLon)

x,y = XY_from_LatLon(sg, originalLat, originalLon)
print(x,y)

lat,lon = LatLon_from_XY(sg, x, y)
print(lat,lon)

Таким образом, я получаю почти одинаковые значения (13000 ~ 13006)

...