Геополитика и экогеодинамика регионов. Том 8 (18). Вып. 2. 2022 г. С. 78-84.
УДК 574.42 А. Дж. Алиева
Вопросы оценки площади растительного покрова на дистанционно зондируемых участках_
Национальное аэрокосмическое агентство, г. Баку, Азербайджанская Республика e-mail: asadzade@rambler.ru
Аннотация. Сформулирована и решена задача оптимальной оценки частичного вегетационного покрова, сформированного на базе нормализованного вегетационного индекса, с учетом известного эффекта насыщения NDVI при высоких значениях LAI. Предложена осуществление коррекции только текущей величины NDVIT, и на этой основе сформирован новый показатель частичного растительного покрова состоящий из скорректированного NDVI. Также предложен интегральный вариант введенного показателя. Введена функция связи LAI = <pk(NDVIT) и при допущении постоянства определенного интеграла этой функции сформирована и решена оптимизационная задача поиска оптимальной функции <p(NDVIT)opt при которой вновь введенный интегральный показатель частичного покрова достигает максимальной величины. Определена взаимосвязь вновь предложенного и известного показателей частичного вегетационного покрова.
Ключевые слова: растительный покров, эффект насыщения, индекс листовой поверхности, нормализованный разностный растительный индекс, оптимизация
Введение
Чтобы охарактеризовать вертикальную проекцию части ландшафта занятой растительностью используется такой показатель как частичный растительный покров, введенный в работе [1]. Этот показатель является важным фенотипическим фактором, характеризующим пространственную характеристику типа растительности. Путем количественной оценки этого показателя может быть осуществлен мониторинг роста растительности [2], а также состояние выращенного урожая [3]. Частичный растительный покров может быть использован для оценки процессов, происходящих на земле [4,5,6,7], а также связи антропогенной деятельности с состоянием окружающей среды [8,9,10]. Известны различные алгоритмы для оценки частичного растительного покрова, такие как алгоритмы оценки относительного преобладания с использованием максимальных и минимальных значений вегетационных индексов [11], алгоритмы спектральной смеси [12], алгоритмы спектральной обучаемой классификации [13], алгоритмы на базе физических моделей [14] и др. Особое место среди вышеуказанных алгоритмов занимают первая группа алгоритмов, основанных на оценках вегетационных индексов. Несколько подробно рассмотрим эти алгоритмы.
Частичный вегетационный покров, вычисляемый на основе индекса листовой поверхности (LAI), определяемый как [15];
fCl = 1- exp(—kp • LAI) (1)
где: fc - обозначение частичного растительного покрова; kp - коэффициент ослабления;
Частичный вегетационный покров, вычисляемый на основе
нормализованного разностного растительного индекса [16]:
f _ NDVIT-NDVIS fC2 = NDVIX-NDVIS
(2)
Частичный вегетационный покров, вычисляемый квадратично по отношению f2 [17].
f _ (ndvit-ndvis\2
fC3 = \NDVI,X-NDVIS)
(3)
В формулах (2) и (3) приняты:
NDVIT - текущая величина этого индекса;
NDVIS - значение этого индекса для незаросшей почвы
NDVIm - максимально возможное значение этого индекса.
Вместе с тем, известный эффект насыщения NDVI при увеличение LAI привело к необходимости корректировки значения NDVI следующим образом [18].
NDVIK = NDVIm + (NDVIS — NDVIm) • exp[-kpLAl] (4)
С учетом (2) и (4) получаем, что корректируется не текущее значение, а разность NDVIS — NDVIm, т.е. достаточно широкий диапазон изменения NDVIT . Как нам представляется, имеет смысл корректировать не разность (NDVIS — NDVIm), а только текущую величину этого индекса, т.е. NDVIT.
Предлагаемый метод
С учетом вышесказанного предлагается новое определение частичного растительного покрова, по сути являющуюся скорректированной величиной NDVIT.
fC4 = NDVIt • fCi = NDVIt • [l — exp[—kp • LAI]] (5)
Введем на рассмотрение оптимизируемую функцию связи LAI и NDVIT, т.е.
LAI = <k(NDVIT) , (6)
Далее, рассмотрим следующий функционал, являющийся предлагаемым интегральным показателем частичного растительного покрова F1
р1 = CvlTr* NDVIt [l — exp[—kp<k(NDVIT)]] dNDVIT
0 (7)
В отношении функции <k(NDVIT) введем следующее ограничительное условие
CVVlTr <k(NDVIT) dNDVIT = С, С = const (8)
где: NDVIDmaT°= 1 NDVITo = 0. (9)
С - показатель, характеризующий ограничительное условие, налагаемое на искомый оптимальный вид функции (6).
С учетом (7) и (8) составим задачу оптимизации, заключающуюся по сути в выборе наилучшей функции связи в смысле достижения максимума интегрального показателя частичного растительного покрова F1
F± = NDVIT [1 - exp[-kp^k(NDVIT)]\ dNDVIT - X уk(NDVIT) dNDVIT -
С]
(10)
где: X - множитель Лагранжа.
Вычислим оптимальную функцию коррекции
LAI = Vk(NDVlT)opt Согласно [19], решение оптимизационной задачи (10) должно удовлетворить условию
d{NDVIT[l-exp[-kppk(NDVIT)]\-X-cpk(NDVIT)} _ ^
dpk(NDVIT) = ( )
Из выражения (11) получаем:
NDVIT • кр • exp[-kpyk(NDVIT)] -X = 0 (12)
Из выражения (12) находим:
ЛС^М^^Ц«] (13)
Для вычисления значения множителя Лагранжа достаточно вставить <Pk(NDVIT)opt под интеграл в выражении (8) и осуществив интегрирование определить величину X. Не вдаваясь математические детали этот операции, обозначим вычисленное значение множителя Лагранжа как X0.
Таким образом, решение оптимизационной задачи (10) имеет вид
<P,<mviT)oPt = 1rin[NDDVp] (14)
Легко проверить, что при решении (14) функционал (10) достигает максимума т.к. вторая производная интегранта в (10) по искомой функции оказывается отрицательной величиной. С учетом (1) и (14) выражения для fc приобретает следующий вид
Гс1 = 1-ехр^п[^\ = 1-< ь
Следовательно, вновь введенный показатель fcопределится как
fc = NDVIT • 1---—] = NDVIT - — (16)
Jc4 ' NDVIT-kp\ ' kp v '
На основе выражений (2) и (16) определим связь между показателями fC2 и Из (16) находим
NDVIt = fc+Yp
Р (17)
Из (2) находим
NDVIT = fC2 (NDVIm - NDVIS) + NDVIS (18)
Из выражений (17) и (18) получаем
fc4 + fc2 (NDVIs - NDVI J = -(гр + NDVIs) (19)
X^fc,+X2^fC2 = 1 Р (20)
NDVIT-kp
(15)
fc4 •
— 1
где ¿1 = к—т (21)
£°+NDVIs
Кр
= ndvi^NDVs (22)
2 tp+NDV>S ( )
Как видно из выражений (21) и (22), Л1 и X2 разнополярны, а следовательно и разность X1fC и X2fC2 равно единице.
Аналогичным образом можно вычислить взаимосвязь fC и fC3 используя выражения (16) и (3).
Таким образом, предлагаемый алгоритм максимального учета эффекта насыщения N D V в предложенном показателе частичного вегетационного покрова, основанного на N DV заключается в следующем:
Коррекция текущей величины NDVIT по выражению (5). Введение функции связи LAI и NDVIk
LAI = <k(NDVIk) Поиск вида скорректированной функции связи при которой целевой функционал F1 достиг бы максимальной величины путем решения оптимизационной задачи (10). Определения максимального значения fCи на этой основе взаимосвязи fC и fC .
JC2 JC4max
Выводы
Таким образом, сформулирована и решена задача уточнения показателя частичного вегетационного покрова, сформированного на базе нормализованного вегетационного индекса, с учетом известного эффекта насыщения N DVI при высоких значениях LAI. В отличие от известного метода учета указанного насыщения предложена осуществление коррекции только текущей величины NDVI, а не разности (NDVIm — NDVI). На этой основе предложен новый показатель частичного растительного покрова состоящий из скорректированного NDVI. Введение функции связи LAI = <k(NDVIT) и допущение о постоянстве определенного интеграла этой функции позволило сформировать и решить оптимизационную задачу поиска оптимальной функции <(NDVI)0pt при которой вновь введенный интегральный показатель частичного покрова достигает максимальной величины. Полученные результаты при этом позволяют определить взаимосвязь вновь предложенного и известного показателей частичного вегетационного покрова.
Литература
1. Deardorff J. W. Efficient Prediction of Ground Surface Temperature and Moisture, with Inclusion of a Layer of Vegetation // Journal of Geophysical Research. 1978. № 83. рр.1889-1903.
2. Allen R. G., Pereira L. S. Estimating crop coefficients from fraction of ground cover and height // Irrigation Science, Springer Science+Business Media. 2009. Volume 28, Issue 1. pp. 17-34
3. Casa A. de la, Ovando G., Bressanini L., Martinez J., Diaz G., Miranda C. Soybean crop coverage estimation from NDVI images with different spatial
Ameea A. jlw.
resolution to evaluate yield variability in a plot. ISPRS Journal of Photogrammetry and Remote Sensing. 2018. Vol.146. pp.531-547
4. Foley J. A., Levis S., Costa M. H., Cramer W., Pollard D. Incorporating Dynamic Vegetation Cover within Global Climate Models // Ecological Applications. Vol. 10, №. 6. 2000. pp. 1620-1632.
5. Li F., Kustas W. P., Prueger J. H., Neale C. M. U., Jackson T. J. Utility of remote sensing-based two-source energy balance model under low- and high-vegetation cover conditions // Journal of Hydrometeorology. Vol. 6. 2005. pp. 878-891.
6. Wang G. X., Liu G. S., Li C. J., Yang Y. The variability of soil thermal and hydrological dynamics with vegetation cover in a permafrost region // Agricultural and Forest Meteorology. Vol. 162-163, 15. 2012. pp. 44-57.
7. Wei X. H., Li, Q., Zhang M. F., Giles-Hansen K., Liu W. F., Fan H. B., Wang Y., Zhou G. Y., Piao S. L., Liu S. R. Vegetation cover-another dominant factor in determining global water resources in forested regions // Global Change Biology. 2018. pp. 786-795.
8. Jiang M., Tian S., Zheng Z., Zhan Q., He Y. Human Activity Influences on Vegetation Cover Changes in Beijing, China, from 2000 to 2015. Remote Sens. 2017, 9(3), 271. https://doi.org/10.3390/rs9030271
9. Tong L., Zhou J., Li X., Qian Y. T., Gao Y. S. Region-Based Structure Preserving Nonnegative Matrix Factorization for Hyperspectral Unmixing // IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing. 2017. Vol. 10, Issue: 4. pp. 1575-1588.
10. Xin Z. B., Xu, J. X., Zheng W. Spatiotemporal variations of vegetation cover on the Chinese Loess Plateau (1981-2006): Impacts of climate changes and human activities. Science in China Series D, Earth Sciences. 2008. № 51(1), pp.67-78.
11. Gutman G., Ignatov A. Satellite-derived green vegetation fraction for the use in numerical weather prediction models // Advances in Space Research. 1997. Vol. 19, Issue 3. pp. 477-480.
12. Roberts D. A., Gardner M., Church R., Ustin S., Scheer G., Green R. O. Mapping Chaparral in the Santa Monica Mountains Using Multiple Endmember Spectral Mixture Models. 1998. Vol. 65, Issue 3. pp. 267-279 .
13. Friedl M. A., McIver D. K., Hodges J. C. F., Zhang X. Y. ; Muchoney D., Strahler A. H. ; Woodcock, C. E. ; Gopal, S. ; Schneider, A. ; Cooper, A. ; Baccini, A., Gao F., Schaaf C. Global land cover mapping from MODIS: algorithms and early results // Remote Sensing of Environment. 2002. vol. 83, issue 1-2. pp. 287-302.
14. Chopping M., Su L., Rango A., Martonchik J., Peters D., Laliberte A. Remote sensing of woody shrub cover in desert grasslands using MISR with a geometric-optical canopy reflectance model // Remote Sensing of Environment. 2008 vol. 112, issue 1. pp. 19-34.
15. Nilson T. A theoretical analysis of the frequency of gaps in plant stands // Agricultural Meteorology. 1971. № 8. pp. 25-38.
16. Qi J., Marsett R. C., Moran M. S., Goodrich D. C., et al. Spatial and temporal dynamics of vegetation in the San Pedro River basin area // Agricultural and Forest Meteorology. 2000. vol. 105, issue 1-3. pp. 55-68.
17. Carlson T. N., Ripley D. A. On the relation between NDVI, fractional vegetation cover, and leaf area index. Remote Sensing of Environment. 1997. vol. 62, issue 3. pp. 241-252 .
18. Baret F., Guyot G. Potentials and limits of vegetation indices for LAI and APAR assessment // Remote Sensing of Environment. 1991. vol. 35, issue 2-3. pp. 161-173.
19. Эльсгольц Л. Э. Дифференциальные уравнения на вариационное исчисление. М. Наука. 1974. 432 с.
A. J. Aliyeva
Issues of estimating the area of vegetation cover in remotely sensed areas_
National Aerospace Agency, Baku, Republic of Azerbaijan e-mail: asadzade@rambler.ru
Abstract. The problem of optimal estimation of partial vegetation cover formed on the basis of a normalized vegetation index is formulated and solved, taking into account the known effect of NDVI saturation at high LAI values. The correction of only the current value of [NDVI]_t is proposed, and on this basis a new indicator of partial vegetation cover consisting of adjusted NDVI is formed. An integral version of the introduced indicator is also proposed. The coupling function LAI=y_k ([NDVI]_t) is introduced and, assuming the constancy of a certain integral of this function, the optimization problem of finding the optimal function p([NDVI]_t )_opt is formed and solved at which the newly introduced integral indicator of partial cover reaches the maximum value. The interrelation of the newly proposed and known indicators of partial vegetation cover is determined.
Keywords: vegetation cover, saturation effect, leaf surface index, normalized difference vegetation index, optimization
References
1. Deardorff J. W. (1978) Efficient Prediction of Ground Surface Temperature and Moisture, with Inclusion of a Layer of Vegetation // Journal of Geophysical Research. 1978. № 83. pp.1889-1903. (in English)
2. Allen R. G., Pereira L. S. Estimating crop coefficients from fraction of ground cover and height // Irrigation Science, Springer Science+Business Media. 2009. Volume 28, Issue 1. pp. 17-34 (in English)
3. Casa A. de la, Ovando G., Bressanini L., Martinez J., Diaz G., Miranda C. Soybean crop coverage estimation from NDVI images with different spatial resolution to evaluate yield variability in a plot. ISPRS Journal of Photogrammetry and Remote Sensing. 2018. Vol.146. pp.531-547. (in English)
4. Foley J. A., Levis S., Costa M. H., Cramer W., Pollard D. Incorporating Dynamic Vegetation Cover within Global Climate Models // Ecological Applications. Vol. 10, №. 6. 2000. pp. 1620-1632. (in English)
5. Li F., Kustas W. P., Prueger J. H., Neale C. M. U., Jackson T. J. Utility of remote sensing-based two-source energy balance model under low- and high-vegetation cover conditions // Journal of Hydrometeorology. Vol. 6. 2005. pp. 878-891. (in English)
6. Wang G. X., Liu G. S., Li C. J., Yang Y. The variability of soil thermal and hydrological dynamics with vegetation cover in a permafrost region // Agricultural and Forest Meteorology. Vol. 162-163, 15. 2012. pp. 44-57. (in English)
7. Wei X. H., Li, Q., Zhang M. F., Giles-Hansen K., Liu W. F., Fan H. B., Wang Y., Zhou G. Y., Piao S. L., Liu S. R. Vegetation cover-another dominant factor in determining global water resources in forested regions // Global Change Biology. 2018. pp. 786-795. (in English)
8. Jiang M., Tian S., Zheng Z., Zhan Q., He Y. Human Activity Influences on Vegetation Cover Changes in Beijing, China, from 2000 to 2015. Remote Sens. 2017, 9(3), 271. https://doi.org/10.3390/rs9030271. (in English)
9. Tong L., Zhou J., Li X., Qian Y. T., Gao Y. S. Region-Based Structure Preserving Nonnegative Matrix Factorization for Hyperspectral Unmixing // IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing. 2017. Vol. 10, Issue: 4. pp. 1575-1588. (in English)
10. Xin Z. B., Xu, J. X., Zheng W. Spatiotemporal variations of vegetation cover on the Chinese Loess Plateau (1981-2006): Impacts of climate changes and human activities. Science in China Series D, Earth Sciences. 2008. № 51(1), pp.67-78. (in English)
11. Gutman G., Ignatov A. Satellite-derived green vegetation fraction for the use in numerical weather prediction models // Advances in Space Research. 1997. Vol. 19, Issue 3. pp. 477-480. (in English)
12. Roberts D. A., Gardner M., Church R., Ustin S., Scheer G., Green R. O. Mapping Chaparral in the Santa Monica Mountains Using Multiple Endmember Spectral Mixture Models. 1998. Vol. 65, Issue 3. pp. 267-279. (in English)
13. Friedl M. A., McIver D. K., Hodges J. C. F., Zhang X. Y. ; Muchoney D., Strahler A. H. ; Woodcock, C. E. ; Gopal, S. ; Schneider, A. ; Cooper, A. ; Baccini, A., Gao F., Schaaf C. Global land cover mapping from MODIS: algorithms and early results // Remote Sensing of Environment. 2002. vol. 83, issue 1-2. pp. 287-302. (in English)
14. Chopping M., Su L., Rango A., Martonchik J., Peters D., Laliberte A. Remote sensing of woody shrub cover in desert grasslands using MISR with a geometric-optical canopy reflectance model // Remote Sensing of Environment. 2008 vol. 112, issue 1. pp. 19-34. (in English)
15. Nilson T. A theoretical analysis of the frequency of gaps in plant stands // Agricultural Meteorology. 1971. № 8. pp. 25-38. (in English)
16. Qi J., Marsett R. C., Moran M. S., Goodrich D. C., et al. Spatial and temporal dynamics of vegetation in the San Pedro River basin area // Agricultural and Forest Meteorology. 2000. vol. 105, issue 1-3. pp. 55-68. (in English)
17. Carlson T. N., Ripley D. A. On the relation between NDVI, fractional vegetation cover, and leaf area index. Remote Sensing of Environment. 1997. vol. 62, issue 3. pp. 241-252 . (in English)
18. Baret F., Guyot G. Potentials and limits of vegetation indices for LAI and APAR assessment // Remote Sensing of Environment. 1991. vol. 35, issue 2-3. pp. 161-173. (in English)
19. El'sgol'c L. E. Differencial'nye uravneniya na variacionnoe ischislenie. M. Nauka. 1974. 432 s. (in Russian)
nocmynuna epeda^uw 01.04.2022 г.