Страница 3 из 6
Re: Привязка растровой карты Арктики
Добавлено: 26 дек 2015, 09:36
KotAlex
Denis Rykov писал(а): ...
Огромное спасибо за информацию ! Приятно сознавать, что с вами можно иметь серьезные дела. Подобная отзывчивость - к сожалению, такая редкость в нынешнем мире ...
Позвольте уточнить всю последовательность действий :
1. Сперва делаем то, что в моем предыдущем посте под п.2, то есть, для заданной проекции ( proj=... )
для пары координат каждой точки, заданных в долготе-широте, получаем пары координат в этой проекции
при помощи утилиты cs2cs.
2. Затем вводим изображение карты, смотрим пиксельные координаты тех точек, географические координаты которых были использованы для расчета в утилите cs2cs ( предыдущий пункт этого сообщения ) и передаем
эти пикселные координаты и рассчитанные для них проекционные координаты в утилиту gdal_translate. Также передаем в эту утилиту имя файла изображения карты.
3. В результате получаем файл изображения карты, на котором наблюдаются координаты точек в проекционных координатах.
4. Нужно еще дополнительно озаботиться тем, чтобы эти проекционные координаты из п. 3 преобразовать в географические.
Правильно я изложил ?
Re: Привязка растровой карты Арктики
Добавлено: 26 дек 2015, 12:50
Denis Rykov
Не совсем.
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 преобразовать в географические.
Не понял вопрос.
Re: Привязка растровой карты Арктики
Добавлено: 26 дек 2015, 13:45
KotAlex
Еще раз большое спасибо ! Кажется, уже можно сказать - слава тебе господи, ухлопали проблему ...
Denis Rykov писал(а):
KotAlex писал(а):4. Нужно еще дополнительно озаботиться тем, чтобы эти проекционные координаты из п. 3 преобразовать в географические.
Не понял вопрос.
Это не вопрос. Я неудачно выразился. Я хотел сказать следующее :
В результате всех описанных действий мы получаем файл изображения карты, на котором наблюдаются координаты точек в проекционных координатах ( например, в статусной строке при перемещении курсора мыши - как вы сами написали ). Однако, если я хочу наблюдать не проекционные, а географические координаты ( широта - долгота ), то надо еще дополнительно преобразовать проекционные координаты в географические. А для этого надо или где-то найти подходящий класс, выполняющий подобное преобразование, или написать этот класс самому.
Re: Привязка растровой карты Арктики
Добавлено: 26 дек 2015, 14:09
trir
для этого есть
proj и куча его портов
Re: Привязка растровой карты Арктики
Добавлено: 26 дек 2015, 17:58
KotAlex
trir писал(а):для этого есть
proj и куча его портов
Опять не угадал.
Нужно что-нибудь для Delphi / C++ Builder
Re: Привязка растровой карты Арктики
Добавлено: 26 дек 2015, 18:41
Denis Rykov
Re: Привязка растровой карты Арктики
Добавлено: 26 дек 2015, 18:45
gamm
KotAlex писал(а):trir писал(а):для этого есть
proj и куча его портов
Опять не угадал.
Нужно что-нибудь для Delphi / C++ Builder
и Дельфи, и Буилдер, уже лет двадцать как, умеют подключать dll, в том числе и PROJ4
Re: Привязка растровой карты Арктики
Добавлено: 26 дек 2015, 19:49
KotAlex
Ну что ж, вы заслужили себе еще одно спасибо ...
И еще одно. Вас не затруднит привести здесь определение строки proj ( proj=... )
для использованной вами выше в привязке растра полярной равнопромежуточной проекции ?
Я хочу в понедельник поэкспериментировать с этим в QGIS. Просто не помню, имеется там такая проекция или нет.
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
Re: Привязка растровой карты Арктики
Добавлено: 26 дек 2015, 22:30
KotAlex
Александр Мурый писал(а):
В QGIS эта проекция есть, называется "North_Pole_Azimuthal_Equidistant", код EPSG:102016
Спасибо и вам также. Уже начинаю руку набивать на раздаче репутаций ...
Сейчас КуГиса под рукой нет, посмотреть не могу. Но где-то там у него в каком-то окошке огромный список всяких проекций. Функция поиска при этом реализована столь своеобразно, что скорее удавишься, чем найдешь то, что надо ...
Re: Привязка растровой карты Арктики
Добавлено: 28 дек 2015, 13:46
KotAlex
Что нужно сделать для того, чтобы увидеть географические координаты для файлов геотифф,
выложенных в постах :
Александр Мурый » 24 дек 2015, 18:03
Denis Rykov » 25 дек 2015, 22:40 ?
Я ввожу оба файла в QGIS через "Слой" -> "Добавить растровый слой", в нижней строке окна видны координаты, очевидно, в проекции. Затем кнопкой "Преобразование координат" в правом нижнем углу вызываю окно "Свойства проекта / Система координат", включаю "Автоматическое перепроецирование координат" и ставлю Систему координат "Pulkovo 1942" или "WGS 84" - изображение или улетает за пределы экрана, или меняется весьма интересным образом, или вообще ничего не происходит ...
Re: Привязка растровой карты Арктики
Добавлено: 28 дек 2015, 13:50
trir
Вы хоть представляете какой формы должно стать ваше изображение в этих проекциях? Лично у меня воображения не хватает
Re: Привязка растровой карты Арктики
Добавлено: 28 дек 2015, 13:56
Александр Мурый
KotAlex писал(а):Что нужно сделать для того, чтобы увидеть географические координаты для файлов геотифф,
Лучше всего подключить модуль "Захват координат", выбрать там желаемую СК и работать с ней.
А вот вид растров в "полярных" координатах с перепроецированием "на лету" на глобус реально не для слабонервных, это да.
Re: Привязка растровой карты Арктики
Добавлено: 28 дек 2015, 14:22
KotAlex
Александр Мурый писал(а):
А вот вид растров в "полярных" координатах с перепроецированием "на лету" на глобус реально не для слабонервных, это да.
Вы имеете в виду такую картинку, как пониже ? Вид действительно шикарно-потрясный. И не для слабонервных, хотя мои нервы выдержали :

- CCC.jpg (196.45 КБ) 9960 просмотров
А можно ли как-то увидеть сетку координат, которая на вашем скриншоте из поста ( 24 дек 2015, 18:03 ) нанесена красными линиями ?
Re: Привязка растровой карты Арктики
Добавлено: 28 дек 2015, 14:29
trir
кажется пропала чутли не половина картинки