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，效果如下：

效果图.png

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