Arcpy处理修改shapefile FeatureClass 线要素坐标
参考资料:
1. 关于坐标转换计算: 在另一篇文章中有参考资料 http://www.cnblogs.com/beautifulplanet/p/4309222.html
2. 如何用arcpy进行处理呢,先找到这个链接
http://gis.stackexchange.com/questions/17096/edit-end-points-in-polyline-python-arcmap10,
从def offsetFirstPointInLine(line_geom,X_distance,Y_distance)函数
中找到如何提取每一条线记录的每个点坐标:
- geom = r.getValue("SHAPE")
- array = geom.getPart(0)
那array 是什么呢?
3. 上面那个例子只处理了每条线的一个点,而我需要对每个点进行处理,那怎么构成其中的new line呢?又去找ArcGIS的帮助,找到PolyLine类
http://help.arcgis.com/zh-cn/arcgisdesktop/10.0/help/index.html#/na/000v000000n2000000/
看看其中polyline是如何由点组装起来的。
解决方案:
经过上面三步梳理,整理代码如下
- >>> import arcpy,math
- ... pi = 3.14159265358979324
- ... a = 6378245.0
- ... ee = 0.00669342162296594323
- ... def transformLat(x,y):
- ... ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * math.sqrt(abs(x))
- ... ret += (20.0 * math.sin(6.0 * x * pi) + 20.0 * math.sin(2.0 * x * pi)) * 2.0 / 3.0
- ... ret += (20.0 * math.sin(y * pi) + 40.0 * math.sin(y / 3.0 * pi)) * 2.0 / 3.0
- ... ret += (160.0 * math.sin(y / 12.0 * pi) + 320 * math.sin(y * pi / 30.0)) * 2.0 / 3.0
- ... return ret
- ... def transformLon(x,y):
- ... ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * math.sqrt(abs(x))
- ... ret += (20.0 * math.sin(6.0 * x * pi) + 20.0 * math.sin(2.0 * x * pi)) * 2.0 / 3.0
- ... ret += (20.0 * math.sin(x * pi) + 40.0 * math.sin(x / 3.0 * pi)) * 2.0 / 3.0
- ... ret += (150.0 * math.sin(x / 12.0 * pi) + 300.0 * math.sin(x / 30.0 * pi)) * 2.0 / 3.0
- ... return ret
- ... def offsetPoint(wgLat,wgLon):
- ... dLat = transformLat(wgLon - 105.0, wgLat - 35.0)
- ... dLon = transformLon(wgLon - 105.0, wgLat - 35.0)
- ... radLat = wgLat / 180.0 * pi
- ... magic = math.sin(radLat)
- ... magic = 1 - ee * magic * magic
- ... sqrtMagic = math.sqrt(magic);
- ... dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi)
- ... dLon = (dLon * 180.0) / (a / sqrtMagic * math.cos(radLat) * pi)
- ... new_point = arcpy.Point(wgLon+dLon, wgLat+dLat)
- ... return new_point
- ... def offsetFirstPointInLine(line_geom):
- ... new_point = arcpy.Point()
- ... new_array = arcpy.Array()
- ... array = line_geom.getPart(0)
- ... for x in range(0,array.count):
- ... old_point = array[x]
- ... new_point = offsetPoint(old_point.Y,old_point.X)
- ... new_array.add(new_point)
- ... new_line = arcpy.Polyline(new_array,SR)
- ... new_array.removeAll()
- ... return new_line
- ... fc = r"E:\movedata-lonlat\guodao.shp"
- ... cur= arcpy.UpdateCursor(fc)
- ... SR= arcpy.Describe(fc).spatialReference
- ... print SR
- ... for r in cur:
- ... geom = r.getValue("SHAPE")
- ... r.setValue("SHAPE",offsetFirstPointInLine(geom))
- ... cur.updateRow(r)
- ... del r,cur
将处理后的数据再进行Project : WGS 1984 Web Mercator.prj,与谷歌地图进行叠加显示位置一致。
有几句代码标红,是遇到数据缺失的问题后修改的,过程如下:
遇到的问题:
但是...后来放大数据,才发现有缺失,如图(紫色是原始数据,红色是处理后得到的数据)

new_array.add(new_point)
new_line = arcpy.Polyline(new_array)
打印new_array也是38个点,但是打印new_line.pointCount就只有3个点了。
重新分析测试数据(WGS-84 地理坐标系/经纬度 下进行),启用编辑状态,原始的线要素上点比较密集
用 Arcpy Polyline lose points 去谷歌 搜到 http://stackoverflow.com/questions/14248618/shape-information-lost-when-using-arcpy-polyline-object-as-dictionary-value
好像是要添加坐标系,另外又进行了实验:
给这条线加上几个点,其中有几个点距离比较近:
再次处理,线变成了:
推测是不是在构成线的时候有什么距离限制?
用ArcPy Polyline accuracy、 ArcPy Polyline tolerance 谷歌,搜到了http://gis.stackexchange.com/questions/86728/minimum-shapelength-using-arcpy-polyline
其中关键的一句 “ I recommend setting the Geometry's Spatial Reference before inputting the coordinates. The SR has its own XYTolerance and XYResolution”
因此对上面的程序进行修改,增加红色部分。
Arcpy处理修改shapefile FeatureClass 线要素坐标的更多相关文章
- arcpy arcgis python实例教程--原点夹角距离定义线(坐标正算)
arcpy arcgis python实例教程--原点夹角距离定义线(坐标正算) 商务合作,科技咨询,版权转让:向日葵,135-4855__4328,xiexiaokui#qq.com 此地理处理工具 ...
- ArcEngine中合并断开的线要素(根据几何判断)
在上一篇ArcEngine环境下合并断开的线要素(根据属性)随笔中介绍了如何通过shp文件属性表中相同的属性字段进行线要素的合并.今天刚把通过几何条件判断的方式连接断开的线要素的ArcGIS 插件完成 ...
- ArcEngine环境下合并断开的线要素(根据属性)
1.遇到的问题: 最近遇到在线要素(矢量数据)中,一条完整的道路.河流等往往是断开的,如下图1所示: 2.思路: 在ArcGIS Desktop中没有相关的工具可以将这些断开的线要素进行自动合并,今天 ...
- geotools修改shapefile 属性名乱码问题
在GeoServer中文社区的讨论地址为:http://opengeo.cn/bbs/read.php?tid=1701&page=e&#a 使用geotools修改shapefile ...
- ArcMAp对线要素进行平滑处(打断)
一:工具简单介绍 -- ArcMAp10.1的高级编辑工具中提供了对线/面要素进行概括/平滑处理的工具. 概括工具.平滑工具分别例如以下:(首先得开启编辑状态 --- 才干够对要素的属性进行更改).选 ...
- ArcGIS 网络分析[1.4] 制作点线要素时需要注意的地方
有很多同学虽然成功做好了网络数据集,但是分析时会出现这样的问题: 这是为什么呢? 这有三个可能的原因: 1. 两个点之间所有的道路没有连通(问题出在点.线数据上). 2. 网络数据集出现了孤立的点位置 ...
- Shapefile文件中的坐标绘制到屏幕时的映射模式设置
pDC->SetMapMode(MM_ANISOTROPIC ); //首先选择MM_ANISOTROPIC映射模式,其它映射模式都不合适 pDC->SetWindowExt( max(a ...
- arcgis 面或线要素类上的搜索游标
import arcpy infc = arcpy.GetParameterAsText(0) # Identify the geometry field # desc = arcpy.Describ ...
- ArcEngine 创建线要素图层
在创建要素图层的时候,默认的几何类型是Polygon: Dim objectClassDescription As IObjectClassDescription = New FeatureClass ...
随机推荐
- with和上下文管理器
with和上下文管理器 如果你有时间阅读源码的习惯,可能会看到一些优秀的代码会出现带有with关键字的语句. 对于系统资源如文件,数据库连接,socket而言,应用程序打开这些资源并执行完业务逻辑之后 ...
- Linux高性能服务器编程:高性能服务器程序框架
服务器有三个主要模块: (1)I/O处理单元 (2)逻辑单元 (3)存储单元 1.服务器模型 C/S模型 逻辑:服务器启动后,首先创建一个或多个监听socket,并调用bind函数将其绑定到服务器感兴 ...
- import torch 报错
1.进入官网 https://pytorch.org/ 2.复制command到anaconda环境,即可
- 吴裕雄--天生自然MySQL学习笔记:MySQL 及 SQL 注入
如果通过网页获取用户输入的数据并将其插入一个MySQL数据库,那么就有可能发生SQL注入安全的问题. 本章节将为大家介绍如何防止SQL注入,并通过脚本来过滤SQL中注入的字符. 所谓SQL注入,就是通 ...
- one_day_one_linuxCmd---netstat命令
<坚持每天学习一个 linux 命令,今天我们来学习 netstat 命令> 前言:netstat 命令用于显示各种网络相关信息,如网络连接,路由表,接口状态等,还可以很方便查看本地机器上 ...
- java后台开发细节记录
1. ResultMap是程序员控制SQL查询结果和实体类的映射关系,而不是sql语句中字段的重命名,所以在sql语句中还是要按照原来字段的格式进行书写.
- APP-计算器
在网课上学习了计算器的制作方法,并自己做了小常识,看完一便后,感觉自己会了,思想也明白了,但是当关掉视频真正开始做的时候会发向许多的问题,自己在前的好多语法用的并不是很熟练,但是反复的查阅资料,观看教 ...
- 1811 06 pygame 的继续开发
早上看了 高数和python 好像系统没有保存 桑心啊 关于游戏背景的制作 游戏背景就是 背景在移动 而主人物还在原位置的 常常用于跑酷游戏类 背景开始绘制两张图像 一张完全 ...
- 基于Token的身份验证
最近了解下基于 Token 的身份验证,跟大伙分享下.很多大型网站也都在用,比如 Facebook,Twitter,Google+,Github 等等,比起传统的身份验证方法,Token 扩展性更强, ...
- LeetCode——456.132模式
给定一个整数序列:a1, a2, ..., an,一个132模式的子序列 ai, aj, ak 被定义为:当 i < j < k 时,ai < ak < aj.设计一个算法,当 ...