PostGIS 入门：PostgreSQL 中的空间查询

在处理海量空间地理数据时，传统的关系型数据库往往力不从心。而 PostgreSQL 配合其著名的空间数据库扩展 PostGIS，则摇身一变，成为了目前业界最强大的开源空间数据库。本文将为您介绍 PostGIS 的核心概念，并展示几个实用的 SQL 空间查询示例。

1. 什么是 PostGIS？

PostGIS 是 PostgreSQL 的空间数据库扩展。它引入了新的空间数据类型（如 Geometry 和 Geography），并提供了数百个用于空间数据分析 and 处理的 SQL 函数。

2. PostGIS 的核心数据类型

• Geometry（几何）：使用平面坐标系。适用于局部的、平面的空间计算。
• Geography（地理）：使用球面坐标系（经纬度）。自动基于地球曲率计算距离，计算精度高，但计算开销大。

3. 常用 SQL 空间查询示例

3.1 创建表

CREATE TABLE points_of_interest (
    id SERIAL PRIMARY KEY,
    name VARCHAR(100),
    geom GEOMETRY(Point, 4326) -- WGS84 坐标系
);

3.2 插入地理空间数据

INSERT INTO points_of_interest (name, geom) 
VALUES ('北京天安门', ST_GeomFromText('POINT(116.4074 39.9042)', 4326));

3.3 计算两个空间要素之间的距离

SELECT ST_Distance(
    ST_GeogFromText('POINT(116.4074 39.9042)'), 
    ST_GeogFromText('POINT(116.3970 40.0156)')
) AS distance_meters;

3.4 范围查询

-- 寻找距离天安门 5000米以内的所有点
SELECT name FROM points_of_interest
WHERE ST_DWithin(geom::geography, ST_GeogFromText('POINT(116.4074 39.9042)'), 5000);