Форумы GIS-Lab.info

Нужно привязать растр по огромному количеству точек.
Пробую через gdal_translate при помощи параметра -gcp, но примерно на 8 тыс. точке вылетает ошибка: Как я понимаю это из-за огромного числа точек (226695 точек), так как при небольшом количестве все нормально. Как обойти ограничение (если проблема в этом)? Может через GDAL API как-то можно?
Заранее спасибо.

а в растре сколько точек?

2048x4445

wowka1319 писал(а):2048x4445

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

Растр со спутника. Там свой формат. Ошибка исключена.

меня смущает, что для такого небольшого растра столько точек привязки... там походу сам скрипт с ума сходит..
вон пацаны пишут много импортных слов.. там у дядьки аналогичная проблема:
http://lists.osgeo.org/pipermail/gdal-d ... 31744.html

wowka1319 писал(а):Растр со спутника. Там свой формат. Ошибка исключена.

это не ответ. У вас что, собственный спутник, или вы сами распаковывали поток? если нет, то у вас Level1 как минимум, и более-менее стандартный формат.

Да, у меня Level1, а точнее формат HRPT.L1B, из которого я и вытащил растр и lat/lon. Там в каждой строке 51 пара lat/lon. В моем примере в снимке 4445 строк, что в итоге дает 51*4445 = 226695 точек. Я конечно же могу проредить массив lat/lon до ~8300 точек (примерно столько gdal_translate переваривает), но очень не хочется пренебрегать исключенными точками - рука не поднимается.

GDAL вроде умеет сразу переводить HDF4=>GeoTIFF, например grasswiki.osgeo.org/wiki/AVHRR

У меня не HDF4, а HRPT.L1B

по моему, автоматизация: порезать на (226695/8300) файлов с полной по пиксельной привязкой с последующей трансформацией 28 файлов в правильную СК решит задачу в течении часа.

GDAL 2.0 умеет работать с HRPT.L1B. Задача решена.
А вообще даже если и не умел, то в GDAL есть API для указания GCP.