Проблема с модулем i.topo.corr

Ответить
Vedmed
Участник
Сообщения: 70
Зарегистрирован: 19 янв 2012, 14:48
Репутация: 6

Проблема с модулем i.topo.corr

Сообщение Vedmed » 16 фев 2015, 08:58

Установлена GRASS 7.0, система win 7 (если это имеет значение).
Суть проблемы. На первом шаге все замечательно, растр иллюминации строится без проблем. На втором применяю команду:
i.topo.corr input=DCELL_Landsat_B2@PERMANENT output=Illum basemap=Illum_SRTM@PERMANENT zenith=43.34897272
В результате имею однотонный растр. Начинаю разбираться, очень сильный разброс значений, небольшой ряд пикселов имеют совершенно дикие величины.
Range of data: min = -32540.0136844285 max = 53433.5563319935
Спойлер
-32540.013684--32202.862429 |1
-7927.972072--7590.820817 |1
-5905.064542--5567.913287 |1
-2196.400737--1859.249482 |1
-1859.249482--1522.098228 |2
-1522.098228--1184.946973 |2
-1184.946973--847.795718 |1
-847.795718--510.644463 |1
-510.644463--173.493208 |3
-173.493208-163.658047 |9791983
163.658047-500.809302 |76933
500.809302-837.960557 |338
837.960557-1175.111812 |42
1175.111812-1512.263067 |18
1512.263067-1849.414322 |5
1849.414322-2186.565577 |7
2186.565577-2523.716832 |1
2523.716832-2860.868087 |1
2860.868087-3198.019342 |1
3872.321852-4209.473107 |1
5220.926872-5558.078127 |1
9266.741931-9603.893186 |1
9603.893186-9941.044441 |1
35564.539819-35901.691074 |1
41296.111153-41633.262408 |1
53096.405077-53433.556332 |1
* |6130651
Так и должно быть или я что-то делаю не так? По какой причине значения пикселов улетают в космос и проваливаются в бездну? Что с этим делать, посоветуйте?

PS
Спойлер
-0.48392--0.478101|1
-0.454823--0.449004|1
-0.443185--0.437366|1
-0.437366--0.431546|3
-0.425727--0.419908|2
-0.408269--0.40245|1
-0.40245--0.39663|1
-0.384992--0.379173|2
-0.379173--0.373353|4
-0.373353--0.367534|4
-0.367534--0.361715|3
-0.361715--0.355895|3
-0.355895--0.350076|5
-0.350076--0.344257|4
-0.344257--0.338438|9
-0.338438--0.332618|10
-0.332618--0.326799|9
-0.326799--0.32098|11
-0.32098--0.31516|19
-0.31516--0.309341|12
-0.309341--0.303522|20
-0.303522--0.297702|18
-0.297702--0.291883|23
-0.291883--0.286064|32
-0.286064--0.280245|26
-0.280245--0.274425|39
-0.274425--0.268606|33
-0.268606--0.262787|38
-0.262787--0.256967|37
-0.256967--0.251148|54
-0.251148--0.245329|64
-0.245329--0.23951|52
-0.23951--0.23369|68
-0.23369--0.227871|58
-0.227871--0.222052|106
-0.222052--0.216232|91
-0.216232--0.210413|112
-0.210413--0.204594|129
-0.204594--0.198774|123
-0.198774--0.192955|138
-0.192955--0.187136|165
-0.187136--0.181317|177
-0.181317--0.175497|175
-0.175497--0.169678|217
-0.169678--0.163859|241
-0.163859--0.158039|244
-0.158039--0.15222|282
-0.15222--0.146401|285
-0.146401--0.140582|299
-0.140582--0.134762|318
-0.134762--0.128943|310
-0.128943--0.123124|355
-0.123124--0.117304|404
-0.117304--0.111485|416
-0.111485--0.105666|443
-0.105666--0.099846|509
-0.099846--0.094027|553
-0.094027--0.088208|546
-0.088208--0.082389|656
-0.082389--0.076569|652
-0.076569--0.07075|640
-0.07075--0.064931|730
-0.064931--0.059111|766
-0.059111--0.053292|804
-0.053292--0.047473|854
-0.047473--0.041654|850
-0.041654--0.035834|998
-0.035834--0.030015|1079
-0.030015--0.024196|1093
-0.024196--0.018376|1218
-0.018376--0.012557|1219
-0.012557--0.006738|1292
-0.006738--0.000919|1405
-0.000919-0.004901|1432
0.004901-0.01072|1577
0.01072-0.016539|1712
0.016539-0.022359|1749
0.022359-0.028178|1857
0.028178-0.033997|2047
0.033997-0.039817|2157
0.039817-0.045636|2277
0.045636-0.051455|2367
0.051455-0.057274|2417
0.057274-0.063094|2643
0.063094-0.068913|2885
0.068913-0.074732|3003
0.074732-0.080552|3240
0.080552-0.086371|3429
0.086371-0.09219|3778
0.09219-0.098009|3783
0.098009-0.103829|4075
0.103829-0.109648|4318
0.109648-0.115467|4669
0.115467-0.121287|4914
0.121287-0.127106|5132
0.127106-0.132925|5367
0.132925-0.138745|5751
0.138745-0.144564|6212
0.144564-0.150383|6469
0.150383-0.156202|6943
0.156202-0.162022|7317
0.162022-0.167841|7485
0.167841-0.17366|8070
0.17366-0.17948|8449
0.17948-0.185299|8945
0.185299-0.191118|9226
0.191118-0.196937|9742
0.196937-0.202757|10012
0.202757-0.208576|10641
0.208576-0.214395|11251
0.214395-0.220215|11622
0.220215-0.226034|12056
0.226034-0.231853|12467
0.231853-0.237673|13256
0.237673-0.243492|13627
0.243492-0.249311|14412
0.249311-0.25513|15000
0.25513-0.26095|15818
0.26095-0.266769|16317
0.266769-0.272588|17362
0.272588-0.278408|18111
0.278408-0.284227|18921
0.284227-0.290046|19542
0.290046-0.295865|20335
0.295865-0.301685|21406
0.301685-0.307504|22445
0.307504-0.313323|23288
0.313323-0.319143|24118
0.319143-0.324962|25294
0.324962-0.330781|26462
0.330781-0.3366|27614
0.3366-0.34242|28493
0.34242-0.348239|30045
0.348239-0.354058|30730
0.354058-0.359878|32343
0.359878-0.365697|33700
0.365697-0.371516|35718
0.371516-0.377336|37013
0.377336-0.383155|38195
0.383155-0.388974|39735
0.388974-0.394793|41360
0.394793-0.400613|42955
0.400613-0.406432|44289
0.406432-0.412251|46337
0.412251-0.418071|47839
0.418071-0.42389|49644
0.42389-0.429709|51383
0.429709-0.435528|53189
0.435528-0.441348|55145
0.441348-0.447167|57172
0.447167-0.452986|59228
0.452986-0.458806|61173
0.458806-0.464625|63501
0.464625-0.470444|65641
0.470444-0.476264|67971
0.476264-0.482083|69999
0.482083-0.487902|72859
0.487902-0.493721|74825
0.493721-0.499541|77071
0.499541-0.50536|79402
0.50536-0.511179|82099
0.511179-0.516999|85129
0.516999-0.522818|87261
0.522818-0.528637|90797
0.528637-0.534456|93195
0.534456-0.540276|96306
0.540276-0.546095|99512
0.546095-0.551914|101955
0.551914-0.557734|104784
0.557734-0.563553|107968
0.563553-0.569372|111529
0.569372-0.575192|114426
0.575192-0.581011|117081
0.581011-0.58683|121211
0.58683-0.592649|123538
0.5926
49-0.598469|127592
0.598469-0.604288|130703
0.604288-0.610107|133522
0.610107-0.615927|136198
0.615927-0.621746|139114
0.621746-0.627565|142478
0.627565-0.633384|145873
0.633384-0.639204|148882
0.639204-0.645023|151963
0.645023-0.650842|155392
0.650842-0.656662|158060
0.656662-0.662481|162388
0.662481-0.6683|165523
0.6683-0.67412|169157
0.67412-0.679939|173691
0.679939-0.685758|177340
0.685758-0.691577|183506
0.691577-0.697397|189127
0.697397-0.703216|198457
0.703216-0.709035|211248
0.709035-0.714855|228804
0.714855-0.720674|254376
0.720674-0.726493|283518
0.726493-0.732312|303594
0.732312-0.738132|274543
0.738132-0.743951|246525
0.743951-0.74977|223243
0.74977-0.75559|208441
0.75559-0.761409|198404
0.761409-0.767228|191461
0.767228-0.773047|186656
0.773047-0.778867|181831
0.778867-0.784686|177789
0.784686-0.790505|173014
0.790505-0.796325|169384
0.796325-0.802144|165929
0.802144-0.807963|161743
0.807963-0.813783|158617
0.813783-0.819602|154706
0.819602-0.825421|150192
0.825421-0.83124|147207
0.83124-0.83706|144216
0.83706-0.842879|139863
0.842879-0.848698|136374
0.848698-0.854518|133038
0.854518-0.860337|130300
0.860337-0.866156|126233
0.866156-0.871975|121811
0.871975-0.877795|118743
0.877795-0.883614|114959
0.883614-0.889433|111975
0.889433-0.895253|108205
0.895253-0.901072|103470
0.901072-0.906891|99994
0.906891-0.912711|96081
0.912711-0.91853|91781
0.91853-0.924349|88423
0.924349-0.930168|84165
0.930168-0.935988|79725
0.935988-0.941807|76334
0.941807-0.947626|72210
0.947626-0.953446|67325
0.953446-0.959265|63292
0.959265-0.965084|58601
0.965084-0.970903|53696
0.970903-0.976723|48201
0.976723-0.982542|41687
0.982542-0.988361|35509
0.988361-0.994181|28390
0.994181-1|20817
*|2342270

Vedmed
Участник
Сообщения: 70
Зарегистрирован: 19 янв 2012, 14:48
Репутация: 6

Re: Проблема с модулем i.topo.corr

Сообщение Vedmed » 17 фев 2015, 00:35

Кажется я начал понимать в чем проблема. Если верить формулам отсюда, то, согласно им, происходит деление на околонулевые величины. Методы cosine и minnaert прямо делят на растр иллюминации (а в нём околонулевых значений полно, см. статистику в прошлом сообщении), с методом c-factor сложнее, если я правильно понял, для его рассчета нужно выразить "m", тогда:
m=(ref_o-a)/cos_i
где ref_o - обрабатываемый растр, cos_i - растр иллюминации, а "а" - азимут Солнца. Вновь деление на растр иллюминации, отчего "m", при околонулевых значениях растра, оказывается устремлена в бесконечность и при расчете "це фактора":
с=a/m
"це" оказывается крайне незначительной величиной. В этом случае, формула расчета исправленной освещенности практически повторяет формулу cosine со всеми ее недостатками.
Единственный выход, получается, метод percent, добавление единицы к растру иллюминации избавляет от деления на ноль. Правда остается загвоздка. Ряд значений все равно уходит далеко за геотиффовскую планку в 255 единиц. Что с ними делать? Я думаю нужно перекодировать растр, но на каких правилах? Все что выше 255 приравнять к 255, или как то иначе?
Спойлер
0.081281-2.398584|14817
2.398584-4.715887|1035
4.715887-7.03319|587
7.03319-9.350493|414
9.350493-11.667796|744
11.667796-13.985099|642
13.985099-16.302402|545
16.302402-18.619704|664
18.619704-20.937007|1276
20.937007-23.25431|4739
23.25431-25.571613|15847
25.571613-27.888916|35160
27.888916-30.206219|53517
30.206219-32.523522|128338
32.523522-34.840825|393508
34.840825-37.158128|867010
37.158128-39.475431|1370578
39.475431-41.792734|1841722
41.792734-44.110037|2242715
44.110037-46.42734|2522707
46.42734-48.744643|2687862
48.744643-51.061946|2675237
51.061946-53.379249|2514154
53.379249-55.696551|2343676
55.696551-58.013854|2289412
58.013854-60.331157|2158377
60.331157-62.64846|1919345
62.64846-64.965763|1820506
64.965763-67.283066|1688025
67.283066-69.600369|1509130
69.600369-71.917672|1323738
71.917672-74.234975|1127648
74.234975-76.552278|906106
76.552278-78.869581|705525
78.869581-81.186884|562182
81.186884-83.504187|454660
83.504187-85.82149|380943
85.82149-88.138793|330273
88.138793-90.456096|298042
90.456096-92.773399|266012
92.773399-95.090701|237114
95.090701-97.408004|213271
97.408004-99.725307|190211
99.725307-102.04261|168446
102.04261-104.359913|148925
104.359913-106.677216|132027
106.677216-108.994519|114684
108.994519-111.311822|100617
111.311822-113.629125|90226
113.629125-115.946428|80287
115.946428-118.263731|71584
118.263731-120.581034|63310
120.581034-122.898337|56744
122.898337-125.21564|50736
125.21564-127.532943|44253
127.532943-129.850246|38958
129.850246-132.167548|34368
132.167548-134.484851|30943
134.484851-136.802154|28167
136.802154-139.119457|25179
139.119457-141.43676|22584
141.43676-143.754063|20444
143.754063-146.071366|18066
146.071366-148.388669|16037
148.388669-150.705972|14506
150.705972-153.023275|13028
153.023275-155.340578|11640
155.340578-157.657881|10722
157.657881-159.975184|9596
159.975184-162.292487|9414
162.292487-164.60979|8438
164.60979-166.927093|8769
166.927093-169.244395|7380
169.244395-171.561698|6443
171.561698-173.879001|5961
173.879001-176.196304|5617
176.196304-178.513607|5268
178.513607-180.83091|5015
180.83091-183.148213|4895
183.148213-185.465516|4632
185.465516-187.782819|4447
187.782819-190.100122|4382
190.100122-192.417425|4177
192.417425-194.734728|4030
194.734728-197.052031|3905
197.052031-199.369334|3802
199.369334-201.686637|3698
201.686637-204.00394|3650
204.00394-206.321242|3623
206.321242-208.638545|3589
208.638545-210.955848|3508
210.955848-213.273151|3720
213.273151-215.590454|3388
215.590454-217.907757|3413
217.907757-220.22506|3402
220.22506-222.542363|3299
222.542363-224.859666|3217
224.859666-227.176969|3212
227.176969-229.494272|3087
229.494272-231.811575|3050
231.811575-234.128878|3121
234.128878-236.446181|3101
236.446181-238.763484|3069
238.763484-241.080787|2986
241.080787-243.398089|2902
243.398089-245.715392|2877
245.715392-248.032695|2821
248.032695-250.349998|2821
250.349998-252.667301|2776
252.667301-254.984604|2780
254.984604-257.301907|4391
257.301907-259.61921|4494
259.61921-261.936513|4518
261.936513-264.253816|4896
264.253816-266.571119|4755
266.571119-268.888422|4839
268.888422-271.205725|5079
271.205725-273.523028|5236
273.523028-275.840331|5239
275.840331-278.157634|5419
278.157634-280.474937|5529
280.474937-282.792239|5701
282.792239-285.109542|6009
285.109542-287.426845|5866
287.426845-289.744148|5551
289.744148-292.061451|5330
292.061451-294.378754|5067
294.378754-296.696057|4896
296.696057-299.01336|4891
299.01336-301.330663|4745
301.330663-303.647966|4437
303.647966-305.965269|4277
305.965269-308.282572|4181
308.282572-310.599875|3766
310.599875-312.917178|3623
312.917178-315.234481|3387
315.234481-317.551784|3030
317.551784-319.869086|2735
319.869086-322.186389|2591
322.186389-324.503692|2509
324.503692-326.820995|2202
326.820995-329.138298|2122
329.138298-331.455601|1987
331.455601-333.772904|1833
333.772904-336.090207|1648
336.090207-338.40751|1581
338.40751-340.724813|1430
340.724813-343.042116|1298
343.042116-345.359419|1198
345.359419-347.676722|1099
347.676722-349.994025|1019
349.994025-352.311328|955
352.311328-354.628631|876
354.628631-356.945933|836
356.945933-359.263236|731
359.263236-361.580539|672
361.580539-363.897842|622
363.897842-366.215145|594
366.215145-368.532448|503
368.532448-370.849751|502
370.849751-373.167054|430
373.167054-375.484357|409
375.484357-377.80166|359
377.80166-380.118963|299
380.118963-382.436266|305
382.436266-384.753569|264
384.753569-387.070872|224
387.070872-389.388175|258
389.388175-391.705478|195
391.705478-394.02278|176
394.02278-396.340083|155
396.340083-398.657386|138
398.657386-400.974689|136
400.974689-403.291992|142
403.291992-405.609295|108
405.609295-407.926598|113
407.926598-410.243901|97
410.243901-412.561204|71
412.561204-414.878507|76
414.878507-417.19581|79
417.19581-419.513113|72
419.513113-421.830416|56
421.830416-424.147719|53
424.147719-426.465022|44
426.465022-428.782325|48
428.782325-431.099627|44
431.099627-433.41693|42
433.41693-435.734233|39
435.734233-438.051536|29
438.051536-440.368839|32
440.368839-442.686142|21
442.686142-445.003445|25
445.003445-447.320748|18
447.320748-449.638051|20
449.638051-451.955354|18
451.955354-454.272657|16
454.272657-456.58996|17
456.58996-458.907263|11
458.907263-461.224566|9
461.224566-463.541869|11
463.541869-465.859172|16
465.859172-468.176475|7
468.176475-470.493777|11
470.493777-472.81108|11
472.81108-475.128383|8
475.128383-477.445686|8
477.445686-479.762989|3
479.762989-482.080292|3
482.080292-484.397595|7
484.397595-486.714898|2
486.714898-489.032201|1
489.032201-491.349504|5
491.349504-493.666807|3
493.666807-495.98411|2
495.98411-498.301413|4
498.301413-500.618716|3
500.618716-502.936019|1
502.936019-505.253322|6
505.253322-507.570624|1
507.570624-509.887927|3
509.887927-512.20523|3
512.20523-514.522533|4
514.522533-516.839836|3
516.839836-519.157139|3
519.157139-521.474442|2
521.474442-523.791745|2
523.791745-526.109048|1
526.109048-528.426351|2
528.426351-530.743654|2
533.060957-535.37826|2
535.37826-537.695563|2
544.647471-546.964774|1
546.964774-549.282077|1
549.282077-551.59938|2
553.916683-556.233986|1
558.551289-560.868592|2
563.185895-565.503198|1
567.820501-570.137804|2
572.455107-574.77241|1
574.77241-577.089713|1
588.676227-590.99353|1
*|21088357
Ну и напоследок, математиком не являюсь ни разу, поэтому мог где-то что-то не так понять. Кто-нибудь вообще делает эту самую топографическую коррекцию? Как боретесь с разбросом значений?

alexandr cherepanov
Гуру
Сообщения: 534
Зарегистрирован: 30 ноя 2006, 13:31
Репутация: 116
Откуда: Moscow

Re: Проблема с модулем i.topo.corr

Сообщение alexandr cherepanov » 17 фев 2015, 11:31

i.topo.corr - Computes topographic correction of reflectance.
Какие 255 значений? Reflectance, то посчитан?

Vedmed
Участник
Сообщения: 70
Зарегистрирован: 19 янв 2012, 14:48
Репутация: 6

Re: Проблема с модулем i.topo.corr

Сообщение Vedmed » 17 фев 2015, 13:57

alexandr cherepanov писал(а):
i.topo.corr - Computes topographic correction of reflectance.
Какие 255 значений? Reflectance, то посчитан?
Эээ, какой reflectance имеется ввиду? Top-of-atmospheric или Surface Reflectance? Если surface и это то, что получается в результате атмосферной коррекции через i.landsat.toar и i.atcorr, то да, посчитан. Собственно на этом снимке и пытаюсь победить тени. Если что-то понимаю неправильно, поясните, с дешифровкой ДДЗ плотно разбираться начал недавно, могу что-то упустить или недопонять.
Последний раз редактировалось Vedmed 17 фев 2015, 14:19, всего редактировалось 1 раз.

alexandr cherepanov
Гуру
Сообщения: 534
Зарегистрирован: 30 ноя 2006, 13:31
Репутация: 116
Откуда: Moscow

Re: Проблема с модулем i.topo.corr

Сообщение alexandr cherepanov » 17 фев 2015, 14:17

Это смотря для чего Вы делаете коррекцию. В принципе чтобы просто посчитать хватит и TOA.

Vedmed
Участник
Сообщения: 70
Зарегистрирован: 19 янв 2012, 14:48
Репутация: 6

Re: Проблема с модулем i.topo.corr

Сообщение Vedmed » 17 фев 2015, 14:21

alexandr cherepanov писал(а):Это смотря для чего Вы делаете коррекцию. В принципе чтобы просто посчитать хватит и TOA.
Не, мне нужно именно на атмосферноскорректированном снимке тени выправить. Хочу со всей логикой и последовательностью этого дела разобраться, чтобы иметь возможность сопоставлять результаты дешифровки разновременных космоснимков.

Ответить

Вернуться в «GRASS»

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

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 4 гостя