Зональная статистика в QGIS – не показывает часть значений растра.

[Natalia Novoselova]

Не могу понять

Загружаю растр – продукт MODIS LAI_500m, MCD15A2H https://lpdaac.usgs.gov/dataset_discove ... 15a2h_v006
Загружаю его тайлами, которые сначала склеиваю в QGIS (Raster- Miscellaneous-Merge) в один растр, а затем пересохраняю, чтобы присвоить нереальным значениям (от 10.01 до 25.4) NoData.
Итоговый растр имеет значения от 0 до 9.991.
Что похоже на правду (т.к. для параметра LAI это число нужно поделить на scale factor 0.1 , т.е. умножить на 10 = 99,9%).
Считаю по этому растру зональную статистику в QGIS (Raster - Zonal Statistics) на наложенных полигональных объектах. И в ней наибольшее значение параметра max (т.е. наибольшее вообще значение растра, рассчитанное зональной статистикой для всех полигональных объектов) равно только 7! majotiry тоже равно 7!
Думала,может быть мои объекты случайно не захватили области растра, где он от 7 до 9.991. Чтобы проверить я сделала векторный слой из 1 объекта, в который включила весь Merged растр. То же самое. tmax= 7. majotiry на этом пробном оно почему-то равно 0, мб из-за его обрисовки.

Главное, куда в обоих случаях зональная статистика дела значения растра от 7 до 9.991?!

[Natalia Novoselova]

Тыкая идентификатором в разные пиксели растра, там где LAI должен бвть наибольшим (т.е. внутри больших фрагментов лесов) - я и правда не вижу значений выше 7.
Хотя я не могу проверить так каждый пиксель. Может и есть где-то.

То есть - может быть глюк не от зональной статистики, а от самого продукта? Который показывает максимум до 9.991..
не понимаю..

[alexandr cherepanov]

Посмотрите гистограмму, например у этой гранулы MCD15A2H.A2017241.h13v09.006.2017250135218.hdf нет значений больше 70. 250-255 - это заливка.

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

Filename: C:\Users\alex\Downloads\MCD15A2H.A2017241.h13v09.006.2017250135218.hdf
Dims: Full Scene (5,760,000 points)

Basic Stats    Min    Max       Mean      Stdev
     Band 1      0    255  63.753157  90.345868

Histogram       DN     Npts    Total   Percent   Acc Pct
Band 1           0      140      140    0.0024    0.0024
                 1     3509     3649    0.0609    0.0634
                 2    14607    18256    0.2536    0.3169
                 3    51732    69988    0.8981    1.2151
                 4    59815   129803    1.0385    2.2535
                 5   109966   239769    1.9091    4.1627
                 6   211492   451261    3.6717    7.8344
                 7   181574   632835    3.1523   10.9867
                 8   220178   853013    3.8225   14.8093
                 9   154172  1007185    2.6766   17.4859
                10   206907  1214092    3.5921   21.0780
                11   259314  1473406    4.5020   25.5800
                12   222516  1695922    3.8631   29.4431
                13   216766  1912688    3.7633   33.2064
                14   258273  2170961    4.4839   37.6903
                15   218874  2389835    3.7999   41.4902
                16   193829  2583664    3.3651   44.8553
                17   170628  2754292    2.9623   47.8176
                18   120868  2875160    2.0984   49.9160
                19   140891  3016051    2.4460   52.3620
                20   134555  3150606    2.3360   54.6980
                21    84462  3235068    1.4664   56.1644
                22   123085  3358153    2.1369   58.3013
                23    70922  3429075    1.2313   59.5326
                24    71772  3500847    1.2460   60.7786
                25    77355  3578202    1.3430   62.1216
                26    52708  3630910    0.9151   63.0366
                27    59324  3690234    1.0299   64.0666
                28    44067  3734301    0.7651   64.8316
                29    33760  3768061    0.5861   65.4177
                30    39118  3807179    0.6791   66.0969
                31    29401  3836580    0.5104   66.6073
                32    23634  3860214    0.4103   67.0176
                33    20253  3880467    0.3516   67.3692
                34    20433  3900900    0.3547   67.7240
                35    15009  3915909    0.2606   67.9845
                36    14086  3929995    0.2445   68.2291
                37    12366  3942361    0.2147   68.4438
                38    10322  3952683    0.1792   68.6230
                39     8320  3961003    0.1444   68.7674
                40     9344  3970347    0.1622   68.9296
                41     8760  3979107    0.1521   69.0817
                42     7483  3986590    0.1299   69.2116
                43     7316  3993906    0.1270   69.3386
                44     7044  4000950    0.1223   69.4609
                45     7761  4008711    0.1347   69.5957
                46     8584  4017295    0.1490   69.7447
                47     8346  4025641    0.1449   69.8896
                48    12286  4037927    0.2133   70.1029
                49    10443  4048370    0.1813   70.2842
                50     9687  4058057    0.1682   70.4524
                51     8061  4066118    0.1399   70.5923
                52    10845  4076963    0.1883   70.7806
                53    10152  4087115    0.1762   70.9569
                54     9398  4096513    0.1632   71.1200
                55    12877  4109390    0.2236   71.3436
                56    21774  4131164    0.3780   71.7216
                57    18856  4150020    0.3274   72.0490
                58    11999  4162019    0.2083   72.2573
                59    24274  4186293    0.4214   72.6787
                60    18414  4204707    0.3197   72.9984
                61    22284  4226991    0.3869   73.3853
                62   103059  4330050    1.7892   75.1745
                63    15726  4345776    0.2730   75.4475
                64    14468  4360244    0.2512   75.6987
                65   110795  4471039    1.9235   77.6222
                66   245107  4716146    4.2553   81.8775
                67     2478  4718624    0.0430   81.9206
                68     2309  4720933    0.0401   81.9606
                69      846  4721779    0.0147   81.9753
                70     7754  4729533    0.1346   82.1099
                71        0  4729533    0.0000   82.1099
                72        0  4729533    0.0000   82.1099
                73        0  4729533    0.0000   82.1099
                74        0  4729533    0.0000   82.1099
                75        0  4729533    0.0000   82.1099
                76        0  4729533    0.0000   82.1099
                77        0  4729533    0.0000   82.1099
                78        0  4729533    0.0000   82.1099
                79        0  4729533    0.0000   82.1099
                80        0  4729533    0.0000   82.1099
                81        0  4729533    0.0000   82.1099
                82        0  4729533    0.0000   82.1099
                83        0  4729533    0.0000   82.1099
                84        0  4729533    0.0000   82.1099
                85        0  4729533    0.0000   82.1099
                86        0  4729533    0.0000   82.1099
                87        0  4729533    0.0000   82.1099
                88        0  4729533    0.0000   82.1099
                89        0  4729533    0.0000   82.1099
                90        0  4729533    0.0000   82.1099
                91        0  4729533    0.0000   82.1099
                92        0  4729533    0.0000   82.1099
                93        0  4729533    0.0000   82.1099
                94        0  4729533    0.0000   82.1099
                95        0  4729533    0.0000   82.1099
                96        0  4729533    0.0000   82.1099
                97        0  4729533    0.0000   82.1099
                98        0  4729533    0.0000   82.1099
                99        0  4729533    0.0000   82.1099
               100        0  4729533    0.0000   82.1099
               101        0  4729533    0.0000   82.1099
               102        0  4729533    0.0000   82.1099
               103        0  4729533    0.0000   82.1099
               104        0  4729533    0.0000   82.1099
               105        0  4729533    0.0000   82.1099
               106        0  4729533    0.0000   82.1099
               107        0  4729533    0.0000   82.1099
               108        0  4729533    0.0000   82.1099
               109        0  4729533    0.0000   82.1099
               110        0  4729533    0.0000   82.1099
               111        0  4729533    0.0000   82.1099
               112        0  4729533    0.0000   82.1099
               113        0  4729533    0.0000   82.1099
               114        0  4729533    0.0000   82.1099
               115        0  4729533    0.0000   82.1099
               116        0  4729533    0.0000   82.1099
               117        0  4729533    0.0000   82.1099
               118        0  4729533    0.0000   82.1099
               119        0  4729533    0.0000   82.1099
               120        0  4729533    0.0000   82.1099
               121        0  4729533    0.0000   82.1099
               122        0  4729533    0.0000   82.1099
               123        0  4729533    0.0000   82.1099
               124        0  4729533    0.0000   82.1099
               125        0  4729533    0.0000   82.1099
               126        0  4729533    0.0000   82.1099
               127        0  4729533    0.0000   82.1099
               128        0  4729533    0.0000   82.1099
               129        0  4729533    0.0000   82.1099
               130        0  4729533    0.0000   82.1099
               131        0  4729533    0.0000   82.1099
               132        0  4729533    0.0000   82.1099
               133        0  4729533    0.0000   82.1099
               134        0  4729533    0.0000   82.1099
               135        0  4729533    0.0000   82.1099
               136        0  4729533    0.0000   82.1099
               137        0  4729533    0.0000   82.1099
               138        0  4729533    0.0000   82.1099
               139        0  4729533    0.0000   82.1099
               140        0  4729533    0.0000   82.1099
               141        0  4729533    0.0000   82.1099
               142        0  4729533    0.0000   82.1099
               143        0  4729533    0.0000   82.1099
               144        0  4729533    0.0000   82.1099
               145        0  4729533    0.0000   82.1099
               146        0  4729533    0.0000   82.1099
               147        0  4729533    0.0000   82.1099
               148        0  4729533    0.0000   82.1099
               149        0  4729533    0.0000   82.1099
               150        0  4729533    0.0000   82.1099
               151        0  4729533    0.0000   82.1099
               152        0  4729533    0.0000   82.1099
               153        0  4729533    0.0000   82.1099
               154        0  4729533    0.0000   82.1099
               155        0  4729533    0.0000   82.1099
               156        0  4729533    0.0000   82.1099
               157        0  4729533    0.0000   82.1099
               158        0  4729533    0.0000   82.1099
               159        0  4729533    0.0000   82.1099
               160        0  4729533    0.0000   82.1099
               161        0  4729533    0.0000   82.1099
               162        0  4729533    0.0000   82.1099
               163        0  4729533    0.0000   82.1099
               164        0  4729533    0.0000   82.1099
               165        0  4729533    0.0000   82.1099
               166        0  4729533    0.0000   82.1099
               167        0  4729533    0.0000   82.1099
               168        0  4729533    0.0000   82.1099
               169        0  4729533    0.0000   82.1099
               170        0  4729533    0.0000   82.1099
               171        0  4729533    0.0000   82.1099
               172        0  4729533    0.0000   82.1099
               173        0  4729533    0.0000   82.1099
               174        0  4729533    0.0000   82.1099
               175        0  4729533    0.0000   82.1099
               176        0  4729533    0.0000   82.1099
               177        0  4729533    0.0000   82.1099
               178        0  4729533    0.0000   82.1099
               179        0  4729533    0.0000   82.1099
               180        0  4729533    0.0000   82.1099
               181        0  4729533    0.0000   82.1099
               182        0  4729533    0.0000   82.1099
               183        0  4729533    0.0000   82.1099
               184        0  4729533    0.0000   82.1099
               185        0  4729533    0.0000   82.1099
               186        0  4729533    0.0000   82.1099
               187        0  4729533    0.0000   82.1099
               188        0  4729533    0.0000   82.1099
               189        0  4729533    0.0000   82.1099
               190        0  4729533    0.0000   82.1099
               191        0  4729533    0.0000   82.1099
               192        0  4729533    0.0000   82.1099
               193        0  4729533    0.0000   82.1099
               194        0  4729533    0.0000   82.1099
               195        0  4729533    0.0000   82.1099
               196        0  4729533    0.0000   82.1099
               197        0  4729533    0.0000   82.1099
               198        0  4729533    0.0000   82.1099
               199        0  4729533    0.0000   82.1099
               200        0  4729533    0.0000   82.1099
               201        0  4729533    0.0000   82.1099
               202        0  4729533    0.0000   82.1099
               203        0  4729533    0.0000   82.1099
               204        0  4729533    0.0000   82.1099
               205        0  4729533    0.0000   82.1099
               206        0  4729533    0.0000   82.1099
               207        0  4729533    0.0000   82.1099
               208        0  4729533    0.0000   82.1099
               209        0  4729533    0.0000   82.1099
               210        0  4729533    0.0000   82.1099
               211        0  4729533    0.0000   82.1099
               212        0  4729533    0.0000   82.1099
               213        0  4729533    0.0000   82.1099
               214        0  4729533    0.0000   82.1099
               215        0  4729533    0.0000   82.1099
               216        0  4729533    0.0000   82.1099
               217        0  4729533    0.0000   82.1099
               218        0  4729533    0.0000   82.1099
               219        0  4729533    0.0000   82.1099
               220        0  4729533    0.0000   82.1099
               221        0  4729533    0.0000   82.1099
               222        0  4729533    0.0000   82.1099
               223        0  4729533    0.0000   82.1099
               224        0  4729533    0.0000   82.1099
               225        0  4729533    0.0000   82.1099
               226        0  4729533    0.0000   82.1099
               227        0  4729533    0.0000   82.1099
               228        0  4729533    0.0000   82.1099
               229        0  4729533    0.0000   82.1099
               230        0  4729533    0.0000   82.1099
               231        0  4729533    0.0000   82.1099
               232        0  4729533    0.0000   82.1099
               233        0  4729533    0.0000   82.1099
               234        0  4729533    0.0000   82.1099
               235        0  4729533    0.0000   82.1099
               236        0  4729533    0.0000   82.1099
               237        0  4729533    0.0000   82.1099
               238        0  4729533    0.0000   82.1099
               239        0  4729533    0.0000   82.1099
               240        0  4729533    0.0000   82.1099
               241        0  4729533    0.0000   82.1099
               242        0  4729533    0.0000   82.1099
               243        0  4729533    0.0000   82.1099
               244        0  4729533    0.0000   82.1099
               245        0  4729533    0.0000   82.1099
               246        0  4729533    0.0000   82.1099
               247        0  4729533    0.0000   82.1099
               248        0  4729533    0.0000   82.1099
               249        0  4729533    0.0000   82.1099
               250     4689  4734222    0.0814   82.1914
               251        0  4734222    0.0000   82.1914
               252        0  4734222    0.0000   82.1914
               253     6264  4740486    0.1087   82.3001
               254  1018733  5759219   17.6863   99.9864
               255      781  5760000    0.0136  100.0000

[Natalia Novoselova]

осмотрите гистограмму, например у этой гранулы MCD15A2H.A2017241.h13v09.006.2017250135218.hdf нет значений больше 70. 250-255 - это заливка.

Не поняла. Где смотреть это на гистограмме? В какой колонке там идут значения?
И как эту гистограмму посмотреть - в QGIS?

И в любом случае, если в продукте нет значений выше 7, тогда почему растр показывает значения до 9.991?

[Natalia Novoselova]

А здесь (с другим продуктом - GPP_500m, MOD17A2H) обратная странность. Растр показывает макс значение 0.067.. а гистограмма (и зональная статика) дают максимум больше 0.1 (з.с. его показывает 0.113)
Все в QGIS 2.18

наверное, из рода ошибки..

[alexandr cherepanov]

В этом поле не показывается реальное максимальное значение, Вы там можете вбить любое число хоть 1000000. Реальные значения и их распределение есть на гистограмме.

Выберите Стиль->Значения мин/макс.
Мин/макс, Точность - фактические(медленнее) ->Загрузить -> Применить.

[Natalia Novoselova]

Выберите Стиль->Значения мин/макс.
Мин/макс, Точность - фактические(медленнее) ->Загрузить -> Применить.

Так получилось, спасибо!

И относительно продутктов MODIS странно, что они вместо nodata поставили значения. И если в продуктах GPP, FPAR, LAI nodata имеет действительно нереальные значения, то в продуктах NDVI и EVI (MOD13Q1) за пределами продукта и внутри него в прогалах с отсутствием данных это значение "0"!! Число, которое входит в диапазон РЕЛЬАНЫХ значений индекса!

Только остается, все удалять, что ниже нуля для NDVI. Но зачем они так сделали..