Не работает: космоснимки GeoTIFF и метафайл (QGIS)

[gimran]

Здравствуйте!
Столкнулся со следующей проблемой: скачал файлы Landsat5-TM.
Решил объединить поканально:
Растр — Объединить.
Объединение прошло. Открыл выходной файл.
Вместо координат вида ХХ.000 «улетело» в сторону ХХ0000.
И так все космоснимки.
Решил проверить.
Добавил один из каналов в «чистый проект»: та же картина.
Открываю «Растр — Информация»:

Driver: GTiff/GeoTIFF
Files: C:/LT51670232010160MOR00/L5167023_02320100609_B70.TIF
Size is 8051, 7271
Coordinate System is:
PROJCS["WGS 84 / UTM zone 39N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",51],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32639"]]
Origin = (553800.000000000000000,5994000.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Metadata:
AREA_OR_POINT=Point
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 553800.000, 5994000.000) ( 51d49'21.23"E, 54d 5'28.92"N)
Lower Left ( 553800.000, 5775870.000) ( 51d47'9.52"E, 52d 7'50.63"N)
Upper Right ( 795330.000, 5994000.000) ( 55d30'28.53"E, 54d 0'34.61"N)
Lower Right ( 795330.000, 5775870.000) ( 55d18'29.38"E, 52d 3'16.42"N)
Center ( 674565.000, 5884935.000) ( 53d36'23.59"E, 53d 5'7.07"N)
Band 1 Block=8051x1 Type=Byte, ColorInterp=Gray

Открываю метафайл (отрывок):

GROUP = L1_METADATA_FILE
GROUP = METADATA_FILE_INFO
ORIGIN = "Image courtesy of the U.S. Geological Survey"
.......
BAND_COMBINATION = "1234567"
PRODUCT_UL_CORNER_LAT = 54.0913663
PRODUCT_UL_CORNER_LON = 51.8225643
PRODUCT_UR_CORNER_LAT = 54.0096324
PRODUCT_UR_CORNER_LON = 55.5074688
PRODUCT_LL_CORNER_LAT = 52.1310011
PRODUCT_LL_CORNER_LON = 51.7859815
PRODUCT_LR_CORNER_LAT = 52.0548471
PRODUCT_LR_CORNER_LON = 55.3077499
PRODUCT_UL_CORNER_MAPX = 553800.000
PRODUCT_UL_CORNER_MAPY = 5994000.000
PRODUCT_UR_CORNER_MAPX = 795300.000
PRODUCT_UR_CORNER_MAPY = 5994000.000
PRODUCT_LL_CORNER_MAPX = 553800.000
PRODUCT_LL_CORNER_MAPY = 5775900.000
PRODUCT_LR_CORNER_MAPX = 795300.000
PRODUCT_LR_CORNER_MAPY = 5775900.000
.....
GROUP = PROJECTION_PARAMETERS
REFERENCE_DATUM = "WGS84"
REFERENCE_ELLIPSOID = "WGS84"
GRID_CELL_SIZE_THM = 30.000
GRID_CELL_SIZE_REF = 30.000
ORIENTATION = "NUP"
RESAMPLING_OPTION = "CC"
MAP_PROJECTION = "UTM"
END_GROUP = PROJECTION_PARAMETERS
GROUP = UTM_PARAMETERS
ZONE_NUMBER = 39
END_GROUP = UTM_PARAMETERS
END_GROUP = L1_METADATA_FILE
END

То же самое с другими космоснимками (более десяти штук).

То есть вместо чтения из PRODUCT_ХХ_CORNER_LAT/_LON = 51.8225643 идёт чтение из PRODUCT_ХХ_CORNER_MAPX/_MAPY = 5994000.000

Прошу помочь в решении ошибки.

Проверено на версиях 1.7.0-Trunk ревизия 15028 и 1.7.0-Trunk ревизия 14989

[Voltron]

Боюсь, проблема не QGIS, а в GDAL.
Если выполнять склейку из командной строки результат такой же?

[gimran]

Не пытался, так как забыл то, как производить склейку из командной строки.
Если не сложно, напомните, пожалуйста, процедуру и ключи.
Помню, что, кажется, нужно выбрать «Растр → Объединить», выбрать нужные слои, а затем сгенерированную команду скопировать, вставить в OSGeo4W Shell и добавить какие-то ключи...

[Voltron]

Ключи добавляются только если тех, что есть в сгенерированной команде не хватает. А так можно копировать один в один и выполнять

[gimran]

А какие вообще ключи для данной команды есть?

[Voltron]

gdal_merge

[gimran]

Решил проблему:

По шагам:
вызываю OSGeo4W Shell
ввожу:
gdal17
нажимаю Enter
ввожу:
gdal_merge.bat -of GTiff -o E:/cosmo/osen/save/L5164022_02220100604.tif E:/cosmo/osen/save/LT51640222010155MOR00/L5164022_02220100604_B10.TIF E:/cosmo/osen/save/LT51640222010155MOR00/L5164022_02220100604_B20.TIF E:/cosmo/osen/save/LT51640222010155MOR00/L5164022_02220100604_B30.TIF E:/cosmo/osen/save/LT51640222010155MOR00/L5164022_02220100604_B40.TIF E:/cosmo/osen/save/LT51640222010155MOR00/L5164022_02220100604_B50.TIF E:/cosmo/osen/save/LT51640222010155MOR00/L5164022_02220100604_B60.TIF E:/cosmo/osen/save/LT51640222010155MOR00/L5164022_02220100604_B70.TIF
нажимаю Enter

на выходе файл:
L5164022_02220100604.tif

QGIS:
добавляю этот слой на карту.
«Растр — Информация»

Driver: GTiff/GeoTIFF
Files: E:/cosmo/osen/save/L5164022_02220100604.tif

PROJCS["WGS 84 / UTM zone 40N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],

«Установки — Свойства проекта»: ставлю WGS 84 / UTM zone 40N.

Перепроецирование:
«Растр — перепроецировать»
Указываю выходной файл, target SRS: указываю WGS84, Memory caching — ставлю галочку (20 Mb).
Копирую информацию внизу окна.

Командная строка:
ввожу:
gdal17
нажимаю Enter
ввожу:
gdalwarp -s_srs EPSG:32640 -t_srs EPSG:4326 -wm 20 -of GTiff E:/cosmo/osen/save/L5164022_02220100604.tif E:/cosmo/osen/save/L5164022_02220100604_.tif
нажимаю Enter

На выходе: файл L5164022_02220100604_.tif

Добавляю в QGIS и для проверки слой, заведомо «перекрывающий» территорию, охватываемую снимком.

Теперь снимок «попал» на нужное место (в не зависимости от зон WGS84 выбранной проекции слоя).

[gimran]

Проверил. Эта операция без командной строки, в самом QGIS, тоже прекрасно работает с указанными выше параметрами

[gimran]

Чуть не забыл, в «Объединении» нужно ставить галочку на «склеить поканально».