你好，arcgis中如何将多个shp文件，比如点文件进行合并，怎么用python编程啊

你是想学PYTHON编程还是只是想只得到这个PYTHON文件。可以给你提供一条简洁的途径用modelbuilder来实现，将多个SHP文件拖入进去，再把MERGE工具拖进去，双击modelbuilder中的merge工具框设置，再双击output

dataset框设置输出。然后将这些shp文件用倒数第二个按钮添加链接的功能将他们一个个与merge工具框链接起来。最后点击model-export-to

srcipt-python

就会输出一个python文件，可以用记事本打开查看里面的代码。

要计算ArcGIS面要素在某一方向上的最长距离，可以使用以下步骤：

使用arcpy库中的Polygon对象读取面要素。

计算要截取的平行线的起点和终点坐标，可以使用math库中的sin和cos函数计算。

使用arcpy库中的Polyline对象创建平行线。

使用arcpy库中的Intersect_analysis函数计算面要素和平行线的交点。

使用math库中的sqrt函数计算每个交点到起点的距离。

找到最长的距离并返回。

下面是一个示例代码：

import arcpyimport math

# 读取面要素polygon = arcpy.Polygon(arcpy.Array([arcpy.Point(*coords) for coords in polygon_coords]))

# 角度

angle = 112# 起点坐标

start_x, start_y = polygon.centroid.X, polygon.centroid.Y

# 计算终点坐标我爱编程网

end_x = start_x + math.cos(math.radians(angle)) * 1000 # 1000是平行线长度

end_y = start_y + math.sin(math.radians(angle)) * 1000# 创建平行线

line_array = arcpy.Array([arcpy.Point(start_x, start_y), arcpy.Point(end_x, end_y)])line = arcpy.Polyline(line_array)

# 计算交点

intersect_points = arcpy.Intersect_analysis([polygon, line], arcpy.Geometry())

# 计算每个交点到起点的距离

distances = [math.sqrt((point.X - start_x) ** 2 + (point.Y - start_y) ** 2) for point in intersect_points]

# 找到最长距离

max_distance = max(distances)

print('最长距离:', max_distance)

注意，上述代码中的polygon_coords是一个包含面要素坐标的列表。另外，代码中使用了ArcGIS的arcpy库，需要安装ArcGIS软件并配置Python环境才能正常运行。

1、点类型的转栅格：

PointToRaster example 1 (Python

window)

Converts point features to a raster dataset.

import arcpy

from arcpy import env

env.workspace = "c:/data"

arcpy.PointToRaster_conversion("ca_ozone_pts.shp", "ELEVATION",

"c:/output/ca_elev", "MAXIMUM", "", 2000)

PointToRaster example 2

(stand-alone script)

Converts point features to a raster dataset.

# Name: PointToRaster_Ex_02.py

# Description: Converts point features to a raster dataset.

# Requirements: ArcInfo

# Import system modules

import arcpy

from arcpy import env

# Set environment settings

env.workspace = "C:/data"

# Set local variables

inFeatures = "ca_ozone_pts.shp"

valField = "ELEVATION"

outRaster = "c:/output/ca_elev02"

assignmentType = "MAXIMUM"

priorityField = ""

cellSize = 2000

# Execute PointToRaster

arcpy.PointToRaster_conversion(inFeatures, valField, outRaster,

assignmentType, priorityField, cellSize)

2、面类型的转栅格：

PolygonToRaster example 1

(Python window)

Converts polygon features to a raster dataset.

import arcpy

from arcpy import env

env.workspace = "c:/data"

arcpy.PolygonToRaster_conversion("ca_counties.shp", "NAME",

"c:/output/ca_counties.img",

"MAXIMUM_AREA", "MALES", 0.25)

PolygonToRaster example 2

(stand-alone script)

Converts polygon features to a raster dataset.

# Name: PolygonToRaster_Ex_02.py

# Description: Converts polygon features to a raster dataset.

# Requirements: ArcInfo

# Import system modules

import arcpy

from arcpy import env

# Set environment settings

env.workspace = "C:/data"

# Set local variables

inFeatures = "ca_counties.shp"

valField = "NAME"

outRaster = "c:/output/ca_counties"

assignmentType = "MAXIMUM_AREA"

priorityField = "MALES"

cellSize = 0.5

# Execute PolygonToRaster

arcpy.PolygonToRaster_conversion(inFeatures, valField, outRaster,

assignmentType, priorityField, cellSize)

3、线类型转栅格

PolylineToRaster example 1

(Python window)

Converts polyline features to a raster dataset.

import arcpy

from arcpy import env

env.workspace = "c:/data"

arcpy.PolylineToRaster_conversion("roads.shp", "CLASS", "c:/output/roads.img",

"MAXIMUM_COMBINED_LENGTH", "LENGTH", 30)

PolylineToRaster example 2

(stand-alone script)

Converts polyline features to a raster dataset.

# Name: PolylineToRaster_Ex_02.py

# Description: Converts polyline features to a raster dataset.

# Requirements: ArcInfo

# Import system modules

import arcpy

from arcpy import env

# Set environment settings

env.workspace = "C:/data"

# Set local variables

inFeatures = "roads.shp"

valField = "CLASS"

outRaster = "c:/output/roads.tif"

assignmentType = "MAXIMUM_COMBINED_LENGTH"

priorityField = "LENGTH"

cellSize = 30

# Execute PolylineToRaster

arcpy.PolylineToRaster_conversion(inFeatures, valField, outRaster,

assignmentType, priorityField, cellSize)