Страница 1 из 1

датасет с геопривязкой в geotiff

Добавлено: 14 апр 2015, 22:14
mogis
Всем здоровья!
Помогите, пожалуйста, решить проблему с привязкой, тк сам весьма далек от геодезии и картографии, но оч вас всех люблю.
Есть 3 одинаковых набора данных NxN: датасет с долготой, датасет с широтой и датасет со значением параметра в точке. Расстояние между двумя соседними точками - 1 км, соответственно с увеличением широты шаг в градусах между соседними пикселями возрастает.
Нужно это перегнать в gtiff используя gdal инструментарий (gdal_translate). Проблема в том, что не могу понять, что делать с координатами, тк gdal, как я полагаю, хочет перед преобразованием знать от меня угловые координаты рассматриваемого набора данных. И на выходе (если отрисовать на меркаторе) получаются квадраты, а хотелось бы увидеть трапеции.
Пишу с надеждой, что есть какое-то простое решение, если же его нет, то хоть подскажите куда копать.
Премного благодарен!

Re: датасет с геопривязкой в geotiff

Добавлено: 14 апр 2015, 22:42
Александр Мурый
А посмотреть на сами "датасеты" можно?

Re: датасет с геопривязкой в geotiff

Добавлено: 15 апр 2015, 01:05
ericsson
А на сколько велико значение N? Потому как если там совсем немного с севера на юг, плюс - средняя полоса России, то расхождение в километр между масштабами у южной и северной границы набежит не так уж скоро.

Re: датасет с геопривязкой в geotiff

Добавлено: 15 апр 2015, 07:44
mogis
Александр Мурый писал(а):А посмотреть на сами "датасеты" можно?
да,конечно,залью пример как приду на работу
ericsson писал(а):А на сколько велико значение N? Потому как если там совсем немного с севера на юг, плюс - средняя полоса России, то расхождение в километр между масштабами у южной и северной границы набежит не так уж скоро.
500х500, но хотелось бы, чтобы расхождений не было вовсе :D
И я видимо ошибся разделом, мне надо было скорее в gdal/ogr, извиняюсь,не сориентировался.
В принципе решение придумал, почему-то не обратил внимание, что у gdal_translate есть -gcp опция, те я могу задать помимо ul(x,y) lr(x,y) еще и ll(x,y) ur(x,y) и получить что хочу. Или задать 1000 точек по 2 на каждую строчку.

Re: датасет с геопривязкой в geotiff

Добавлено: 15 апр 2015, 07:59
gamm
mogis писал(а):gdal, как я полагаю, хочет перед преобразованием знать от меня угловые координаты рассматриваемого набора данных.
вы неправильно полагаете. Координаты должны быть в проекции, в той, в которой у вас расстояние 1км (или в любой равнопромежуточной). Получите нормальный растр в той же проекции, который затем можете корежить в соответствии с пожеланиями.