第三章 基于 GIS 的马山生态环境与旅游空间分析
3.1 生态环境时空演变分析
3.1.1 基于 NDVI 的马山植被覆盖度时序变化分析（2000—2025）
在 GEE 平台上，用 Landsat 影像（30 m 分辨率）每隔 5 年提取一次生长季 NDVI 中值合成影像。NDVI 的算法是：
式中，NIR 为近红外波段反射率，RED 为红光波段反射率。NDVI 取值范围为 -1 至 1，值越高表示植被覆盖越好。
图3.1 画了马山半岛 6 个年份的 NDVI 分布。从图上看，西边和南边的山区植被一直比较好（NDVI 多在 0.3—0.6 之间），东北方向靠近城区的地方就差一些（0.1—0.3），说明城镇化确实挤压了植被空间。
对研究区各年 NDVI 影像进行统计分析，结果如表3.1 所示。
马山半岛 NDVI 均值从 2000 年的 0.2878 经波动后上升至 2025 年的 0.3132，整体植被覆盖持续改善。
又在 GEE 上把 2000—2025 年每一年的 NDVI 均值都提出来画了个折线图（图3.2）。总体上看是往上走的，拟合了一条线性趋势线，斜率是 +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 差值图（图3.3），看看哪些地方植被变好了、哪些变差了。统计下来，75.2% 的区域 NDVI 是增加的（植被在恢复），24.8% 是下降的（植被在退化），平均变化是 +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 的土地利用变化分析
基于 GEE 平台和 Landsat 影像，采用 NDVI、MNDWI、NDBI 光谱指数决策树分类，将土地利用划分为林地、耕地/草地、建设用地、水域和裸地五类（图3.4）。
整体来看马山还是以林地为主，25 年里林地占比一直在 80%—93% 之间。建设用地虽然总量不大，但变化很明显，从 2000 年的 3.7% 涨到 2008 年的 8.6%（那时候大搞建设），后来政策上开始控制了就慢慢降回到 2023 年的 2.3%。耕地和草地面积则是一路缩小，从 8%—12% 降到了 4%—5%，主要是退耕还林和部分土地被转成了旅游用地。
结合旅游开发进程，土地利用变化可分为三个阶段：
结合旅游开发进程，LUCC 变化分为三阶段：旅游开发期（2000—2008）建设用地从 3.7% 升至 8.6%；建设高峰期（2008—2016）建设用地维持 5%—9% 高位，林地逐步恢复；生态修复期（2016—2025）建设用地降至约 3%，林地恢复至 91%—93%，与 NDVI 趋势高度一致。
Sentinel-2（10 m 分辨率）交叉验证表明（图3.5），两组数据在建设用地（r=0.854, p<0.01）和林地（r=0.719, p<0.05）趋势上高度一致，验证了本研究 LUCC 演变趋势的科学稳健性。
3.1.3 旅游开发核心区的生态环境变化特征
对灵山胜境与拈花湾两大核心区进行微观诊断：
- 灵山核心区：持续演进与局部生态压力固化。灵山大佛自 1997 年落成后经历多期外围扩建（如 2008 年梵宫、2009 年五印坛城及大型露天停车场）。NDVI 差值分析表明，虽然大佛主坡面等早期建设区域在生态保护下展现了显著的植被绿化恢复（ΔNDVI > +0.090），但梵宫展馆及东侧的大型停车场因地表重度硬化，导致该局部区块植被指数出现剧烈下跌（ΔNDVI = -0.190），构成了局部不可逆的生态斑块割裂。
- 拈花湾核心区：高强度空间重塑与生态替代。拈花湾小镇（2015年底开园）的建设使约 1.2 km² 的区域发生了剧烈土地转变。LUCC 变化检测显示，该区域由原生的“山林-荒草-坑塘”自然地表迅速重构为高密度的“商业小镇-人造水景-观赏园林”。虽然内部补种了高密度景观绿化，但相较于原生自然林地，整体生物量仍大幅减少（ΔNDVI = -0.114），这属于高经济附加值的人工环境对原生生态底色的置换，降低了局部的生态调节与涵养功能。从社会经济驱动力来看，王昌静对拈花湾运营模式的研究指出，该项目存在"地产化趋势严重"的问题，规划中拟配套小批量住宅且过度追求纯粹商业地产文旅项目以实现收益最大化[6]。这种资本逐利与过度地产化倾向，正是导致本研究在遥感影像上观察到大面积地表硬化与生态指数剧烈下跌的核心社会经济驱动因素。
3.2 旅游资源空间分布分析
3.2.1 马山旅游资源的空间分布特征
基于高德地图 API 采集的旅游景点 POI 共 129 处。空间分布呈显著的双核集聚、外围稀疏格局：第一核心为灵山胜境—拈花湾片区，集中了全岛近 60% 的景点 POI；第二核心为龙头渚—环湖公园带。半岛东南部和北部景点分布稀疏，反映出旅游开发高度依赖两大核心 IP，南部优质岸线尚未得到有效开发。
3.2.2 旅游资源的核密度分析
采用核密度估计（KDE）方法对景点 POI 进行空间聚集分析（搜索半径 1000 m），结果如图3.6所示：灵山—拈花湾形成高密度"双子核心"，呈典型的单中心极化聚集格局；龙头渚形成次级密度区，峰值仅为核心区的约 1/3；半岛东部和南部沿湖带核密度趋近于零，构成大面积"冷区"。
3.2.3 资源缓冲区叠加分析
为定量评估核心景区对周边旅游服务设施的空间辐射效应，本研究以灵山胜境和拈花湾禅意小镇两个核心景区入口为中心，分别建立 500 m、1000 m 和 2000 m 的同心圆缓冲区，与住宿设施（434 处）和餐饮设施（172 处）的 POI 数据进行空间叠加分析，统计各缓冲区内的设施数量和覆盖比例（表3.2）。
表3.2 核心景区缓冲区内旅游服务设施统计
| 缓冲中心 | 缓冲半径 | 景点数 | 住宿设施 | 占比 | 餐饮设施 | 占比 |
|---|---|---|---|---|---|---|
| 灵山胜境 | 500 m | 47 | 1 | 0.2% | 2 | 1.2% |
| 灵山胜境 | 1000 m | 79 | 132 | 30.4% | 35 | 20.3% |
| 灵山胜境 | 2000 m | 86 | 269 | 62.0% | 86 | 50.0% |
| 拈花湾 | 500 m | 16 | 68 | 15.7% | 54 | 31.4% |
| 拈花湾 | 1000 m | 25 | 84 | 19.4% | 62 | 36.0% |
| 拈花湾 | 2000 m | 30 | 100 | 23.0% | 67 | 39.0% |
| 两核合并 | 2000 m | 115 | 367 | 84.6% | 151 | 87.8% |
分析表明：灵山呈"景区外围集聚"模式（500m 内几乎无设施，1000m 内住宿达 132 处）；拈花湾呈"景区内融合"模式（500m 内即有 68 处住宿和 54 处餐饮）。两核合并 2000m 缓冲区覆盖了 84.6% 的住宿和 87.8% 的餐饮，外围形成显著的"设施真空带"。
3.3 旅游设施空间配套分析
3.3.1 住宿设施的空间分布与等级结构
住宿 POI 共 434 处（民宿客栈 255 处、旅馆 150 处、星级酒店 29 处），呈金字塔型供给与强景区依附特征。灵山—群丰社区为最高密度集聚区，拈花湾—嶂青社区为次高值，南部沿湖带面临"设施真空"。
3.3.2 餐饮设施的空间分布与核密度特征
餐饮 POI 共 172 处，以地方中餐为主导（58.1%）。核密度与住宿高度耦合，南部环湖区域餐饮匮乏，存在"游在南部、吃在北部"的空间错位。
3.3.3 交通设施可达性分析
基于 OSM 路网（506 条路段、179.8 km）评估交通可达性：路网密度高但高等级进出通道单一；500m 公交覆盖率仅 33.2%（表3.3），南部环湖带和东部山区几乎处于公交盲区；77 处停车场集中于核心景区，固化了自驾依赖。
表3.3 公交站点步行等时圈覆盖率
| 步行半径 | 覆盖面积 | 覆盖率 | 未覆盖面积 |
|---|---|---|---|
| 300 m | 5.86 km² | 14.3% | 35.22 km² |
| 500 m | 13.64 km² | 33.2% | 27.44 km² |
| 1000 m | 30.40 km² | 74.0% | 10.68 km² |
3.3.4 旅游设施配套综合评价
综合评价表明：设施空间高度极化（两核 2000m 覆盖 84.6% 住宿和 87.8% 餐饮）；公交覆盖偏低（500m 仅 33.2%）；业态结构单一（民宿占 93.3%，中餐占 58.1%）。