Зональная статистика в QGIS – не показывает часть значений растра.
- Natalia Novoselova
- Гуру
- Сообщения: 3020
- Зарегистрирован: 15 янв 2013, 20:14
- Репутация: 69
- Ваше звание: Лиса
- Откуда: **
- Контактная информация:
Зональная статистика в QGIS – не показывает часть значений растра.
Не могу понять
Загружаю растр – продукт 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?!
Загружаю растр – продукт 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?!
- Вложения
-
- zonal-stat_ex2.jpg (280.71 КБ) 3458 просмотров
- Natalia Novoselova
- Гуру
- Сообщения: 3020
- Зарегистрирован: 15 янв 2013, 20:14
- Репутация: 69
- Ваше звание: Лиса
- Откуда: **
- Контактная информация:
Re: Зональная статистика в QGIS – не показывает часть значений растра.
Тыкая идентификатором в разные пиксели растра, там где LAI должен бвть наибольшим (т.е. внутри больших фрагментов лесов) - я и правда не вижу значений выше 7.
Хотя я не могу проверить так каждый пиксель. Может и есть где-то.
То есть - может быть глюк не от зональной статистики, а от самого продукта? Который показывает максимум до 9.991..
не понимаю..
Хотя я не могу проверить так каждый пиксель. Может и есть где-то.
То есть - может быть глюк не от зональной статистики, а от самого продукта? Который показывает максимум до 9.991..
не понимаю..
-
- Гуру
- Сообщения: 534
- Зарегистрирован: 30 ноя 2006, 13:31
- Репутация: 116
- Откуда: Moscow
Re: Зональная статистика в QGIS – не показывает часть значений растра.
Посмотрите гистограмму, например у этой гранулы 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
- Гуру
- Сообщения: 3020
- Зарегистрирован: 15 янв 2013, 20:14
- Репутация: 69
- Ваше звание: Лиса
- Откуда: **
- Контактная информация:
Re: Зональная статистика в QGIS – не показывает часть значений растра.
Не поняла. Где смотреть это на гистограмме? В какой колонке там идут значения?alexandr cherepanov писал(а): ↑05 ноя 2017, 10:37осмотрите гистограмму, например у этой гранулы MCD15A2H.A2017241.h13v09.006.2017250135218.hdf нет значений больше 70. 250-255 - это заливка.
И как эту гистограмму посмотреть - в QGIS?
И в любом случае, если в продукте нет значений выше 7, тогда почему растр показывает значения до 9.991?
- Natalia Novoselova
- Гуру
- Сообщения: 3020
- Зарегистрирован: 15 янв 2013, 20:14
- Репутация: 69
- Ваше звание: Лиса
- Откуда: **
- Контактная информация:
Re: Зональная статистика в QGIS – не показывает часть значений растра.
А здесь (с другим продуктом - GPP_500m, MOD17A2H) обратная странность. Растр показывает макс значение 0.067.. а гистограмма (и зональная статика) дают максимум больше 0.1 (з.с. его показывает 0.113)
Все в QGIS 2.18
наверное, из рода ошибки..
Все в QGIS 2.18
наверное, из рода ошибки..
- Вложения
-
- gluk1.jpg (60.21 КБ) 3396 просмотров
-
- Гуру
- Сообщения: 534
- Зарегистрирован: 30 ноя 2006, 13:31
- Репутация: 116
- Откуда: Moscow
Re: Зональная статистика в QGIS – не показывает часть значений растра.
В этом поле не показывается реальное максимальное значение, Вы там можете вбить любое число хоть 1000000. Реальные значения и их распределение есть на гистограмме.
Выберите Стиль->Значения мин/макс.
Мин/макс, Точность - фактические(медленнее) ->Загрузить -> Применить.
Выберите Стиль->Значения мин/макс.
Мин/макс, Точность - фактические(медленнее) ->Загрузить -> Применить.
- Natalia Novoselova
- Гуру
- Сообщения: 3020
- Зарегистрирован: 15 янв 2013, 20:14
- Репутация: 69
- Ваше звание: Лиса
- Откуда: **
- Контактная информация:
Re: Зональная статистика в QGIS – не показывает часть значений растра.
Так получилось, спасибо!alexandr cherepanov писал(а): ↑06 ноя 2017, 09:25Выберите Стиль->Значения мин/макс.
Мин/макс, Точность - фактические(медленнее) ->Загрузить -> Применить.
И относительно продутктов MODIS странно, что они вместо nodata поставили значения. И если в продуктах GPP, FPAR, LAI nodata имеет действительно нереальные значения, то в продуктах NDVI и EVI (MOD13Q1) за пределами продукта и внутри него в прогалах с отсутствием данных это значение "0"!! Число, которое входит в диапазон РЕЛЬАНЫХ значений индекса!
Только остается, все удалять, что ниже нуля для NDVI. Но зачем они так сделали..
Кто сейчас на конференции
Сейчас этот форум просматривают: Ahrefs [Bot] и 18 гостей