Автоматическая оценка облачности по методике ACCA в QGIS

[karapuz2009]

Всем доброго времени суток. После запуска алгоритма выводиться прогресс выполнения, после чего программа завершается. В папке которую указываю для промежуточных результатов появляются файлы каналов с 2 по 6 и формируется файл маски, но получаемые спектральные отражательные способности получаются с отрицательными значениями, а файл маски пустой. Входящие данные взяты с http://glovis.usgs.gov/, папки названы латиницей, пробелов и недопустимых символов не имеют. В чем может быть проблема?

[KolesovDmitry]

В чем может быть проблема?

Пока не понятно, нужно смотреть. Можете сбросить ссылку на номер сцены?

[karapuz2009]

ID: LT51730182007226MOR00
http://earthexplorer.usgs.gov/download/ ... 00?node=GV

[sergei-kzn]

Добрый вечер,спасибо за отзыв, сегодня вечером скачаю снимки, посмотрю. в чем проблема пока не могу сказать.

[karapuz2009]

Может быть проблема в разделителях?

[sergei-kzn]

каких разделителях?

[karapuz2009]

В файле метаданных, вместо точки запятая или наоборот.

[sergei-kzn]

мм, они разве не везде одинаковые? я как-то этот вопрос не пробивал

[karapuz2009]

Просто откуда могут взяться отрицательные значения? Может при расчете cos(x) при определении отражательной способности?

[sergei-kzn]

отрицательные получаются все значения? или только часть?
зы я только через пару часов доберусь до кода, и смогу посмотреть что к чему.

[karapuz2009]

Тепловой канал норма.

[sergei-kzn]

хорошо, постараюсь разобраться что к чему.

[sergei-kzn]

мде, снимки не смогу скачать, у них сайт "на обслуживании"

[karapuz2009]

Сайт уже открыт для скачивания.

[KolesovDmitry]

После запуска алгоритма выводиться прогресс выполнения, после чего программа завершается. В папке которую указываю для промежуточных результатов появляются файлы каналов с 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 и плагина?