статистика изображения - баг или фича

[manisha]

Итак, берем любой снимок, для определенности монитор me1_2337_101_1_fusion_nevinnomys_frag_nat.tif из этого архива http://gis-lab.info/data/samples/monito ... r-nevin.7z

Не загружая его в QGIS, смотрим, рассчитана ли для него статистика изображения (утилитой Information из GDALTools). Ага, статистика рассчитана

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

  Minimum=3.000, Maximum=255.000, Mean=101.370, StdDev=47.443

Теперь загружаем его в QGIS, открываем свойства слоя и настраиваем отображение как нам нужно (скажем, указываем 2 стандартных отклонения и растяжение до мин/максимума). Нажимаем ОК.
QGIS радостно рассчитывает статистику (ладно, допустим, ему чем-то не нравится уже имеющаяся).

Ну теперь то статистика есть, верно? На всякий случай проверяем - ага, есть, и она идентичная тому, что было.

Ладно. Удаляем слой из QGIS. Добавляем его заново, снова открываем свойства, указываем нужные настройки, нажимаем ок - и? Снова QGIS рассчитывает статистику!

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

Чем обусловлено такое поведение? Есть причины или это серьезный баг?

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

сходу такое замечание, а где именно хранится статистика? Я думал что в тэгах геотиф, однако нет.
У меня сейчас сломан GDALTools, но просто gdalinfo ничего похожего на то про что ты говоришь вроде не показывает.

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

C:\temp>gdalinfo me1_2337_101_1_fusion_nevinnomys_frag_nat.tif
Driver: GTiff/GeoTIFF
Files: me1_2337_101_1_fusion_nevinnomys_frag_nat.tif
Size is 2505, 3545
Coordinate System is `'
Metadata:
  TIFFTAG_XRESOLUTION=96
  TIFFTAG_YRESOLUTION=96
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 3545.0)
Upper Right ( 2505.0,    0.0)
Lower Right ( 2505.0, 3545.0)
Center      ( 1252.5, 1772.5)
Band 1 Block=2505x1 Type=Byte, ColorInterp=Red
Band 2 Block=2505x1 Type=Byte, ColorInterp=Green
Band 3 Block=2505x1 Type=Byte, ColorInterp=Blue

[manisha]

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

Driver: GTiff/GeoTIFF
Files: /media/DATA/tmp/monitor2/me1_2337_101_1_fusion_nevinnomys_frag_nat.tif
Size is 2505, 3545
Coordinate System is `'
Metadata:
  TIFFTAG_XRESOLUTION=96
  TIFFTAG_YRESOLUTION=96
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 3545.0)
Upper Right ( 2505.0,    0.0)
Lower Right ( 2505.0, 3545.0)
Center      ( 1252.5, 1772.5)
Band 1 Block=2505x1 Type=Byte, ColorInterp=Red
  Min=3.000 Max=255.000 
  Minimum=3.000, Maximum=255.000, Mean=101.370, StdDev=47.443
  Metadata:
    STATISTICS_MINIMUM=3
    STATISTICS_MAXIMUM=255
    STATISTICS_MEAN=101.37018821527
    STATISTICS_MEDIAN=99
    STATISTICS_MODE=57
    STATISTICS_STDDEV=47.44344372184
    STATISTICS_HISTONUMBINS=256
    STATISTICS_HISTOMIN=0
    STATISTICS_HISTOMAX=255
    LAYER_TYPE=athematic
    STATISTICS_HISTOBINVALUES=0|0|0|11774|20184|10092|10092|25230|5887|2523|11774|30276|13456|20184|32799|9251|2523|15138|33640|21025|19343|37004|13456|4205|21025|41209|22707|25230|36163|9251|1682|24389|47096|33640|35322|47096|20184|5887|38686|48778|22707|5046|24389|74849|43732|49619|71485|26071|7569|36163|93351|50460|50460|105125|40368|10933|40368|149698|65598|66439|119422|45414|11774|31117|115217|71485|54665|95874|47937|8410|31117|105966|55506|70644|104284|31117|18502|43732|100079|81577|70644|124468|50460|21025|53824|101761|59711|18502|37004|104284|68962|77372|112694|47937|15138|50460|105125|76531|69803|107648|50460|14297|50460|126150|73167|67280|122786|47096|16820|49619|106807|90828|67280|126991|59711|17661|38686|119422|74849|76531|106807|48778|10092|44573|75690|72326|25230|47096|116058|74849|63916|103443|58870|17661|49619|80736|56347|57188|74849|39527|15979|24389|80736|40368|45414|82418|44573|12615|30276|71485|35322|40368|52983|38686|14297|19343|68121|50460|36163|55506|29435|13456|18502|44573|43732|26071|46255|22707|8410|17661|31958|23548|10092|10933|24389|29435|23548|26071|20184|9251|5887|21025|20184|16820|14297|13456|4205|8410|21866|17661|20184|10092|15979|4205|9251|10092|9251|8410|11774|8410|2523|6728|9251|7569|5887|5046|5887|5046|2523|4205|1682|2523|841|5046|6728|9251|5046|841|2523|1682|2523|4205|1682|5887|2523|2523|4205|841|3364|1682|1682|2523|1682|1682|4205|1682|841|2523|5046|3364|1682|2523|2523|3364|3364|1682|1682|0|0|0|841|2523|841|0|1682|52142|

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

это у тебя GDALTools\Info показывает?

[manisha]

Да.
gdalinfo в командной строке тоже самое показывает:

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

$ gdalinfo /media/DATA/tmp/monitor2/me1_2337_101_1_fusion_nevinnomys_frag_nat.tif
Driver: GTiff/GeoTIFF
Files: /media/DATA/tmp/monitor2/me1_2337_101_1_fusion_nevinnomys_frag_nat.tif
Size is 2505, 3545
Coordinate System is `'
Metadata:
  TIFFTAG_XRESOLUTION=96
  TIFFTAG_YRESOLUTION=96
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 3545.0)
Upper Right ( 2505.0,    0.0)
Lower Right ( 2505.0, 3545.0)
Center      ( 1252.5, 1772.5)
Band 1 Block=2505x1 Type=Byte, ColorInterp=Red
  Min=3.000 Max=255.000 
  Minimum=3.000, Maximum=255.000, Mean=101.370, StdDev=47.443
  Metadata:
    STATISTICS_MINIMUM=3
    STATISTICS_MAXIMUM=255
    STATISTICS_MEAN=101.37018821527
    STATISTICS_MEDIAN=99
    STATISTICS_MODE=57
    STATISTICS_STDDEV=47.44344372184
    STATISTICS_HISTONUMBINS=256
    STATISTICS_HISTOMIN=0
    STATISTICS_HISTOMAX=255
    LAYER_TYPE=athematic
    STATISTICS_HISTOBINVALUES=0|0|0|11774|20184|10092|10092|25230|5887|2523|11774|30276|13456|20184|32799|9251|2523|15138|33640|21025|19343|37004|13456|4205|21025|41209|22707|25230|36163|9251|1682|24389|47096|33640|35322|47096|20184|5887|38686|48778|22707|5046|24389|74849|43732|49619|71485|26071|7569|36163|93351|50460|50460|105125|40368|10933|40368|149698|65598|66439|119422|45414|11774|31117|115217|71485|54665|95874|47937|8410|31117|105966|55506|70644|104284|31117|18502|43732|100079|81577|70644|124468|50460|21025|53824|101761|59711|18502|37004|104284|68962|77372|112694|47937|15138|50460|105125|76531|69803|107648|50460|14297|50460|126150|73167|67280|122786|47096|16820|49619|106807|90828|67280|126991|59711|17661|38686|119422|74849|76531|106807|48778|10092|44573|75690|72326|25230|47096|116058|74849|63916|103443|58870|17661|49619|80736|56347|57188|74849|39527|15979|24389|80736|40368|45414|82418|44573|12615|30276|71485|35322|40368|52983|38686|14297|19343|68121|50460|36163|55506|29435|13456|18502|44573|43732|26071|46255|22707|8410|17661|31958|23548|10092|10933|24389|29435|23548|26071|20184|9251|5887|21025|20184|16820|14297|13456|4205|8410|21866|17661|20184|10092|15979|4205|9251|10092|9251|8410|11774|8410|2523|6728|9251|7569|5887|5046|5887|5046|2523|4205|1682|2523|841|5046|6728|9251|5046|841|2523|1682|2523|4205|1682|5887|2523|2523|4205|841|3364|1682|1682|2523|1682|1682|4205|1682|841|2523|5046|3364|1682|2523|2523|3364|3364|1682|1682|0|0|0|841|2523|841|0|1682|52142|

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

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

C:\temp>gdalinfo --version
GDAL 1.7.2, released 2010/04/23

у тебя?

[manisha]

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

$ gdalinfo --version
GDAL 1.6.3, released 2009/11/19

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

Пожалуйста, перезагрузи тестовый снимок (или перераспакуй из архива) и посмотри gdalinfo ДО загрузки в кугис.

[manisha]

Посмотрел - тоже самое. Этот файл и НЕ был загружен в QGIS, загружал я другую копию. Которую тоже смотрел через gdalinfo ДО и ПОСЛЕ - ничего в статистике не изменилось.

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

Для начала, как мы объясним, что у нас для одного файла разное показывается? Точнее, что у меня статистика вообще не показывается. Проверил 1.6 - тоже что и для 1.7

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

кстати, почему у нас с тобой разные файлы, у меня 3х канальный, у тебя 1?

[manisha]

Другой файл из того же набора

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

$ gdalinfo /media/DATA/tmp/monitor/ME1_PSA_2337_101_2_RG.tif
Driver: GTiff/GeoTIFF
Files: /media/DATA/tmp/monitor/ME1_PSA_2337_101_2_RG.tif
Size is 6959, 6647
Coordinate System is:
PROJCS["WGS 84 / UTM zone 37N",
    GEOGCS["WGS 84",
        DATUM["unknown",
            SPHEROID["unretrievable - using WGS84",6378137,298.257223563,
                AUTHORITY["EPSG","0"]],
            AUTHORITY["EPSG","14824"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",39],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32637"]]
Origin = (696975.477493361919187,4981335.531404845416546)
Pixel Size = (7.999999999999998,-8.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  696975.477, 4981335.531) ( 41d29'50.58"E, 44d57'29.34"N)
Lower Left  (  696975.477, 4928159.531) ( 41d28'36.78"E, 44d28'47.61"N)
Upper Right (  752647.477, 4981335.531) ( 42d12'8.09"E, 44d56'25.96"N)
Lower Right (  752647.477, 4928159.531) ( 42d10'33.54"E, 44d27'45.27"N)
Center      (  724811.477, 4954747.531) ( 41d50'17.27"E, 44d42'39.01"N)
Band 1 Block=6959x1 Type=UInt16, ColorInterp=Gray

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

Ну вот, тоже самое. Теперь рассказывай, как ты статистику получаешь.

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

C:\temp>gdalinfo ME1_PSA_2337_101_2_RG.tif
Driver: GTiff/GeoTIFF
Files: ME1_PSA_2337_101_2_RG.tif
Size is 6959, 6647
Coordinate System is:
PROJCS["WGS 84 / UTM zone 37N",
    GEOGCS["WGS 84",
        DATUM["unknown",
            SPHEROID["unretrievable - using WGS84",6378137,298.257223563
                AUTHORITY["EPSG","0"]],
            AUTHORITY["EPSG","14824"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",39],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32637"]]
Origin = (696975.477493361920000,4981335.531404845400000)
Pixel Size = (7.999999999999998,-8.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  696975.477, 4981335.531) ( 41d29'50.58"E, 44d57'29.34"N)
Lower Left  (  696975.477, 4928159.531) ( 41d28'36.78"E, 44d28'47.61"N)
Upper Right (  752647.477, 4981335.531) ( 42d12'8.09"E, 44d56'25.96"N)
Lower Right (  752647.477, 4928159.531) ( 42d10'33.54"E, 44d27'45.27"N)
Center      (  724811.477, 4954747.531) ( 41d50'17.27"E, 44d42'39.01"N)
Band 1 Block=6959x1 Type=UInt16, ColorInterp=Gray

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

Ок. Предпроблема решена, дело в том, что я смотрел чистый TIF, а если распаковать еще и AUX, то тогда все как у тебя. Теперь говори, что дальше.

[manisha]

Дальше сама проблема: каждый раз при добавлении файла в QGIS статистика рассчитывается заново при измененнии параметров отображения.

Да, под тифом выше я имел ввиду пару tif + aux.