GDALWarp (latitude or longitude exceeded limits)

Ответить
Voin
Интересующийся
Сообщения: 42
Зарегистрирован: 13 сен 2013, 12:05
Репутация: 3
Откуда: Berlin

GDALWarp (latitude or longitude exceeded limits)

Сообщение Voin » 30 янв 2017, 12:38

Добрый день,

пытаюсь трансформировать снимок Sentinel 1 в какую-нибудь проекцию (не важно какую, например UTM или вэб-меркатор). Обычно с этим не было никаких проблем, пока не столкнулся со следующей ситуацией. Снимок пересекает -/+180 меридиан, опорные точки к снимку лежат по обе строны меридиана, например:

GCP[ 15]: Id=16, Info= (7815,0) -> (175.175867378252,-66.8310013895681,-2.82004475593567e-06)

GCP[ 21]: Id=22, Info= (0,507) -> (-178.840974107832,-68.4810317707879,-8.65384936332703e-06)

при попытке спроецировать получаю следующую ошибку "ERROR 1: latitude or longitude exceeded limits"

Подозреваю что проблема/ограничения в самой библиотеке, GDAL/OGR/PROJ4 ...??? Кто-нибудь сталкивался с этой проблемой и есть ли пути решения?

Пользуюсь GDAL С/С++ API, версия 2.1.3 (на машинах с 2.1 и 2.1.2 то же самое)

Заранее благодарен

Аватара пользователя
Евгений Брусникин
Новоприбывший
Сообщения: 14
Зарегистрирован: 31 май 2014, 14:21
Репутация: 0

Re: GDALWarp (latitude or longitude exceeded limits)

Сообщение Евгений Брусникин » 31 янв 2017, 17:15

Пересчитайте исходные долготы в 0 +360,

andverh
Новоприбывший
Сообщения: 8
Зарегистрирован: 29 сен 2014, 04:51
Репутация: 3

Re: GDALWarp (latitude or longitude exceeded limits)

Сообщение andverh » 01 фев 2017, 05:55


Voin
Интересующийся
Сообщения: 42
Зарегистрирован: 13 сен 2013, 12:05
Репутация: 3
Откуда: Berlin

Re: GDALWarp (latitude or longitude exceeded limits)

Сообщение Voin » 01 фев 2017, 11:13

Евгений Брусникин писал(а):Пересчитайте исходные долготы в 0 +360,
такой подход работает не со всеми целевыми проекциями. С EPSG4326 проблем нет, все работает. А вот UTM обрезается на 180 меридиане. С псевдо-меркатором (EPSG3857) вообще ничего не выходит.


Спасибо. Судя по всему, это единственный выход для меня - резать изображения на две части.

gamm
Гуру
Сообщения: 4044
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1049
Ваше звание: программист
Откуда: Казань

Re: GDALWarp (latitude or longitude exceeded limits)

Сообщение gamm » 01 фев 2017, 17:20

Voin писал(а):Спасибо. Судя по всему, это единственный выход для меня - резать изображения на две части.
можно сделать проще - просто убрать в исходниках проверку, и перекомпилировать. Под VC все собирается минут за 20 ...

Voin
Интересующийся
Сообщения: 42
Зарегистрирован: 13 сен 2013, 12:05
Репутация: 3
Откуда: Berlin

Re: GDALWarp (latitude or longitude exceeded limits)

Сообщение Voin » 01 фев 2017, 17:29

gamm писал(а):
Voin писал(а):Спасибо. Судя по всему, это единственный выход для меня - резать изображения на две части.
можно сделать проще - просто убрать в исходниках проверку, и перекомпилировать. Под VC все собирается минут за 20 ...
в моем случае это не проще. Программа работает на нескольких машинах под разными ОС - Ubuntu, CentOS, Solaris. У меня нет прав на перекомпиляцию билиотек. Разве что себе отдельно GDAL собирать, под свой проект..

alexandr cherepanov
Гуру
Сообщения: 534
Зарегистрирован: 30 ноя 2006, 13:31
Репутация: 116
Откуда: Moscow

Re: GDALWarp (latitude or longitude exceeded limits)

Сообщение alexandr cherepanov » 01 фев 2017, 21:05

Если UTM подходит, что мешает перевести опорные точки в метры в 59 или 1 зону и трансформировать?
Если нужно отображать в 3857, тогда конечно резать.

gamm
Гуру
Сообщения: 4044
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1049
Ваше звание: программист
Откуда: Казань

Re: GDALWarp (latitude or longitude exceeded limits)

Сообщение gamm » 02 фев 2017, 18:49

Voin писал(а):в моем случае это не проще.
тогда Plan-B, для UTM. Выгружаете GCP из TIF, каким-нибыдь пакетным редактором (под Унихом их хватает) отнимаете (прибавляете) к долготе 6 градусов (или 12, если 6 мало), трансформируете с указанием соответствующей зоны, и назначаете результаты правильную зону. Немного через назад, но не требует перекомпиляции.

С остальными проекциями подобный фокус может пройти тоже (если соответствующим образом поменять центральный меридиан, а вместо загадочных номеров использовать нормальное задание проекции +proj)

Ответить

Вернуться в «GDAL/OGR»

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

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