Географическая привязка данных с помощью QGIS

[updates-bot]

Обсуждение статьи "Географическая привязка данных с помощью QGIS"

http://gis-lab.info/qa/georef-qgis.html

Обновления:
03.01.2010 Статья переделана с учетом новой функциональности привязчика
21.08.2010 Статья еще раз переделана с учетом еще более новой функциональности привязчика
21.08.2010 Добавлен файл использованный в статье для примера.

[p2004r]

Делаю по шагам как написано в статье и такое впечатление, что создаваемый wld файл просто не применяется при загрузке растрового слоя. Результат должен быть повернутое градусов на 15 изображение исходного плана. Может qgis принципиально не может показывать растровые слои повернутые относительно друг друга?

[Максим Дубинин]

версия QGIS?

[Максим Дубинин]

Статья несколько устарела, переделал под новый привязчик, также приложен файл использованный для примеров.

[p2004r]

версия QGIS?

Из дебиан-тестинг репозитория поставил.

[Максим Дубинин]

очень старая версия, попробуйте 1.5 найти или из свн собрать
можете также попробовать использовать полиномиальное вместо линейного преобразования

[Evgenia]

Дорого времени!
Подскажите, пожалуйста, как можно в QGIS трансформировать уже привязанный растр? В идеале хотелось бы перейти от Пулково42 сразу к UTM/WGS84.

QGIS: 1.5.0
OS: Win Vista

[Максим Дубинин]

Пока нельзя, привязчик не понимает исходные координаты не в локальной системе.

Если у вас уже привязанный растр и вы хотите его перепроецировать используйте Raster\Reproject
Либо удалить привязку растра и привязывайте по новой.

[Evgenia]

Большое спасибо!

[Rayna]

Проблема - какой "тяжести" растр может быть привязан? Склеила картинку в фотошопе, сохранила в джипеге - 150 Мб. После набора точек ругается -

[Максим Дубинин]

размер растра тут непричем, перепроверяйте точки, попробуйте отключать-включать по одной

[Rayna]

Прошу прощения - ответ на вопрос был прост - я не проверила путь к растру, а в него закрались кириллистические символы, вот гишка и возмутилась...

[bolotoved]

1) В инструкции не указан последний шаг
Нажать кнопку "Начать привязку". Для некоторых блондинок (типа меня) это фатально.

2) Примечание о том, что привязанные растры перепривязать нельзя, лучше внести в статью. Я пару часов мучился, пока понял в чем дело, а сейчас вот и прочитал, что с такой проблемой уже сталкивались.
Кстати, как правильно убрать привязку из GeoTiff? Я поступил богомерзки: пересохранил geotiff в фотошопе, при этом привязка грохнулась, но подозреваю, что есть более правильный метод.

[Voltron]

Кстати, как правильно убрать привязку из GeoTiff? Я поступил богомерзки: пересохранил geotiff в фотошопе, при этом привязка грохнулась, но подозреваю, что есть более правильный метод.

ИМХО, это единственный вариант, впрочем, вместо Photoshop можно использовать любой граф. редактор.

[Александр Мурый]

Есть ещё вариант через gdal_translate (-co PROFILE=BASELINE) и listgeo:
~$ gdalinfo -noct -nomd ASTGTM_N43E039_dem.tif

Код: Выделить всё

Driver: GTiff/GeoTIFF
Files: ASTGTM_N43E039_dem.tif
       ASTGTM_N43E039_dem.tif.aux.xml
Size is 3601, 3601
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
Origin = (38.999861111111109,44.000138888888891)
Pixel Size = (0.000277777777778,-0.000277777777778)
Corner Coordinates:
Upper Left  (  38.9998611,  44.0001389) ( 38d59'59.50"E, 44d 0'0.50"N)
Lower Left  (  38.9998611,  42.9998611) ( 38d59'59.50"E, 42d59'59.50"N)
Upper Right (  40.0001389,  44.0001389) ( 40d 0'0.50"E, 44d 0'0.50"N)
Lower Right (  40.0001389,  42.9998611) ( 40d 0'0.50"E, 42d59'59.50"N)
Center      (  39.5000000,  43.5000000) ( 39d30'0.00"E, 43d30'0.00"N)
Band 1 Block=3601x1 Type=Int16, ColorInterp=Gray
  Min=0.000 Max=2834.000 
  Minimum=0.000, Maximum=2834.000, Mean=160.715, StdDev=378.829

~$ listgeo -tfw ASTGTM_N43E039_dem.tif

ASTGTM_N43E039_dem.tfw

Код: Выделить всё

0.0002777778
0.0000000000
0.0000000000
-0.0002777778
39.0000000000
44.0000000000

~$ gdal_translate -of GTiff -co "PROFILE=BASELINE" ASTGTM_N43E039_dem.tif ASTGTM_N43E039_dem__nongeo.tif

~$ gdalinfo -noct -nomd ASTGTM_N43E039_dem__nongeo.tif

Код: Выделить всё

Driver: GTiff/GeoTIFF
Files: ASTGTM_N43E039_dem__nongeo.tif
       ASTGTM_N43E039_dem__nongeo.tif.aux.xml
Size is 3601, 3601
Coordinate System is `'
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 3601.0)
Upper Right ( 3601.0,    0.0)
Lower Right ( 3601.0, 3601.0)
Center      ( 1800.5, 1800.5)
Band 1 Block=3601x1 Type=Int16, ColorInterp=Gray
  Min=0.000 Max=2834.000 
  Minimum=0.000, Maximum=2834.000, Mean=160.715, StdDev=378.829

Теперь можем делать с растром, чего нам надо (не меняя размер), а потом восстановим привязку с помощью geotifcp.

P.S. Опции PROFILE=BASELINE и TFW=YES вместе не работают (привязка сохраняется).