node-gdal实践-读写shp数据

gdal用于读写矢量和栅格数据集，gdal常常使用c，c++做开发，也有部分其他语言扩展如py，node。但由于其他语言的官方例子不多，且存在很多问题，所以入门非常难。本文主要面向使用node语言做工具且对GIS软件有二次开发基础的同学，介绍使用gdal读写经典shp的入门教程。

一 node-gdal安装

Github地址：https://github.com/naturalatlas/node-gdal 使用nodejs安装： npm install gdal --save。 注意--save不能漏，笔者一开始未加，导致安装失败。

二 HelloWorld

api文档地址：http://naturalatlas.github.io/node-gdal/classes/gdal.html 由于资料很少，一切唯有官网api可供参考。本篇是个高级的helloword，用于初步描述如何使用node-gdal。下面贴出node的测试脚本gdal.js代码说明如何操作：

    var gdal=require('gdal'); 
	//r只读  r+读取修改增加  w可写
	var dataset = gdal.open("test.shp",'w');
	//获取数据集中第一个图层
	var layer = dataset.layers.get(0);
	//添加一个字段叫name，字符型
	layer.fields.add(new gdal.FieldDefn('name', gdal.OFTString));
	//构造一个feature
	var feature = new gdal.Feature(layer);
	//给字段赋值
	feature.fields.set('name', '测试');
	//构造图形并赋值
	var lineString = new gdal.LineString();
	lineString.points.add(new gdal.Point(118.5,32.1));
	lineString.points.add(new gdal.Point(118.6,32.4));
	feature.setGeometry(lineString);
	//图层添加构造的feature
	layer.features.add(feature);
	//将对shp图层的更改写入磁盘文件
	layer.flush();
	console.log("图形json串:"+layer.features.next().getGeometry().toJSON());
	console.log("图层中要素数量: " + layer.features.count());
	console.log("图层字段名称: " + layer.fields.getNames());
	console.log("边界范围: " + JSON.stringify(layer.getExtent()));
	console.log("图层坐标系: " + (layer.srs ? layer.srs.toWKT() : 'null'));

node执行输出如下：

E:\gdal>node gdal.js
图形json串:{ "type": "LineString", "coordinates": [ [ 118.5, 32.1 ], [ 118.6, 32.4 ] ] }
图层中要素数量: 1
图层字段名称: Id,name
边界范围: {"minX":0,"minY":0,"maxX":118.6,"maxY":32.4}
图层坐标系: GEOGCS["GCS_WGS_1984",DATUM["WGS_1984",SPHEROID["WGS_84",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]

在ArcMap中打开shp，效果如下： 注意： 如上文：var dataset = gdal.open("test.shp",'w');，w模式代表可写，否则增加字段，增加Feature会失败，提示 当前数据源是只读数据。

转载于:https://my.oschina.net/freegis/blog/912284