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中,大多数矢量数据支持编辑，有两种编辑模式：

1. 通过下面的函数来查看矢量图层是否支持编辑操作：

    caps = layer.dataProvider().capabilities()
   

   使用dataProvider所提供的函数可以使编辑直接生效。具体的方法详见PyQGISDeveloperCookbook的4.2节。

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。这里需要注意的是：

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

2. "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

参考资料

1. http://geospatialdesktop.com/2009/02/creating_a_standalone_gis_application_1/
2. http://docs.qgis.org/2.2/en/docs/pyqgis_developer_cookbook/
3. http://qgis.org/api/2.2/
4. http://gpx2shp.sourceforge.jp/
5. https://github.com/tkrajina/gpxpy
6. https://github.com/GeospatialPython/pyshp
7. http://trac.osgeo.org/proj/
8. http://all-geo.org/volcan01010/2012/11/change-coordinates-with-pyproj/
9. http://www.kaij.org/blog/?p=105