Мой опыт работы: опытный программист, историк-любитель, полный новичок в ГИС (поэтому я использую QGIS). Хотя я научился использовать инструмент геореференции QGIS.
Ввод: предположим, у меня есть изображение, подобное этому Википедия: Файл: VoyagesOfRabbanBarSauma.jpg . Предположим, что это «голый» файл JPG без каких-либо географических метаданных, и автор не может предоставить такую информацию. Входное изображение лицензировано под GNU FDL, под которым я с удовольствием выпускаю свои производные работы.
Цель: импортировать путь, показанный на карте Википедии, в QGIS, чтобы я мог изменить проекцию, наложить высоту, добавить метки для физической географии и т. Д.
Предостережение: одним из обходных путей может быть просто ввести широту / длину местоположений на карте в новый проект ГИС. Я не хочу идти по этому пути, потому что я хотел бы иметь запись оценки этого конкретного автора географической информации.
Что я пробовал: я применил QGIS Georeferencer к растровому файлу JPG с кучей точек по краям континентов. Расстановка точек внутри континента так себе.
Что я хотел бы знать: искажения этой конкретной карты довольно заметны вокруг северного края Евразии, и кажется, что это довольно стандартная проекция. Мне кажется, что должен быть алгоритмический способ определения проекции и параметров проекции, которые сгенерировали эту карту, и восстановления точного неинтерполяционного преобразования между пикселями изображения и широтой / долготой?
Я вижу два вопроса по этому вопросу:
- Определение системы координат шейп-файла, когда неизвестно?
- Неизвестная система координат на старом чертеже
Кажется, что они говорят о методе грубой силы с участием GDAL или QGIS, а новичку в ГИС они кажутся несколько трудоемкими (в основном это геореференция и подгонка их к известным проекциям), но я хочу попробовать их и посмотреть, смогу ли я автоматизировать их.
Инструкции для обоих этих GIS.SE немного выше, чем позволяет мой опыт, и любые советы о том, как именно сделать этот метод грубой силы, будут очень полезны. Например, вопрос 2 выше показывает, как использовать gdaltransform для преобразования списка координат в широту / долготу. Входные координаты xy расположены в пикселях?
источник
Карта использует широту и долготу, как если бы они были декартовыми координатами. Это версия равноугольной или "Plate Carree" проекции. Таким образом, все, что вам нужно сделать, это сместить и изменить масштаб координат x и y, чтобы получить долготу и широту соответственно.
Инструмент географической привязки в ArcGIS 10.0 быстро создает ссылки для оценки преобразования. Хотя он не сообщает о преобразовании, можно экспортировать таблицу ссылок
R
и оценить там коэффициенты. Поместив центральный пиксель в верхний левый угол версии изображения с высоким разрешением в (0,0), получимЧерные контуры страны эталонного слоя (поставляемого ESRI с предыдущими версиями ArcView) показывают, насколько хорошо слои совпадают после применения такого преобразования. Среднеквадратичные ошибки составляют 6-7 км вдоль каждой оси: около одного пикселя на изображении, почти столько, сколько можно ожидать.
источник
Карта, скорее всего, будет схематичной, и расположение городов может быть не идеальным, что делает их рискованными для использования при географической привязке.
Наилучшим подходом может быть захват карты на глаз, используя соответствующую базу в QGIS. Он не будет полностью соответствовать оригиналу, но маршруты, вероятно, будут ориентировочными, поэтому вы вряд ли потеряете точность.
источник