osr.TransformPoint переворачивает x и y, если предоставлен файл в формате ENVI RAW
Я хочу преобразовать угловые координаты проецируемого файла UTM в longlat. Если я читаю координаты из GeoTIFF, функция работает, как и ожидалось. Но после преобразования необработанного файла GeoTiff в ENVI (с помощью gdal_translate) координаты меняются местами:
from osgeo import gdal
from osgeo import osr
def corner2longlat(fname):
dataset = gdal.Open(fname, gdal.GA_ReadOnly)
ds_geotrans = dataset.GetGeoTransform()
X = ds_geotrans[0]
Y = ds_geotrans[3] + dataset.RasterYSize * ds_geotrans[5]
srs_projection = dataset.GetProjectionRef()
srs = osr.SpatialReference()
srs.ImportFromWkt(srs_projection)
srsLatLong = srs.CloneGeogCS()
ct = osr.CoordinateTransformation(srs, srsLatLong)
latlon = ct.TransformPoint(X, Y)
print(srs_projection)
print([latlon[0], latlon[1]])
image_tif = "all_bands.tif"
image_raw = "all_bands.raw"
PROJCS["WGS 84 / UTM zone 43N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY[" EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["степень",0,0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY ["EPSG","4326"]],ПРОЕКЦИЯ["Поперек_Меркатора"],ПАРАМЕТР["широта_начала",0],ПАРАМЕТР["центральный_меридиан",75],ПАРАМЕТР["масштабный_коэффициент",0,9996],ПАРАМЕТР["false_easting" ,500000],PARAMETER["false_northing",0],UNIT["meter",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Neverthing",NORTH] ,ОРГАН ВЛАСТИ["EPSG","32643"]]
[41.463751886708124, 74.99976050569691]
PROJCS["WGS 84 / UTM zone 43N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY[" EPSG","6326"]],PRIMEM["Гринвич",0],UNIT["Градус",0,0174532925199433]],PROJECTION["Поперек_Меркатора"],PARAMETER["широта_начала",0],PARAMETER["central_meridian", 75],ПАРАМЕТР["масштабный_коэффициент",0,9996],ПАРАМЕТР["false_easting",500000],ПАРАМЕТР["false_northing",0],ЕДИНИЦА["метр",1,АВТОРИЗАЦИЯ["EPSG","9001"]], ОСЬ["Восток",ВОСТОК],ОСЬ["Север",СЕВЕР]]
[74.99976050569691, 41.463751886708124]
Кто-нибудь знает, как это могло произойти?
заранее спасибо
Лукас