Привязка растровой карты Арктики

Системы координат, проекции, преобразования, привязка
KotAlex
Участник
Сообщения: 79
Зарегистрирован: 24 дек 2015, 11:28
Репутация: 1

Re: Привязка растровой карты Арктики

Сообщение KotAlex » 26 дек 2015, 09:36

Denis Rykov писал(а): ...
Огромное спасибо за информацию ! Приятно сознавать, что с вами можно иметь серьезные дела. Подобная отзывчивость - к сожалению, такая редкость в нынешнем мире ...
Позвольте уточнить всю последовательность действий :

1. Сперва делаем то, что в моем предыдущем посте под п.2, то есть, для заданной проекции ( proj=... )
для пары координат каждой точки, заданных в долготе-широте, получаем пары координат в этой проекции
при помощи утилиты cs2cs.

2. Затем вводим изображение карты, смотрим пиксельные координаты тех точек, географические координаты которых были использованы для расчета в утилите cs2cs ( предыдущий пункт этого сообщения ) и передаем
эти пикселные координаты и рассчитанные для них проекционные координаты в утилиту gdal_translate. Также передаем в эту утилиту имя файла изображения карты.

3. В результате получаем файл изображения карты, на котором наблюдаются координаты точек в проекционных координатах.

4. Нужно еще дополнительно озаботиться тем, чтобы эти проекционные координаты из п. 3 преобразовать в географические.

Правильно я изложил ?

Аватара пользователя
Denis Rykov
Гуру
Сообщения: 3376
Зарегистрирован: 11 апр 2008, 21:09
Репутация: 529
Ваше звание: Author
Контактная информация:

Re: Привязка растровой карты Арктики

Сообщение Denis Rykov » 26 дек 2015, 12:50

Не совсем.
KotAlex писал(а):3. В результате получаем файл изображения карты, на котором наблюдаются координаты точек в проекционных координатах.
В результате получаем файл c GCP, по которым и будет выполняться трансформация на следующем этапе с помощью утилиты gdalwarp:

Код: Выделить всё

gdalwarp -t_srs "+proj=aeqd +lat_0=90 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs" -r cubicspline -order 1 -co COMPRESS=NONE "/tmp/arctic.jpg" "/home/rykov/tmp/arctic.tif"
Параметр -t_srs используется для того, чтобы назначить выходному растру указанную СК после трансформации.
KotAlex писал(а):4. Нужно еще дополнительно озаботиться тем, чтобы эти проекционные координаты из п. 3 преобразовать в географические.
Не понял вопрос.
Spatial is now, more than ever, just another column- The Geometry Column.

KotAlex
Участник
Сообщения: 79
Зарегистрирован: 24 дек 2015, 11:28
Репутация: 1

Re: Привязка растровой карты Арктики

Сообщение KotAlex » 26 дек 2015, 13:45

Еще раз большое спасибо ! Кажется, уже можно сказать - слава тебе господи, ухлопали проблему ...
Denis Rykov писал(а):
KotAlex писал(а):4. Нужно еще дополнительно озаботиться тем, чтобы эти проекционные координаты из п. 3 преобразовать в географические.
Не понял вопрос.
Это не вопрос. Я неудачно выразился. Я хотел сказать следующее :
В результате всех описанных действий мы получаем файл изображения карты, на котором наблюдаются координаты точек в проекционных координатах ( например, в статусной строке при перемещении курсора мыши - как вы сами написали ). Однако, если я хочу наблюдать не проекционные, а географические координаты ( широта - долгота ), то надо еще дополнительно преобразовать проекционные координаты в географические. А для этого надо или где-то найти подходящий класс, выполняющий подобное преобразование, или написать этот класс самому.

trir
Гуру
Сообщения: 5286
Зарегистрирован: 09 апр 2010, 19:30
Репутация: 1015
Ваше звание: просто мимо прохожу
Откуда: Ё-бург

Re: Привязка растровой карты Арктики

Сообщение trir » 26 дек 2015, 14:09

для этого есть proj и куча его портов

KotAlex
Участник
Сообщения: 79
Зарегистрирован: 24 дек 2015, 11:28
Репутация: 1

Re: Привязка растровой карты Арктики

Сообщение KotAlex » 26 дек 2015, 17:58

trir писал(а):для этого есть proj и куча его портов
Опять не угадал.
Нужно что-нибудь для Delphi / C++ Builder

Аватара пользователя
Denis Rykov
Гуру
Сообщения: 3376
Зарегистрирован: 11 апр 2008, 21:09
Репутация: 529
Ваше звание: Author
Контактная информация:

Re: Привязка растровой карты Арктики

Сообщение Denis Rykov » 26 дек 2015, 18:41

Spatial is now, more than ever, just another column- The Geometry Column.

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

Re: Привязка растровой карты Арктики

Сообщение gamm » 26 дек 2015, 18:45

KotAlex писал(а):
trir писал(а):для этого есть proj и куча его портов
Опять не угадал.
Нужно что-нибудь для Delphi / C++ Builder
и Дельфи, и Буилдер, уже лет двадцать как, умеют подключать dll, в том числе и PROJ4

KotAlex
Участник
Сообщения: 79
Зарегистрирован: 24 дек 2015, 11:28
Репутация: 1

Re: Привязка растровой карты Арктики

Сообщение KotAlex » 26 дек 2015, 19:49

Denis Rykov писал(а):http://www.naadsm.org/iproj4-delphii
Ну что ж, вы заслужили себе еще одно спасибо ...
И еще одно. Вас не затруднит привести здесь определение строки proj ( proj=... )
для использованной вами выше в привязке растра полярной равнопромежуточной проекции ?
Я хочу в понедельник поэкспериментировать с этим в QGIS. Просто не помню, имеется там такая проекция или нет.

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

Re: Привязка растровой карты Арктики

Сообщение Александр Мурый » 26 дек 2015, 20:10

KotAlex писал(а):[Вас не затруднит привести здесь определение строки proj ( proj=... )
для использованной вами выше в привязке растра полярной равнопромежуточной проекции ?
Выше уже было:

Код: Выделить всё

+proj=aeqd +lat_0=90 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
В QGIS эта проекция есть, называется "North_Pole_Azimuthal_Equidistant", код EPSG:102016
Редактор материалов, модератор форума

KotAlex
Участник
Сообщения: 79
Зарегистрирован: 24 дек 2015, 11:28
Репутация: 1

Re: Привязка растровой карты Арктики

Сообщение KotAlex » 26 дек 2015, 22:30

Александр Мурый писал(а): В QGIS эта проекция есть, называется "North_Pole_Azimuthal_Equidistant", код EPSG:102016
Спасибо и вам также. Уже начинаю руку набивать на раздаче репутаций ...
Сейчас КуГиса под рукой нет, посмотреть не могу. Но где-то там у него в каком-то окошке огромный список всяких проекций. Функция поиска при этом реализована столь своеобразно, что скорее удавишься, чем найдешь то, что надо ...

KotAlex
Участник
Сообщения: 79
Зарегистрирован: 24 дек 2015, 11:28
Репутация: 1

Re: Привязка растровой карты Арктики

Сообщение KotAlex » 28 дек 2015, 13:46

Что нужно сделать для того, чтобы увидеть географические координаты для файлов геотифф,
выложенных в постах :
Александр Мурый » 24 дек 2015, 18:03
Denis Rykov » 25 дек 2015, 22:40 ?
Я ввожу оба файла в QGIS через "Слой" -> "Добавить растровый слой", в нижней строке окна видны координаты, очевидно, в проекции. Затем кнопкой "Преобразование координат" в правом нижнем углу вызываю окно "Свойства проекта / Система координат", включаю "Автоматическое перепроецирование координат" и ставлю Систему координат "Pulkovo 1942" или "WGS 84" - изображение или улетает за пределы экрана, или меняется весьма интересным образом, или вообще ничего не происходит ...

trir
Гуру
Сообщения: 5286
Зарегистрирован: 09 апр 2010, 19:30
Репутация: 1015
Ваше звание: просто мимо прохожу
Откуда: Ё-бург

Re: Привязка растровой карты Арктики

Сообщение trir » 28 дек 2015, 13:50

Вы хоть представляете какой формы должно стать ваше изображение в этих проекциях? Лично у меня воображения не хватает

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

Re: Привязка растровой карты Арктики

Сообщение Александр Мурый » 28 дек 2015, 13:56

KotAlex писал(а):Что нужно сделать для того, чтобы увидеть географические координаты для файлов геотифф,
Лучше всего подключить модуль "Захват координат", выбрать там желаемую СК и работать с ней.
А вот вид растров в "полярных" координатах с перепроецированием "на лету" на глобус реально не для слабонервных, это да.
Редактор материалов, модератор форума

KotAlex
Участник
Сообщения: 79
Зарегистрирован: 24 дек 2015, 11:28
Репутация: 1

Re: Привязка растровой карты Арктики

Сообщение KotAlex » 28 дек 2015, 14:22

Александр Мурый писал(а): А вот вид растров в "полярных" координатах с перепроецированием "на лету" на глобус реально не для слабонервных, это да.
Вы имеете в виду такую картинку, как пониже ? Вид действительно шикарно-потрясный. И не для слабонервных, хотя мои нервы выдержали :
CCC.jpg
CCC.jpg (196.45 КБ) 8767 просмотров
А можно ли как-то увидеть сетку координат, которая на вашем скриншоте из поста ( 24 дек 2015, 18:03 ) нанесена красными линиями ?

trir
Гуру
Сообщения: 5286
Зарегистрирован: 09 апр 2010, 19:30
Репутация: 1015
Ваше звание: просто мимо прохожу
Откуда: Ё-бург

Re: Привязка растровой карты Арктики

Сообщение trir » 28 дек 2015, 14:29

кажется пропала чутли не половина картинки

Ответить

Вернуться в «Координаты и привязка»

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

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