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 ...
随机推荐
- 牛客小白月赛18——Forsaken的三维数点
这个是一个简单题,不过因为想到比标程时间复杂度更低的方法就尝试了一下. 思路:虽然加点是三维数点,但是我们要求的是半径的大小,这样的话,就可以转变为一维的问题. 标程的解法是,用树状数组维护,然后二分 ...
- dp--P1439 最长公共子序列(LCS)
题目描述 给出1-n的两个排列P1和P2,求它们的最长公共子序列. 输入格式 第一行是一个数n, 接下来两行,每行为n个数,为自然数1-n的一个排列. 输出格式 一个数,即最长公共子序列的长度 找出两 ...
- Windows下MariaDB数据搬家问题
背景:公司买了一台服务器要将原来老服务器数据库数据导入新的服务器上,两台服务器环境如下 做好准备后开始实施,老数据库进行停库 找到MariaDB的安装目录下的整data目录进行拷贝,然后到新服务器进行 ...
- 吴裕雄--天生自然 JAVASCRIPT开发学习:Number 对象
<!DOCTYPE html> <html> <head> <meta charset="utf-8"> <title> ...
- ZJNU 2356 - 六学家
“选出来三个六学家,他们的编号是i,j,k,满足i<j<k,且a[k]=a[j]-a[i]” 所以输入第i个数a[i]时,直接让答案加上前i-1个数中能构成差值为a[i]的数量即可 然后让 ...
- [转载]matlab视频读取函数VideoReader
看到以前matlab中读取视频多 使用mmreader等(参考<matlab读取/播放视频的函数>),而现在matlab有一个专门的视频读取类VideoReader完成视频读取的功能. 相 ...
- CMake变量(提供信息的变量)
目录 CMAKE_VERSION CMAKE_MAJOR_VERSION CMAKE_MINOR_VERSION CMAKE_PATCH_VERSION CMAKE_TWEAK_VERSION CMA ...
- Python7DAY-基础语法.md
# python简介 python交互器的作用:主要是是为了调试代码.
- selector的使用,android:clickable="true"
<ImageView android:id="@+id/patrol_buzzer_btn" android:layout_width="80dp" an ...
- webgis笔记
3.8(02) .特点:由服务端进行数据管理 开源的GO sever WMS/WCS/WTS 1sever/2engine/3database/4standard 扩展的空间数据库,存矢量.栅格.直接 ...