Geo question: Как сгенерировать файл .wld с некоторыми наземными контрольными точками?

Question

ОК, поэтому у меня есть jpeg и куча контрольных точек земли (x/y пикселей и соответствующие ссылки lat/lon).

Я хотел бы создать глобальный файл .wld для сопровождения jpeg из командной строки. Моя система координат - это Карты Google, то есть EPSG: 900913.

Я знаю, что могу использовать gdal_transform для генерации .vrt, учитывая gcps, но мне нужен файл .wld. (Не совсем ясно, в чем разница, но это определенно то, что мне нужно!)

У кого-нибудь есть идеи, как это сделать?

Благодаря Richard

Answer 1

Мир файл в основном 6 строки ASCII текстовый файл определения вашего геопривязку. Если у вас есть набор GCP, вам нужно сопоставить их (используя некоторый инструмент, например gdal), с одним аффинным преобразованием.

Я не верю, что утилиты командной строки gdal дают вам возможность напрямую создавать файл мира, хотя некоторые драйверы в GDAL будут делать это для вас, когда вы пишете изображение, если вы установите WORLDFILE = yes в драйвере , Вам нужно будет проверить драйвер для вашего конкретного формата, чтобы узнать, поддерживает ли он это.

Если это не так, вы можете сделать это легко вручную. Просто создайте файл .VRT, используя GCP, и посмотрите на него в текстовом редакторе. Он будет иметь раздел, как это:

<GeoTransform>440720.0, 60, 0.0, 3751320.0, 0.0, -60.0</GeoTransform> 

Это «геокоординаты» является аффинным преобразованием используется мировым файлом. Все, что вам нужно сделать, это создать файл ASCII, который помещает, что с одним значением в каждой строке, например, так:

60 
0.0 
0.0 
-60.0 
440720.0 
3751320.0 

Это будет правильный файл .WLD для вашего приложения.

FYI - 6 чисел - это размер пикселя x, сдвиг y на х значение, сдвиг x на пиксель y, происхождение x, затем y начало координат. (Сдвиги обеспечивают возможности поворота/сдвига в аффинном преобразовании. Обычно они будут равны 0/0, так как обычно вы хотите ортотрансформированные изображения).

Подробнее см. Wikipedia's entry on Worldfiles.

Answer 2

Существует Python скрипт gcps2wld.py в хранилище GDAL, который переводит множество опорных точек на файл в первое приближение порядка в мировом формате файл.

Я хотел бы добавить к ответу Рида, что можно дать в качестве опорных точек параметров командной строки для gdal_translate:

gdal_translate -gcp 1 2 3 4 -gcp 6 7 8 9 [-gcp ...] -of GTiff inp.img out.tif