Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
.
meta.txtGoogle переводчиком с китайского на english перевел, там вроде эта формула упоминается как старая модель, затем приводятся мои формулы как новая модель на замену. При этом про коэффициенты a не слова. Попробую поискать еще инфы, если будут результаты - отпишусь. Хотя с примерами туго, на китайском это что-ли вбить...gamm писал(а):только там немного другое, Ref=(a2×EV^2+a1×EV+a0)×d^2/100/cos(SolZ) и т.д.
Этот софт мне не нужен. Смотрите стартовое сообщение темы. Я работаю с l1a данными, а указанный софт преобразует бортовые данные в l0 и l1a, но l1a я могу скачать прямо с их сайта (по квотам, но все же).Pilot писал(а):у китайцев есть два модуля для первичной обработки данных FY-3 MPT/HRPT Data Preprocessing software packages
FY3L0pp 1.1 version, FY3L1pp 1.1 version http://www.nsmc.cma.gov.cn/en/NSMC/Cont ... ils01.html
У меня, к сожалению, нет денег, чтобы покупать сейчас проприетарное ПО. Да и наверняка их модуль тоже требует какой-нибудь проги вроде ENVI. Пиратским софтом мне пользоваться никак нельзя (разве что как вспомогательным средством, не более, которое потом откинем).Pilot писал(а):Также можете обратиться в Сканэкс, у них есть модуль обработки данных FY-3 от raw до l1b уровня.
Там идут 14 чисел в строчку, я так понимаю там gain и bias вперемежку, они там чередуются? СПС за помощь, попробую в ближайшее время, отпишусь о результатах. Значения reflectance >1 это я так понимаю "сбойные" пиксели?pendduduk писал(а):Добрый день,
Если Вам нужно получить значения отражательной способности на кромке атмосферы (TOA Reflectance), то:
1. Читайте из HDF массив с коэффициентами (множитель и сдвиг, подряд по каждому каналу) "RefSB_Cal_Coefficients" - для VIRR
2. Читаете массив с зенитным углом Солнца "SolarZenith".
3. Приводите массив углов к градусам применив slope и intercept.
4. Переводите массив углов в радианы.
5. Читаете массив DN яркости и для каждого пикселя выполняете Pn=(DNn*gain[n]+DNn*bias[n])/cos(Z). Z - зенитный угол в радианах.
значения будут примерно 0-3, в среднем 0.1-0.4
gain,bias,gain,bias...pendduduk писал(а):Там идут 14 чисел в строчку, я так понимаю там gain и bias вперемежку, они там чередуются?
Не факт, может быть и блик.pendduduk писал(а):Значения reflectance >1 это я так понимаю "сбойные" пиксели?
Честно говоря не смотрел, но подозреваю что все по аналогии.pendduduk писал(а):А по mersi что-нибудь знаете?
На счет DNn - это значение яркости для n-ого канала? Кроме того, в VIRR нет DN, есть EV_RefSB, я так понимаю, нужно использовать его, да? У меня получились значения в диапазоне 0-18000, среднее 20-40,pendduduk писал(а):Добрый день,
Если Вам нужно получить значения отражательной способности на кромке атмосферы (TOA Reflectance), то:
1. Читайте из HDF массив с коэффициентами (множитель и сдвиг, подряд по каждому каналу) "RefSB_Cal_Coefficients" - для VIRR
2. Читаете массив с зенитным углом Солнца "SolarZenith".
3. Приводите массив углов к градусам применив slope и intercept.
4. Переводите массив углов в радианы.
5. Читаете массив DN яркости и для каждого пикселя выполняете Pn=(DNn*gain[n]+DNn*bias[n])/cos(Z). Z - зенитный угол в радианах.
значения будут примерно 0-3, в среднем 0.1-0.4
Да конечно, нужно использовать яркость каналаcmath10 писал(а):На счет DNn - это значение яркости для n-ого канала? Кроме того, в VIRR нет DN, есть EV_RefSB, я так понимаю, нужно использовать его, да? У меня получились значения в диапазоне 0-18000, среднее 20-40,
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость
© GIS-Lab и авторы, 2002-2017. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов (подробнее).