✅ 博主简介：擅长数据搜集与处理、建模仿真、程序设计、仿真代码、论文写作与指导，毕业论文、期刊论文经验交流。
✅ 具体问题可以私信或扫描文章底部二维码。
(1) 本文基于3S技术（遥感、地理信息系统和全球定位系统）对怒江流域岩桑树水电站近坝区的滑坡敏感性进行了全面评估，旨在构建有效的滑坡风险管理体系。首先，利用高分辨率遥感影像解译与详实的野外调查数据，本文在ARCGIS平台上建立了滑坡灾害的空间数据库，并依照地质信息元数据标准对区域内的滑坡灾害进行整理与分析。通过遥感解译和实地考察，研究区的滑坡发育类型、规模、物质组成和成因机制得到了系统的识别与归纳，形成了一套成熟的区域地质灾害调查方法。该方法结合了遥感影像处理、现场实地验证与数据统计的综合手段，为准确、高效地掌握研究区滑坡灾害提供了科学依据。地质、水文地质、地形地貌等资料的结合使得滑坡数据库更加全面，能够为进一步的滑坡敏感性评价提供可靠的基础数据。
(2) 在滑坡敏感性评价中，选择适宜的制图单元是影响结果精度的重要因素。本文对滑坡敏感性评价中使用的制图单元进行了详细探讨，将其分为网格单元、地域单元、均一条件单元、自然斜坡单元和地形单元，并着重探索了自然斜坡单元的智能划分方式。在研究中，为找出最适合的滑坡敏感性评价制图单元，选择了自然斜坡单元和网格单元两种形式进行对比。通过对比发现，制图单元的选择对敏感性评价结果影响显著。自然斜坡单元具有较好的自然边界表现，能够有效地反映滑坡的发育特征，而网格单元则具有规则性，适合于数据的系统处理。最终的敏感性评价结果表明，不同制图单元的选择对滑坡风险的识别和预测存在一定的差异，建议在具体应用中根据研究目的选择合适的单元类型。
(3) 在滑坡敏感性影响因子的选取和分析中，本文总结了前人研究的经验，构建了较为完善的滑坡敏感性影响因子体系。基于自然斜坡单元和网格单元两种制图单元，详细阐述了对影响因子的提取步骤，并对因子提取过程中的地形起伏度因子的尺度效应进行了讨论。研究中，地形、降雨、地质条件和植被覆盖度等多个因素被纳入影响因子体系，并利用确定性系数和滑坡相对密度对不同因子的敏感性进行定量和定性分析。在分析中，采用了变维分形理论对滑坡分布特征与各影响因子的关系进行探讨，发现不同因子的组合条件对于滑坡的发育具有显著的影响。结合滑坡的实际分布数据，研究总结出最适合滑坡发育的条件组合，并基于此为滑坡敏感性评价模型的建立提供了科学的依据。
(4) 在滑坡敏感性评价模型的建立过程中，本文采用了组合赋权法与Logistic回归法对不同制图单元的影响因子进行了敏感性计算。组合赋权法结合了主观赋权与客观赋权的优势，能够有效地解决各因子之间的权重问题；而Logistic回归法则以其对多因子组合影响的敏感性适用于滑坡敏感性分析。研究结果显示，基于规则网格单元的Logistic回归模型在本研究区的滑坡敏感性评价中具有较高的适用性和精度。通过对模型结果的准确度进行比较，本文采用了经验概率和ROC曲线对各滑坡敏感性模型进行了分析，结果表明，基于网格单元的Logistic模型最适合岩桑树水电站近坝区的滑坡敏感性评价。这一模型能够有效地预测滑坡发生的可能性，为水电站的安全运行与滑坡风险管理提供了技术支撑。
(5) 通过上述研究，本文建立了一套从基础数据采集、遥感图像处理、滑坡灾害信息提取、滑坡影响因子分析到滑坡敏感性区划的综合评价体系。在此过程中，利用遥感影像与GIS技术有效地提取了滑坡灾害的空间特征，并通过多种理论方法对滑坡影响因子进行了系统的分析，最终形成了适用于岩桑树水电站库区的滑坡敏感性评价模型。研究表明，3S技术在滑坡敏感性评价中的应用能够显著提高工作效率和结果精度，为类似区域的地质灾害防治提供了有益的技术参考。
import geopandas as gpd
import rasterio
import numpy as np
import matplotlib.pyplot as plt
from rasterio.plot import show
from sklearn.linear_model import LogisticRegression
# 读取地形因子栅格数据
slope_raster = 'slope_factor.tif'
with rasterio.open(slope_raster) as src:
slope_data = src.read(1)
profile = src.profile
# 设置坡度阈值，提取高风险区域
slope_threshold = 30
high_risk_slope = np.where(slope_data > slope_threshold, 1, 0)
# 计算高风险区域面积
high_risk_area = np.sum(high_risk_slope)
total_area = slope_data.size
high_risk_percentage = (high_risk_area / total_area) * 100
print(f"高坡度风险区域占比: {high_risk_percentage:.2f}%")
# 显示地形坡度栅格图与高风险掩膜
fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(14, 7))
ax1.set_title('Slope Factor Raster')
show(slope_data, ax=ax1, cmap='terrain')
ax2.set_title('High Slope Risk Mask')
show(high_risk_slope, ax=ax2, cmap='Reds')
plt.show()
# 保存高风险区域掩膜为新的栅格文件
output_path = 'high_slope_risk_mask.tif'
with rasterio.open(output_path, 'w', **profile) as dst:
dst.write(high_risk_slope.astype(rasterio.uint8), 1)
# 使用Logistic回归进行滑坡敏感性分析（示例）
X = slope_data.reshape(-1, 1)
y = high_risk_slope.reshape(-1)
model = LogisticRegression().fit(X, y)
print(f"Logistic回归模型系数: {model.coef_}")