GDAL Geolocation Array?

Программы и алгоритмы для обработки данных дистанционного зондирования: ERDAS, ENVI и другие.
Ответить
CMorgan
Новоприбывший
Сообщения: 5
Зарегистрирован: 03 сен 2008, 15:59
Репутация: 0

GDAL Geolocation Array?

Сообщение CMorgan » 04 сен 2008, 13:16

Доброго дня.

Пытаюсь использовать библиотеку GDAL в проекте на Visual Studio.
Конкретно нужно через механизм Geolocation Array трансформировать изображение в проекцию, за основу взял описание RFC 4 в Gdal и туториал по варпу:
http://gdal.org/warptut.html.

Для примера взял принятые данные MODIS, прогнал их через IMAPP, получил MOD01,MOD02,MOD03. Из Mod021KM вытащил датасэты EV_1KM_RefSB, Latitude, Longitude, раскидал это по трем отдельным файлам в формате ENVI. Пытаюсь теперь преобразовать в проекцию данные из EV_1KM_RefSB, в качестве Geoloacation Array's используя файлы с Latitude, Longitude. Перед открытием всех файлов привожу их к виртуальному формату vrt и работаю уже с ним.

Никак не могу получить какое-либо изображение:(. на выходе получается файл с заголовком и файл с данными (размером 2 байта :shock: ).

Помогите кто-чем может, скажите что я делаю не так?

CMorgan
Новоприбывший
Сообщения: 5
Зарегистрирован: 03 сен 2008, 15:59
Репутация: 0

Re: GDAL Geolocation Array?

Сообщение CMorgan » 05 сен 2008, 11:20

Ну с размером в 2байта более менее справился (проблема была в пропущенных значениях в полях Latitude, Longitude). Картинка получается, но кривая - все сжато в полоску - что-то я криво задаю в параметрах проекции.

Поля Geolocation и SRS в vrt исходного файла (MODIS):

<SRS>GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],
TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,AUTHORITY["EPSG",
"9108"]],AXIS["Lat",NORTH],AXIS["Long",EAST],
AUTHORITY["EPSG","4326"]]</SRS>
<Metadata/>
<Metadata domain="GEOLOCATION">
<MDI key="X_DATASET">C:\T\Lat.vrt</MDI>
<MDI key="X_BAND">1</MDI>
<MDI key="Y_DATASET">C:\T\Lat.vrt</MDI>
<MDI key="Y_BAND">1</MDI>
<MDI key="PIXEL_OFFSET">0</MDI>
<MDI key="PIXEL_STEP">4</MDI>
<MDI key="LINE_OFFSET">0</MDI>
<MDI key="LINE_STEP">4</MDI>
</Metadata>

Вывожу в географическую

Никто не привязывал через GDAL API? :(

Ответить

Вернуться в «Обработка ДДЗ»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость