JSTS介绍和功能简单示例

JSTS是一个符合OGC标准的，包含空间拓扑功能的JavaScript类库。和turf类似，区别在于，JSTS是通过原始JTS Java 源代码翻译来的，保留了JTSAPI；而turf来自MapBox。个人感觉JSTS的文档不如turf明了(一些网址还404了)，但是可以直接适用OpenLayers（以下简称OL）的几何对象，对JTS API熟悉的话在OL工程里可能会更方便。

JSTS的参考文档比较分散，GitHub库里有比较全面的简介，这个地址里有一些简单的示例，JTS API直接用的JavaAPI，看起来并不直观，另外使用了一些空间数据类型读取类代替JavaAPI里的IO部分。

所以若想比较完整的了解JSTS类库和基础使用，建议按照顺序，先浏览Github库的README.md，直接用谷歌浏览器的机翻大概看一遍，然后看一下简单示例，到具体的代码编写时，数据读写部分参考http://bjornharrtell.github.io/jsts/1.6.1/doc/index.html，空间功能部分参考JavaAPI。

注，以下个人应用示例用的是Vue2+TS，若你的环境不一样请自行修改

安装和调用

没啥好说的，使用npm/yarn啥的自行安装，或者引用在线地址
引用时，若使用npm包，由于jsts没有类型声明文件，固然有一些比较好的方法，我这里还是图省事直接用的require。调用代码如下

// 引入OL几何类
import Geometry from "ol/geom/Geometry";
import { LineString, LinearRing, MultiLineString, MultiPoint, MultiPolygon, Point, Polygon } from "ol/geom";// 引入jsts
const jsts = require("jsts/dist/jsts.min.js");// 实例化OL解析类
const OLParser = new jsts.io.OL3Parser();
// 注入OL几何对象
OLParser.inject(Point, LineString, LinearRing, Polygon, MultiPoint, MultiLineString, MultiPolygon);

功能示例

先定义了一个函数类型

type TopoFunc = (...params: any) => Geometry | null;

缓冲区

/*** 缓冲区* @param geom* @param len 缓冲区距离，单位米* @returns*/
const getBuffer: TopoFunc = (geom: Geometry, len: number) => {const jstsGeom = OLParser.read(geom);if (!jstsGeom.isValid()) {console.error("几何对象出现拓扑错误,请检查修复");return null;}const buffered = jstsGeom.buffer(len);return OLParser.write(buffered);
};

取交

/*** 取交* @param geom* @param geomB* @returns*/
const intersects: TopoFunc = (geom: Geometry, geomB: Geometry) => {const jstsGeom = OLParser.read(geom);const jstsGeomB = OLParser.read(geomB);if (!jstsGeom.isValid() || !jstsGeomB.isValid()) {console.error("几何对象出现拓扑错误,请检查修复");return null;}const difference = jstsGeom.intersection(jstsGeomB);return OLParser.write(difference);
};

取补

/*** 取geom中geomB的补集* @param geom* @param geomB* @returns*/
const getDifference: TopoFunc = (geom: Geometry, geomB: Geometry) => {const jstsGeom = OLParser.read(geom);const jstsGeomB = OLParser.read(geomB);if (!jstsGeom.isValid() || !jstsGeomB.isValid()) {console.error("几何对象出现拓扑错误,请检查修复");return null;}const difference = jstsGeom.difference(jstsGeomB);return OLParser.write(difference);
};

融合

/*** 融合* @param geom* @param geomB* @returns*/
const union: TopoFunc = (geom: Geometry, geomB: Geometry) => {const jstsGeom = OLParser.read(geom);const jstsGeomB = OLParser.read(geomB);if (!jstsGeom.isValid() || !jstsGeomB.isValid()) {console.error("几何对象出现拓扑错误,请检查修复");return null;}const difference = jstsGeom.union(jstsGeomB);return OLParser.write(difference);
};

对等差分

/*** 对等差分* @param geom* @param geomB* @returns*/
const symDifference: TopoFunc = (geom: Geometry, geomB: Geometry) => {const jstsGeom = OLParser.read(geom);const jstsGeomB = OLParser.read(geomB);if (!jstsGeom.isValid() || !jstsGeomB.isValid()) {console.error("几何对象出现拓扑错误,请检查修复");return null;}const difference = jstsGeom.symDifference(jstsGeomB);return OLParser.write(difference);
};

备注

工作中遇到个拓扑错误，多边形自相交

尝试使用GeometryPrecisionReducer降低精度，但还是有问题，就索性代码里直接抛出，手动修改数据来修复吧。GeometryPrecisionReducer的代码也附下，执行没啥问题，但是结果依然不好，感兴趣可以研究下

function topoFix(jstsGeom: any) {if (jstsGeom.isValid()) {return jstsGeom;}const reducer = new jsts.precision.GeometryPrecisionReducer();// 这里直接使用reduce方法会直接报错// 若使用reducer.fixPolygonalTopology会正常执行，但返回的几何不大对，后面使用会报一样的错return reducer.reduce(jstsGeom);
}

错误