AgsSAT Multiannual (2017-2021) Sentinel-2 Vegetation Indices Composites

AgsSAT Multiannual (2017-2021) Sentinel-2 Vegetation Indices Composites The 2,815 images available for the state of Aguascalientes, Mexico for the years 2017 to 2021 were processed using the Open Data Cube (ODC) platform [Lewis et al. (2017), Gavin et al. (2018), https://www.opendatacube.org/]. These images correspond to multiple coverages of the region of interest. The images were then used to generate cloud-free annual composites by applying geometric median (geomedian) algorithm, as defined in [Roberts et al. (2017)]. Geomedian algorithm produces a pixel-level summary for every pixel, in this case this means that each summary corresponds to a 10m x 10m region in the territory and its observations throughout a calendar year. All these summary pixels form a 12-band (coastal aerosol, blue, green, red, vegetation red edge 5, vegetation red edge 6, vegetation red edge 7, near-infrared, narrow nir, water vapor, swir1 and swir2) composite of the state of Aguascalientes. Another product called GeoMad was generated, which calculates the robust dispersion statistic called MAD, as defined in [Roberts, D., Dunn, B., & Mueller, N. (2018)]. In the resulting image composite, each of the three-pixel bands represents the variation over three distances: Spectral Distance (smad), Euclidean Distance (emad) and the Bray-Curtis Distance (bcmad). More bands were generated to represent different environmental conditions during the study years (2017-2021), these conditions can be captured by analyzing various combinations of bands, these combinations are also called spectral indices, which allow detecting vegetation, presence of water, urbanization, etc., Finally, 28 indices divided into 4 categories were calculated: Vegetation Indices (Atmospherically Resistant Vegetation Index, Kaufman 1972) (Enhanced Vegetation Index, Huete 2002): (Modified Soil Adjusted Vegetation Index, Qi Et Al. 1994) (Normalized Difference Chlorophyll Index, Mishra & Mishra, 2012) (Normalised Difference Moisture Index, Gao 1996) (Normalized Difference Vegetation Index, Rouse 1973) (Optimized Soil Adjusted Vegetation Index, Rondeaux. 1996) (Simple Ratio Vegetation Index Jordan, C.F.1 969) (Soil Adjusted Vegetation Index, Huete 1988) (Visible Atmospherically Resistant Index, Gittleson 2002) Built-up Indexes (Band Ration For Built-up Area, Waqar 2012) (Built-up Area Extraction Index, Bouzekri 2015) (Built-up Index, He Et Al. 2010) (Index-based Built-up Index, Xu 2008) (New Built-up Index, Jieli Et Al. 2010) (Normalized Difference Built-up Index, Zha 2003) (Normalized Built-up Area Index, Waqar 2012) (Urban Index, Kawamura 1996) Water Indices (Modified Normalized Difference Water Index, Xu 1996) (Normalized Difference Water Index, Mcfeeters 1996) (Water Index, Fisher 2016) Other Indices (Bare Soil Index, Rikimaru Et Al. 2002) (Bare Soil Index, Wanhui 2004) (Burn Area Index, Martin 1998) (Clay Minerals Ratio, Drury 1987) (Ferrous Minerals Ratio, Segal 1982) (Iron Oxide Ratio, Segal 1982) (Normalized Burn Ratio, Lopez Garcia 1991) (Normalised Difference Snow Index, Hall 1995).AgsSAT多年份（2017-2021）哨兵二号（Sentinel-2）植被指数合成数据集。针对墨西哥阿瓜斯卡连特斯州2017至2021年期间可用的2815幅影像，采用开放数据立方体（Open Data Cube，ODC）平台进行处理[Lewis等，2017；Gavin等，2018；https://www.opendatacube.org/]。此类影像对应研究区域的多次覆盖采集。随后，依据[Roberts等，2017]中定义的几何中值（geomedian）算法，生成无云年度合成影像。几何中值算法可针对每个像素生成像素级汇总结果，本研究中，每个汇总结果对应辖区内10米×10米的区域及其在一个日历年中的观测数据。所有此类汇总像素共同构成阿瓜斯卡连特斯州的12波段合成影像，波段依次为：海岸气溶胶波段、蓝光波段、绿光波段、红光波段、植被红边5波段、植被红边6波段、植被红边7波段、近红外波段、窄近红外波段、水汽波段、短波红外1波段（swir1）与短波红外2波段（swir2）。 此外，本研究还生成了另一款名为GeoMad的产品，其可计算稳健离散统计量MAD，依据[Roberts, D., Dunn, B., & Mueller, N., 2018]中的定义。在最终生成的影像合成产品中，三个像素波段分别表征三种距离的变异程度：光谱距离（smad）、欧氏距离（emad）与布莱-柯蒂斯距离（bcmad）。 研究期间（2017-2021）还生成了更多用于表征不同环境条件的波段，此类环境条件可通过分析波段的多种组合加以捕捉，此类组合亦被称为光谱指数，可用于检测植被、水体存在情况、城市化进程等。最终，共计算得到分为4大类的28项指数： 1. 植被指数： - 抗大气植被指数（Atmospherically Resistant Vegetation Index，Kaufman，1972） - 增强型植被指数（Enhanced Vegetation Index，Huete，2002） - 修正土壤调整植被指数（Modified Soil Adjusted Vegetation Index，Qi等，1994） - 归一化差异叶绿素指数（Normalized Difference Chlorophyll Index，Mishra & Mishra，2012） - 归一化差异水分指数（Normalised Difference Moisture Index，Gao，1996） - 归一化差异植被指数（Normalized Difference Vegetation Index，Rouse，1973） - 优化土壤调整植被指数（Optimized Soil Adjusted Vegetation Index，Rondeaux，1996） - 简单比值植被指数（Simple Ratio Vegetation Index，Jordan, C.F.，1969） - 土壤调整植被指数（Soil Adjusted Vegetation Index，Huete，1988） - 可见光抗大气指数（Visible Atmospherically Resistant Index，Gittleson，2002） 2. 建筑指数： - 建筑区波段比值指数（Band Ratio For Built-up Area，Waqar，2012） - 建筑区提取指数（Built-up Area Extraction Index，Bouzekri，2015） - 建筑指数（Built-up Index，He等，2010） - 基于指数的建筑指数（Index-based Built-up Index，Xu，2008） - 新型建筑指数（New Built-up Index，Jieli等，2010） - 归一化差异建筑指数（Normalized Difference Built-up Index，Zha，2003） - 归一化建筑区指数（Normalized Built-up Area Index，Waqar，2012） - 城市指数（Urban Index，Kawamura，1996） 3. 水体指数： - 修正归一化差异水体指数（Modified Normalized Difference Water Index，Xu，1996） - 归一化差异水体指数（Normalized Difference Water Index，Mcfeeters，1996） - 水体指数（Water Index，Fisher，2016） 4. 其他指数： - 裸土指数（Bare Soil Index，Rikimaru等，2002） - 裸土指数（Bare Soil Index，Wanhui，2004） - 火烧面积指数（Burn Area Index，Martin，1998） - 黏土矿物比值指数（Clay Minerals Ratio，Drury，1987） - 铁质矿物比值指数（Ferrous Minerals Ratio，Segal，1982） - 氧化铁比值指数（Iron Oxide Ratio，Segal，1982） - 归一化燃烧比率（Normalized Burn Ratio，Lopez Garcia，1991） - 归一化差异积雪指数（Normalised Difference Snow Index，Hall，1995）