Атмосферная коррекция i.atcorr
Добавлено: 05 ноя 2014, 17:04
Провожу атмосферную коррекцию данных Landsat 8 OLI (каналы 1-7), GRASS 7.0 beta 3, ОС Windows 7 64 bit.
Растры каналов предварительно переведены в radiance, тип 32float.
Вопросы:
1) нужно ли корректировать input range (и если нужно, то до каких значений), который по умолчанию [0,255]
2) похоже, что в справке неверно указаны параметры файла атмосферной коррекции для Landsat 8, а именно
определила это при чтении результата команды - там пишется канал при считывании файла с параметрами. Проблема возникает при коррекции 7 канала с номером 124 - информация о канале интерпретируется верно, но модуль выдает предупреждение и работа прекращается
3) при коррекции канала 5 все значения становятся равными 1 - никак не могу понять в чем дело (geotiff во вложении).
В общем и целом, изменений в композитах на основе атмосферно скорректированных данных по сравнению с DN/radiance не заметила, хотя по идее они же должны стать ярче/контрастнее?
Растры каналов предварительно переведены в radiance, тип 32float.
Вопросы:
1) нужно ли корректировать input range (и если нужно, то до каких значений), который по умолчанию [0,255]
2) похоже, что в справке неверно указаны параметры файла атмосферной коррекции для Landsat 8, а именно
Спойлер
17 Landsat 8, а нужно 18 Landsat 8
Спойлер
115 Landsat 8 Coastal Aerosol Band (0.427nm - 0.459nm), а нужно 116
116 Landsat 8 Blue Band (436nm - 527nm), нужно 117
117 Landsat 8 Green Band (512nm-610nm), нужно 118
118 Landsat 8 Red Band (625nm-691nm), нужно 119
119 Landsat 8 Panchromatic Band (488nm-692nm), нужно 120
120 Landsat 8 NIR Band (829nm-900nm), нужно 121
121 Landsat 8 Cirrus Band (1340nm-1409nm), нужно 122
122 Landsat 8 SWIR1 Band (1515nm - 1697nm), нужно 123
123 Landsat 8 SWIR2 Band (2037nm - 2355nm), нужно 124
116 Landsat 8 Blue Band (436nm - 527nm), нужно 117
117 Landsat 8 Green Band (512nm-610nm), нужно 118
118 Landsat 8 Red Band (625nm-691nm), нужно 119
119 Landsat 8 Panchromatic Band (488nm-692nm), нужно 120
120 Landsat 8 NIR Band (829nm-900nm), нужно 121
121 Landsat 8 Cirrus Band (1340nm-1409nm), нужно 122
122 Landsat 8 SWIR1 Band (1515nm - 1697nm), нужно 123
123 Landsat 8 SWIR2 Band (2037nm - 2355nm), нужно 124
Спойлер
i.atcorr --overwrite --verbose input=LC8_2014_162_band7_radiance@Landsat elevation=srtm_bt_aoi_utm35@Landsat parameters=G:\grassdata/bt_transition_potential/Landsat/.tmp/unknown/5080.8 output=LC8_2014_162_band6_atcor rescale=0,1
wavelength less than 0.25 micron:
let's take s(l)=s(0.25)
WARNING: Unsupported iwave value: 124
wavelength less than 0.25 micron:
let's take s(l)=s(0.25)
WARNING: Unsupported iwave value: 124
Спойлер
18 - geometrical conditions=Landsat 8
6 11 9.14 24.345 48.234 - month day hh.ddd longitude latitude
2 - atmospheric mode=midaltitude summer
1 - aerosols model=continental
0 - visibility [km] (aerosol model concentration)
0.022 - aerosol optical depth 550 nm (from MODIS data)
-1.067 - mean target elevation above sea level [km]
-1000 - sensor height (here, sensor on board a satellite)
121 - 5th band of OLI Landsat 8
6 11 9.14 24.345 48.234 - month day hh.ddd longitude latitude
2 - atmospheric mode=midaltitude summer
1 - aerosols model=continental
0 - visibility [km] (aerosol model concentration)
0.022 - aerosol optical depth 550 nm (from MODIS data)
-1.067 - mean target elevation above sea level [km]
-1000 - sensor height (here, sensor on board a satellite)
121 - 5th band of OLI Landsat 8