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


13

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

Завантажуючи файл .csv, ви можете вибрати, якою має бути широта і яка довгота, чи є спосіб зробити те ж саме для растерів?

Оскільки у мене є растр 193, було б чудово, якщо ви можете запропонувати команду, яку можна виконувати ітеративно.

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

На зображенні показана карта Індії та перша смуга файлу 3B31, тобто поверхневі опади. Дві чорні лінії - паралелі 40 і -40. Я хотів би, щоб вони були горизонтальними, а отже, дані про дощі збігаються з Індією.


Тестове зображення було б дуже корисним.
користувач30184

@ user30184 gis.stackexchange.com/questions/89261/trmm-data-download дає безліч посилань.
AndreJ

Відповіді:


7

Я спробував hb-файл 3b31 від 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.


Це саме те, що я шукав, хоча я не в змозі належним чином використовувати команди, які ви написали. Чи можете ви пояснити їх використання більш докладно?
Паоло Сантіні

Помістіть їх у текстовий файл, названий 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

Здається, що файли 3b42 tif мають правильну геореференцію, в той час як файли в форматі hd 3b31 не мають.
AndreJ

Схоже, HDF є складним і не ідеально підтримується GDAL trac.osgeo.org/gdal/wiki/HDF .
користувач30184

Саме деякі з тих, кого ви знайдете в Mirador, є геореференційованими, вибираючи властивості NetCDF у вікні налаштування, інші - ні.
Паоло Сантіні

@ user30184 постачальник даних повинен включати належні теги геореференцій метаданих. Ви отримуєте те, за що заплатили ...
Анджей
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.