Вопрос по MODIS Reprojection Tool

[roooto]

Всем доброго времени суток. Пытаюсь решить задачу при помощи утилиты MODIS Reprojection Tool ссылка . Возникает вопрос, а может ли этот программный продукт сделать то, что мне нужно .
Постараюсь понятно описать задачу: у меня есть некоторый набор данных в синусоидальной проекции MODIS с 1км разрешением (назовем его продукт "А").
Для примера ячейки:

primer.jpg (37.48 КБ) 9305 просмотров

А так же есть файл MOD03*.hdf с матрицей геолокации.
На выходе нужно получить продукт (назовем его "B") с данными из продукта А, но для координат MOD03*.hdf (т.к. в продукте MOD03 размерность 2030x1354 то размерность матрицы продукта "B" так же должна быть размерностью 2030x1354).
Надеюсь я понятно объяснил.
Всем спасибо за ответы.

З.Ы. если MODIS Reprojection Tool не способна выполнить эту задачу, то хотелось бы узнать название утилиты, которая мне подошла бы.

[Максим Дубинин]

если у вас уже есть данные в проекции, то зачем вам MOD03 и MRT?

Что значит "получить продукт с данными из А, для координат..."? У вас уже данные в координатах. Вам другие нужны?

[roooto]

Что значит "получить продукт с данными из А, для координат..."? У вас уже данные в координатах. Вам другие нужны?

Да, мне нужны другие. Мне нужно, чтобы данные, которые привязаны к синусоидальной проекции MODIS стали бы привязаны к координатам в файле MOD03*.hdf

[Максим Дубинин]

Либо вы странного хотите, либо странно выражаете мысли.

Если вы хотите перевести данные из одной системы координат в другую, то так и скажите. Вам нужно чтобы данные были в единицах измерения широта/долгота?

[roooto]

Либо вы странного хотите, либо странно выражаете мысли.

Скорее хочу странного)

Если вы хотите перевести данные из одной системы координат в другую, то так и скажите. Вам нужно чтобы данные были в единицах измерения широта/долгота?

Давайте остановимся пока на этом. Хочу узнать долготу и широту для каждого пикселя данных 1km, которые изначально привязаны к синусоидальной проекции MODIS.

[Максим Дубинин]

gdal

или оно же в QGIS/gdaltools если нужен GUI

[roooto]

gdal

или оно же в QGIS/gdaltools если нужен GUI

Спасибо. Решил.

У меня возникла другая проблема с данным ПО - поэтому решил не создавать новую тему, а спросить тут.
Теперь задача другая. У меня имеются данные MOD09 и файл MOD03 (для MOD09). Мне нужно конвертировать MOD09 в синусоидальную проекцию, чтобы получились данные для определенной ячейки MODIS. Использую для этого ModisSwathTool.
Все бы хорошо, но на выходе (если я все правильно понимаю) должна получиться матрица 1200:1200, а получается 1200:2113. Ячейка h23v03.
Параметры OUTPUT_SPACE_UPPER_LEFT_CORNER и OUTPUT_SPACE_LOWER_RIGHT_CORNER я брал из файла MOD09GA соответствующей ячейки (h23v03).

Ниже привожу файл *.prm, и результат. Спасибо.

Код: Выделить всё

INPUT_FILENAME = D:\convert\MOD09.A2011185.0430.005.2012162143545.hdf
GEOLOCATION_FILENAME = D:\convert\MOD03.A2011185.0430.005.2011207143613.hdf
INPUT_SDS_NAME = 1km_Surface_Reflectance_Band_6
OUTPUT_SPATIAL_SUBSET_TYPE = LAT_LONG
OUTPUT_SPACE_UPPER_LEFT_CORNER (LONG LAT) = 99.999999975 59.999999995
OUTPUT_SPACE_LOWER_RIGHT_CORNER (LONG LAT) = 93.343429595 49.999999996
OUTPUT_FILENAME = D:\convert_out\out
OUTPUT_FILE_FORMAT = HDF_FMT
KERNEL_TYPE (CC/BI/NN) = NN
OUTPUT_PROJECTION_NUMBER = SNSOID
OUTPUT_PROJECTION_PARAMETER = 6371007.181 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
OUTPUT_PROJECTION_SPHERE = 8
OUTPUT_PIXEL_SIZE = 926.6

Результат работы:

im.png (71.78 КБ) 9224 просмотра

[Максим Дубинин]

в чем проблема непонятно

результат правильный или нет? то что охват несколько другой - не обязательно проблема, простым указанием OUTPUT_SPACE углов не получится получить точно заданный охват. Данные результирующие правильные? Если да - дообрезайте как вам нужно постпроцессингом в том же гдале.

PS: Пожалуйста, не пишите "спасибо, решил", лучше ничего не пишите. Решили по другому - пишите как именно, решили как советовали - тоже пишите. Здесь не клуб благотворительный, здесь лаборатория, люди меняются опытом.

[roooto]

в чем проблема непонятно

результат правильный или нет? то что охват несколько другой - не обязательно проблема, простым указанием OUTPUT_SPACE углов не получится получить точно заданный охват. Данные результирующие правильные? Если да - дообрезайте как вам нужно постпроцессингом в том же гдале.

PS: Пожалуйста, не пишите "спасибо, решил", лучше ничего не пишите. Решили по другому - пишите как именно, решили как советовали - тоже пишите. Здесь не клуб благотворительный, здесь лаборатория, люди меняются опытом.

Хорошо, в будущем учту. Решил я проблему с помощью утилиты cs2cs. Вот тут есть немного о ней, только команды я использовал другие: ссылка.

По поводу вопроса... Я не знаю, правильный ли результат. У меня сама задача сводится к тому чтобы получить заданный охват и только после этого я могу узнать правильность результата. Проще говоря, мне нужно из продукта MOD09 получить продукт MOD09GA. Но я не знаю как это сделать

[roooto]

Получилось. Вот что поменял:

Код: Выделить всё

OUTPUT_SPATIAL_SUBSET_TYPE = PROJ_COORDS
OUTPUT_SPACE_UPPER_LEFT_CORNER (X Y) = 5559752.5983333355 6671703.118000003
OUTPUT_SPACE_LOWER_RIGHT_CORNER (X Y) = 6671703.118000003 5559752.5983333355

Получается 45% совпадение перепроектированных данных продукта MOD09 и продукта MOD09GA. Думаю, у меня не верно задан

Код: Выделить всё

OUTPUT_PROJECTION_NUMBER = SNSOID
OUTPUT_PROJECTION_PARAMETER = 6371007.181 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0

Подскажите, как настраивать OUTPUT_PROJECTION_PARAMETER для SNSOID ?