QGIS 笔记三：利用pyQGIS进行shapefile坐标转换

QGIS 笔记三：利用pyQGIS进行shapefile坐标转换

源码github地址在此，记得点星：
https://github.com/brandonxiang/example-pyQGIS

前言

不断有人问我关于Arcpy的问题，但是我却真的一点都不了解Arcpy。最近有个task是将一个shapefile中的所有点的坐标从wgs84转换成为gcj。我决心试用一下Arcpy，结果发现非常难使用。接口非常不人性化，而且，文档也不够详细。使用基本上是靠模仿内置的python文件去写自己的脚本，路径是path\to\ArcGIS\Desktop10.2\ArcToolbox\Scripts。

然后，我又考虑使用fiona进行数据处理，但是python的编码问题又难到我了，python2.7出现乱码，3.4的也有问题。然后，我又考虑到使用gdal进行处理转换，转成geojson后，再做处理，但是编码的问题又出现了。这告诫我们不要将属性表的名字命名成中文的。

最终，我选择使用我比较熟悉的pyQGIS。

思路

QGIS的优点：

• 开源数据和插件的海量支持
• pyQGIS

参考pyqgis developer cookbook，它有个pdf版本，可以方便搜索。

流程图

WGS84的shapefile文件转换成为geojson，然后利用geojson-python-utils把坐标转换成为国测局坐标。最后转为对应的shapefile。

实现

那么坐标转换使用的是接口QgsVectorFileWriter中的writeAsVectorFormat输入的几个参数分别是

• vlayer 输入的矢量图层
• input_json输出的矢量对象
• utf-8 编码
• GeoJSON 输出的数据类型

error=QW.writeAsVectorFormat(vlayer, input_json,'utf-8',None,'GeoJSON')

函数将会返回error是一个int标识。如果它等于QW.NoError，即为0，输出正确。然后对数据进行转换。读取对应的geojson数据将geometry取出进行转换，保存成为另一个geojson文件。

def convert(input, output):
    with open(input, 'r') as file:
        data = file.read()
        geojson = json.loads(data)
        features = geojson['features']
        for feature in features:
            feature['geometry'] = wgs2gcj(feature['geometry'])
        with open(output, 'w') as wfile:
            wfile.write(json.dumps(geojson))

最终将geojson转回shapefile即可。