2021-09-09

最新推荐文章于 2022-06-01 11:35:32 发布

肥宝欣

最新推荐文章于 2022-06-01 11:35:32 发布

阅读量397

点赞数 3

文章标签： python

本文链接：https://blog.youkuaiyun.com/m0_46313004/article/details/120196787

版权

该博客介绍了如何通过Python的QGIS库（PyQGIS）解决大规模矢量图层与底图之间的偏移问题。作者通过循环遍历图层中的每个要素并调整其顶点坐标，实现图层的整体平移，避免了顶点编辑导致的死机问题。代码中运用了QgsVectorLayerEditUtils的moveVertex方法进行坐标修改，确保了效率和准确性。最终展示了平移前后对比效果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

@运用pyQgis实现矢量图层平移

原因

之前在网上下载了矢量图层。但系下载的矢量图层跟我的底图存在一定的偏移，如果运用顶点编辑的话，逻辑上是可以的，奈何文件太大，顶点太多，一不留神就死机了，基于上述原因就运用pyQgis脚本文件来实现。

开始吧

#1、运用iface.activeLayer方法获取当前的活动图层。
activelayer = iface.activeLayer()
print(activelayer.name())
#2、使用循环获取图层里面的每一个要素，我这里的要素是MultiPolygon
for f in activelayer.getFeatures():
	#3、设置A局部变量是为了遍历顶点的编号，顶点的编号根据帮助文档说是从0开始的。
    a=0
    #4、运用f.geometry().vertices()遍历每个要素的顶点，因为到时要调用x()\y()来获取顶点修改后的坐标。
    for i in f.geometry().vertices():
   #5、moveVertex方法输入4个要素：要修改为怎样的X坐标。要修改为怎样的Y坐标。要编辑的要素ID。要编辑的顶点ID。	
   		QgsVectorLayerEditUtils(activelayer).moveVertex(i.x()+10000,i.y()-352,f.id(),a)
		a=a+1

代码1-3就不说了，主要是代码4、5：
其中代码4中： f.geometry().vertices()。f的类型是QgsFeature。然后运用geometry() 方法，返回了一个QgsGeometry方法，再调用 vertices()方法，得出 QgsVertexIterator对象，然后就可以使用for循环遍历了。遍历可以返回每个顶点的QgsPoint对象，就可以调用其中的x()\y()方法。

代码5中：class qgis.core.QgsVectorLayerEditUtils(layer: QgsVectorLayer)，主要调用moveVertex(self, x: float, y: float, atFeatureId: int, atVertex: int)→ bool这个方法。
其实可以运用下QgsGeometry对象中的”moveVertex(self, x: float, y: float, atVertex: int)→ bool“方法，按道理都是可以的。 x(),y()需要加多少，看自己偏移的需要了。

脚本怎么用不说了。
上效果：
前在这里插入图片描述
后