[JavaScript知识库]node.js读写shp文件转换数据的坐标系

接上文：IDEA中配置运行node.js

一、安装所需的npm包

1.读写shp文件的js库gdal

安装命令

npm install gdal -g --save

node-gdal库地址：https://github.com/naturalatlas/node-gdal
文档地址：http://naturalatlas.github.io/node-gdal/classes/gdal.html

2.修改shp文件数据坐标系的js库gcoord

安装命令

npm install gcoord -g

gcoord库地址：https://github.com/hujiulong/gcoord

二、转换代码

两个包安装完成后就可以开始写代码啦，完整的转换代码如下：

const gcoord = require('gcoord');
const gdal = require("gdal");

//r只读  r+读取修改增加  w创建, ESRI Shapefile即shp格式的全称
var dataset = gdal.open("C:\\Users\\syc\\Desktop\\my_road_data\\sz_main_road_test\\mainRoad.shp", "r+", "ESRI Shapefile");
//获取shp文件的第一个图层，我们的这个地图只有一个图层所以layers.get(0)即可
var layer = dataset.layers.get(0);
//获取图层的的数据记录数，看看有没有正确读取到
console.log("number of features: " + layer.features.count());
var i = 0;
var sum = 0;
//遍历数据对每一条记录的坐标数据进行坐标系转换
layer.features.forEach((feature) => {

    //转换前
    // console.log(feature.fid);
    // console.log(feature.getGeometry().toJSON());
    // console.log(feature.fields.toJSON());

    //转换为JSON对象，用gcoord进行geom数据的坐标系转换
    var geojson = JSON.parse(feature.getGeometry().toJSON());
    gcoord.transform(geojson, gcoord.GCJ02, gcoord.WGS84);
    //转换后
    // console.log(geojson);

    //将转换后的数据保存回feature
    feature.setGeometry(gdal.Geometry.fromGeoJson(geojson));

    //将修改保存到shp文件中
    layer.features.set(feature.fid, feature);

    //对比发现保存后的数据点的经度的精度多出了一位，可能是保存浮点数过程中的一个小bug吧，不过没什么影响，尾数精度都十几位了
    // console.log(layer.features.get(feature.fid).getGeometry().toJSON());

    i++;
    sum++;
    if (i === 1000) {
        console.log("已完成数据转换" + sum + "条...............");
        i = 0;
    }
})
console.log("<-----------------转换完成，共转换完成"+ sum + "条数据------------------------>");

//测试时也可用迭代器方法进行迭代，更多操作请参考node-gdal使用文档
// var feature = layer.features.next();