调用arcpy，计算每一年的每一个住房shp数据到上一年污染点shp数据的最短距离

在ArcGIS Pro或Python环境中，你可以使用 arcpy 和 spatial analysis 工具集来计算每个住房Shapefile (SHP) 到最近的污染点Shapefile (Pollution Points SHP) 的最短距离。这里是一个简单的步骤概述：

1. 导入所需的模块：

import arcpy
from arcpy import env

2. 设置环境变量（如果你需要指定工作空间、路径等）：

env.workspace = "your_workspace_path"  # 替换为你的工作目录

3. 加载两个Shapefile数据集：

houses_shp = "Houses.shp"
pollution_points_shp = "PollutionPoints.shp"

houses = arcpy.mp.FeatureClassToFeatureClass_management(houses_shp, "", "temp_houses")
pollution_points = arcpy.mp.FeatureClassToFeatureClass_management(pollution_points_shp, "", "temp_pollution_points")

4. 创建空间参考以确保两者在同一投影下：

if houses.spatialReference != pollution_points.spatialReference:
    arcpy.Project_management(houses, "temp_houses_projected", pollution_points.spatialReference)
    houses = "temp_houses_projected"

5. 计算最短距离（例如使用近邻分析工具）：

shortest_distances = "shortest_distances.shp"
arcpy.Near_analysis(houses, pollution_points, shortest_distances, "MIN", "NEW_GEOMETRY")

6. 结果将在`shortest_distances.shp`中，每一行对应住房点及其到最近污染点的距离。

记得在操作完成后删除临时文件：

arcpy.Delete_management("temp_houses")
arcpy.Delete_management("temp_pollution_points")
arcpy.Delete_management("temp_houses_projected")