arcpy使用SelectLayerByLocation进行空间位置选择、SelectLayerByAttribute进行属性选择，使用CopyFeatures_management进行选择要素导出。

最新推荐文章于 2024-06-11 14:58:46 发布
原创 最新推荐文章于 2024-06-11 14:58:46 发布 · 3.5k 阅读 · 17 ·
CC 4.0 BY-SA版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
文章标签：

#python #数据库 #算法

该博客介绍了如何利用arcpy库在Python中进行地理国情监测的数据处理。首先，通过MakeFeatureLayer_management创建图层并使用SelectLayerByLocation选择位于指定区域内的要素。接着，用SelectLayerByAttribute进行属性选择，并反转选择。最后，当有匹配的要素时，使用CopyFeatures_management将选择结果导出到新的文件夹。这是一个GIS数据处理和分析的示例。

（地理国情监测）arcpy使用SelectLayerByLocation进行空间位置选择、SelectLayerByAttribute进行属性选择，使用CopyFeatures_management进行选择要素导出。

coding = 'utf-8'
import os #加载os库
import arcpy #加载arcpy库

file_dir_Source ="E:\\001.yangdebin\\001.集水区\\国情监测水网数据\\test.gdb\\CHydDataset\\UV_HYDL" 
file_dir_BenDi1 ="E:\\001.yangdebin\\001.集水区\\JC19_340000备份.gdb\\CHydDataset\\CV_HYDL" 

# Make a layer and select cities which overlap the chihuahua polygon
arcpy.MakeFeatureLayer_management(file_dir_Source
最低0.47元/天 解锁文章
利用python (Arcpy) 实现ArcGIS中SelectLayerByLocation的批量裁剪
midori99的博客
01-29 2163
在ArcGIS中我们经常会使用Select By Location功能来裁剪矢量。 例如，要实现根据美国各个州的图层去裁剪一个全美的图层，得到各个州分别的结果。然而美国有50个州，常规做法是使用Select By Location五十次才能实现，操作非常麻烦。 可以用ArcGIS自带的python窗口编写一段简单的代码实现批量裁剪，解放双手。 这是官方对SelectLayerByLocation的解释说明：https://desktop.arcgis.com/zh-cn/arcmap/10.3
在ArcEngine实现ArcMap中的SelectByLocation拓扑选择功能
03-16
在ArcMap中有Select by location这一功能，即根据feature之间 位置的拓扑关系 来选择 要素，如相交，相切，包含等，本程序在ArcEngine中实现了这一功能，本人调试过，能够运行。跟大家分享一下！
Python之Arcpy按位置选择批处理(SelectLayerByLocation)
toto+的博客
03-31 7619
工作中，往常使用FME批量压盖选取不同城市的Link、Node数据（存储于GDB地理数据库），能否通过Python引用Arcpy顺利实现上述功能？ 已知Arcgis中的手动操作步骤： 1、使用城市范围tile.shp按位置选择Link.shp，导出城市范围选择出的Link.shp 2、使用城市范围选择出的Link.shp，按位置选择导出本城市的Node.shp 记录下尝试过程及思考： ...
arcpy SelectByLocation_management 按位置选择
A873054267的博客
12-11 2758
import arcpy from arcpy import env sourceGDB=r"d:\polygon.gdb" env.workspace=sourceGDB target="merge" out_feature_class="p2" arcpy.MakeFeatureLayer_management(out_feature_class, 'newlinks_lyr') ar...
关于arcpy如何针对每一条数据做Select By Location
weixin_48719457的博客
06-11 326
注意：在searchcursor中需要加["SHAPE@"]获得图形，然后再用遍历到的图形row[0] 做按位置选择。目标：要实现遍历feature1图层里的数据，并用按位置选择来筛选出与他100米以内的其他点位。
arcgis中Python脚本的使用
02-22
在处理矢量数据时，`arcpy`也提供了丰富的函数，如`CopyFeatures`用于复制图层，`Append`用于合并多个图层，`SelectLayerByAttribute`和`SelectLayerByLocation`用于选择满足特定条件或位置的要素。 除了基本的地理...
ARCGIS如何导出特定要素的shp文件
08-09
ArcGIS 提供了多个地理处理工具，如 `Select Layer By Attribute` 和 `Select Layer By Location`，可用于创建临时选择集，随后使用 `Copy Features` 或 `Feature Class to Shapefile` 工具将选中要素导出为新的 ...
没有达到出图效果：只显示索引图层中当前出图要素和当前出图要素范围里面（根据用户选择的重叠度或者缓冲距离）的要素和底图图层，其余的全部设为不可见，底图图层里面的要素始终显示（即图层开关自动打开显示）
最新发布
11-28
因此，我们采用另一种方法：在创建临时图层时，我们使用MakeFeatureLayer创建一个图层，然后使用SelectLayerByLocation选择要素，最后通过设置图层的显示为“仅显示所选要素”来实现。但是，ArcPy中并没有直接设置...
出现旧的问题：索引要素当前出图要素范围里面的要素又不显示了
11-26
arcpy.SelectLayerByLocation_management( temp_lyr_name, "INTERSECT", index_geom ) # 获取选中的OID oid_field = arcpy.Describe(temp_lyr_name).OIDFieldName oids = [row[0] for row in arcpy.da....
arcpy程序最好不要依赖于arcmap运行程序MakeFeatureLayer_management，而要用CopyFeatures_management
Never Settle
03-21 1050
可以选择使用CopyFeatures_management函数先将生成的文件保存到本地，然后在执行其他函数时，函数输入参数不要写为arcmap里现有的图层，而是指定为已经在本地的文件的绝对路径，这样就不会在gis中执行函数，大大提高了运行速率。arcpy程序最好不要MakeFeatureLayer_management，这样会在arcmap中产生图层layer，再去调用类似相交分析的功能时，就会在gis中执行，减慢运行速率。
Arcpy中arcpy.management.SelectLayerByAttribute的用法
ZRong_H的博客
05-07 1311
在根据指定的属性条件从图层中选择要素时，我们需要用到arcpy.management.SelectLayer这个函数。该代码将从要素图层中选择满足人口数量大于等于10000的要素，并将选择结果保存为新的要素类。
arcpy实现空间查询_调用ArcPy查询和选择数据
weixin_39556590的博客
02-15 936
1. 按位置选择要素SelectLayerByLocation_management语法SelectLayerByLocation_management (in_layer, {overlap_type}, {select_features}, {search_distance}, {selection_type})测试主要代码如下：Exp1：import arcpyarcpy.env.works...
arcpy.SelectLayerByAttribute_management
darer49的博客
04-08 6184
用python在图层里建立一个新的要素选择：对于字段类型是数字：arcpy.SelectLayerByAttribute_management("pp","NEW_SELECTION",'"ID" = 1')对于字段类型是字符串：arcpy.SelectLayerByAttribute_management("国道","NEW_SELECTION",'"ROADCODE" = \'G107
arcpy 中SelectLayerByAttribute_management使用注意事项
weixin_30410999的博客
03-29 1133
http://blog.sina.com.cn/s/blog_54b09dc90102wgpi.html 转载于:https://www.cnblogs.com/aa123www456/p/10619812.html
arcpy学习随记（八）--- SelectLayerByAttribute函数的选择类型参数选项
nxy_918的专栏
02-04 786
ADD_TO_SELECTION：基于查询从当前选择集添加到当前的记录中（当存在一个选择内容时，会将生成的选择内容添加到现有选择内容中。SWITCH_SELECTION（反向选择）：切换选择内容，选择当前没有被选择的记录（将所选的所有记录从选择内容中移除，将未选取的所有记录添加到选择内容中，当指定该选项时将忽略表达式）（将生成的选择内容从现有选择内容中移除。SUBSET_SELECTION（两个选择的公共内容选择，交集）：结合现有选择集，选择两者共同记录（选取生成的选择内容与现有选择内容交集。
【ArcPy】根据空间条件选择要素，与创建图层问题
猴扯
05-20 3869
我的场景：有一个点要素集，我想知道距离每个点最近的几个点要素是什么 以前数据量不大的时候，这种情况我就直接暴力二层循环两两计算距离，但是随着自己对代码要求的提高，不能接受如此笨重的做法，遂查到空间选择SelectLayerByLocation_management。 这一方法输入的第一个参数（即查询对象）必须是图层，以往我都是要素集操作，不熟悉arcpy中的图层，因此也遇到了问题。 首先贴出...
根据属性选择要素(SelectLayerByAttribute)：使用Arcpy进行批量选择
baidu_41846932的博客
04-15 6063
根据属性选择要素(SelectLayerByAttribute)：使用Arcpy进行批量选择 在ArcMap中可以通过属性选择(SelectLayerByAttribute)根据某一/某些特定字段选择数据，但当数据较多时，手动选择查询条件较为复杂，此时可以在python窗口通过编写python脚本来进行批量选择： import arcpy query_List = [] for i in range(len(object_id_List)): if(i==0): query_List
python gdal实现arcpy的SelectLayerByAttribute_management函数
qq_20373723的博客
03-30 682
import os from osgeo import ogr, gdal shp_path = './1.shp' ds = ogr.Open(shp_path,0) layer = ds.GetLayer() layer.SetAttributeFilter("majority = 1 or majority = 2") driver = ogr.GetDriverByName('ESRI Shapefile') out_ds = driver.CreateDataSource('./shp/t.s
Arcpy 复制带有要素数据集的数据库
GIS on the road happy_h
09-08 4895
复制带有要素数据集的数据库，在复制要素数据集的时候要更改工作空间，xiande
C#中使用ArcPy选择要素导出为新图层
10-07
<think>我们注意到用户的问题是在C#中使用ArcPy选择要素并导出为新图层。但是，ArcPy是Python库，通常用于ArcGIS的Python脚本环境。在C#中，我们通常使用ArcObjects或者ArcGIS Runtime SDK for .NET来实现GIS功能。然而，用户要求使用ArcPy，这可能意味着需要在C#中调用Python脚本（包含ArcPy代码）来实现功能。 因此，我们可以这样解决： 1. 编写一个Python脚本，该脚本使用ArcPy选择要素并导出为新图层（例如，使用Select_analysis工具）。 2. 在C#中调用这个Python脚本。 步骤详细说明： 第一步：编写Python脚本（例如：select_and_export.py） 该脚本接收参数：输入要素路径、选择条件、输出要素路径。 第二步：在C#中通过Process类调用Python解释器执行该脚本，并传递参数。 但是请注意，这种方法需要在运行环境中安装Python和ArcPy（即安装ArcGIS Pro或ArcGIS Desktop），并且C#程序能够访问到Python环境。 另外，用户也可以考虑使用ArcGIS Engine或ArcGIS Runtime SDK for .NET中的选择功能，直接在C#中实现，而不通过ArcPy。但既然用户指定了ArcPy，我们就按照上述方法。 参考引用[1]提到Python窗口的优势，但这里我们是在C#中调用Python脚本，所以需要确保环境配置正确。 下面我们分别给出Python脚本示例和C#调用示例。 ### Python脚本示例（select_and_export.py） ```python import arcpy import sys # 获取参数 input_feature = sys.argv[1] # 输入要素路径 where_clause = sys.argv[2] # 选择条件 output_feature = sys.argv[3] # 输出要素路径 # 执行选择 arcpy.Select_analysis(input_feature, output_feature, where_clause) ``` 注意：在ArcPy中，选择工具为`Select_analysis`，它位于`arcpy.analysis`模块下，但通常可以直接使用`arcpy.Select_analysis`。 ### C#调用示例 在C#中，我们可以使用`System.Diagnostics.Process`来启动Python进程并执行脚本。 ```csharp using System.Diagnostics; public void ExecutePythonScript(string pythonExePath, string scriptPath, string arguments) { ProcessStartInfo start = new ProcessStartInfo(); start.FileName = pythonExePath; // 例如：@"C:\Python27\ArcGISPro\python.exe" start.Arguments = string.Format("{0} {1}", scriptPath, arguments); start.UseShellExecute = false; start.RedirectStandardOutput = true; start.RedirectStandardError = true; using (Process process = Process.Start(start)) { using (StreamReader reader = process.StandardOutput) { string result = reader.ReadToEnd(); Console.Write(result); } using (StreamReader reader = process.StandardError)) { string error = reader.ReadToEnd(); Console.Write(error); } } } // 调用示例 string pythonPath = @"C:\Program Files\ArcGIS\Pro\bin\Python\envs\arcgispro-py3\python.exe"; // 根据实际安装路径修改 string scriptPath = @"C:\Scripts\select_and_export.py"; string inputFeature = @"C:\data\input.shp"; string whereClause = "\"POPULATION\" > 1000000"; // 注意：在命令行传递参数时，引号需要转义 string outputFeature = @"C:\data\output.shp"; string arguments = string.Format("\"{0}\" \"{1}\" \"{2}\"", inputFeature, whereClause, outputFeature); ExecutePythonScript(pythonPath, scriptPath, arguments); ``` 注意： 1. 上面的C#代码中，`whereClause`参数传递时使用了转义引号。这是因为在命令行中，如果条件包含空格或特殊字符，需要用引号括起来。 2. 确保Python脚本的路径和Python解释器的路径正确。 3. 如果ArcPy模块不在Python的系统路径中，可能需要先设置环境变量。通常在ArcGIS Pro安装后，其自带的Python环境已经包含了ArcPy模块。 另外，引用[2]是关于获取要素图层的别名，引用[3]是一系列Arcpy实验，但与本问题关联不大。 如果用户希望更深入地集成，也可以考虑使用IronPython（在.NET环境中运行Python），但ArcPy依赖于ArcGIS的Python环境，而ArcGIS Pro当前使用Python3，IronPython支持Python2.7，且不支持ArcPy，因此不推荐。 所以，在C#中调用外部Python进程是较为可行的方案。 相关问题： 1. 如何在C#中调用Python脚本并传递参数？ 2. 除了使用ArcPy，在C#中如何直接使用ArcGIS Engine或ArcGIS Runtime SDK选择要素并导出？ 3. 如何确保C#程序调用的Python环境包含ArcPy？ </think>基于您的查询，以下是C#中调用ArcPy实现要素选择并导出为新图层的完整方案： ### 实现原理 在C#中需通过**进程调用**执行Python脚本（包含ArcPy代码），核心步骤分为： 1. 创建Python脚本（使用`Select_analysis`工具） 2. C#通过`System.Diagnostics.Process`启动Python进程 3. 参数传递与结果处理 --- ### 步骤1：编写Python脚本（SaveAsLayer.py） ```python import arcpy, sys # 接收C#传递的参数 input_fc = sys.argv[1] # 输入要素类路径 output_fc = sys.argv[2] # 输出图层路径 sql_query = sys.argv[3] # 选择条件 # 执行选择并导出 arcpy.Select_analysis( in_features=input_fc, out_feature_class=output_fc, where_clause=sql_query ) ``` --- ### 步骤2：C#调用代码 ```csharp using System.Diagnostics; public void ExportSelectedFeatures() { // 1. 设置路径参数 string pythonExe = @"C:\Python27\ArcGISPro\python.exe"; // ArcGIS Pro Python路径 string scriptPath = @"C:\GIS\SaveAsLayer.py"; // Python脚本路径 string inputFeature = @"C:\Data\Cities.shp"; // 输入要素 string outputLayer = @"C:\Output\SelectedCities.shp"; // 输出图层 string whereClause = "\"Population\" > 1000000"; // 选择条件（注意转义引号） // 2. 配置进程参数 ProcessStartInfo startInfo = new() { FileName = pythonExe, Arguments = $"\"{scriptPath}\" \"{inputFeature}\" \"{outputLayer}\" \"{whereClause}\"", UseShellExecute = false, RedirectStandardOutput = true, CreateNoWindow = true }; // 3. 执行并获取结果 using (Process process = Process.Start(startInfo)) { using (StreamReader reader = process.StandardOutput) { string result = reader.ReadToEnd(); Console.WriteLine("执行结果: " + result); } process.WaitForExit(); } } ``` --- ### 关键说明 1. **参数传递**： - `whereClause`需遵守ArcGIS SQL语法（如字段名带空格需用双引号） - C#字符串中的引号需转义（`\"`） 2. **路径要求**： - Python路径通常为：`C:\Program Files\ArcGIS\Pro\bin\Python\envs\arcgispro-py3\python.exe` - 文件路径避免包含空格或特殊字符 3. **错误处理**： - 添加`try-catch`捕获异常 - 检查`process.ExitCode`（0表示成功） --- ### 备选方案：直接使用ArcObjects 若需避免Python依赖，可用ArcObjects的`IFeatureClass`接口实现： ```csharp // 需引用ESRI.ArcGIS.Geodatabase IFeatureClass featureClass = ...; IQueryFilter queryFilter = new QueryFilterClass(); queryFilter.WhereClause = "Population > 1000000"; IFeatureClass outputClass = featureClass.Select(queryFilter); outputClass.ExportToShapefile(outputPath); ``` 此方法无需Python环境但需ArcEngine许可[^2]。 ---
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
表情包 代码片
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

木易GIS

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值