Неправильная трансформация в Georeferencer

Здесь нужно сообщать об ошибках в QGIS, с подробным описанием шагов и версии ПО.
Ответить
Аватара пользователя
Aleksandr Dezhin
Активный участник
Сообщения: 244
Зарегистрирован: 03 фев 2010, 10:47
Репутация: 14
Откуда: Москва

Неправильная трансформация в Georeferencer

Сообщение Aleksandr Dezhin » 15 янв 2012, 04:26

Всем привет,

Что-то непонятное получается у меня с georeferencer. Проблема явно связана с какой-то особенностью исходного файла - с другими исходными файлами все получается нормально. Но вот в чем эта особенность я понять не могу.

Итак:

1) Беру файл http://dl.dropbox.com/u/725384/georef-problem/main.jpg

2) Конвертирую в Geotif:

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

$ gdal_translate -of GTiff -co COMPRESS=LZW -co TILED=YES main.jpg main.tif
Результат - http://dl.dropbox.com/u/725384/georef-problem/main.tif

3) Привязываю в georeferencer в EPSG:900913 к проекту отсюда.
Points-файл, Скриншот 1, Скриншот 2.

4) Запускаю трансформацию, в результате получается довольно странный результат.

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

Раньше мне приходилось сталкиваться с похожей проблемой при привязке gif файлов, но раньше эта проблема решалась конвертацией gif в tiff. Сейчас я попробовал разные способы конвертации jpg в tiff, но это не приносит никаких результатов. Так же я попробовал почистить EXIF-заголовки при помощи:

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

$ jhead -purejpg main.jpg
Попробовал трансформировать JPEG напрямую, тот же результат.

Версия QGIS - сегодняшняя версия из git.
Версия GDAL - 1.6, попробовал с GDAL-1.8.1 - тоже не влияет на результат.

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

$ gdalinfo main.jpg 
Driver: JPEG/JPEG JFIF
Files: main.jpg
Size is 11891, 11904
Coordinate System is `'
Metadata:
  EXIF_Orientation=1
  EXIF_XResolution=(72)
  EXIF_YResolution=(72)
  EXIF_ResolutionUnit=2
  EXIF_Software=Adobe Photoshop CS3 Windows
  EXIF_DateTime=2011:08:25 16:02:32
  EXIF_ColorSpace=65535
  EXIF_PixelXDimension=11891
  EXIF_PixelYDimension=11904
Image Structure Metadata:
  SOURCE_COLOR_SPACE=YCbCr
  INTERLEAVE=PIXEL
  COMPRESSION=JPEG
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,11904.0)
Upper Right (11891.0,    0.0)
Lower Right (11891.0,11904.0)
Center      ( 5945.5, 5952.0)
Band 1 Block=11891x1 Type=Byte, ColorInterp=Red
  Image Structure Metadata:
    COMPRESSION=JPEG
Band 2 Block=11891x1 Type=Byte, ColorInterp=Green
  Image Structure Metadata:
    COMPRESSION=JPEG
Band 3 Block=11891x1 Type=Byte, ColorInterp=Blue
  Image Structure Metadata:
    COMPRESSION=JPEG

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

$ gdalinfo main.tif
Driver: GTiff/GeoTIFF
Files: main.tif
Size is 11891, 11904
Coordinate System is `'
Metadata:
  EXIF_Orientation=1
  EXIF_XResolution=(72)
  EXIF_YResolution=(72)
  EXIF_ResolutionUnit=2
  EXIF_Software=Adobe Photoshop CS3 Windows
  EXIF_DateTime=2011:08:25 16:02:32
  EXIF_ColorSpace=65535
  EXIF_PixelXDimension=11891
  EXIF_PixelYDimension=11904
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,11904.0)
Upper Right (11891.0,    0.0)
Lower Right (11891.0,11904.0)
Center      ( 5945.5, 5952.0)
Band 1 Block=256x256 Type=Byte, ColorInterp=Red
Band 2 Block=256x256 Type=Byte, ColorInterp=Green
Band 3 Block=256x256 Type=Byte, ColorInterp=Blue

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

$ identify -verbose main.jpg
Image: main.jpg
  Format: JPEG (Joint Photographic Experts Group JFIF format)
  Class: DirectClass
  Geometry: 11891x11904+0+0
  Resolution: 72x72
  Print size: 165.153x165.333
  Units: PixelsPerInch
  Type: TrueColor
  Endianess: Undefined
  Colorspace: RGB
  Depth: 8-bit
  Channel depth:
    red: 8-bit
    green: 8-bit
    blue: 8-bit
  Channel statistics:
    Red:
      min: 0 (0)
      max: 255 (1)
      mean: 207.194 (0.812527)
      standard deviation: 78.664 (0.308486)
      kurtosis: 1.59561
      skewness: -1.68845
    Green:
      min: 0 (0)
      max: 255 (1)
      mean: 217.078 (0.851284)
      standard deviation: 58.3912 (0.228985)
      kurtosis: 2.66807
      skewness: -1.74044
    Blue:
      min: 0 (0)
      max: 255 (1)
      mean: 215.56 (0.845333)
      standard deviation: 57.1598 (0.224156)
      kurtosis: 2.13567
      skewness: -1.58958
  Image statistics:
    Overall:
      min: 0 (0)
      max: 255 (1)
      mean: 213.277 (0.836381)
      standard deviation: 65.4848 (0.256803)
      kurtosis: 2.62794
      skewness: -1.80081
  Rendering intent: Undefined
  Interlace: None
  Background color: white
  Border color: rgb(223,223,223)
  Matte color: grey74
  Transparent color: black
  Compose: Over
  Page geometry: 11891x11904+0+0
  Dispose: Undefined
  Iterations: 0
  Compression: JPEG
  Quality: 78
  Orientation: TopLeft
  Properties:
    date:create: 2012-01-13T20:30:28+04:00
    date:modify: 2011-08-25T16:03:44+04:00
    exif:ColorSpace: 65535
    exif:Compression: 6
    exif:DateTime: 2011:08:25 16:02:32
    exif:ExifImageLength: 11904
    exif:ExifImageWidth: 11891
    exif:ExifOffset: 164
    exif:JPEGInterchangeFormat: 302
    exif:JPEGInterchangeFormatLength: 11687
    exif:Orientation: 1
    exif:ResolutionUnit: 2
    exif:Software: Adobe Photoshop CS3 Windows
    exif:XResolution: 720000/10000
    exif:YResolution: 720000/10000
    jpeg:colorspace: 2
    jpeg:sampling-factor: 2x2,1x1,1x1
    signature: abc5ae9f9deadb77262363c6bb2e31f94424ba0c6cf4b5fb62451b0b952e6702
    xapMM:DerivedFrom: 
  Profiles:
    Profile-exif: 11995 bytes
    Profile-iptc: 7 bytes
      unknown[2,0]: 
    Profile-xmp: 3829 bytes
  Artifacts:
    verbose: true
  Tainted: False
  Filesize: 11.72MBB
  Number pixels: 141.6MB
  Pixels per second: 51.29MB
  User time: 2.140u
  Elapsed time: 0:03.759
  Version: ImageMagick 6.6.5-6 2011-08-20 Q16 http://www.imagemagick.org

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

Re: Неправильная трансформация в Georeferencer

Сообщение Александр Мурый » 15 янв 2012, 10:06

Может не в тему, но непонятен пункт №2: зачем из JPG делать геотифф, если результат привязки -- тот же геотифф? Не проще ли привязывать исходный JPG?
Редактор материалов, модератор форума

Аватара пользователя
Aleksandr Dezhin
Активный участник
Сообщения: 244
Зарегистрирован: 03 фев 2010, 10:47
Репутация: 14
Откуда: Москва

Re: Неправильная трансформация в Georeferencer

Сообщение Aleksandr Dezhin » 15 янв 2012, 10:23

amuriy писал(а):Может не в тему, но непонятен пункт №2: зачем из JPG делать геотифф, если результат привязки -- тот же геотифф? Не проще ли привязывать исходный JPG?
У меня есть ощущение, что так оно на больших картинках пошустрее на этапе расстановки контрольных точек работает. Хотя точных замеров не производил. В любом случае на результат это не вляет, т.е. если трансформировать JPEG та же "петрушка" получается.

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

Re: Неправильная трансформация в Georeferencer

Сообщение Александр Мурый » 15 янв 2012, 11:30

Попробовал привязать к EPSG:3857 вместо EPSG:900913 (это ж одно и то же по сути, так?)
Похоже, что растр сел на место:
QGIS_georect_EPSG_3857.png
QGIS_georect_EPSG_3857.png (359.25 КБ) 6310 просмотров
Редактор материалов, модератор форума

Аватара пользователя
Aleksandr Dezhin
Активный участник
Сообщения: 244
Зарегистрирован: 03 фев 2010, 10:47
Репутация: 14
Откуда: Москва

Re: Неправильная трансформация в Georeferencer

Сообщение Aleksandr Dezhin » 15 янв 2012, 19:20

Да, с EPSG:3857 все нормально. Странно как-то :)

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

Re: Неправильная трансформация в Georeferencer

Сообщение Александр Мурый » 15 янв 2012, 22:46

Вот, например, в попавшейся статейке годовой давности делается вывод, что единственный кошерный вариант сферического меркатора --- EPSG:3857 (как самый свежий и действительно утверждённый код EPSG).
Редактор материалов, модератор форума

Ответить

Вернуться в «Ошибки QGIS»

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

Сейчас этот форум просматривают: Semrush [Bot] и 5 гостей