AMMI Automatic Mangrove Map and Index: Novelty for Efficiently Monitoring Mangrove Changes with the Case Study in Musi Delta, South Sumatra, Indonesia

,


Introduction
Mangroves are woody plants that grow inland and at the marine boundary (coastal), especially in tropical and subtropical areas [1,2]. Mangrove ecosystems are characterized by high environmental dynamics, e.g., temperature, sedimentation, and tidal currents [3]. Mangroves are highly benefcial to coastal ecosystems and shallow waters due to their contribution to the coastal zone, productivity, and biologically essential ecosystems [4]. Te environment has specifc characteristics that are generally infuenced by freshwater from the land via rivers and saltwater from the sea [5,6]. Te mangrove ecosystem provides services such as nursery areas for many marine fsheries and nutrient cycling [7], habitat for wildlife species, the landing site for thousands of migratory birds [3], and biodiversity [8]. In the context of climate change, mangroves play an important role in carbon sequestration as they can sequester carbon in the atmosphere through photosynthetic processes, and most of them are also stored in the soil [9][10][11][12][13][14].
Mangroves are vulnerable to anthropogenic and natural disturbances. Te Food and Agriculture Organization of the United Nations (FAO) [15] reports that 20% of the world's mangroves have been lost to deforestation since the 1980s. Anthropogenic factors that play a signifcant role include urban expansion [16][17][18][19][20][21], functional change to aquaculture and shrimp farming [22], illegal exploitation for fuelwood and construction materials [23], or damage from natural factors such as storms, hurricanes, and tsunamis [24,25].
Te world's attention is now focused on climate change, where the physical manifestations of threats are often referred to as hazards or climate hazards [26] and rising sea levels signifcantly impact coastal environments. Most mangroves do not keep pace with sea level rise because sediments are not as high and a limited area is available for landward migration [27]. In tropical countries where mangrove forests thrive, Sierra-Correa and Kintz [28] emphasized that the long-term threat of sea level rise requires coastal planning to avoid much more signifcant losses.
Te Indonesian Minister of Environment and Forestry's opening speech at the International Conference on Sustainable Mangrove Ecosystems in BAli on April 18, 2017 [29] states that Indonesia has a mangrove ecosystem area of 3.5 million hectares. Te government manages 2.2 million hectares, and communities manage the remaining 1.3 million hectares, located in 257 districts/cities, most of which are degraded.
For decades, remote sensing via satellite imagery has been widely used to monitor the condition of mangroves. Its ability to cover broader areas and its temporally make it ideal from time to time [30][31][32], and more than 1300 scientifc remote sensing papers on mangroves have been published on various topics [33]. Mapping mangroves with multispectral, medium-resolution images is the most popular data source. Landsat 8, launched on February 11, 2013, has 11 bands, with the spatial resolution of the panchromatic band being 15 m. Sentinel-2A was launched on June 23, 2015, and Sentinel-2B was launched on March 7, 2017. Both provide 13 multispectral bands with a spatial resolution of 10 m for vegetation detection. Compared to previous Landsat satellites, the beam resolution was increased to 16 bits, and the signal-to-noise ratio was signifcantly improved. Tese advances improved Landsat 8's ability to discriminate vegetation [34]. Tese satellite images are easy to collect and open to access, and free software is now available to support image processing. However, the problem is that the results of mangrove research using remote sensing have stagnated at only the journal level for academics and researchers. It is difcult for interested practitioners and local communities to apply the results.
Tis research aims to develop a simple formula to automatically trace, capture, and map mangroves and inform the index based on satellite imagery. Te goal is to monitor, assess, and manage the mangrove condition, and this can be performed by anyone interested in mangroves, including the central government, local government, and local communities. With the requirements of the next 10-15 years, the mapping of mangroves will cover a larger area. Google Earth Engine (GEE) is a platform that provides multitemporal satellite image data archives. Te platform provides applications for various user-created formulas, allowing the possibility of using an algorithm to access mangroves [35,36]. Collecting satellite imagery through the GEE helps users remove clouds, is free from mosaic processes, and is not limited in its spatial coverage.
Technically, AMMI creates an efective and efcient mangrove map that eliminates traditional work such as manual digitization or other classical classifcation methods.
Scientifcally, automatic mangrove mapping can cover a much broader area more quickly. Terefore, this method will provide information not only on the presence and spatial extent of mangroves but also inform the relationship of mangroves to coastal geomorphology [37] and their ecological conditions [31], which can eventually be used to determine and assess the health status of mangroves.
Te case study is near the Musi River Delta in South Sumatra Province, Indonesia ( Figure 1). Te area is administratively divided into two parts. Te northwestern part is a small part of the Sembilang National Park (SNP) under the central government's control. Te other southeastern part is the Air Telang Protected Forest (ATPF) under the control of the South Sumatra provincial government.

Material and Methods
where, ρ ƛ : real refectance on the earth surface, Mρ: bandspecifc multiplicative rescaling factor, Q cal : Quantized and calibrated standard product pixel values (DN), Aρ: bandspecifc additive rescaling factor, and θ SUN : local sun elevation angle All the above parameters are stored in the MTL fle of each image packet. Te real refections as results of the processing of L and sat and Sentinel images are annotated as ρ Blue, ρ Green, ρ Red, ρ NIR , ρ SWIR1 , and ρ SWIR2 for Blue, Green, Red, NIR, SWIR1, and SWIR2, respectively. In Landsat 8 OLI, a panchromatic band increases the spatial resolution to 15 m. Te result of the fusion, due to the higher spatial resolution, makes the image visually clearer [38]. However, the color pattern still obtained images with a resolution of 30 m × 30 m [39]. A panchromatic band increases spatial resolution, and spatial sharpening will be visually helpful in interpreting the mangrove [34]. Tis work converts the digital number of images into surface refectance for radiometric correction using DOS1 in the semiautomatic classifcation plugin (SCP) module in QGIS vers. 3.4 Madeira.  Figure 1: Te case study is located in the Musi River delta, South Sumatra, Indonesia.  International Journal of Forestry Research maximum slope in the increase of refectance from ρ Red to ρ NIR has been a good indicator of the leaf level's chlorophyll content and at the canopy level [40,41]. Kuenzer et al. [31] and [42] found that the spectral response of mangroves in the wavelength range of 380-750nm (ρ Blue − ρ Red ) is feeble, while it is strong in the spectrum of 750-2500 nm (ρ NIR − ρ SWIR2 ), especially concerning to leaf structure and properties, water content, and mangrove biochemistry. An illustration of the spectral response in mangroves is shown in Figure 2(a). Te visible refectance is lower than in ρ NIR due to chlorophyll absorption and leaf cell wall scattering. Further explained, the refectance of ρ NIR is principally controlled by the walls of the spongy mesophyll cells, with healthier leaves tending to have more substantial refectance in ρ NIR as they refect excessive amounts of incoming energy from the electromagnetic spectrum. In contrast, stressed leaves will have lower refectance due to cell structure changes. Te leaf water content is the primary determinant of refectance in the ρ SWIR region.
From the reviews above, it is generally accepted as the basic principles of the optical-remotely sensed to detect vegetation: chlorophyll-based, sensitive in the ρ Red − ρ NIR wavelength range and water content-based in leaves which are sensitively in the ρ NIR − ρ SWIR1 wavelength range. Many formulas have been available in mangrove research, probably almost a hundred, as revealed by Xue and Su [43] and Kobayashi et al. [44]. Even though the revealed formulas cannot automatically delineate mangroves, manual digitization is still used in mangrove research to separate mangroves and other features. As is known, the work is very tedious, tiring, and time-consuming.
Te existing vegetation indices (VIs) can be broadly grouped into two main streams: (1) the chlorophyll-based content of the vegetation with all its variants, and (2) the water/moisture-based content with all its variants.

Te Existing Vegetation Indices as References
(1) Formula Based on Chlorophyll Content. Healthy vegetation based on chlorophyll content refects more ρ NIR and ρ Green light than other wavelengths and absorbs more red and blue light. As an indicator of health, chlorophyll strongly absorbs visible light, and the cell structure of leaves strongly refects ρ NIR . When the plant is dehydrated, diseased, or afected by a disease, the sponge layer deteriorates and the plant absorbs more near-infraRed light instead of refecting it. Tus, observing how ρ NIR changes compared to ρ Red provides an accurate indication of the presence of chlorophyll, which correlates with plant health. Stressed mangroves have greater refectance in visible light, particularly in the ρ Red regions, likely due to a decrease in chlorophyll content and an increase in carotenoids in the leaves [42].
Tere are formulas, the oldest and simplest of which, the ratio vegetation index (RVI) using ρ Red and ρ NIR , was formulated as (ρ Red /ρ NIR ) proposed by Jordan [45]. Te RVI is widely used for estimating and monitoring green biomass, especially when there is a dense vegetation cover. Tis index is sensitive to vegetation and correlates with plant biomass. However, when vegetation cover is sparse (less than 50%), the RVI is more sensitive to atmospheric efects, and its representation of biomass is weak [44].
Te most popular is the normalized diference vegetation index (NDVI) proposed by Rouse et al. [46], which has been the most common and widely used in mangrove remote sensing for more than two decades [41,47]. NDVI quantifes vegetation by measuring the diference between ρ NIR (which is strongly refected by vegetation) and ρ Red (which is absorbed by vegetation), formulated as (ρ NIR − ρ Red )/ (ρ NIR + ρ Red ). NDVI was initially developed to monitor plant growth in plantation environments. However, this formula is adopted and applied in mangrove research. It is estimated that this formula has been used and applied in hundreds of mangrove research papers.
Several derivatives of NDVI have also been proposed to address the limitations, including the perpendicular vegetation index (PVI) [48], the soil-adjusted vegetation index (SAVI) [49], the atmospherically resistant vegetation index (ARVI) [50], and the global environment monitoring index (GEMI) [51]. Tese attempted to incorporate intrinsic corrections for one or more confounding factors. Several new generation algorithms are proposed for estimating biogeophysical variables to take advantage of modern sensors' improved performance and characteristics and eliminate confounding factors. Despite these factors, NDVI remains a valuable tool for quantitatively monitoring vegetation in terms of photosynthetic capacity at a spatial scale appropriate for various phenomena.
(2) Formula Based on Water/Moisture Content. Te Infrared Index (II) was proposed by Hardisky et al. [52] and was perhaps the frst to propose a moisture-based index using ρ NIR and ρ SWIR1 , formulated as (ρ NIR − ρ SWIR1 )/ (ρ NIR + ρ SWIR1 ). Te results showed that the decrease of II in the canopy is correlated with the increasing salinity of soil in salt forests. Tis plot of water content shows a signifcant decrease in canopy moisture with increasing soil salinity. In summary, a combination of NDVI and II can detect morphological and physiological changes associated with moisture stress. However, using longer wavelengths is a more direct indicator of water content. Te normalized diference water index (NDWI) was proposed by Gao [47] based on research using moderate resolution imaging spectrometer (MODIS) and airborne visible infrared imaging spectrometer (AVIRIS) images to monitor vegetation changing based on the liquid water content in the canopy. Te NDWI is originally formulated as (ρ(0.86 μm) − ρ(1.24 μm))/(ρ(0.86 μm) + ρ(1.24 I μm)). Te vegetation index has been widely applied to Landsat images, especially in mangrove research, thus the index, universally, is formulated as (ρ NIR − ρ SWIR1 )/(ρ NIR + ρ SWIR1 ). Te liquid water absorption in ρ NIR is negligible and presents absorption at ρ SWIR . Vegetation canopy scattering enhances water absorption. As a result, NDWI is sensitive to changes in the liquid water content of vegetation canopies. Atmospheric aerosol scattering efects from the ρ NIR to ρ SWIR1 wavelength region are weak, so NDWI is less sensitive to the atmospheric efect than NDVI. However, NDWI does not completely remove the efect of soil background, but by increasing the vegetation fraction, the NDWI value increases. Tucker [40] frst suggested that the ρ SWIR1 wavelength was the best-suited band for monitoring the water status at the plant canopy from space.
Normalized diference moisture index (NDMI) [53] is formulated (ρ NIR − ρ SWIR1 )/(ρ NIR + ρ SWIR1 ). However, retaining the term moisture, there is no better term, but the point is related to the wetness that includes the water content of vegetation, water absorbance in the fresh leaf, and soil wetness that will afect the sensitivity to soil and plant moisture. Te diference between ρ NIR and ρ SWIR1 appears in the ability of ρ SWIR1 wavelengths to absorb water, so that the index value can be used to estimate water content in the vegetation [54]. In the green leaves, the ρ NIR band has more refectance than the other bands, and the reduction in ρ SWIR1 refectance compared to ρ NIR is due to water absorption. Wetness change is a good indicator and the single most consistent indicator of forest change, including lighter disturbance/partial cuts, because it captures changes in ρ SWIR1 .
Mangrove discrimination indices (MDI) [34] intended to separate mangrove and nonmangrove vegetation using ρ NIR and ρ SWIR1 or ρ SWIR2 and formulated as (ρ NIR − ρ SWIR1 )/ (ρ SWIR1 ). In his fndings, ρ SWIR1 or ρ SWIR2 can increase the diference between mangrove and non-mangrove vegetation. In the application, when using MDI1 (ρ NIR and ρ SWIR1 ) to separate between mangroves and other vegetation, it is not yet clear whether it is best to move to MDI2 with the replacement from ρ SWIR1 to ρ SWIR2 .
All existing vegetation indices (VIs) and the formulas proposed by previous researchers are presented in Table 3.

Spectral Response of Segara Anakan Mangrove
Vegetations as a Reference. Something is quite interesting in the mangrove environment of Segara Anakan, Cilacap, Central Java, Indonesia (Figure 1), reported by Winarso and Purwanto [55]. In the logged mangrove areas, shrubs such as Derris trifoliata and Acanthus illicifolius close the felled mangrove area so that the existing mangrove seedlings and saplings cannot develop. It is known that these shrubs are closely related and included in the mangrove association [56]. Using NDVI analysis, these shrubs show a very high index and even exceed the true mangrove due to strongly refecting the ρ NIR . Anyone who is not careful will be deceived, as if it is like dense mangroves.
Te authors identifed spectral responses in seven different land cover types surrounding the Segara Anakan mangroves, i.e., shrubs, Nypa, rice felds, land forest, mangroves, settlements, and waters. Each consists of ten plots recorded in the corrected Landsat 8 OLI image. Te average spectral response of each land cover is shown in Figure 2(b).
Te spectral response was collected from various features in Segara Anakan mangrove forests in the Landsat 8 OLI corrected images to create the automatic formula to separate mangrove and nonmangrove vegetation and other features.

Index Accuracy Assessment.
Te accuracy assessment of the research results consists of 2 stages: the spatial extent accuracy and the canopy density index assessments. Te extent and boundary of mangroves with other vegetation will be visually clear using the RGB composite image of ρ NIR-ρ SWIR1-ρ Red . Te problem is how to automatically trace and capture what visually looks like a mangrove through an algorithm. In this research, the mangrove map is a map of mangroves as a mangrove community in the form of a tree community, from very sparse, which can be recorded in the image, to dense canopies. It cannot capture shrubs that are usually a part of the mangrove associated with the ecosystem. Similarly, the canopy index is a relative index that does not use the number of trees within a certain area. Classifed as a dense index/the high index is the high value captured in the pixels as a spectral response from the execution of an algorithm.
Te accurate assessment of the canopy density index in the study uses 200 randomly distributed points. It uses simple statistical methods to determine the linear relationship between the index created by this algorithm versus the indexes created by several previous vegetation indices. calculated from other green vegetation using all the information attached to the satellite images. In the paper, the boundary of the mangrove vegetation index is examined using the canopy closure approach. Vegetation, soil, water, and seasonal and diurnal intertidal interactions are essential features that contribute to the pixel composition of mangroves in satellite remotely sensed images [31]. Mangrove presence becomes sharper in the Landsat image when displayed through an RGB (Red-Green-Blue) composite of ρ NIR − ρ SWR1 − ρ Red , as shown in Figure 3(a). However, the problem is how to automatically capture spatial extent as shown in the visual appearance. Te properties of the vegetation will refect ρ NIR strongly, while ρ Red plays an essential role in determining vegetation due to photosynthetic activities and ρ SWIR1 sensitivity to evaporation, the liquid water content in the leaf, and tidal inundation [31,47,55,[57][58][59]. In dense mangroves, ρ SWIR1 will be slightly higher than ρ Red but will change twice as high in in nonmangroves (Figure 2(b)).

Results and Discussion
Similarly, in the nonmangrove, ρ NIR is high but gradually decreases in mangroves and is absorbed in the water. Based on these spectral characteristics, it is possible to separate Table 3: Summaries of the existing vegetation indices (Vis) and the formulas.  mangroves from other nonmangroves using a combination of ρ NIR , ρ SWIR1, and ρ Red refectances. Te normalized difference water index (NDWI) [60] is used to delineate and sharpen the water features (inland and open waters) and reduce the spectral refectances of feature elements in the land. Tis index is applicable in identifying mangroves because of the liquid water component in the leaf canopy. Likewise, NDVI [46] can also be used as a basis for delineating land and eliminating/weakening the spectral of marine features.

VI's Formula based on chlorophyll content
Based on Figure 2, therefore, to create the mangrove maps, there are two steps. Te frst is delineating the land boundary by increasing the spectral forest vegetation and weakening/eliminating the spectral response of other nonforest areas. Te second is tracing and capturing the mangroves on the land.

Improve the Spectral Response of Vegetation.
Te frst step aims to improve the spectral response of vegetation (land forests and mangroves) while reducing/weakening the spectral response of marine objects (water, coral, and mud), shrubs, settlements, and open land. Although the NDVI formula can capture land and eliminate water features, it is less accurate in distinguishing land vegetation. At this step, it is necessary to preserve vegetation features while ignoring other features, hence the slight modifcation of the NDVI formula by replacing the ρ NIR in the denominator with ρ SWIR1 . Tis will efectively improve the spectral response for both land and mangrove forests. To trace the land vegetation using ρ Red , ρ NIR , and ρ SWIR1 in the following equation: As known, ρ NIR is sensitive to all vegetation types, but the greater part will be absorbed in the water environment. Likewise, ρ Red will show lower refectance in the vegetation environment due to chlorophyll absorption and higher refectance in the water environment.
Te combination of ρ Red and ρ NIR , referring to the spectral response in Figure 3(b), cannot separate mangrove forests from terrestrial forests because they have the same spectral value. It would be more efective to use ρ SWIR1 which has diferent index ranges and is longer. Figure 3(c) compares the sharpness of separating the forest canopy from other vegetation, including shrubs, by replacing ρ SWIR1 in the denominator with ρ Red in the NDVI formula.

Identifcation and Tracing of Mangroves.
Te diference in ρ SWIR1 index between the mangrove and nonmangrove forest is thought to be a diference in water content or moisture in the forest canopy. To trace the mangrove extent is by using ρ NIR , ρ SWIR1, and ρ Red in the following equation: Te constant of 0.65 in equation (3) reduces the ρ Red refectance to avoid the infnite value in the mangrove edge bordering the sea. In the outer mangrove, when mangroves   International Journal of Forestry Research border the sea, the index of 1 pixel will be overestimated because ρ SWIR1 is lower than ρ Red (Figure 3(d) and zoomed in Figure 3(e)). the intensity of the ρ Red to obtain a more precise boundary is reduced and an index that is appropriate for the actual conditions. Based on the two equations, the automatic mangrove map and index (AMMI) in this study can be written as follows: Te results of AMMI execution using the combination of ρ Red, ρ NIR, and ρ SWIR1 in radiometrically corrected Landsat 8 OLI, September 9, 2019, trace and capture the spatial extent and present the relative canopy density of the mangrove in one band of the grayscale image is shown in Figure 4(a). Te magnitude of the spectral response in the spectral color is shown in Figure 4(b). Te fgures show that mangroves will refect a stronger spectral response, while nonmangrove features will be weaker.

Accuracy Assessment.
Spatially, the AMMI captures the mangroves from sparse mangroves, indicated by the low spectral sensitivity with an index of about 5, to dense mangroves (>20), as shown in Figure 5(a), and the index below 5 is classifed as nonmangrove. Te relationship of the index to the NDVI, using 200 randomly distributed points (Figure 3(a)), has a correlation value (R 2 ) of 0.62 ( Figure 5(b)), and the relationship to NDWI/NDMI shows a very high correlation value, even reaching 0.99 ( Figure 5(c)).
As is well known, NDWI/NDMI were initially created and applied for the terrestrial forest, but it was also explained that forest canopy density from satellite imagery is closely related to the water content in the canopy.

Monitoring Mangrove Changes Using Multitemporal
Images. To study the evolution of mangroves chronologically and monitor their condition, the AMMI performs well in detecting mangrove changes. Various satellite images can be used, namely Landsat 5 TM 1989 (Figure 6 Figure 6(e) shows this was an inset, a small clip from a highspatial-resolution Google satellite to evaluate the algorithm's performance. Te accuracy of the mangrove extent is comparable to the high resolution of Google's satellite 2022 covering the research area. AMMI only showed the spectral response of mangroves and eliminated other spectral responses such as shrubs, settlements, dry land, or forest land.
Te spatial accuracy of AMMI in capturing and trapping the mangrove changes will be clearer, as shown in Figure 7. Figure 7(a) depicts the initial conditions before the SNP mangroves were damaged. Changes in mangroves in the SNP area are generally caused by logging/felling to create ponds/aquaculture and typically bloom between 1989 and 2002, as illustrated in Figure 7(b). A symmetrical square felling pattern determined the conversion to ponds. However, in the following imagery (Figure 7(c)), many of these ponds have been abandoned and converted back to mangrove forests, and Figure 4(d) shows mangroves restored to their pre-establishment condition. Some active ponds are found outside the mangrove forest.
Te correlation value between AMMI and NDWI/ NDMI reached 0.99, indicating that the AMMI's performance is comparable to that of the NDWI/NDMI index. Based on these results, AMMI is a breakthrough in NDWI/ NDMI automatic innovation that automatically traces and captures mangroves and informs the index.

Learn about Mangrove
Changes from the Mangrove Production Forests. NDWI/NDMI vegetation indices in the mangrove research are less commonly used than NDVI. However, NDWI/NDMI has many advantages, including quickly identifying minor damage and light disturbances. According to Otero et al. [61], mangrove age from sprout to 7 years old will correlate with a sharply increasing NDWI/ NDMI index. More than seven years old are no longer correlated and are frequently reduced. Tis fnding is valuable because their research was in Matang Mangrove Forest Reserve (MMFR), Malaysia, one of the mangrove production forests rarely found in the mangrove forests in the world. Further research revealed that reducing the NDWI index at mature to old age in mangroves was caused by the gaping process of the canopy. Will the index shift in mangroves from maturity to old age occur when using another index, such as NDVI. Goessens et al. [62] reported the status of mangrove forest production in Matang, Malaysia, located at 4.82 N and 100.59 E, providing an exciting research location for studying the evolution of mangroves.
To avoid breaking the code of ethics and infringing on the rights of neighboring countries, the authors describe at a glance a small part of Malaysia's Matang mangrove forest using multi-temporal Landsat images, as shown in Figure 8.  investigated the high NDVI in young vegetation, which has a bright green color, ρ Green increased after leaf emergence and decreased after canopy closure during early growth, while ρ Red continued to decline. According to these fndings, the highest NDVI is not always found in the most densely forested mangroves, nor is it always found in the oldest and healthiest mangroves. NDVI levels are typically highest in mature mangroves aged 7 to 10 years. Furthermore, as seen in NDWI/ NDMI, the NDVI decreases slightly with increasing mangrove age while remaining in the moderately high range. Besides detecting damage due to logging, NDMI is also sensitive to forest disturbances caused by diseases related to forest health. Reference [64] and vegetation damage due to drought efects [65]; it is also benefcial to monitor water status as an early warning against drought [66].

Application of AMMI in a Broader Area Using the GEE Facilities.
In applying automatic mangrove mapping to a larger area, the GEE facility was used to reveal the Sundarban mangrove forest, stretching from eastern India to Bangladesh, using surface refectance of Landsat 8, 2021 imagery. Based on visual interpretation informs, the Sundarbans mangrove habitat area of 621925 ha is divided into India (214769 ha) and Bangladesh (407156 ha) (Figure 9(a)).
Analysis using the AMMI algorithm shows that the true mangrove in Sundarban India is 183296 ha (85% of habitat), while that in Sundarban Bangladesh is only 189733 ha (47%). Sundarban Bangladesh mangroves are sparsely distributed and concentrated only around rivers, with a very low canopy density and are inhabited mainly by shrubs and open land (Figure 9(b)).
Te Sundarban mangrove has become an icon of the world's largest mangroves. Is there any donor agency that intends to replant and reforest Sundarban mangroves in Bangladesh?

Conclusion
Te AMMI application in mangrove mapping in Musi Delta, South Sumatra, Indonesia, based on the spatial distribution and the index of canopy density using Landsat 5 TM, Landsat 7 ETM, Landsat 8 OLI, and Sentinel 2 performs well. Te mangrove extent has been spatially traced and captured, corresponding to and matching the visual satellite images. Te canopy density index of AMMI versus NDVI has a correlation value (R 2 ) of 0.62 in correlation diagrams, and AMMI versus NDWI/NDMI has a signifcant value with an R 2 of 0.99.
Te paper revealed how to automatically trace and capture mangroves for future research and eliminate manual digitizing, which is exhausting and time-consuming and frequently results in inaccuracy due to misinterpretation. Operationally, it is simple to use, produces mangrove maps quickly, and efciently monitors the mangrove condition from time to time. Te instantaneous application may be sufcient for monitoring mangrove conditions in protected mangrove forests by local communities, practitioners, and conservationists.
Scientifcally, for future research, mapping of mangroves over a larger area and describing the surrounding physical environment is the baseline data for the mangrove condition. Other supporting data, such as tree density and their age per unit area, biota content statistics, and geo-biochemical data, will determine the mangrove health index.

Data Availability
Data are available in Supplementary.

Conflicts of Interest
Te authors declare that they have no conficts of interest.

Mangrove habitat
India Bangladesh Ponds Ba dy ad ha ri R .
A ip a n g a s ia R .