MRTSwath не могу победить определение углов при проецирован

[CMorgan]

Ситуация следующая.
Есть изображения MODIS прошедшие через IMAPP. Изображения с тарелки, соответственно есть изредка полосы и отсутствие данных в начале и в конце. Я написал прожку, которая определяет углы непрерывного изображения и генерит prm файл. Файл генерится и поворачивается в нужную проекцию, но почему-то оказывается сдвинуть в растре вбок (если поменять углы UL LR, на UR LL та же картина только файл сдвинут вверх).

Т.е. углы UL, LR, LL в растр попадают, UR обрезан, а слева здоровенная белая полоса. При этом при просмотре в Arcview сгенеренное изображение ложится правильно, т.е. какая-то ошибка с обсчетом расположения изображения на растровом поле.

Пробовал поворачивать в Альберса и в Географическую, в обоих случаях растр оказывается сдвинутым.
Пример PRM Файла

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

INPUT_FILENAME = D:/Modis_Test/Terra1323jul14/MOD02QKM.AM0907140921.hdf
GEOLOCATION_FILENAME = D:/Modis_Test/Terra1323jul14/MOD03.AM0907140921.hdf
INPUT_SDS_NAME = EV_250_RefSB, 1,0
OUTPUT_SPATIAL_SUBSET_TYPE = LAT_LONG
OUTPUT_SPACE_UPPER_LEFT_CORNER (LONG LAT) = 49.1711 56.777
OUTPUT_SPACE_LOWER_RIGHT_CORNER (LONG LAT) = 6.96207 39.7162
OUTPUT_FILENAME = D:\Modis_Test\Terra1323jul14\Georeferenced
OUTPUT_FILE_FORMAT = GEOTIFF_FMT
KERNEL_TYPE (CC/BI/NN) = CC
OUTPUT_PROJECTION_NUMBER = ALBERS
OUTPUT_PROJECTION_PARAMETER = 0.0 0.0 65.0 52.0 100.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 = 500.0

Помогите советом? В какую сторону копать?

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

не совсем понятно, растр оказывается сдвинут или неправильно обрезан или и то и другое?
вы проверяли соответствие углов из PRM и углов растра после перепроецирования?
Какая версия MRTSwath? Не старая ли, я помню у меня тоже были проблемы с тем, что он неправильно обрезал, но это было довольно давно.

[CMorgan]

Растр неправильно обрезан.

Версия: MODIS Reprojection Tool Swath (v2.1 January 2006)

Углы после перепроецирования соответствуют тому что задано в файле PRM (попробовал менять координаты UL, DR между собой - ничего не изменилось, картинка такая же)

Во вложении воткнул файл с тем что у меня получается, это еще более менее случай только слегка правый верхний угол обрезан. При перепроецировании в географическую у меня вообще вся северо-восточная часть снимка рубится .

PRM файл для одного из каналов вложенного файла

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

INPUT_FILENAME = D:/Modis_Test/Aqua1336jul14/MOD02QKM.PM0907140936.hdf
GEOLOCATION_FILENAME = D:/Modis_Test/Aqua1336jul14/MOD03.PM0907140936.hdf
INPUT_SDS_NAME = EV_250_RefSB, 1,0
OUTPUT_SPATIAL_SUBSET_TYPE = LAT_LONG
OUTPUT_SPACE_UPPER_LEFT_CORNER (LONG LAT) = 26.4883 55.5386
OUTPUT_SPACE_LOWER_RIGHT_CORNER (LONG LAT) = 67.5729 35.3332
OUTPUT_FILENAME = D:\Modis_Test\Aqua1336jul14\Georeferenced
OUTPUT_FILE_FORMAT = GEOTIFF_FMT
KERNEL_TYPE (CC/BI/NN) = CC
OUTPUT_PROJECTION_NUMBER = ALBERS
OUTPUT_PROJECTION_PARAMETER = 0.0 0.0 65.0 52.0 100.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 = 500.0

Jgw:

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

 2500.00000000000000
                   0.00000000000000
                   0.00000000000000
               -2500.00000000000000
            -5714999.48082194570000
             7868075.23227590880000

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

можете выложить пример?
ftp://gis-lab.info/incoming

[CMorgan]

Выложил.

MODIS_SWATH_Test.zip