Заставить QGIS интерпретировать координаты как long-lat вместо lat-long?

13

Я работаю с файлами TRMM 3B31 по осадкам, предоставленным НАСА. Когда я загружаю их в QGIS, выбирая проекцию WGS84 вместе с другими файлами форм, карты отображаются вертикально, то есть QGIS интерпретирует как широту, а не долготу, и наоборот. Я пробовал gdal_translate option -a_ullr, выбирая правильные координаты, но это не решает проблему: изображения только сжимаются и растягиваются, чтобы удовлетворить мои требования, но это явно не то, что я хотел бы иметь.

Когда вы загружаете CSV-файл, вы можете выбрать широту и долготу. Есть ли способ сделать то же самое для растров?

Поскольку у меня есть 193 растра, было бы здорово, если бы вы могли предложить команду, которая может выполняться итеративно.

введите описание изображения здесь

На изображении показана карта Индии и первая полоса файла 3B31, то есть поверхностные осадки. Две черные линии - это параллели 40 и -40. Я хотел бы, чтобы они были горизонтальными, и поэтому данные о дожде перекрываются с Индией.

Паоло Сантини
источник
Тестовое изображение было бы очень полезно.
user30184
@ user30184 gis.stackexchange.com/questions/89261/trmm-data-download дает много ссылок.
AndreJ

Ответы:

7

Я попробовал 3b31 hdf файл от mirador и нашел ту же ориентацию, с которой вы столкнулись. Поэтому я создал GCP для вращения против часовой стрелки:

gdal_translate -of VRT HDF4_SDS:UNKNOWN:"3b31.hdf":0 0.vrt
gdal_translate -of VRT -gcp 0 0 -180 -40 -gcp 160 0 -180 40 -gcp 0 720 180 -40 -gcp 160 720 180 40 0.vrt 1.vrt
gdalwarp -r bilinear -t_srs EPSG:4326 1.vrt 0.tif

с этим результатом:

введите описание изображения здесь

Мой поворот должен следовать тегу метаданных Origin = SOUTHWEST.

Andrej
источник
Это именно то, что я искал, хотя я не могу правильно использовать команды, которые вы написали. Можете ли вы объяснить их использование более подробно?
Паоло Сантини
Поместите их в текстовый файл с именем rotate.batв той же папке, что и файл HDF. Проверьте имя исходного файла, запустите оболочку OSGEO4W, перейдите в папку с помощью cd ... и запустите пакет. Старая добрая навигация по DOS всегда полезна ;-)
AndreJ
3

Я скачал 3B42RT.2016031806.7.7day.tif и .tfw, но не вижу ничего странного в геопривязке.

введите описание изображения здесь

Gdalinfo отчет об этом изображении

gdalinfo 3B42RT.2016031806.7.7day.tif
Driver: GTiff/GeoTIFF
Files: 3B42RT.2016031806.7.7day.tif
       3B42RT.2016031806.7.7day.tfw
Size is 1440, 480
Coordinate System is `'
Origin = (-180.000000000000000,60.000000000000000)
Pixel Size = (0.250000000000000,-0.250000000000000)
Metadata:
  TIFFTAG_DATETIME=2016:03:18 10:25:21
  TIFFTAG_DOCUMENTNAME=./image/3B42RT.2016031806.7.7day.tif
  TIFFTAG_IMAGEDESCRIPTION=IDL TIFF file
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_SOFTWARE=IDL 7.1.1, ITT Visual Information Solutions
  TIFFTAG_XRESOLUTION=100
  TIFFTAG_YRESOLUTION=100
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (-180.0000000,  60.0000000)
Lower Left  (-180.0000000, -60.0000000)
Upper Right ( 180.0000000,  60.0000000)
Lower Right ( 180.0000000, -60.0000000)
Center      (   0.0000000,   0.0000000)
Band 1 Block=1440x2 Type=UInt16, ColorInterp=Gray
user30184
источник
Кажется, что файлы 3b42 tif имеют правильную географическую привязку, а файлы 3b31 hdf - нет.
AndreJ
Скорее всего, HDF сложен и не полностью поддерживается GDAL trac.osgeo.org/gdal/wiki/HDF .
user30184
Именно некоторые из тех, что вы найдете в Mirador, имеют географическую привязку, выбирая свойства NetCDF в окне настройки, другие - нет.
Паоло Сантини
@ user30184 поставщик данных должен включить правильные теги метаданных с географической привязкой. Вы получаете то, за что заплатили ...
AndreJ