实现代码

import arcpy

from arcpy import env

import os

import ConversionUtils

infc1= ConversionUtils。gp。GetParameterAsText(0)

infc2= ConversionUtils。gp。GetParameterAsText(1)

out_workspace=ConversionUtils。gp。GetParameterAsText(2)

file_name=ConversionUtils。gp。GetParameterAsText(3)

arcpy。CreateFeatureclass_management(out_workspace,file_name,"POLYLINE")

cur3 =arcpy。da。InsertCursor(out_workspace+"/"+file_name,("FID","SHAPE@","ID"))

array = arcpy。Array()

i=0

cur1=arcpy。da。SearchCursor(infc1,["FID","SHAPE@XY"])

for row1 in cur1:来`自+优-尔^论:文,网www.chuibin.com +QQ752018766-

array。removeAll()

fid1=row1[0]

x1,y1=row1[1]

cur2=arcpy。da。SearchCursor(infc2,["FID","SHAPE@XY"])

for row2 in cur2:

fid2=row2[0]

x2,y2=row2[1]

if fid1==fid2:

array。add(arcpy。Point(x1,y1,0,0))

array。add(arcpy。Point(x2,y2,0,0))

cur3。insertRow((fid1,arcpy。Polyline(array),0))

if cur2:

del cur2

然后，将QuickBird影像、TM遥感影像和SRTM高程数据进行集处理，将坐标系统统一转换为高斯克吕格面坐标系统（带号为Beijing 1954 GK Zone 20N，单位为：米）。然后在ArcGIS里面的GeoReFerence模块中以TM遥感影像数据为目标栅格，通过设置控制点将SRTM高程数据、QuickBird影像和道路图层进行相对配准。根据研究目标区域建立矢量图层面，利用数据分析工具中的掩膜提取工具对QuickBird影像、TM遥感影像和SRTM高程数据进行裁剪，得到研究目标区域的删格集成数据。利用Buffer工具对北山公路矢量数据进行还原并对北山公路的滑坡影响距离（600米）进行缓冲区分析得出影响范围。利用Conversion工具箱里的栅格转点工具对遥感数据的七个波段进行数据转换，接着用Create TIN创建TIN，再用TIN to Raster工具转为删格数据，最后用Raster to Point工具进行栅格转点。