И снова об импорте растра из OZI

[leon_ti]

Дорогие товарищи, прошу Вашей помощи!
Я, согласно рекомендациям из QA, с помощью GDAL сконвертировал привязанный .gif в GeoTIFF.
Полученный файл по gdalinfo выдает следующее:

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

Driver: GTiff/GeoTIFF
Files: 111.tiff
Size is 4494, 3223
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["Pulkovo 1942",
        DATUM["Pulkovo_1942",
            SPHEROID["Krassowsky 1940",6378245,298.2999999999985,
                AUTHORITY["EPSG","7024"]],
            TOWGS84[23.92,-141.27,-80.9,-0,0.35,0.82,-0.12],
            AUTHORITY["EPSG","6284"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4284"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",27],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (631804.224248144659214,6843929.279902421869338)
Pixel Size = (6.349476833277608,-6.349476833277608)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  631804.224, 6843929.280) ( 29d29'23.16"E, 61d40'48.62"N)
Lower Left  (  631804.224, 6823464.916) ( 29d28'30.31"E, 61d29'48.15"N)
Upper Right (  660338.773, 6843929.280) ( 30d 1'40.69"E, 61d40' 9.56"N)
Lower Right (  660338.773, 6823464.916) ( 30d 0'36.45"E, 61d29' 9.39"N)
Center      (  646071.499, 6833697.098) ( 29d45' 2.74"E, 61d34'59.88"N)

При импорте его в QGIS я получаю слой, координаты в котором имеют непонятный для меня вид:


Соответственно, когда я пытаюсь состыковать этот слой с векторным (из OSM), я испытываю фиаско.

Подскажите, пожалуйста, как это исправить?

P.S.
Версия QGIS (Ubuntu Linux):

1.9.0-Master
Ревизия
exported
Версия Qt при сборке
4.8.3
Текущая версия Qt
4.8.3
Версия GDAL/OGR при сборке
1.9.0
Текущая версия GDAL/OGR
1.9.0
Версия GEOS
3.3.3
Версия PostgreSQL
9.1.6
Версия SpatiaLite
3.1.0-RC2
Версия QWT
6.0.0-svn
PROJ.4 Version
470
QScintilla2 Version

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

При импорте его в QGIS я получаю слой, координаты в котором имеют непонятный для меня вид

А какой вид координат является понятным для вас? На первый взгляд, всё выглядит вполне нормально.

Соответственно, когда я пытаюсь состыковать этот слой с векторным (из OSM), я испытываю фиаско.

Видимо, слой из OSM в другой СК/проекции (WGS84?). Попробуйте включить перепроецирование "на лету" (флажок внизу справа).

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

Прочитать это:
http://gis-lab.info/qa/proj-sk-faq.html

Сделать так:

[leon_ti]

Спасибо! Действительно, я не ожидал увидеть метровые координаты. Думал, что растр криво привязан
Кстати, Вы не подскажете - в моей инсталляции не работают плагины OpenLayers и плагины, работающие с проекциями растров (меню Растр-Проекции).
QGis выдает синтаксические ошибки в питоновых скриптах типа такой:

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

Ошибка при выполнении сценария Python:

Traceback (most recent call last):
  File "/usr/share/qgis/python/plugins/GdalTools/GdalTools.py", line 317, in doWarp
    from tools.doWarp import GdalToolsDialog as Warp
  File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 438, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
  File "/usr/share/qgis/python/plugins/GdalTools/tools/doWarp.py", line 33, in 
    from dialogSRS import GdalToolsSRSDialog as SRSDialog
  File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 438, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
  File "/usr/share/qgis/python/plugins/GdalTools/tools/dialogSRS.py", line 50
    def proj4string(self):
      ^
SyntaxError: invalid syntax

В чем тут может быть дело?

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

если вы хотите увеличить шансы на получение ответа на ваши вопросы - не валите их в кучу.

Отделите вопросы не относящиеся к теме в отдельные новые темы и, возможно, кто-нибудь посмотрит в чем проблемы.

2.6

[leon_ti]

Максим, спасибо за ненавязчивое напоминание о культуре форумного общения
Понимаю, что вопрос об ошибках в плагинах выходит за рамки темы в заголовке.
На всякий случай - вдруг кто-нибудь найдет эту тему поиском - сообщаю о решении питоновских проблем, которое сработало у меня.
Решение тут: https://github.com/slarosa/qgis-openlay ... 51154b2cb6
Найдено тут: http://hub.qgis.org/issues/6916

[Sofia_S]

Здравствуйте! Имеется карта генштабовская, привязанная в Ozi Explorer. Растр в .gif. файл привязки .map Конвертировала согласно инструкции с сайта gis-lab.info/qa/qgis-ozi.html
в Geotiff. получился файл .tiff без данных о привязке.
Посмотрела в QGis - 'этот слой не содержит данных о проекции.

Что я делаю не так? Может файлы привязки у меня неправильные? Стоит Quantum GIS 1.8.0 Lisboa на Windows XP
Растр и файл привязки прилагаются.

050k--n37-038-4.map

(3.76 КБ) 702 скачивания

050k--n37-038-4.rar

(1.41 МБ) 397 скачиваний

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

Попробовал с вашим файлом

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

gdalwarp -r near 050k--n37-038-4.gif out.tif

Результат такой, какой должен быть:

c:\temp\1>gdalinfo out.tif
Driver: GTiff/GeoTIFF
Files: out.tif
Size is 2572, 2890
Coordinate System is:
PROJCS["unnamed",
GEOGCS["Pulkovo 1942",
DATUM["Pulkovo_1942",
SPHEROID["Krassowsky 1940",6378245,298.2999999999985,
AUTHORITY["EPSG","7024"]],
AUTHORITY["EPSG","6284"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4284"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",39],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (354167.632581056790000,6081772.511196699900000)
Pixel Size = (7.069508894554487,-7.069508894554487)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 354167.633, 6081772.511) ( 36d43'48.29"E, 54d50'20.28"N)
Lower Left ( 354167.633, 6061341.630) ( 36d44'25.15"E, 54d39'20.08"N)
Upper Right ( 372350.409, 6081772.511) ( 37d 0'46.43"E, 54d50'38.13"N)
Lower Right ( 372350.409, 6061341.630) ( 37d 1'18.70"E, 54d39'37.82"N)
Center ( 363259.021, 6071557.071) ( 36d52'34.65"E, 54d44'59.38"N)

Версия GDAL:

c:\temp\1>gdalinfo --version
GDAL 1.9.2, released 2012/10/08

[Sofia_S]

Вот странно: я делаю то же самое, информация:

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

Driver: GTiff/GeoTIFF
Files: C:/out.tif
       C:/out.tif.aux.xml
Size is 2572, 2890
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["Pulkovo 1942",
        DATUM["Pulkovo_1942",
            SPHEROID["Krassowsky 1940",6378245,298.2999999999985,
                AUTHORITY["EPSG","7024"]],
            AUTHORITY["EPSG","6284"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4284"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",39],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (354167.632581056790000,6081772.511196699900000)
Pixel Size = (7.069508894554487,-7.069508894554487)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  354167.633, 6081772.511) ( 36d43'48.29"E, 54d50'20.28"N)
Lower Left  (  354167.633, 6061341.630) ( 36d44'25.15"E, 54d39'20.08"N)
Upper Right (  372350.409, 6081772.511) ( 37d 0'46.43"E, 54d50'38.13"N)
Lower Right (  372350.409, 6061341.630) ( 37d 1'18.70"E, 54d39'37.82"N)
Center      (  363259.021, 6071557.071) ( 36d52'34.65"E, 54d44'59.38"N)
Band 1 Block=2572x3 Type=Byte, ColorInterp=Palette

Но в ArcMap (собственно для чего затевался весь сыр-бор) полученный Геотифф не отображается с привязкой. Пишет, что отсутствуют данные о привязке.

[Донецков]

Под OZI существует утилита, которая из самого OZI сохраняет для растра world файл...

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

Но в ArcMap (собственно для чего затевался весь сыр-бор) полученный Геотифф не отображается с привязкой. Пишет, что отсутствуют данные о привязке.

Видимо, из-за этого: PROJCS["unnamed"
Может, стоит попробовать через TIF+TFW?

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

gdalwarp -r near -co "PROFILE=BASELINE" -co "TFW=YES" 050k--n37-038-4.gif out.tif

Ну или "насильно" задать нужную проекцию в ArcGIS.

[Sofia_S]

Всем спасибо за ответы, наверное, да, была задана неспроецированная система координат,
можно задать спроецированую, например, таким образом:

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

gdalwarp -t_srs EPSG:32637 -r near -of GTiff -co "PROFILE=BASELINE" -co "TFW=YES" C:/temp/100k--n37-030.gif C:/temp/100k--n37-030.tif

была задана WGS84/UTM zone 37

получившийся файл

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

Driver: GTiff/GeoTIFF
Files: C:/temp/100k--n37-030.tif
       C:/temp/100k--n37-030.tfw
       C:/temp/100k--n37-030.aux
Size is 2661, 3086
Coordinate System is `'
Origin = (466818.659217173760000,6133200.597179904600000)
Pixel Size = (12.866578165900000,-12.866578165900000)
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  466818.659, 6133200.597) 
Lower Left  (  466818.659, 6093494.337) 
Upper Right (  501056.624, 6133200.597) 
Lower Right (  501056.624, 6093494.337) 
Center      (  483937.641, 6113347.467) 
Band 1 Block=2661x3 Type=Byte, ColorInterp=Palette

Несмотря на это, ArcМap выдает сообщение о неопределенной системе координат, но, оказывается, проецирует растры совершенно правильно.

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

Вы сами себе усложняете жизнь. Разумеется Arcmap выдает сообщение о неопределенной СК, ведь у вас:

Coordinate System is `'

Я бы попробовал взять результат выше с PROJCS["unnamed"] и задать ему в каталоге СК из описания. Но все равно странно, потому что описание выглядит правильным и все должно распознаваться как надо.

[Sofia_S]

Я бы попробовал взять результат выше с PROJCS["unnamed"] и задать ему в каталоге СК из описания. Но все равно странно, потому что описание выглядит правильным и все должно распознаваться как надо.

Думаю, это какая-то ошибка ArcGis. Потому что сейчас проверяю с геотиффом, сделанным без всяких файлов привязки, тоже правильно отображает, хотя предупреждение выдает.

[Игорь Белов]

К. О. напоминает, что GeoTIFF потому и Geo, что содержит информацию о привязке внутри TIFF, в тегах.