Растровые таблицы MapInfo неверно открываются

[Constantine]

Здравствуйте!
Столкнулся вот с какой проблемой в новой версии QGIS 2.8. Растровые таблицы, созданные в Mapinfo 10.0 стали открываться где то в экваториальной Африке вместо северо-западного региона России. В более ранних версиях QGIS всё было на своих местах. Не подскажете как справиться с проблемой?
Система координат проекта WGS84 UTM 36N. Таблица той же проекции. Перепроецирование "на лету" включено.

P.S. Изображение я искривил внутри намеренно, чтобы "секретные люди" не нажаловались.

[Донецков]

Странно, система координат WGS84 UTM 36N подхватывается автоматом, но координаты отображаются растра, а не привязки...
Возможное предложение - есть несколько утилит, которые привязку tab конвертируют в tfw, а с ней должно открыться нормально, особенно если добавить prj...
А вообще-то зачем "наш" лесной планшет было привязывать к UTM, если он сделан на Пулково 42?

[Constantine]

Здравствуйте!
Надо понимать, что с версии 2.8 QGIS потеряла совместимость с растровыми таблицами MapInfo. Кстати сказать, растры полученные экспортом окна ложатся на место. Есть просьба. Не подскажете пример(ссылку) утилит - конверторов работающих в linux. Что касается выбора системы координат. Мы используем в работе космическую съемку в сочетании с GPS приемниками и лесные карты "посажены" на снимки. Поэтому такая система координат.
И все таки плохо, что QGIS потеряла совместимость с MapInfo. для меня это серьезная проблема в работе. Видимо придется делать откат версии. Выхода похоже нет.

[SergeyRyzhkov]

Не уверен что поможет, но попробуйте через GDAL сделать геотиф на основе привязки.
viewtopic.php?f=30&t=18732

У меня была похожая проблема (причем КуГис 2.6.1 !). Так вот пришлось сначала убирать привязку из тиф,
а потом ее опять "делать" уже на основе информации из tab

[Constantine]

К норме растр приводит перепроецирование в ту же систему координат. Жаль только, что придется создавать дополнительную библиотеку растров специально для работы в QGIS. Мне кажется то, что QGIS неверно открывает растровые таблицы MapInfo это "bug" программы. Может стоит сообщить разработчикам, что скажете?

[SergeyRyzhkov]

Не уверен на счет QGIS, но все возможно.
Может быть косвенно "виноват" QGIS ,
У Вас GDAL какой версии? (с которой работает QGIS, понимаете что их (GDAL) может быть несколько)
Что говорит gdalinfo Example.tif ?

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

не факт, что это ошибка QGIS. У вас противоречащая информация в TIF и TAB.

В GeoTIFF:

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

gdalinfo Example.tif
Warning 1: TIFFFetchNormalTag:Incompatible type for "RichTIFFIPTC"; tag ignored
Warning 1: TIFFFetchNormalTag:Incompatible type for "RichTIFFIPTC"; tag ignored
Driver: GTiff/GeoTIFF
Files: Example.tif
Size is 4808, 4808
Coordinate System is `'
Metadata:
  TIFFTAG_DATETIME=2015:04:15 13:22:45
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_SOFTWARE=Adobe Photoshop CS6 (Windows)
  TIFFTAG_XRESOLUTION=200
  TIFFTAG_YRESOLUTION=200
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 4808.0)
Upper Right ( 4808.0,    0.0)
Lower Right ( 4808.0, 4808.0)
Center      ( 2404.0, 2404.0)
Band 1 Block=4808x1 Type=Byte, ColorInterp=Red
Band 2 Block=4808x1 Type=Byte, ColorInterp=Green
Band 3 Block=4808x1 Type=Byte, ColorInterp=Blue
sim@sim-think:~/work/Example$ gdalinfo Example.tif
Warning 1: TIFFFetchNormalTag:Incompatible type for "RichTIFFIPTC"; tag ignored
Warning 1: TIFFFetchNormalTag:Incompatible type for "RichTIFFIPTC"; tag ignored
Driver: GTiff/GeoTIFF
Files: Example.tif
Size is 4808, 4808
Coordinate System is `'
Metadata:
  TIFFTAG_DATETIME=2015:04:15 13:22:45
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_SOFTWARE=Adobe Photoshop CS6 (Windows)
  TIFFTAG_XRESOLUTION=200
  TIFFTAG_YRESOLUTION=200
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 4808.0)
Upper Right ( 4808.0,    0.0)
Lower Right ( 4808.0, 4808.0)
Center      ( 2404.0, 2404.0)
Band 1 Block=4808x1 Type=Byte, ColorInterp=Red
Band 2 Block=4808x1 Type=Byte, ColorInterp=Green
Band 3 Block=4808x1 Type=Byte, ColorInterp=Blue

С TAB:

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

gdalinfo Example.tif
Warning 1: TIFFFetchNormalTag:Incompatible type for "RichTIFFIPTC"; tag ignored
Warning 1: TIFFFetchNormalTag:Incompatible type for "RichTIFFIPTC"; tag ignored
Driver: GTiff/GeoTIFF
Files: Example.tif
Size is 4808, 4808
Coordinate System is `'
Metadata:
  TIFFTAG_DATETIME=2015:04:15 13:22:45
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_SOFTWARE=Adobe Photoshop CS6 (Windows)
  TIFFTAG_XRESOLUTION=200
  TIFFTAG_YRESOLUTION=200
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 4808.0)
Upper Right ( 4808.0,    0.0)
Lower Right ( 4808.0, 4808.0)
Center      ( 2404.0, 2404.0)
Band 1 Block=4808x1 Type=Byte, ColorInterp=Red
Band 2 Block=4808x1 Type=Byte, ColorInterp=Green
Band 3 Block=4808x1 Type=Byte, ColorInterp=Blue
sim@sim-think:~/work/Example$ gdalinfo Example.tif
Warning 1: TIFFFetchNormalTag:Incompatible type for "RichTIFFIPTC"; tag ignored
Warning 1: TIFFFetchNormalTag:Incompatible type for "RichTIFFIPTC"; tag ignored
Driver: GTiff/GeoTIFF
Files: Example.tif
       Example.TAB
Size is 4808, 4808
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["unnamed",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563],
            TOWGS84[0,0,0,-0,-0,-0,0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",33],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]
GeoTransform =
  585828.5819285229, 1.276317582846247, 0.01719112283963887
  6787703.501272072, 0.01835195661816134, -1.262473791044528
Metadata:
  TIFFTAG_DATETIME=2015:04:15 13:22:45
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_SOFTWARE=Adobe Photoshop CS6 (Windows)
  TIFFTAG_XRESOLUTION=200
  TIFFTAG_YRESOLUTION=200
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  585828.582, 6787703.501) ( 34d35'52.01"E, 61d12'51.39"N)
Lower Left  (  585911.237, 6781633.527) ( 34d35'47.62"E, 61d 9'35.22"N)
Upper Right (  591965.117, 6787791.737) ( 34d42'43.21"E, 61d12'49.22"N)
Lower Right (  592047.772, 6781721.763) ( 34d42'38.11"E, 61d 9'33.06"N)
Center      (  588938.177, 6784712.632) ( 34d39'15.23"E, 61d11'12.26"N)
Band 1 Block=4808x1 Type=Byte, ColorInterp=Red
Band 2 Block=4808x1 Type=Byte, ColorInterp=Green
Band 3 Block=4808x1 Type=Byte, ColorInterp=Blue

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

Дмитрия Барышникова надо спросить для начала про приоритет привязки из внешнего файла и внутренних тегов в GDAL. В нем при живом табе берется информация из него.

Если интересно, можете оформить репорт вот тут, мы посмотрим тогда в кугисе.

[SergeyRyzhkov]

Насколько я помню,когда разбирался, то порядок таков,
сначала определяется привязка в самом файле (через теги), если ее нет, то определяется из файла привязки.
Поэтому и приходится использовать метод "очищения" файла, а потом заново привязывать (как правильно написал Максим, когда привязка в самом тифе и в табе имеет противоречие)

Очищение делается так:
gdal_translate -co "PROFILE=BASELINE" ...

UPD
Хотя, что-то я не увидел противоречие. Тиф - чистый и так, без привязки внутри...

[Донецков]

При аналогии с ecw+tab, в ArcGIS и QGIS привязка бралась из ecw, а в MapInfo из tab... Это было заметно, когда СК в ecw и tab были различны...
Вот-вот, и я для начала открыл и пересохранил растр в Paint, после чего вся привязка "убивается" гарантированно, и все равно QGIS показывает чепуху...

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

Наверное, скажу очевидную вещь, но проблема лечится gdalwarp

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

gdalwarp Example.tif Example_warp.tif

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

Саша, не то и было уже.

Я что-то действительно поспешил с выводом, тиф чистый

Кто-то может вывод gdalinfo сделать и показать на 1.11? А то у меня наша версия 2.0 и в ней всё ок

[SergeyRyzhkov]

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

C:\_dataGIS>gdalinfo --version
GDAL 1.11.1, released 2014/09/24

C:\_dataGIS>gdalinfo Example.tif
Warning 1: TIFFFetchNormalTag:Incompatible type for "RichTIFFIPTC"; tag ignored
Warning 1: TIFFFetchNormalTag:Incompatible type for "RichTIFFIPTC"; tag ignored
Driver: GTiff/GeoTIFF
Files: Example.tif
       Example.TAB
Size is 4808, 4808
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["unnamed",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563],
            TOWGS84[0,0,0,-0,-0,-0,0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",33],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]
GeoTransform =
  585828.5819285229, 1.276317582846244, 0.01719112283964264
  6787703.501272073, 0.01835195661808237, -1.262473791044608
Metadata:
  TIFFTAG_DATETIME=2015:04:15 13:22:45
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_SOFTWARE=Adobe Photoshop CS6 (Windows)
  TIFFTAG_XRESOLUTION=200
  TIFFTAG_YRESOLUTION=200
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  585828.582, 6787703.501) ( 34d35'52.01"E, 61d12'51.39"N)
Lower Left  (  585911.237, 6781633.527) ( 34d35'47.62"E, 61d 9'35.22"N)
Upper Right (  591965.117, 6787791.737) ( 34d42'43.21"E, 61d12'49.22"N)
Lower Right (  592047.772, 6781721.763) ( 34d42'38.11"E, 61d 9'33.06"N)
Center      (  588938.177, 6784712.632) ( 34d39'15.23"E, 61d11'12.26"N)
Band 1 Block=4808x1 Type=Byte, ColorInterp=Red
Band 2 Block=4808x1 Type=Byte, ColorInterp=Green
Band 3 Block=4808x1 Type=Byte, ColorInterp=Blue

[SergeyRyzhkov]

Для интереса открыл в QGIS Desktop 2.8.1 - действительно косяк ...
Сделал, как правильно предложил Александр Мурый, (gdalwarp), то есть "внедрили" в тиф привязку - стало все OK

[Constantine]

Здравствуйте.
Растр у меня гарантировано чистый, потому как все привязывалось в MapInfo своими руками и сажалось на место ими же по "космосу". И, конечно, перепроецирование в ту же проекцию все нормализует с образованием нового растра GeoTif.
У меня в руках библиотека этих изображений ( около 5000 растров), которые до версии QGIS 2.0, кажется, открывались как есть, а теперь чтобы воспользоваться каждым файлом надо создать новый растр. Сейчас меня интересует существует ли способ не удваивая библиотеку растров, и не портя исходный растр как изображение (так как этими изображениями пользуются для построения чертежей), пользоваться имеющейся у меня в распоряжении библиотекой сканированных карт. Есть ли решение?