Привязка растра по 200+ тыс. точкам

[wowka1319]

Нужно привязать растр по огромному количеству точек.
Пробую через gdal_translate при помощи параметра -gcp, но примерно на 8 тыс. точке вылетает ошибка:

Скрин консоли с ошибкой

mapping_error.png (5.21 КБ) 9619 просмотров

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

[trir]

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

[wowka1319]

2048x4445

[gamm]

2048x4445

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

[wowka1319]

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

[jerry-maori]

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

[gamm]

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

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

[wowka1319]

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

[gamm]

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

[wowka1319]

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

[Boris]

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

[wowka1319]

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