Date: 2014-12-18 21:09
Summary: 本文記錄了使用Python進(jìn)行QGIS二次開發(fā)的過(guò)程，以及過(guò)程中遇到的問(wèn)題。由于內(nèi)容較多，將文章分成2個(gè)部分。這是第一部分。
使用Python進(jìn)行QGIS二次開發(fā)2
3. 自定義地圖工具
在做GPX文件編輯工具的時(shí)候，主要存在幾個(gè)難點(diǎn)：
- gpx文件是不能直接編輯的，必須先轉(zhuǎn)成shp格式才可以編輯。
- gpx中存儲(chǔ)的只是經(jīng)緯度、高度、時(shí)間戳信息，沒(méi)有投影信息。
解決方法：
之前本人曾做過(guò)GPS軌跡可視化的小系統(tǒng)，其中有一個(gè)地方就是gpx轉(zhuǎn)shp格式。當(dāng)時(shí)使用了一個(gè)gpx2shp的工具（主頁(yè)在這里），但是試了一下發(fā)現(xiàn)，這個(gè)工具把GPX點(diǎn)轉(zhuǎn)成了線類型，我們需要的是點(diǎn)類型，所以不能使用了。干脆直接解析gpx文件提取點(diǎn)，然后轉(zhuǎn)成shp，難度也不大。github上找到了gpxpy這個(gè)庫(kù)，試了一下還是很好用的。用這個(gè)庫(kù)解析gpx里面的軌跡點(diǎn)之后，用pyshp（https://github.com/GeospatialPython/pyshp）把點(diǎn)寫成shp文件。
用Proj4（http://trac.osgeo.org/proj/）這個(gè)庫(kù)把經(jīng)緯度點(diǎn)轉(zhuǎn)換成需要的投影點(diǎn)，然后才能和柵格圖匹配上。可以使用pip安裝pyproj，當(dāng)前是1.9.3的版本。具體的用法很簡(jiǎn)單，可以參考這篇帖子，具體步驟就是：
-
先定義一個(gè)需要的投影
# 從proj4字符串創(chuàng)建投影 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代碼創(chuàng)建投影 wgs84=pyproj.Proj("+init=EPSG:4326") -
對(duì)經(jīng)緯度進(jìn)行轉(zhuǎn)換
x, y = wgs84(p.longitude, p.latitude)
QGIS中編輯矢量
在QGIS中,大多數(shù)矢量數(shù)據(jù)支持編輯，有兩種編輯模式：
-
通過(guò)下面的函數(shù)來(lái)查看矢量圖層是否支持編輯操作：
caps = layer.dataProvider().capabilities()使用dataProvider所提供的函數(shù)可以使編輯直接生效。具體的方法詳見PyQGISDeveloperCookbook的4.2節(jié)。
-
使用編輯緩沖區(qū)
在編輯矢量圖層的時(shí)候，所做的改變可以到需要提交的時(shí)候一起提交，方便進(jìn)行撤銷等操作。具體的方法詳見PyQGISDeveloperCookbook的4.3節(jié)。
MapCanvas中渲染臨時(shí)點(diǎn)線面對(duì)象
要在canvas中顯示一些額外的數(shù)據(jù)可以使用canvas item，有兩個(gè)很有用的canvas item，詳見PyQGISDeveloperCookbook的7.3節(jié)：
-
QgsRubberBand
用來(lái)畫線： r = QgsRubberBand(canvas, False) # False = not a polygon points = [ QgsPoint(-1, -1), QgsPoint(0, 1), QgsPoint(1, -1) ] r. setToGeometry(QgsGeometry. fromPolyline(points), None) 用來(lái)畫多邊形： r = QgsRubberBand(canvas, True) # True = a polygon points = [ [ QgsPoint(-1, -1), QgsPoint(0, 1), QgsPoint(1, -1) ] ] r. setToGeometry(QgsGeometry. fromPolygon(points), None) -
QgsVertexMarker
用來(lái)畫點(diǎn) m = QgsVertexMarker(canvas) m. setCenter(QgsPoint(0, 0))
遇到的坑：
- 在把圖層注冊(cè)到QgsMapLayerRegistry的時(shí)候，這個(gè)管理器會(huì)給圖層重命名！這可是個(gè)大坑。本來(lái)以為加入到里面的圖層能夠根據(jù)名稱找到的，結(jié)果調(diào)試了半天才發(fā)現(xiàn)是這個(gè)問(wèn)題。它重命名的規(guī)則是“圖層名+時(shí)間戳”，在圖層名部分還會(huì)把小數(shù)點(diǎn)(.)、短橫（-）、空格等字符給轉(zhuǎn)換成下劃線（_）。
4. 使用py2exe打包
下面代碼用來(lái)將QGIS和PyQt打包，在命令行中運(yùn)行py2exe，打開命令行到工程文件目錄，運(yùn)行python setup.py py2exe。這里需要注意的是：
一定要保證當(dāng)前目錄中有msvcm90.dll、msvcp90.dll、msvcr90.dll這三個(gè)dll和VC80.CRT的msvcm80.dll、msvcp80.dll、msvcr80.dll三個(gè)文件！
-
"PyQt4.QtNetwork", "PyQt4.QtXml"這兩項(xiàng)要加進(jìn)去，不然會(huì)在運(yùn)行的時(shí)候說(shuō)找不到庫(kù)！
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 } )
代碼
參考資料
- http://geospatialdesktop.com/2009/02/creating_a_standalone_gis_application_1/
- http://docs.qgis.org/2.2/en/docs/pyqgis_developer_cookbook/
- http://qgis.org/api/2.2/
- http://gpx2shp.sourceforge.jp/
- https://github.com/tkrajina/gpxpy
- https://github.com/GeospatialPython/pyshp
- http://trac.osgeo.org/proj/
- http://all-geo.org/volcan01010/2012/11/change-coordinates-with-pyproj/
- http://www.kaij.org/blog/?p=105
