Форумы GIS-Lab.info

Доброго дня.

Пытаюсь использовать библиотеку 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 байта ).

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

Ну с размером в 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?