Преобразование изображений из одной геопространственной системы координат в другую

Question

Я хотел бы сфотографировать изображения с этой карты: http://earthquake.usgs.gov/regional/nca/soiltype/map/ от линейной проекции на карту листовки. Исходные фрагменты находятся в известных, но нестандартных уровнях масштабирования, а карты лифтов нуждаются в плитке мерцателя Mercator XYZ. В принципе, я знаю, как это сделать - у меня есть функции для изменения координат XY в координатах lat-lng на двух картах, и мне просто нужно отображать пиксели для целевой карты в терминах пикселей на исходной карте.

Это, к сожалению, нетривиально, поскольку исходные пиксели распространяются на сотни разных файлов изображений, и я пытаюсь помещать их в сотни изображений. Есть ли программный пакет, который делает это немного более простым? Если нет библиотеки для работы с данными такого типа, кажется, что действительно должно быть ...

Answer 1

Postgis имеет метод RT_ST_Transform, который под капотом использует GdalWarp. Итак, у вас есть, по крайней мере, эти два варианта. Если вы используете Postgis, вам нужно будет фактически зарегистрировать/импортировать изображения в Postgis, используя raster2pgsql, а затем вызывать RT_ST_Transform на каждом из них, а затем выгружать их снова - это может быть написано в некоторой степени с использованием plpgsql (язык сценариев Postgres). Существует какая-то кривая обучения, связанная с использованием растра Postgis, что может оказаться целесообразным, если вы планируете провести любой другой анализ обработки изображений. Вы также можете написать сценарий оболочки (или аналогичный) для автоматизации gdalwarp, если вы не хотите идти по маршруту Postgis.

Answer 2

Для менее формального метода, чем gdalwarp (отличная программа), вы можете проверить плагин Leaflet Leaflet.imageTransform, который может трансформироваться и изображение «на лету» в браузере.