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即可。
