第三章 基于 GIS 的马山生态环境与旅游空间分析
3.1 生态环境时空演变分析
3.1.1 基于 NDVI 的马山植被覆盖度时序变化分析（2000—2025）
为揭示马山半岛近 25 年来植被覆盖的时空演变规律，本研究基于 Google Earth Engine 平台，利用 Landsat 系列卫星影像（Landsat 5 TM / Landsat 8-9 OLI，空间分辨率 30 m），以 5 年为间隔提取 2000、2005、2010、2015、2020、2025 年生长季（5—10 月）的归一化植被指数（NDVI）中值合成影像，以消除云雾和季节波动的干扰。计算公式为：
式中，NIR 为近红外波段反射率，RED 为红光波段反射率。NDVI 取值范围为 -1 至 1，值越高表示植被覆盖越好。
图 3-1 展示了马山半岛 6 个时间断面的 NDVI 空间分布。整体来看，半岛西部和南部山体的植被覆盖持续保持较好（NDVI 值多在 0.3—0.6 之间），而东北部靠近城区的区域 NDVI 值相对较低（0.1—0.3），反映出城镇化对植被的压缩效应。
对研究区各年 NDVI 影像进行统计分析，结果如表 3-1 所示。
表 3-1 马山半岛各年 NDVI 统计特征
| 年份 | 均值 | 标准差 | 最小值 | 最大值 | 中位数 |
|---|---|---|---|---|---|
| 2000 | 0.2878 | 0.0705 | -0.0273 | 0.4247 | 0.3052 |
| 2005 | 0.2503 | 0.0683 | -0.0370 | 0.4172 | 0.2690 |
| 2010 | 0.2858 | 0.0924 | -0.0604 | 0.5049 | 0.3069 |
| 2015 | 0.2611 | 0.0733 | -0.0759 | 0.4370 | 0.2793 |
| 2020 | 0.2928 | 0.0739 | -0.0367 | 0.4456 | 0.3164 |
| 2025 | 0.3132 | 0.0802 | -0.0509 | 0.4899 | 0.3377 |
为进一步揭示研究区植被覆盖的年际演变规律，本研究在 GEE 平台上提取了 2000—2025 年连续 26 年的生长季 NDVI 均值（图 3-2）。整体来看，马山半岛 NDVI 年均值呈显著上升趋势，线性回归斜率为 +0.0027/yr（R² = 0.389），表明 25 年间植被覆盖持续改善。
结合马山半岛旅游开发历程，可将 NDVI 变化划分为三个阶段：
（1）快速开发期（2000—2006 年）：该时期马山半岛处于旅游基础设施集中建设阶段，灵山大佛景区扩建、拈花湾项目启动等大型工程导致大量地表被扰动，NDVI 均值由 2000 年的 0.2902 波动下降至 2006 年的 0.2264，植被覆盖明显退化。
（2）过渡调整期（2007—2015 年）：随着大规模建设活动趋于尾声，NDVI 均值逐步企稳，在 0.25—0.29 区间波动。该时期虽有局部开发活动，但整体破坏力度减弱，植被开始初步恢复。
（3）生态恢复期（2016—2025 年）：NDVI 均值显著上升并稳定在 0.30 以上，2024 年达到观测期最高值 0.3288。这一阶段与无锡市推行的生态修复政策、景区绿化工程以及"绿水青山就是金山银山"发展理念的深入实践密切相关。
为进一步识别植被变化的空间分异特征，本研究利用 QGIS 栅格计算器生成了 2025 年与 2000 年的 NDVI 差值图，并结合 OpenStreetMap 地理标注数据对关键区域进行标注（图 3-3）。结果显示，研究区 75.2% 的面积 NDVI 呈正向变化（植被改善），24.8% 的面积 NDVI 下降（植被退化），整体平均变化量为 +0.025。
从空间分布来看，植被退化区主要与旅游开发和城镇建设高度吻合：拈花湾禅意小镇（ΔNDVI = -0.114）、太悦度假酒店（ΔNDVI = -0.127）以及弘阳三万顷住宅区（ΔNDVI = -0.074）等新建项目所在区域植被退化显著。在灵山景区内部（图 3-3 右下角详情图），梵宫（ΔNDVI = -0.190）因大型建筑占用植被用地退化最为严重，五印坛城和九龙灌浴区域（ΔNDVI = -0.044）同样呈现退化趋势。
与此形成对比的是，西部山区（ΔNDVI = +0.102）和灵山大佛周边（ΔNDVI = +0.096）植被覆盖显著改善，前者得益于自然植被恢复，后者则体现了成熟景区绿化管理的成效。这一空间格局表明，旅游开发对植被的影响具有明显的空间分化效应：新建开发项目导致的植被退化集中而剧烈，而成熟景区和远离开发区的山体植被则总体呈恢复趋势。此外，需要注意的是部分红色区域（如月湾水库附近）实为新建水塘或湖泊，水体 NDVI 天然呈负值，并非植被退化。
3.1.2 基于 LUCC 的土地利用变化分析
本研究利用 Google Earth Engine 平台，基于 2000—2025 年 Landsat 系列影像，采用 NDVI、MNDWI、NDBI 光谱指数的决策树分类方法，将马山半岛土地利用类型划分为林地、耕地/草地、建设用地、水域和裸地五类，逐年提取各类面积占比（图 3-4）。
从整体结构来看，马山半岛以林地为绝对主体，25 年间林地面积占比始终维持在 80%—93% 之间。建设用地占比较低但变化趋势明显，从 2000 年的 3.7% 波动上升至 2008 年的 8.6%（旅游开发高峰），随后在生态修复政策推动下逐步回落至 2023 年的 2.3%。耕地/草地面积呈下降趋势，从早期的 8%—12% 缩减至近年的 4%—5%，反映了退耕还林和部分农用地向旅游服务用地的转化。
结合旅游开发进程，土地利用变化可分为三个阶段：
（1）旅游开发期（2000—2008）： 灵山大佛景区扩建、拈花湾项目启动，建设用地从 3.7% 上升至 8.6%，林地面积出现阶段性下降（2008 年降至 78.9%），耕地/草地波动较大（5%—12%），体现了土地利用的快速转型。
（2）建设高峰期（2008—2016）： 拈花湾禅意小镇建成运营、梵宫修建、太悦度假酒店等项目集中落地，建设用地维持在 5%—9% 的较高水平。林地在该阶段有所恢复（2015 年回升至 87.1%），说明景区绿化工程开始发挥作用。
（3）生态修复期（2016—2025）： 在生态文明建设背景下，建设用地占比持续下降至 3% 左右，林地恢复至 91%—93% 的历史高位，耕地/草地稳定在 4%—5%。这一阶段植被恢复效果显著，与 NDVI 分析所揭示的趋势高度一致。
需要说明的是，基于光谱指数的像素级分类方法在 Landsat 30 m 分辨率下对混合像素的识别存在一定局限性，个别年份（如 2004、2013 年）因可用影像质量较低导致分类精度偏低，2007、2009、2011 年因有效影像不足未纳入分析。为减少影像质量波动对趋势判断的干扰，本研究对分类结果进行了 3 年滑动平均处理。
为验证 Landsat 分类结果的可靠性，本研究引入 Sentinel-2（10 m 分辨率）影像对 2016—2025 年的分类结果进行交叉验证（图 3-5）。采用与 Landsat 相同的光谱指数阈值和决策树规则，对两颗卫星的分类结果进行逐年对比。结果表明，两组数据在建设用地（r=0.854, p<0.01）和林地（r=0.719, p<0.05）的变化趋势上具有高度一致性，均呈现建设用地持续下降、林地面积稳步恢复的总体格局。
在绝对值上，Sentinel-2 识别的建设用地占比系统性高于 Landsat 约 4.2 个百分点，这主要是因为 10 m 分辨率能够更精确地捕捉景区步道、停车场、小型服务设施等细碎建设用地，而 Landsat 30 m 像元将这些区域与周围植被混合，多被归入耕地或林地类别。据此推断，马山半岛近年实际建设用地占比约为 7%—8%，高于 Landsat 直接分类结果的 3%—4%。两组独立数据的趋势一致性证实了本研究 LUCC 分析结论的稳健性。
3.1.3 旅游开发核心区的生态环境变化特征
本研究进一步选取马山半岛最具代表性的两大旅游开发核心区——灵山胜境（含梵宫及周边配套）与拈花湾禅意小镇作为典型对象，结合 3.1.1 和 3.1.2 的分析结果进行微观诊断。
（1）灵山核心区："持续演进与局部生态压力" 灵山大佛自 1997 年落成后，不断进行外围扩展（如 2008 年灵山梵宫落成、2009 年五印坛城落成）。从 NDVI 差值图（图 3-3）可以看出，灵山景区整体在早期经历了明显的植被替换，但在 2016 年至今的生态恢复期，大佛主坡面及早期景区内部呈现出强烈的正向绿化改善（ΔNDVI > +0.090）。然而，梵宫所在位置及景区东侧的大型露天停车场，由于高度硬化的地表，不仅取代了原有林草地，且使得该局部区块的植被指数出现断崖式下跌（ΔNDVI = -0.190），反映了大型人造景观及硬化场地对局部微生态不可逆的改变。
（2）拈花湾核心区："高强度重塑与生态替代" 拈花湾自 2015 年底开园，其前期建设（2010—2015年）位于耿湾及周边近湖丘陵地带。LUCC 变化检测显示，拈花湾的建设使得该片区约 1.2 km² 的土地发生了剧烈的性质转变——由原本的"山林-草荒-坑塘"混合地带，迅速重构为高密度的"建设用地-人造水景-观赏园林"。这一土地整理过程在 NDVI 曲线上留下了深深的"退化凹谷"。虽然小镇内部植入了大量园林绿化，但相较于开发前原生的连片自然林地，地表生物量仍大幅减少（ΔNDVI = -0.114）。这种生态变化属于典型的"用高经济附加值的人工环境置换原生低效自然环境"，在带来巨大旅游收益的同时，客观上也降低了该片区生态涵养功能。
3.1.4 小结
本节通过长时序 NDVI 植被覆盖度分析和 LUCC 土地利用结构变化分析，定量揭示了马山半岛在过去四分之一个世纪中的生态环境演变脉络：
第一，马山半岛的生态演变历程与其旅游开发轨迹具有高度的耦合性，呈现出清晰的 "快速开发破坏（2000-2006） → 局部震荡平衡（2007-2015） → 整体生态恢复（2016至今）" 的阶段性特征。目前，半岛林地覆盖率稳定在 90% 以上，整体生态底色优良。
第二，旅游开发导致生态环境变化的"空间分化效应"极其显著。广大的中高海拔山林地带因受到生态红线保护，植被持续向好；而滨湖低地、山谷平缓地带（特别是灵山、拈花湾及部分高端地产项目区）因承载了高强度的旅游基础设施建设，建设用地扩张对原生植被和耕地产生了显著的局部蚕食效应。
第三，旅游开发的生态压力从"增量扩张"转向"存量固化"。近年来马山新增建设用地已临界停滞，但已建成的大型景区和硬化场地（如超大面积停车场、人造场馆）构成了永久性的生态斑块割裂。如何在严格控制开发边界的前提下，通过对现有硬化区域引入微生态改造和绿色基础设施补充，将是马山旅游后续发展的核心环境挑战。
3.2 旅游资源空间分布分析
3.2.1 马山旅游资源的空间分布特征
（待撰写）
3.2.2 旅游资源的核密度分析
（待撰写）
3.2.3 资源缓冲区叠加分析
（待撰写）
3.3 旅游设施空间配套分析
3.3.1 住宿设施的空间分布与等级结构
（待撰写）
3.3.2 餐饮设施的空间分布与核密度特征
（待撰写）
3.3.3 交通设施可达性分析
（待撰写）
3.3.4 旅游设施配套综合评价
（待撰写）