使用Python进行QGIS二次开发2

2016-08-17 本文已影响0人浪尖儿

Date: 2014-12-18 21:09
Summary: 本文记录了使用Python进行QGIS二次开发的过程，以及过程中遇到的问题。由于内容较多，将文章分成2个部分。这是第一部分。

使用Python进行QGIS二次开发2

3. 自定义地图工具

在做GPX文件编辑工具的时候，主要存在几个难点：

gpx文件是不能直接编辑的，必须先转成shp格式才可以编辑。
gpx中存储的只是经纬度、高度、时间戳信息，没有投影信息。

解决方法：

之前本人曾做过GPS轨迹可视化的小系统，其中有一个地方就是gpx转shp格式。当时使用了一个gpx2shp的工具（主页在这里），但是试了一下发现，这个工具把GPX点转成了线类型，我们需要的是点类型，所以不能使用了。干脆直接解析gpx文件提取点，然后转成shp，难度也不大。github上找到了gpxpy这个库，试了一下还是很好用的。用这个库解析gpx里面的轨迹点之后，用pyshp（https://github.com/GeospatialPython/pyshp）把点写成shp文件。
用Proj4（http://trac.osgeo.org/proj/）这个库把经纬度点转换成需要的投影点，然后才能和栅格图匹配上。可以使用pip安装pyproj，当前是1.9.3的版本。具体的用法很简单，可以参考这篇帖子，具体步骤就是：

先定义一个需要的投影

  # 从proj4字符串创建投影
  isn2004=pyproj.Proj("+proj=lcc +lat_1=64.25 +lat_2=65.75 +lat_0=65 +lon_0=-19 +x_0=1700000 +y_0=300000 +no_defs +a=6378137 +rf=298.257222101 +to_meter=1")
  
  # 使用EPSG代码创建投影
  wgs84=pyproj.Proj("+init=EPSG:4326")

对经纬度进行转换

  x, y = wgs84(p.longitude, p.latitude)

QGIS中编辑矢量

在QGIS中,大多数矢量数据支持编辑，有两种编辑模式：

通过下面的函数来查看矢量图层是否支持编辑操作：
```
 caps = layer.dataProvider().capabilities()
```
使用dataProvider所提供的函数可以使编辑直接生效。具体的方法详见PyQGISDeveloperCookbook的4.2节。
使用编辑缓冲区

在编辑矢量图层的时候，所做的改变可以到需要提交的时候一起提交，方便进行撤销等操作。具体的方法详见PyQGISDeveloperCookbook的4.3节。

MapCanvas中渲染临时点线面对象

要在canvas中显示一些额外的数据可以使用canvas item，有两个很有用的canvas item，详见PyQGISDeveloperCookbook的7.3节：

QgsRubberBand

  用来画线：      
  r = QgsRubberBand(canvas, False) # False = not a polygon
  points = [ QgsPoint(-1, -1), QgsPoint(0, 1), QgsPoint(1, -1) ]
  r. setToGeometry(QgsGeometry. fromPolyline(points), None)
  
  用来画多边形：
  r = QgsRubberBand(canvas, True) # True = a polygon
  points = [ [ QgsPoint(-1, -1), QgsPoint(0, 1), QgsPoint(1, -1) ] ]
  r. setToGeometry(QgsGeometry. fromPolygon(points), None)

QgsVertexMarker

  用来画点
  m = QgsVertexMarker(canvas)
  m. setCenter(QgsPoint(0, 0))

遇到的坑：

在把图层注册到QgsMapLayerRegistry的时候，这个管理器会给图层重命名！这可是个大坑。本来以为加入到里面的图层能够根据名称找到的，结果调试了半天才发现是这个问题。它重命名的规则是“图层名+时间戳”，在图层名部分还会把小数点(.)、短横（-）、空格等字符给转换成下划线（_）。

4. 使用py2exe打包

下面代码用来将QGIS和PyQt打包，在命令行中运行py2exe，打开命令行到工程文件目录，运行python setup.py py2exe。这里需要注意的是：

一定要保证当前目录中有msvcm90.dll、msvcp90.dll、msvcr90.dll这三个dll和VC80.CRT的msvcm80.dll、msvcp80.dll、msvcr80.dll三个文件！

"PyQt4.QtNetwork", "PyQt4.QtXml"这两项要加进去，不然会在运行的时候说找不到库！

 from distutils.core import setup
 import py2exe
 
 py2exe_options = {
     "includes":["sip", "PyQt4.QtCore", "PyQt4.QtGui",
                 "PyQt4.QtNetwork", "PyQt4.QtXml"],
 }
 
 setup(
     windows=["TrackEditer.py"],
     options={
         'py2exe':py2exe_options
     }
 )

代码

ShapeViewer
GPX文件编辑：https://github.com/dslwz2008/GPSTrackEditer

使用Python进行QGIS二次开发2

使用Python进行QGIS二次开发2

QGIS中编辑矢量

MapCanvas中渲染临时点线面对象

遇到的坑：

代码

参考资料

猜你喜欢

热点阅读