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

Вопросы по свободной ГИС QGIS. Сообщения об ошибках, предложения по улучшению, локализация.
Правила форума
Ответить
leon_ti
Новоприбывший
Сообщения: 5
Зарегистрирован: 10 янв 2013, 20:37
Репутация: 0

И снова об импорте растра из 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
Вернуться к началу
Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 792
Ваше звание: званий не имею
Откуда: Москва

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

Сообщение Александр Мурый »

leon_ti писал(а): При импорте его в QGIS я получаю слой, координаты в котором имеют непонятный для меня вид
А какой вид координат является понятным для вас? На первый взгляд, всё выглядит вполне нормально.
leon_ti писал(а): Соответственно, когда я пытаюсь состыковать этот слой с векторным (из OSM), я испытываю фиаско.
Видимо, слой из OSM в другой СК/проекции (WGS84?). Попробуйте включить перепроецирование "на лету" (флажок внизу справа).
Редактор материалов, модератор форума
Вернуться к началу
Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:
Контактная информация пользователя Максим Дубинин

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

Сообщение Максим Дубинин »

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

Сделать так:
Изображение
пристегивайтесь, турбулентность прямо по курсу
Вернуться к началу
leon_ti
Новоприбывший
Сообщения: 5
Зарегистрирован: 10 янв 2013, 20:37
Репутация: 0

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

Сообщение 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
В чем тут может быть дело?
Вернуться к началу
Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:
Контактная информация пользователя Максим Дубинин

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

Сообщение Максим Дубинин »

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

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

2.6
пристегивайтесь, турбулентность прямо по курсу
Вернуться к началу
leon_ti
Новоприбывший
Сообщения: 5
Зарегистрирован: 10 янв 2013, 20:37
Репутация: 0

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

Сообщение leon_ti »

Максим, спасибо за ненавязчивое напоминание о культуре форумного общения :)
Понимаю, что вопрос об ошибках в плагинах выходит за рамки темы в заголовке.
На всякий случай - вдруг кто-нибудь найдет эту тему поиском - сообщаю о решении питоновских проблем, которое сработало у меня.
Решение тут: https://github.com/slarosa/qgis-openlay ... 51154b2cb6
Найдено тут: http://hub.qgis.org/issues/6916
Вернуться к началу
Sofia_S
Интересующийся
Сообщения: 37
Зарегистрирован: 06 фев 2013, 09:10
Репутация: 5
Откуда: Тула

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

Сообщение 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 КБ) 714 скачиваний
050k--n37-038-4.rar
(1.41 МБ) 407 скачиваний
Вернуться к началу
Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:
Контактная информация пользователя Максим Дубинин

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

Сообщение Максим Дубинин »

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

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

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
Интересующийся
Сообщения: 37
Зарегистрирован: 06 фев 2013, 09:10
Репутация: 5
Откуда: Тула

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

Сообщение 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 (собственно для чего затевался весь сыр-бор) полученный Геотифф не отображается с привязкой. Пишет, что отсутствуют данные о привязке.
Вернуться к началу
Донецков
Гуру
Сообщения: 3058
Зарегистрирован: 19 май 2010, 19:44
Репутация: 189

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

Сообщение Донецков »

Под OZI существует утилита, которая из самого OZI сохраняет для растра world файл...
Вернуться к началу
Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 792
Ваше звание: званий не имею
Откуда: Москва

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

Сообщение Александр Мурый »

Sofia_S писал(а): Но в ArcMap (собственно для чего затевался весь сыр-бор) полученный Геотифф не отображается с привязкой. Пишет, что отсутствуют данные о привязке.
Видимо, из-за этого: PROJCS["unnamed"
Может, стоит попробовать через TIF+TFW?

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

gdalwarp -r near -co "PROFILE=BASELINE" -co "TFW=YES" 050k--n37-038-4.gif out.tif
Ну или "насильно" задать нужную проекцию в ArcGIS.
Редактор материалов, модератор форума
Вернуться к началу
Sofia_S
Интересующийся
Сообщения: 37
Зарегистрирован: 06 фев 2013, 09:10
Репутация: 5
Откуда: Тула

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

Сообщение 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 выдает сообщение о неопределенной системе координат, но, оказывается, проецирует растры совершенно правильно.
Вернуться к началу
Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:
Контактная информация пользователя Максим Дубинин

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

Сообщение Максим Дубинин »

Вы сами себе усложняете жизнь. Разумеется Arcmap выдает сообщение о неопределенной СК, ведь у вас:
Sofia_S писал(а):Coordinate System is `'
Я бы попробовал взять результат выше с PROJCS["unnamed"] и задать ему в каталоге СК из описания. Но все равно странно, потому что описание выглядит правильным и все должно распознаваться как надо.
пристегивайтесь, турбулентность прямо по курсу
Вернуться к началу
Sofia_S
Интересующийся
Сообщения: 37
Зарегистрирован: 06 фев 2013, 09:10
Репутация: 5
Откуда: Тула

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

Сообщение Sofia_S »

Я бы попробовал взять результат выше с PROJCS["unnamed"] и задать ему в каталоге СК из описания. Но все равно странно, потому что описание выглядит правильным и все должно распознаваться как надо.
Думаю, это какая-то ошибка ArcGis. Потому что сейчас проверяю с геотиффом, сделанным без всяких файлов привязки, тоже правильно отображает, хотя предупреждение выдает.
Вернуться к началу
Аватара пользователя
Игорь Белов
Гуру
Сообщения: 2241
Зарегистрирован: 04 янв 2011, 22:00
Репутация: 1514
Откуда: Казань

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

Сообщение Игорь Белов »

К. О. напоминает, что GeoTIFF потому и Geo, что содержит информацию о привязке внутри TIFF, в тегах.
The purpose of computing is insight, not numbers
Вернуться к началу
Ответить

Вернуться в «QGIS»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 3 гостя