postgreSQL+postGIS+NodeJS 路径查询：从一无所知到略知一二

目录

一、软件准备

1.下载

2.安装

二、创建数据库，数据准备

1.创建空间库

2.导入数据

3.对数据进行路网拓扑数据计算处理

4.路径计算测试

三、NodeJS 后端代码实现

1.nodeJS 链接 postgreSQL

2. 后端实现

五、效果展示

1.查询效果

2. 效率测试

---

一、软件准备

1.下载

首先，在官网下载即可：

1) PostgreSQL，我用的windows 版 12.4：https://www.enterprisedb.com/downloads/postgres-postgresql-downloads；

2) PostGIS，记得下配套的版本，pg12:http://download.osgeo.org/postgis/windows/;

2.安装

1) PostgreSQL，按照默认配置安装即可，端口默认5432，自定义用户名/密码；

安装成功后打开pgAdmin 4。

2) PostGIS 

选择 Create spatial database，会初始化一个空间数据库：

选择安装目录，这里选择之前安装的PostgreSQL12 （务必）的目录：

安装成功后打开 PostGIS PostGIS Bundle 3 for PostgreSQL x64 12 Shapefile and DBF Loader Exporter

二、创建数据库，数据准备

1.创建空间库

1) sql 添加空间扩展：
切换刚才创建的数据库下，运行下面的命令即可：
CREATE EXTENSION postgis;
CREATE EXTENSION postgis_topology;
CREATE EXTENSION fuzzystrmatch;
CREATE EXTENSION postgis_tiger_geocoder;
CREATE EXTENSION address_standardizer;

2) 手动添加（二选一即可）：

2.导入数据

1) 打开postgis，连接数据库mapdata

add File ，选择shape 文件；

3.对数据进行路网拓扑数据计算处理

执行成功后，执行成功后会生产一个 vertices_pgr 的表，里面包含路网相交点的空间数据

alter table nanjingroad add column source int;
alter table nanjingroad add column target int;

create index road_source_idx on road("source");
create index road_target_idx on road("target");

ALTER TABLE nanjingroad ADD COLUMN length double precision;  
update road set length =st_length(geom);
ALTER TABLE public.nanjingroad ADD column rev_length numeric;
update public.nanjing_road set rev_length=length;

SELECT pgr_createTopology('road',0.00001, 'geom', 'gid');  

转换坐标系(如有报错)

UnhandledPromiseRejectionWarning: error: ST_Distance: Operation on mixed SRID geometries (MultiLineString, 0) != (Point, 4326)

sql 执行以下代码：

select UpdateGeometrySRID('public','nanjingroad', 'geom', 4326);

4.路径计算测试

例如：source：47461，target：42345

SELECT * FROM pgr_dijkstra('SELECT gid as id,source,target,length::float as cost,rev_length::float as reverse_cost FROM public.nanjingroad',47461,42345,true);

三、NodeJS 后端代码实现

1.nodeJS 链接 postgreSQL

1）创建连接池

var pgsql = require('pg');
var pgconnection =new pgsql.Pool({
  host: 'localhost',
  port: 5432,
  user: 'postgres',
  password: 'postgres',
  database: 'mapdata'
});

2. 后端实现

2）关键代码

/**
 * 求当前点最近的路段
 * @param {*} x 
 * @param {*} y 
 */
function getRoadByPoint(x, y) {

    return new Promise((resolve, reject) => {
        pgConnection.connect(function (err, connection) {
            if (err) {
                reject(err)
            } else {
                connection.query(`SELECT * FROM public.nanjingroad ORDER BY ST_Distance(geom,ST_GeometryFromText('POINT(${x} ${y})',4326)) LIMIT 1`, (err, rows) => {
                    if (err) {
                        reject(err)
                    } else {
                        resolve(rows)
                    }
                    // 结束会话
                    connection.release()
                })
            }
        })
    })
}

/**
 * 使用dijkstra 算法求最短路径
 * @param {*} source 
 * @param {*} target 
 */
function getDijkstraPath(source, target) {
    return new Promise((resolve, reject) => {
        pgConnection.connect(function (err, connection) {
            if (err) {
                reject(err)
            } else {
                let sql = `SELECT * from public.nanjingroad where gid in (SELECT edge FROM pgr_dijkstra('SELECT gid as id,source,target,length::float as cost,rev_length::float reverse_cost FROM public.nanjingroad',${source},${target},true))`;
                connection.query(sql, (err, rows) => {
                    if (err) {
                        reject(err)
                    } else {
                        resolve(rows)
                    }
                    // 结束会话
                    connection.release()
                })
            }
        })
    })
}

五、效果展示

1.查询效果

2. 效率测试

数据量：9w+

单次查询用时:500-600ms 左右

本文地址：https://blog.csdn.net/u012123612/article/details/108573859