Форумы GIS-Lab.info

Считали NDSI в 4 и 6 канале в Erdas по продукту mod09, значения получаются больше 1(252.., ) тип данных выходного файла устанавливали float single. Для уровня 1B использование тех же каналов дает корректные значения(от -1 до 1 ).

не подскажете в чем проблема?

Вы для расчетов используете модель в Modeller? Правильно ли указывается исходный тип? Ведь значения MOD09 - мастабированные.
Приведите пример с реальными данными (пиксельными значениями), что должно и что получается.

А не подскажете, как каналы MOD09L2 (например MOD09_L2.2008002.0827.hdf) преобразовать в немасштабированые в Erdas? Во всех руководствах для вычесления NDSI используется продукт MOD09GHK (может там значания немасштабированы?), у нас его нету. MOD09L2 получили выполнив алгоритм ATBD08, который обсуждался в теме "Калибровочные данные для TERRA". Далее вырезаем каналы 500m Surface Reflectance Band 4 и 6, и исспользуем стандартную модель NDSI в Erdas.

Почитали руководство по mod09
http://modis-sr.ltdri.org/MAIN_SURFACE_ ... rGuide.pdf
там указаны для продукта MOD09GHK допустимые значения (Valid Range) -100 до 16000.
выполнили в модели соттветствующую фильтацию значений данных каналов. результат не изменился. выполнили фильтрацию от 0 до 16000 - получились допустимые значения (-1;1)
прибавление +100 ко всем значениям тоже "сработало" в смысле диапазона, но значения немного другие и это ведь не корректно.
Исходный тип signed 16 bit. учет Scale Factor (divide by) результат не меняло(чтиобы прийти не к масштабированым надо ведь делить получается на значение Scale Factor?)

и еще вопроc

В руководстве по расчету NDSI http://www.ucalgary.ca/bgs/files/bgs/po ... o_2006.pdf
написано
band 2 > 11%
band 4>=10%
в тексте, как поняли, это условие требуется для отделения "темных" пикселей.
вопрос 11% - это от чего?
есть предположение: исходный тип signed integer 16 bit - соотвтственно значения от -32768 до +32768. мы взяли 10% - 3276.8. т.е. условие band 2 > 3276.8.


еще на стр. 5 документа есть выражение SCA= a exp(-b*ADD) (1) это имеет отношение к тому что можно использовать 31 и 32 канал для отделения тех участков где температура выше 277 K(283 for Collection 4)
http://www-modis.bu.edu/brdf/Hall_Boston_Outreach.pdf
и как получить значение температуры их тех значений что хранятся в этих каналах?