Страница 2 из 4
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 30 сен 2012, 17:16
karapuz2009
Всем доброго времени суток. После запуска алгоритма выводиться прогресс выполнения, после чего программа завершается. В папке которую указываю для промежуточных результатов появляются файлы каналов с 2 по 6 и формируется файл маски, но получаемые спектральные отражательные способности получаются с отрицательными значениями, а файл маски пустой. Входящие данные взяты с
http://glovis.usgs.gov/, папки названы латиницей, пробелов и недопустимых символов не имеют. В чем может быть проблема?
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 01 окт 2012, 13:08
KolesovDmitry
karapuz2009 писал(а):В чем может быть проблема?
Пока не понятно, нужно смотреть. Можете сбросить ссылку на номер сцены?
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 01 окт 2012, 13:26
karapuz2009
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 01 окт 2012, 15:50
sergei-kzn
Добрый вечер,спасибо за отзыв, сегодня вечером скачаю снимки, посмотрю. в чем проблема пока не могу сказать.
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 01 окт 2012, 16:09
karapuz2009
Может быть проблема в разделителях?
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 01 окт 2012, 16:20
sergei-kzn
каких разделителях?
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 01 окт 2012, 16:22
karapuz2009
В файле метаданных, вместо точки запятая или наоборот.
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 01 окт 2012, 16:23
sergei-kzn
мм, они разве не везде одинаковые? я как-то этот вопрос не пробивал
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 01 окт 2012, 16:27
karapuz2009
Просто откуда могут взяться отрицательные значения? Может при расчете cos(x) при определении отражательной способности?
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 01 окт 2012, 16:29
sergei-kzn
отрицательные получаются все значения? или только часть?
зы я только через пару часов доберусь до кода, и смогу посмотреть что к чему.
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 01 окт 2012, 16:30
karapuz2009
Тепловой канал норма.
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 01 окт 2012, 16:31
sergei-kzn
хорошо, постараюсь разобраться что к чему.
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 01 окт 2012, 20:29
sergei-kzn
мде, снимки не смогу скачать, у них сайт "на обслуживании"
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 02 окт 2012, 11:37
karapuz2009
Сайт уже открыт для скачивания.
Re: Автоматическая оценка облачности по методике ACCA в QGIS
Добавлено: 02 окт 2012, 11:47
KolesovDmitry
karapuz2009 писал(а):После запуска алгоритма выводиться прогресс выполнения, после чего программа завершается. В папке которую указываю для промежуточных результатов появляются файлы каналов с 2 по 6 и формируется файл маски, но получаемые спектральные отражательные способности получаются с отрицательными значениями, а файл маски пустой.
Частично подтверждаю проблему. На моей системе (qgis 1.7.4, AccaPlugin 0.2, Ubuntu) алгоритм отрабатывает, файл маски создается. НО! При этом действительно в промежуточных файлах, создаваемых модулем toar, появляются отрицательные значения. Как следствие и в файле маски возникают сдвиги (так, например, облака были помечены как тени облаков).
Причина появления отрицательных значений ясна: в метаданных приводятся отрицательные Lmin:
Код: Выделить всё
LMIN_BAND1 = -1.520
LMIN_BAND2 = -2.840
LMIN_BAND3 = -1.170
LMIN_BAND4 = -1.510
LMIN_BAND5 = -0.370
LMIN_BAND6 = 1.238
LMIN_BAND7 = -0.150
и, как следствие, после расчета top of atmosphere radiances (как это по-русски будет, кто мне скажет?) для довольно большой части пикселей получается отрицательная.
Проблема известная,
здесь предлагается такое решение:
Corrected at-sensor values (method=corrected)
At-sensor reflectance values range from zero to one, whereas at-sensor radiance must be greater or equal to zero. However, since Lmin can be a negative number then the at-sensor values can also be negative. To avoid these possible negative values and set the minimum possible values at-sensor to zero, this method corrects the radiance to output a corrected at-sensor values using the equations (not for thermal bands):
* radiance = (uncorrected_radiance - Lmin)
* reflectance = radiance / sun_radiance
Вероятно, нужно производить подобную коррекцию и в AccaPlugin (наверное, по-умолчанию?).
Остается вопрос, почему у вас не был создан файл маски. Какая у вас версия ОС, qgis и плагина?