Привязка растра по 200+ тыс. точкам
-
- Участник
- Сообщения: 93
- Зарегистрирован: 04 дек 2013, 02:14
- Репутация: 1
Привязка растра по 200+ тыс. точкам
Нужно привязать растр по огромному количеству точек.
Пробую через gdal_translate при помощи параметра -gcp, но примерно на 8 тыс. точке вылетает ошибка: Как я понимаю это из-за огромного числа точек (226695 точек), так как при небольшом количестве все нормально. Как обойти ограничение (если проблема в этом)? Может через GDAL API как-то можно?
Заранее спасибо.
Пробую через gdal_translate при помощи параметра -gcp, но примерно на 8 тыс. точке вылетает ошибка: Как я понимаю это из-за огромного числа точек (226695 точек), так как при небольшом количестве все нормально. Как обойти ограничение (если проблема в этом)? Может через GDAL API как-то можно?
Заранее спасибо.
Последний раз редактировалось wowka1319 06 фев 2015, 01:19, всего редактировалось 1 раз.
-
- Гуру
- Сообщения: 5286
- Зарегистрирован: 09 апр 2010, 19:30
- Репутация: 1015
- Ваше звание: просто мимо прохожу
- Откуда: Ё-бург
Re: Привязка растра по 200+ тыс. точкам
а в растре сколько точек?
-
- Участник
- Сообщения: 93
- Зарегистрирован: 04 дек 2013, 02:14
- Репутация: 1
Re: Привязка растра по 200+ тыс. точкам
2048x4445
-
- Гуру
- Сообщения: 4056
- Зарегистрирован: 15 окт 2010, 08:33
- Репутация: 1054
- Ваше звание: программист
- Откуда: Казань
Re: Привязка растра по 200+ тыс. точкам
а растр вы где взяли (и координаты пикселей)? мучают меня смутные подозрения, что вы что-то не то делаете, все должно быть намного проще ...wowka1319 писал(а):2048x4445
-
- Участник
- Сообщения: 93
- Зарегистрирован: 04 дек 2013, 02:14
- Репутация: 1
Re: Привязка растра по 200+ тыс. точкам
Растр со спутника. Там свой формат. Ошибка исключена.
- jerry-maori
- Гуру
- Сообщения: 585
- Зарегистрирован: 22 авг 2012, 17:02
- Репутация: 143
- Откуда: Нижний Новгород
Re: Привязка растра по 200+ тыс. точкам
меня смущает, что для такого небольшого растра столько точек привязки... там походу сам скрипт с ума сходит..
вон пацаны пишут много импортных слов.. там у дядьки аналогичная проблема:
http://lists.osgeo.org/pipermail/gdal-d ... 31744.html
вон пацаны пишут много импортных слов.. там у дядьки аналогичная проблема:
http://lists.osgeo.org/pipermail/gdal-d ... 31744.html
-
- Гуру
- Сообщения: 4056
- Зарегистрирован: 15 окт 2010, 08:33
- Репутация: 1054
- Ваше звание: программист
- Откуда: Казань
Re: Привязка растра по 200+ тыс. точкам
это не ответ. У вас что, собственный спутник, или вы сами распаковывали поток? если нет, то у вас Level1 как минимум, и более-менее стандартный формат.wowka1319 писал(а):Растр со спутника. Там свой формат. Ошибка исключена.
-
- Участник
- Сообщения: 93
- Зарегистрирован: 04 дек 2013, 02:14
- Репутация: 1
Re: Привязка растра по 200+ тыс. точкам
Да, у меня Level1, а точнее формат HRPT.L1B, из которого я и вытащил растр и lat/lon. Там в каждой строке 51 пара lat/lon. В моем примере в снимке 4445 строк, что в итоге дает 51*4445 = 226695 точек. Я конечно же могу проредить массив lat/lon до ~8300 точек (примерно столько gdal_translate переваривает), но очень не хочется пренебрегать исключенными точками - рука не поднимается.
-
- Гуру
- Сообщения: 4056
- Зарегистрирован: 15 окт 2010, 08:33
- Репутация: 1054
- Ваше звание: программист
- Откуда: Казань
Re: Привязка растра по 200+ тыс. точкам
GDAL вроде умеет сразу переводить HDF4=>GeoTIFF, например grasswiki.osgeo.org/wiki/AVHRR
-
- Участник
- Сообщения: 93
- Зарегистрирован: 04 дек 2013, 02:14
- Репутация: 1
Re: Привязка растра по 200+ тыс. точкам
У меня не HDF4, а HRPT.L1B
-
- Гуру
- Сообщения: 4205
- Зарегистрирован: 10 апр 2006, 22:34
- Репутация: 433
- Откуда: Париж
Re: Привязка растра по 200+ тыс. точкам
по моему, автоматизация: порезать на (226695/8300) файлов с полной по пиксельной привязкой с последующей трансформацией 28 файлов в правильную СК решит задачу в течении часа.
-
- Участник
- Сообщения: 93
- Зарегистрирован: 04 дек 2013, 02:14
- Репутация: 1
Re: Привязка растра по 200+ тыс. точкам
GDAL 2.0 умеет работать с HRPT.L1B. Задача решена.
А вообще даже если и не умел, то в GDAL есть API для указания GCP.
А вообще даже если и не умел, то в GDAL есть API для указания GCP.
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 5 гостей