Більше

Геореференції за допомогою квантової ГІС та перепроектування до EPSG: 900913

Геореференції за допомогою квантової ГІС та перепроектування до EPSG: 900913


Я намагаюся повторно спроектувати карту на EPSG: 900913 Зображення, яке я повторюю, таке: http://i.imgur.com/HmfLW.png ">

Отже, Х - довгота в градусах, а Y - широта в градусах.

Я використовую такі налаштування трансформації та цільові srs:

Враховуючи ці налаштування, дія «Створити сценарій GDAL» дає мені такий сценарій:

gdal_translate -of GTiff -gcp 39,084 83,0534 12,9639 55,6032 -gcp 569,16 509,313 24,0161 49,838 109,924 489,771 -gcp 14,3701 50,1065 -gcp 420,153 334,656 20,9839 52,2413 -gcp 426,26 581,374 21,2256 48,7345 -gcp 549,618 142,901 23,8623 54,9019 miasta.bmp miasta_translated.bmp gdalwarp -t_srs EPSG : 900913 -r біля-замовлення 1 -co КОМПРЕС = НІЯКОГО miasta_translated.bmp miasta_modified.tif

Дивно, що це не додавання параметра -t_srs до команди gdalwarp, тому я додав його сам.

Я очікую, що miasta_modified.tif міститиме проекцію оригінальної карти Google Maps, але натомість він містить наступне, що не схоже на проекцію Google Maps, найбільш помітною відмінністю є те, що він стискається у вертикальному напрямку: http: // i.imgur.com/9DMyT.png ">


Це незвичне застосування інструменту геореференцій.

Google Mercartor - це спроектована система координат, тому одиниці вимірювання знаходяться у метрах, і вам фактично не потрібні t_srs, оскільки dstX і dstY повинні бути цільовими координатами srs.

Але в цьому випадку я думаю, що вам також потрібно визначитись вручну-s_srsу gdalwarp як незапроектований лат довжиною, інакше програміст буде думати, що ви працюєте в метрах.