Страница 1 из 2
статистика изображения - баг или фича
Добавлено: 13 июл 2010, 19:45
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 рассчитывает статистику!
На больших проектах, в которых добавлено десятки снимков, открытие файла проекта занимает десятки минут. потому что для каждого снимка статистика рассчитывается заново.
Чем обусловлено такое поведение? Есть причины или это серьезный баг?
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:08
Максим Дубинин
сходу такое замечание, а где именно хранится статистика? Я думал что в тэгах геотиф, однако нет.
У меня сейчас сломан 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
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:11
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|
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:12
Максим Дубинин
это у тебя GDALTools\Info показывает?
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:16
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|
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:19
Максим Дубинин
Код: Выделить всё
C:\temp>gdalinfo --version
GDAL 1.7.2, released 2010/04/23
у тебя?
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:20
manisha
Код: Выделить всё
$ gdalinfo --version
GDAL 1.6.3, released 2009/11/19
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:22
Максим Дубинин
Пожалуйста, перезагрузи тестовый снимок (или перераспакуй из архива) и посмотри gdalinfo ДО загрузки в кугис.
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:28
manisha
Посмотрел - тоже самое. Этот файл и НЕ был загружен в QGIS, загружал я другую копию. Которую тоже смотрел через gdalinfo ДО и ПОСЛЕ - ничего в статистике не изменилось.
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:31
Максим Дубинин
Для начала, как мы объясним, что у нас для одного файла разное показывается? Точнее, что у меня статистика вообще не показывается. Проверил 1.6 - тоже что и для 1.7
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:33
Максим Дубинин
кстати, почему у нас с тобой разные файлы, у меня 3х канальный, у тебя 1?
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:38
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
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:41
Максим Дубинин
Ну вот, тоже самое. Теперь рассказывай, как ты статистику получаешь.
Код: Выделить всё
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
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:58
Максим Дубинин
Ок. Предпроблема решена, дело в том, что я смотрел чистый TIF, а если распаковать еще и AUX, то тогда все как у тебя. Теперь говори, что дальше.
Re: статистика изображения - баг или фича
Добавлено: 14 июл 2010, 00:59
manisha
Дальше сама проблема: каждый раз при добавлении файла в QGIS статистика рассчитывается заново при измененнии параметров отображения.
Да, под тифом выше я имел ввиду пару tif + aux.