Вопрос по освещенности склонов

[MISTIK]

Добрый вечер!

Передо мной стоит задача расчета NDVI в предгорной местности. На имеющихся снимках бросается в глаза разница освещенности склонов. Нашел единственную методу по выравниванию контрастности: www.giscraft.ru/methods/method9.shtml. Но для ее выполнения не могу найти данные снимков: Solar Asimuth и Solar Elevation. Подскажите, пожалуйста где их найти? Снимки РапидАй. Единственные показатели со словом азимут, которые я нашел в метаданных это: <opt:illuminationAzimuthAngle uom="deg">1.800516e+02</opt:illuminationAzimuthAngle><opt:illuminationElevationAngle uom="deg">6.739295e+01</opt:illuminationElevationAngle><re:azimuthAngle uom="deg">1.023800e+02</re:azimuthAngle>. Что это за данные?

2 вопрос: нельзя ли компенсировать разницу освещённости склонов в ENVI? Читал здесь что можно создать маски для склонов различной освещенности, это реальная замена вышеуказанному методу?

Если где-то написал бред - не судите строго,я в этом деле новичок...

[ericsson]

Если дата и время съемки известны, посчитать положение солнца - не проблема совсем.

[MISTIK]

Спасибо! Не могли бы Вы описать поподробней или дать ссылку?
С датой легко. А время думаю взять из этой строки: <gml:validTime><gml:TimePeriod><gml:beginPosition>2010-05-27T09:16:43Z<
gml:beginPosition><gml:endPosition>2010-05-27T09:17:13Z</gml:endPosition>

Получается 9:17 среднее. Я же думаю там секунды не так важны?

[Дмитрий Барышников]

Z говорит что это по гринвичу. Вы в каком часовом поясе? По идее должно быть 11-12 ч. Т.к. обычно орбита солнечно-синхронная.
Я думаю:
opt:illuminationAzimuthAngle - угол Солнца
opt:illuminationElevationAngle - высота Солнца
re:azimuthAngle - угол съемки

[ericsson]

http://www.esrl.noaa.gov/gmd/grad/solcalc/azel.html можно проверить совпадение

[MISTIK]

Снимок Кавказа. Значит время снимка 12:17. Посчитал на калькуляторе, вот что получилось. Правда не понимаю,почему он автоматически при выборе Москвы в графе Offset to UTC ставит -3, ведь по идее должно быть +3?


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

[Александр Мурый]

но в конце выдает вот такую ошибку:
Кто нибудь сталкивался с такой?

Картинка не видна, лучше прикрепить её к сообщению.

[MISTIK]

Картинку вставил. Но какие координаты вбивать при расчете солнечного азимута? Их взять где-то из метаданных?

[Saimon]

Какие данные использовались в качестве исходных??
Одна из возможных причин в том что:

Calculations do not work if angular units (x, y; Lat/Lon) are mixed with distance units (z; meters).

[MISTIK]

Хм. Кажется додумался... Проекция ЦМР возможно не UTM (как у снимка). Правда проверить смогу не ранее чем через неделю. Как думаете,ошибка может быть из-за этого?

И какие же все таки координаты вбивать при расчете азимута солнца? Координаты спутника? Они есть в метаданных?

[Saimon]

Ошибка может быть из-за чего угодно, хорошо бы еще увидеть лог из ERDAS и сгенерированный скрипт.
При расчете азимута солнца, ЕМНИП, указываются координаты центра сцены...

[MISTIK]

Добрый вечер!

Ошибку теперь Erdas не выдает. Видимо что-то было с имеющийся ЦМР. Скачал SRTM,но теперь после выполнения нормализации получается вот такая картина:
От снимка видны лишь маленькие фрагменты,все остальное залито черным. Подскажите, пожалуйста из-за чего такое может происходить? В качестве исходных данных использовались: фрагмент снимка РапидАй на нужную территорию, скачанная модель SRTM (изменил проекцию на UTM), данные о положении Солнца взяты из метаданных (Bishop был прав где их найти). Все данные в формате TIFF. И возможно ли вообще применение ЦМР с разрешением 90 м со снимком с 5м разрешением?