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）算法生成无云年度合成影像。

几何中值算法可对每个像素生成像素级统计汇总，本研究中，该汇总对应辖区内10m×10m的地表区域及其在一个日历年中的全部观测结果。

所有汇总像素共同构成阿瓜斯卡连特斯州的12波段合成影像，波段依次为：海岸气溶胶波段、蓝光波段、绿光波段、红光波段、植被红边5波段、植被红边6波段、植被红边7波段、近红外波段、窄近红外波段、水汽波段、短波红外1（swir1）波段与短波红外2（swir2）波段。

同步生成了名为GeoMad的衍生产品，其计算了文献[Roberts, D., Dunn, B., & Mueller, N. 2018]中定义的稳健离散统计量MAD。在最终的影像合成产品中，三个像素波段分别代表三种距离下的离散程度：光谱距离（smad）、欧氏距离（emad）与Bray-Curtis距离（bcmad）。

研究期间（2017-2021）还生成了额外波段，用于表征不同环境条件，这类条件可通过波段的多种组合加以捕捉，此类组合亦被称为光谱指数（spectral indices），可用于检测植被分布、水体赋存、城市化进程等。最终共计算得到28个指数，分为4大类：

### 植被指数
1. 抗大气影响植被指数（Atmospherically Resistant Vegetation Index，Kaufman 1972）
2. 增强型植被指数（Enhanced Vegetation Index，Huete 2002）
3. 修正土壤调整植被指数（Modified Soil Adjusted Vegetation Index，Qi等，1994）
4. 归一化差异叶绿素指数（Normalized Difference Chlorophyll Index，Mishra & Mishra, 2012）
5. 归一化差异水分指数（Normalised Difference Moisture Index，Gao 1996）
6. 归一化差异植被指数（Normalized Difference Vegetation Index，Rouse 1973）
7. 优化土壤调整植被指数（Optimized Soil Adjusted Vegetation Index，Rondeaux, 1996）
8. 简单比值植被指数（Simple Ratio Vegetation Index，Jordan, C.F. 1969）
9. 土壤调整植被指数（Soil Adjusted Vegetation Index，Huete 1988）
10. 可见光抗大气指数（Visible Atmospherically Resistant Index，Gittleson 2002）

### 建筑用地指数
1. 建筑用地波段比值指数（Band Ration For Built-up Area，Waqar 2012）
2. 建筑用地提取指数（Built-up Area Extraction Index，Bouzekri 2015）
3. 建筑指数（Built-up Index，He等，2010）
4. 基于指数的建筑用地指数（Index-based Built-up Index，Xu 2008）
5. 新型建筑指数（New Built-up Index，Jieli等，2010）
6. 归一化差异建筑指数（Normalized Difference Built-up Index，Zha 2003）
7. 归一化建筑用地指数（Normalized Built-up Area Index，Waqar 2012）
8. 城市指数（Urban Index，Kawamura 1996）

### 水体指数
1. 修正归一化差异水体指数（Modified Normalized Difference Water Index，Xu 1996）
2. 归一化差异水体指数（Normalized Difference Water Index，Mcfeeters 1996）
3. 水体指数（Water Index，Fisher 2016）

### 其他指数
1. 裸土指数（Bare Soil Index，Rikimaru等，2002）
2. 裸土指数（Bare Soil Index，Wanhui 2004）
3. 火烧面积指数（Burn Area Index，Martin 1998）
4. 黏土矿物比率指数（Clay Minerals Ratio，Drury 1987）
5. 铁质矿物比率指数（Ferrous Minerals Ratio，Segal 1982）
6. 氧化铁比率指数（Iron Oxide Ratio，Segal 1982）
7. 归一化燃烧比率指数（Normalized Burn Ratio，Lopez Garcia 1991）
8. 归一化差异积雪指数（Normalised Difference Snow Index，Hall 1995)