PostGIS 小白入门详细指南

最新推荐文章于 2025-03-15 18:30:16 发布
最新推荐文章于 2025-03-15 18:30:16 发布
阅读量6.7k 收藏 46
点赞数 24
文章标签： arcgis postgresql 数据库
版权声明：本文为博主原创文章，遵循 CC 4.0 BY 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43905308/article/details/144523347
版权

PostGIS 是 PostgreSQL 数据库的空间扩展，为其添加了地理信息系统（GIS）的功能。无论您是地理信息系统（GIS）爱好者、开发者，还是数据库管理员，学习 PostGIS 都能极大地提升您处理空间数据的能力。本指南将从基础到实战，详细介绍 PostGIS 的安装、配置、操作及项目实践，帮助您快速掌握并应用 PostGIS。

目录

  1. 什么是 PostGIS？
  2. 安装 PostgreSQL 和 PostGIS
  3. PostGIS 基础概念
  4. 创建和管理空间数据库
  5. 常用的空间函数和操作
  6. 导入和导出空间数据
  7. 空间查询示例与案例分析
  8. 空间索引与性能优化
  9. 空间数据可视化
  10. 项目实践：构建一个简单的地理信息系统
  11. 常见问题与解决方法

什么是 PostGIS？

PostGIS 是 PostgreSQL 的一个扩展，它将 PostgreSQL 转变为一个功能强大的空间数据库，使其能够存储、查询和分析地理空间数据。PostGIS 遵循开放地理空间联盟（OGC）的标准，支持各种空间数据类型和操作，如点、线、多边形、缓冲区、空间连接等。

主要功能包括：

  • 存储各种空间数据类型（如几何类型和地理类型）
  • 空间索引，提高查询性能
  • 丰富的空间函数和操作
  • 支持投影和坐标系统转换

应用场景：

  • 城市规划与管理
  • 环境监测与分析
  • 物流与运输优化
  • 地图制作与可视化
  • 位置服务与导航

安装 PostgreSQL 和 PostGIS

1. 安装 PostgreSQL

PostGIS 依赖于 PostgreSQL，因此首先需要安装 PostgreSQL。

Windows:
  1. 访问 PostgreSQL 官方下载页面。
  2. 下载适用于 Windows 的安装程序（通常是通过 EnterpriseDB 提供的安装包）。
  3. 运行安装程序，按照提示完成安装。安装过程中可以选择安装 pgAdmin 和其他工具。
macOS:

使用 Homebrew 安装 PostgreSQL：

brew update
brew install postgresql

启动 PostgreSQL 服务：

brew services start postgresql
Linux:
Debian/Ubuntu:
sudo apt update
sudo apt install postgresql postgresql-contrib

初始化数据库并启动 PostgreSQL：

sudo systemctl enable postgresql
sudo systemctl start postgresql
CentOS/RHEL:
sudo yum install postgresql-server postgresql-contrib
sudo postgresql-setup initdb
sudo systemctl enable postgresql
sudo systemctl start postgresql

2. 安装 PostGIS

Windows:

在安装 PostgreSQL 时，可以选择安装 PostGIS 插件。确保在安装向导中选中 PostGIS 组件。

macOS:

使用 Homebrew 安装 PostGIS：

brew install postgis
Linux:
Debian/Ubuntu:
sudo apt install postgis postgresql-<version>-postgis-<version>

替换 <version> 为您的 PostgreSQL 版本号，如 12、13 等。

CentOS/RHEL:
sudo yum install postgis postgis-utils

3. 验证安装

  1. 登录 PostgreSQL：
psql -U postgres
  1. 创建测试数据库：
CREATE DATABASE test_postgis;
\c test_postgis
  1. 启用 PostGIS 扩展：
CREATE EXTENSION postgis;
CREATE EXTENSION postgis_topology;
  1. 查看 PostGIS 版本：
SELECT PostGIS_Full_Version();

预期结果：

返回 PostGIS 的详细版本信息，表明安装成功。

PostGIS 基础概念

空间数据类型

PostGIS 支持多种空间数据类型，主要包括：

  • GEOMETRY: 通用几何类型，支持各种几何对象，如点、线、多边形等。
  • GEOGRAPHY: 基于地理坐标系的类型，适用于地球表面的测量。
  • Raster: 栅格数据类型，用于存储栅格图像。

几何对象

  • POINT: 表示单个坐标点。
  • LINESTRING: 由一系列点组成的线。
  • POLYGON: 由封闭的线组成的多边形。
  • MULTIPOINT, MULTILINESTRING, MULTIPOLYGON: 分别表示多个点、线、多边形的集合。
  • GEOMETRYCOLLECTION: 含有不同类型几何对象的集合。

空间参考系统（SRS）

空间数据通常有一个参考坐标系统，用于定义数据的地理位置。常见的 SRS 有：

  • EPSG:4326: WGS 84，地理坐标系，使用经纬度表示。
  • EPSG:3857: Web Mercator，投影坐标系，常用于网页地图。

空间索引

PostGIS 使用 GiST（Generalized Search Tree）索引来提高空间查询的性能。空间索引允许数据库在执行空间查询时更快地定位相关数据。

投影与坐标转换

不同的空间参考系统适用于不同的应用场景。PostGIS 提供了强大的函数来进行坐标系统之间的转换，如 ST_Transform，确保空间数据的一致性和准确性。

创建和管理空间数据库

1. 创建数据库

CREATE DATABASE gis_db;

2. 连接到数据库

psql -U postgres -d gis_db

3. 启用 PostGIS 扩展

CREATE EXTENSION postgis;
CREATE EXTENSION postgis_topology;

4. 创建空间表

以创建一个存储城市位置的表为例：

CREATE TABLE cities (
    id SERIAL PRIMARY KEY,
    name VARCHAR(100) NOT NULL,
    population INTEGER,
    geom GEOMETRY(Point, 4326)
);

解释：

  • id: 自增主键。
  • name: 城市名称，非空。
  • population: 城市人口。
  • geom: 存储城市的地理位置，类型为 Point，使用 EPSG:4326 坐标系统。

5. 创建空间索引

空间索引能显著提高空间查询的性能。

CREATE INDEX idx_cities_geom ON cities USING GIST (geom);

解释：

  • CREATE INDEX: 创建索引。
  • idx_cities_geom: 索引名称。
  • USING GIST: 使用 GiST 索引方法。
  • (geom): 针对 geom 字段创建空间索引。

6. 插入空间数据

INSERT INTO cities (name, population, geom)
VALUES (
    'Beijing',
    21540000,
    ST_GeomFromText('POINT(116.4074 39.9042)', 4326)
);

解释：

  • ST_GeomFromText('POINT(116.4074 39.9042)', 4326): 使用 Well-Known Text (WKT) 创建一个点对象，指定坐标系统为 EPSG:4326。

7. 查看空间表结构

\d cities

解释：

  • \d cities: 在 psql 中查看 cities 表的结构，包括字段类型和索引信息。

常用的空间函数和操作

PostGIS 提供了丰富的空间函数，以下是一些常用的函数及其详细解释和示例。

1. 创建几何对象

ST_GeomFromText

从 WKT（Well-Known Text）创建几何对象。

语法：

ST_GeomFromText(text, srid)

示例：

SELECT ST_GeomFromText('POINT(30 10)', 4326);

解释：

创建一个点对象，坐标为 (30, 10)，使用 EPSG:4326 坐标系统。

ST_Point

创建点对象。

语法：

ST_Point(x, y)

示例：

SELECT ST_Point(30, 10);

解释：

创建一个点对象，坐标为 (30, 10)。注意，此函数返回 GEOMETRY 类型，需要进一步指定 SRID。

2. 空间关系函数

ST_Contains

判断一个几何对象是否包含另一个几何对象。

语法：

ST_Contains(geometry A, geometry B)

示例：

SELECT ST_Contains(
    ST_GeomFromText('POLYGON((0 0, 0 10, 10 10, 10 0, 0 0))', 4326),
    ST_GeomFromText('POINT(5 5)', 4326)
);

解释：

判断点 (5,5) 是否在指定的多边形内，结果为 true。

ST_Within

判断一个几何对象是否在另一个几何对象内部。

语法：

ST_Within(geometry A, geometry B)

示例：

SELECT ST_Within(
    ST_GeomFromText('POINT(5 5)', 4326),
    ST_GeomFromText('POLYGON((0 0, 0 10, 10 10, 10 0, 0 0))', 4326)
);

解释：

判断点 (5,5) 是否在指定的多边形内，结果为 true。

ST_Intersects

判断两个几何对象是否相交。

语法：

ST_Intersects(geometry A, geometry B)

示例：

SELECT ST_Intersects(
    ST_GeomFromText('LINESTRING(0 0, 10 10)', 4326),
    ST_GeomFromText('LINESTRING(0 10, 10 0)', 4326)
);

解释：

判断两条线是否相交，结果为 true（交点为 (5,5)）。

3. 空间分析函数

ST_Distance

计算两个几何对象之间的距离。

语法：

ST_Distance(geometry A, geometry B)

示例：

SELECT ST_Distance(
    ST_GeomFromText('POINT(0 0)', 4326),
    ST_GeomFromText('POINT(3 4)', 4326)
);

解释：

计算点 (0,0) 和点 (3,4) 之间的距离，结果为 5。

注意：

  • 使用 GEOGRAPHY 类型可以计算地球表面的实际距离（米）。
  • 使用 GEOMETRY 类型计算的是平面上的距离。

示例（使用 GEOGRAPHY）：

SELECT ST_Distance(
    ST_GeogFromText('SRID=4326;POINT(116.4074 39.9042)'), -- 北京
    ST_GeogFromText('SRID=4326;POINT(-0.1278 51.5074)')  -- 伦敦
);

解释：

计算北京和伦敦之间的地理距离，结果为大约 8135 公里。

ST_Buffer

生成几何对象的缓冲区。

语法：

ST_Buffer(geometry, radius)

示例：

SELECT ST_Buffer(
    ST_GeomFromText('POINT(0 0)', 4326),
    1000
);

解释：

创建以点 (0,0) 为中心，半径为 1000 单位的缓冲区。注意，单位依赖于坐标系统。

注意：

  • 对于 GEOGRAPHY 类型，半径单位为米。
  • 对于 GEOMETRY 类型，半径单位为坐标系的单位（如度）。

示例（使用 GEOGRAPHY）：

SELECT ST_Buffer(
    ST_GeogFromText('SRID=4326;POINT(116.4074 39.9042)'), -- 北京
    10000  -- 10 公里
);
ST_Area

计算多边形的面积。

语法：

ST_Area(geometry)

示例：

SELECT ST_Area(
    ST_GeomFromText('POLYGON((0 0, 0 10, 10 10, 10 0, 0 0))', 4326)
);

解释：

计算指定多边形的面积，结果单位依赖于坐标系统。

注意：

  • 对于 GEOGRAPHY 类型，面积单位为平方米。
  • 对于 GEOMETRY 类型，面积单位为坐标系的单位的平方（如度的平方）。

4. 空间转换函数

ST_Transform

转换几何对象的空间参考系统。

语法：

ST_Transform(geometry, target_srid)

示例：

SELECT ST_Transform(
    ST_GeomFromText('POINT(116.4074 39.9042)', 4326),
    3857
);

解释：

将点 (116.4074, 39.9042) 从 EPSG:4326 转换到 EPSG:3857 坐标系统。

ST_Simplify

简化几何对象，减少其复杂度。

语法：

ST_Simplify(geometry, tolerance)

示例：

SELECT ST_Simplify(
    ST_GeomFromText('LINESTRING(0 0, 1 1, 2 2, 3 3)', 4326),
    1
);

解释：

简化折线，允许的最大偏差为 1 单位。

导入和导出空间数据

1. 使用 shp2pgsql 导入 Shapefile

PostGIS 提供了 shp2pgsql 工具，将 Shapefile 转换为 SQL 并导入数据库。

步骤：

  1. 准备 Shapefile 文件

    假设有一个名为 cities.shp 的 Shapefile 文件，存储城市位置数据。

  2. 转换并导入

    shp2pgsql -I -s 4326 cities.shp cities | psql -U postgres -d gis_db

    参数解释：

    • -I: 创建空间索引。
    • -s 4326: 指定空间参考系统（EPSG:4326）。
    • cities.shp: 要导入的 Shapefile 文件。
    • cities: 目标数据库表名。

  3. 验证导入

    登录 PostgreSQL，检查 cities 表是否存在并包含数据。

    \c gis_db
\d cities
SELECT COUNT(*) FROM cities;

2. 使用 ogr2ogr 导入多种格式

ogr2ogr 是 GDAL 工具的一部分，支持多种空间数据格式。

示例：

ogr2ogr -f "PostgreSQL" PG:"dbname=gis_db user=postgres" cities.geojson -nln cities -overwrite

参数解释：

  • -f "PostgreSQL": 输出格式为 PostgreSQL。
  • PG:"dbname=gis_db user=postgres": 目标数据库连接信息。
  • cities.geojson: 要导入的数据文件（GeoJSON 格式）。
  • -nln cities: 指定目标表名为 cities。
  • -overwrite: 覆盖已有表。

3. 导出空间数据

可以使用 pgsql2shp 或 ogr2ogr 导出数据。

使用 pgsql2shp:
pgsql2shp -f output.shp gis_db "SELECT * FROM cities"

解释：

  • -f output.shp: 指定输出文件名。
  • gis_db: 数据库名称。
  • "SELECT * FROM cities": 导出的数据查询语句。
使用 ogr2ogr:
ogr2ogr -f "GeoJSON" output.geojson PG:"dbname=gis_db user=postgres" cities

解释：

  • -f "GeoJSON": 输出格式为 GeoJSON。
  • output.geojson: 输出文件名。
  • PG:"dbname=gis_db user=postgres": 数据库连接信息。
  • cities: 要导出的表名。

空间查询示例与案例分析

通过具体案例，深入理解 PostGIS 的空间查询功能。

案例 1：查找指定范围内的城市

背景：

在 cities 表中存储了多个城市的位置。需要查找位于某一多边形区域内的城市。

步骤：

  1. 定义多边形区域

    以北京周边的一个简单多边形为例：

    SELECT ST_GeomFromText('POLYGON((116 39, 117 39, 117 40, 116 40, 116 39))', 4326);

  2. 执行查询

    SELECT name, population
FROM cities
WHERE ST_Within(
    geom,
    ST_GeomFromText('POLYGON((116 39, 117 39, 117 40, 116 40, 116 39))', 4326)
);

解释：

  • ST_Within(geom, polygon): 判断城市的 geom 是否在指定的多边形内。

示例结果：

namepopulation
Beijing21540000
Tianjin15620000

案例 2：计算两城市之间的距离

背景：

需要计算数据库中所有城市两两之间的距离。

步骤：

  1. 编写查询

    SELECT
    c1.name AS city1,
    c2.name AS city2,
    ST_Distance(c1.geom::geography, c2.geom::geography) AS distance_meters
FROM
    cities c1, cities c2
WHERE
    c1.id < c2.id;

解释：

  • c1.geom::geography: 将 geom 类型转换为 geography 类型，以计算地理距离（单位为米）。
  • c1.id < c2.id: 避免重复计算同一对城市之间的距离。

示例结果：

city1city2distance_meters
BeijingTianjin120000
BeijingShanghai1060000
TianjinShanghai1180000

案例 3：查找距离某点 10 公里内的城市

背景：

以北京天安门广场为中心，查找周边 10 公里内的城市。

步骤：

  1. 定义中心点

    天安门广场坐标：经度 116.4074, 纬度 39.9042。

  2. 执行查询

    SELECT name, population
FROM cities
WHERE ST_DWithin(
    geom::geography,
    ST_GeogFromText('SRID=4326;POINT(116.4074 39.9042)'),
    10000
);

解释：

  • ST_DWithin: 判断两个地理对象是否在指定距离内。
  • 10000: 距离为 10,000 米（10 公里）。

示例结果：

namepopulation
Beijing21540000
Haidian3000000
Chaoyang5000000

案例 4：查找城市的缓冲区内的道路

背景：

在 cities 表中存储城市位置，在 roads 表中存储道路信息。需要查找某城市 5 公里范围内的道路。

步骤：

  1. 创建 roads 表

    CREATE TABLE roads (
    id SERIAL PRIMARY KEY,
    name VARCHAR(100),
    geom GEOMETRY(LineString, 4326)
);

  2. 创建空间索引

    CREATE INDEX idx_roads_geom ON roads USING GIST (geom);

  3. 插入示例道路数据

    INSERT INTO roads (name, geom)
VALUES
('Road A', ST_GeomFromText('LINESTRING(116.4 39.9, 116.5 39.95)', 4326)),
('Road B', ST_GeomFromText('LINESTRING(116.3 39.8, 116.6 39.85)', 4326)),
('Road C', ST_GeomFromText('LINESTRING(116.45 39.9, 116.55 39.95)', 4326));

  4. 查找城市缓冲区内的道路

    假设查找北京周边 5 公里内的道路：

    SELECT r.name
FROM roads r, cities c
WHERE c.name = 'Beijing'
  AND ST_DWithin(
      r.geom::geography,
      c.geom::geography,
      5000
  );

解释：

  • ST_DWithin: 查找距离城市 geom 不超过 5000 米的道路。
  • r.geom::geography: 将道路几何转换为地理类型，以便计算实际距离。

示例结果：

name
Road A
Road C

空间索引与性能优化

空间索引是提高 PostGIS 空间查询性能的关键。以下是有关空间索引和性能优化的详细介绍。

1. 创建空间索引

空间索引通常使用 GiST（Generalized Search Tree）索引方法。

示例：

CREATE INDEX idx_cities_geom ON cities USING GIST (geom);

解释：

  • CREATE INDEX: 创建索引。
  • idx_cities_geom: 索引名称。
  • USING GIST: 使用 GiST 索引方法。
  • (geom): 针对 geom 字段创建空间索引。

2. 使用 EXPLAIN 分析查询计划

通过 EXPLAIN 可以查看查询是否使用了空间索引，从而判断查询性能。

示例：

EXPLAIN ANALYZE
SELECT name, population
FROM cities
WHERE ST_Within(
    geom,
    ST_GeomFromText('POLYGON((116 39, 117 39, 117 40, 116 40, 116 39))', 4326)
);

解释：

  • EXPLAIN ANALYZE: 显示查询执行计划和实际执行时间。
  • 检查输出中是否包含 Bitmap Heap Scan 或 Index Scan，表明空间索引被使用。

3. 优化查询性能的技巧

  • 使用合适的空间函数：例如，尽量使用 ST_DWithin 而不是 ST_Intersects，因为 ST_DWithin 可以更好地利用空间索引。

    示例：

    SELECT name, population
FROM cities
WHERE ST_DWithin(
    geom::geography,
    ST_GeogFromText('SRID=4326;POINT(116.4074 39.9042)'),
    10000
);

  • 减少转换开销：尽量避免在查询中频繁转换几何类型（如从 GEOMETRY 到 GEOGRAPHY），因为这会增加计算开销。可以考虑在表中直接存储 GEOGRAPHY 类型的字段。

  • 限制查询范围：在空间查询中，先使用边界框（Bounding Box）过滤，再进行精确的空间计算。例如，使用 && 操作符结合边界框。

    示例：

    SELECT name, population
FROM cities
WHERE geom && ST_MakeEnvelope(116, 39, 117, 40, 4326)
  AND ST_Within(
      geom,
      ST_GeomFromText('POLYGON((116 39, 117 39, 117 40, 116 40, 116 39))', 4326)
  );

    解释：

    • geom && ST_MakeEnvelope(...): 首先使用边界框进行粗筛，减少需要精确计算的几何对象数量。

4. 定期维护索引

空间索引可能会随着数据的插入、更新和删除而变得低效。定期进行索引维护，如重建索引，可以保持查询性能。

示例：

REINDEX INDEX idx_cities_geom;

解释：

  • REINDEX INDEX: 重建指定的索引。
  • idx_cities_geom: 要重建的空间索引名称。

空间数据可视化

虽然 PostGIS 本身不提供可视化功能，但可以结合其他工具进行空间数据的可视化。

1. 使用 QGIS

QGIS 是一个开源的 GIS 应用程序，支持与 PostGIS 数据库连接。

步骤：

  1. 安装 QGIS

    下载并安装 QGIS：QGIS 官方下载

  2. 连接到 PostGIS 数据库

    • 打开 QGIS。
    • 在左侧面板中，右键点击 “PostgreSQL” 并选择 “新建连接”。
    • 输入连接名称、主机、端口、数据库名称、用户名和密码。
    • 点击 “测试连接” 确认连接成功，然后点击 “确定”。

  3. 加载空间数据

    • 在 “浏览器” 面板中，展开新建的 PostgreSQL 连接。
    • 展开数据库，选择要加载的表（如 cities）。
    • 右键点击表名，选择 “添加选中的图层”。

  4. 可视化与样式

    • 选择图层，右键点击选择 “属性”。
    • 在 “样式” 选项卡中，可以设置点、线、多边形的样式，如颜色、大小、透明度等。

示例：

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

2. 使用 Web 地图框架

结合 PostGIS 和 Web 地图框架（如 Leaflet、OpenLayers），可以构建动态的在线地图应用。

步骤：

  1. 后端：使用 Node.js 和 Express 查询 PostGIS 数据

    • 安装依赖：

      npm install express pg

    • 编写服务器代码：

      const express = require('express');
const { Pool } = require('pg');
const app = express();
const port = 3000;

const pool = new Pool({
  user: 'postgres',
  host: 'localhost',
  database: 'gis_db',
  password: 'your_password',
  port: 5432,
});

app.get('/cities', async (req, res) => {
  try {
    const result = await pool.query('SELECT name, population, ST_AsGeoJSON(geom) AS geojson FROM cities');
    res.json(result.rows.map(row => ({
      name: row.name,
      population: row.population,
      geometry: JSON.parse(row.geojson),
    })));
  } catch (err) {
    console.error(err);
    res.status(500).send('Server Error');
  }
});

app.listen(port, () => {
  console.log(`Server running at http://localhost:${port}/`);
});

  2. 前端：使用 Leaflet 显示 GeoJSON 数据

    • 创建 index.html 文件：

      <!DOCTYPE html>
<html>
<head>
  <title>PostGIS + Leaflet</title>
  <link rel="stylesheet" href="https://unpkg.com/leaflet/dist/leaflet.css" />
  <style>
    #map { height: 600px; }
  </style>
</head>
<body>
  <h1>城市地图</h1>
  <div id="map"></div>

  <script src="https://unpkg.com/leaflet/dist/leaflet.js"></script>
  <script>
    const map = L.map('map').setView([39.9042, 116.4074], 10);

    L.tileLayer('https://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png', {
      attribution: '&copy; OpenStreetMap contributors'
    }).addTo(map);

    fetch('/cities')
      .then(response => response.json())
      .then(data => {
        L.geoJSON(data, {
          onEachFeature: function (feature, layer) {
            layer.bindPopup(`<b>${feature.name}</b><br>Population: ${feature.population}`);
          }
        }).addTo(map);
      });
  </script>
</body>
</html>

  3. 运行项目

    • 启动后端服务器：

      node server.js

    • 打开浏览器，访问 http://localhost:3000/，即可看到城市地图。

3. 使用 GIS 客户端工具

除 QGIS 外，其他 GIS 客户端如 GeoServer 也可以与 PostGIS 集成，提供 Web 服务接口。

步骤：

  1. 安装 GeoServer

    下载并安装 GeoServer：GeoServer 官方下载

  2. 配置 GeoServer 与 PostGIS 连接

    • 打开 GeoServer 管理界面。
    • 导航到 “数据” -> “工作区”，创建新的工作区。
    • 导航到 “数据” -> “存储”，添加新的 PostGIS 存储。
    • 输入数据库连接信息，测试连接，保存配置。

  3. 发布图层

    • 在新建的存储下，选择要发布的表（如 cities）。
    • 配置图层样式（使用 SLD 或预设样式）。
    • 发布图层后，可以通过 WMS、WFS 等协议进行访问。

示例：

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

项目实践：构建一个简单的地理信息系统

通过一个实际项目，综合应用所学的 PostGIS 知识。

项目背景

构建一个简单的城市信息管理系统，包括以下功能：

  • 存储和管理城市信息（名称、人口、位置）。
  • 查询指定区域内的城市。
  • 计算城市之间的距离。
  • 可视化城市位置和相关数据。

步骤详解

1. 设计数据库

创建数据库和启用 PostGIS：

CREATE DATABASE city_info;
\c city_info
CREATE EXTENSION postgis;
CREATE EXTENSION postgis_topology;

创建 cities 表：

CREATE TABLE cities (
    id SERIAL PRIMARY KEY,
    name VARCHAR(100) NOT NULL,
    population INTEGER,
    geom GEOGRAPHY(Point, 4326)
);

创建空间索引：

CREATE INDEX idx_cities_geom ON cities USING GIST (geom);
2. 导入示例数据

插入示例城市数据：

INSERT INTO cities (name, population, geom)
VALUES
('Beijing', 21540000, ST_GeogFromText('SRID=4326;POINT(116.4074 39.9042)')),
('Shanghai', 24240000, ST_GeogFromText('SRID=4326;POINT(121.4737 31.2304)')),
('Guangzhou', 13500000, ST_GeogFromText('SRID=4326;POINT(113.2644 23.1291)')),
('Shenzhen', 12500000, ST_GeogFromText('SRID=4326;POINT(114.0579 22.5431)')),
('Chengdu', 16330000, ST_GeogFromText('SRID=4326;POINT(104.0668 30.5728)'));
3. 实现功能
功能 1：查询指定区域内的城市

需求：

查找位于 POLYGON 区域内的城市。

实现：

SELECT name, population
FROM cities
WHERE ST_Within(
    geom,
    ST_GeogFromText('SRID=4326;POLYGON((115 30, 115 40, 125 40, 125 30, 115 30))')
);

解释：

查找位于经度 115°至125°，纬度 30°至40°之间的城市。

功能 2：计算城市之间的距离

需求：

计算所有城市两两之间的距离，并找出距离最近的两个城市。

实现：

SELECT
    c1.name AS city1,
    c2.name AS city2,
    ST_Distance(c1.geom, c2.geom) AS distance_meters
FROM
    cities c1, cities c2
WHERE
    c1.id < c2.id
ORDER BY
    distance_meters ASC
LIMIT 1;

解释：

  • 计算所有不同城市对之间的距离。
  • 按距离升序排列，取最短的一对。

示例结果：

city1city2distance_meters
ShenzhenGuangzhou138100
功能 3：可视化城市位置

步骤：

  1. 使用 QGIS 或 Web 地图框架（如 Leaflet）加载 cities 表数据。
  2. 在地图上显示城市位置，并标注名称和人口。

示例（使用 Leaflet）：

参考前述的 空间数据可视化 - 使用 Leaflet 部分。

功能 4：添加城市的缓冲区

需求：

为每个城市创建一个 50 公里的缓冲区，并查询缓冲区内的道路（假设存在 roads 表）。

实现：

  1. 创建缓冲区：

    ALTER TABLE cities ADD COLUMN buffer GEOGRAPHY(Geometry, 4326);

UPDATE cities
SET buffer = ST_Buffer(geom, 50000); -- 50 公里

  2. 查询缓冲区内的道路：

    SELECT r.name
FROM roads r, cities c
WHERE c.name = 'Beijing'
  AND ST_Intersects(r.geom, c.buffer);

解释：

  • ST_Buffer(geom, 50000): 为城市位置创建 50 公里的缓冲区。
  • ST_Intersects: 查找道路与缓冲区相交的记录。
4. 优化与扩展
  • 添加更多属性：如城市的行政区划、经济数据等。
  • 引入更多空间数据：如道路、河流、区域边界等，丰富地理信息系统的功能。
  • 实现高级空间分析：如空间聚类、热点分析、路径优化等。

常见问题与解决方法

1. PostGIS 扩展无法创建

问题原因： PostgreSQL 版本不兼容或 PostGIS 未正确安装。

解决方法:

  • 检查版本兼容性： 确保安装的 PostGIS 版本与 PostgreSQL 版本兼容。参考 PostGIS 支持的 PostgreSQL 版本.

  • 重新安装 PostGIS：

    # 对于 Debian/Ubuntu
sudo apt-get install --reinstall postgis postgresql-<version>-postgis-<version>

# 对于 CentOS/RHEL
sudo yum reinstall postgis postgis-utils

  • 检查库路径： 确保 PostgreSQL 能找到 PostGIS 的共享库文件。检查 shared_preload_libraries 配置项。

2. 空间查询性能低

问题原因： 缺少空间索引或索引未正确使用。

解决方法:

  • 创建空间索引： 确保对几何字段创建了 GiST 空间索引。

    CREATE INDEX idx_cities_geom ON cities USING GIST (geom);

  • 使用 EXPLAIN 分析查询计划： 确认查询是否使用了空间索引。

    EXPLAIN ANALYZE
SELECT name, population
FROM cities
WHERE ST_Within(
    geom,
    ST_GeomFromText('POLYGON((116 39, 117 39, 117 40, 116 40, 116 39))', 4326)
);

  • 优化查询语句： 使用更高效的空间函数和查询逻辑，如 ST_DWithin 替代 ST_Intersects，结合边界框过滤。

3. 坐标系统不一致导致查询失败

问题原因： 不同几何对象使用了不同的空间参考系统。

解决方法:

  • 使用 ST_Transform 函数转换几何对象到相同的空间参考系统。

    SELECT ST_Transform(geom, 4326)
FROM roads
WHERE id = 1;

  • 确保所有空间数据使用统一的 SRS。 在创建表时，明确指定 SRID。

4. 导入数据失败

问题原因： 数据格式不正确或缺少必要的字段。

解决方法:

  • 检查数据文件格式： 确保数据文件（如 Shapefile、GeoJSON）的格式正确，符合 PostGIS 导入工具的要求。

  • 确认目标表的结构与数据文件匹配： 字段名称、数据类型、空间参考系统等应一致。

  • 使用日志信息查找具体错误原因： 检查 PostgreSQL 和导入工具的错误日志，定位问题。

5. 插入空间数据时报错

问题原因： 几何对象不符合指定的 SRID 或格式错误。

解决方法:

  • 检查 SRID 是否正确： 使用 ST_SetSRID 或 ST_GeogFromText 指定正确的 SRID。

    INSERT INTO cities (name, population, geom)
VALUES (
    'Test City',
    1000000,
    ST_SetSRID(ST_Point(120, 30), 4326)
);

  • 验证几何对象的正确性： 使用 ST_IsValid 检查几何对象是否有效。

    SELECT ST_IsValid(geom)
FROM cities
WHERE name = 'Test City';
比松
关注
PostGIS 小白入门指南：行政区划边界 PostgreSQL
XwlyVue的博客
09-17 301
在本文中，我们介绍了如何使用 PostGIS 扩展来获取行政区划边界数据，并进行一些基本的空间分析操作。PostGIS 是一个用于 PostgreSQL 数据库的地理空间扩展，它提供了一套功能强大的工具和函数，用于处理地理空间数据。在本文中，我们将介绍如何使用 PostGIS 扩展来获取行政区划边界数据，并进行一些基本的空间分析操作。在实际使用中，请确保数据的准确性和合法性，并根据需要进行适当的数据清洗和处理。一旦数据导入成功，就可以使用 PostGIS 提供的函数和操作符来查询和分析行政区划边界数据。
PostGis -基础、Springboot 整合、电子围栏处理
何xiao树
07-30 1108
说明：PostGIS是PostgreSQL对象关系数据库的空间数据库扩展。它增加了对地理对象的支持，允许在SQL中运行位置查询。扩展插件，空间数据处理 postgis安装。mybatisplus 自定义类型处理。
Postgis-学习 入门 教程 教学 手册-Postgresql
08-09
PostGIS简介、安装、新建空间数据库、简单SQL语句查询练习、几何数据操作和练习、空间关系操作与练习、空间链接与练习、投影数据与练习、空间索引与练习…… 还包括PostgreSQL for spatial的配置，验证，postgis函数简介等信息。 非常详细。
【职场技术栈学习实录】PostgreSQL 与 PostGIS详解
最新发布
driftingalone98的博客
03-15 1101
本系列是我作为一名后端开发，在日常工作中积累的技术学习与实践笔记。在这里，我将分享工作中后端开发领域的技术栈、架构设计、代码实现以及问题解决的经验与心得。出于保护公司隐私的原因，实际业务场景无法详细展开，但我会通过抽象化的需求和通用的技术问题，来介绍与探讨如何在实际工作中应用这些技术。
PostgreSQL + PostGIS：空间数据存储及管理解决方案
GIS工具
09-02 2011
PostgreSQL 是一个开源关系数据库管理系统(RDBMS)，以其稳健性、灵活性和遵守 SQL 标准而闻名。它于 20 世纪 80 年代诞生于加州大学伯克利分校，现已发展成为世界各地开发人员和数据库管理员的必备工具。巩固 PostgreSQL 受欢迎程度的一些功能包括：可扩展性：用户可以定义自己的数据类型、运算符、聚合函数、索引方法等。并发性和性能：它能够在不影响数据完整性的情况下处理多个并发事务。
PostGIS 快速入门
饮一瓢
09-13 4万+
转自:http://live.osgeo.org/zh/quickstart/postgis_quickstart.html PostGIS 快速入门 　　PostGIS 是 PostgreSQL 关系数据库的空间操作扩展。它为 PostgreSQL 提供了存储、查询和修改空间关系的能力。本文中 ‘PostgreSQL’ 指代基本的关系数据库功能，而 ‘PostGIS’ 指代扩展的空间
POSTGIS教程
10-03 1万+
一 什么是PostgreSQL和PostGIS1.1 什么是PostgreSQL说起数据库，大家耳熟能详的商业数据库产品当推Oracle、微软的SqlServer和IBM的DB2等，而开源数据库中则有两大产品MySQL和PostgreSQL。PostgreSQL 是一种对象-关系型数据库管理系统（ORDBMS），也是目前功能最强大、特性最丰富和最复杂的自由软件数据库系统。它起源于伯克利（BSD）的
PostGIS五分钟入门【空间数据库】
新缸中之脑
05-09 3020
在本文中，我们将介绍 PostGIS 的一些基础知识及其功能，以及一些可用于简化解决方案或提高性能的提示和技巧。推荐：用快速搭建可编程3D场景。简而言之 - PostGIS 是一个 Postgres 扩展，增加了对存储和操作空间数据类型的支持。当我们构建在地图上存储、操作和可视化数据的软件应用程序时，我们通常需要使用空间数据存储。我们可以使用谷歌地图或类似应用程序作为典型地理空间可视化软件功能的一个很好的例子。
PostGIS教程一：PostGIS介绍
热门推荐
不睡觉的怪叔叔的博客
12-22 9万+
一、什么是空间数据库 PostGIS是一个空间数据库，Oracle Spatial和SQL Server(2008和之后版本）也是空间数据库。 但是这意味着什么？是什么使普通数据库变成空间数据库？ 简短的答案是... 空间数据库像存储和操作数据库中其他任何对象一样去存储和操作空间对象。 下面简短介绍了空间数据库的发展，然后回顾了将空间数据与数...
PostGIS学习教程六：几何图形（geometry）
数字人生
11-15 3130
在前面的章节中，我们已经往数据库中加载了数据，现在让我们来先看一些简单的例子。在pgAdmin中，再次选择nyc数据库并打开SQL查询工具。将下面的SQL代码粘贴到pgAdmin SQL Editor窗口中（删除默认情况下可能存在的任何文本），然后执行。上面的示例创建了一个表（geometries），然后向该表中插入5个几何图形数据（geometry）：一个点（POINT）一条线（LINESTRING）一个多边形（POLYGON）
CentOS7下PostgreSQL与PostGIS插件离线部署指南 用于centos7的x86-64 包括postgis3.3
09-23
本文详细介绍了如何在CentOS7环境中离线安装PostgreSQL及其配套PostGIS插件的方法。提供了完整的离线仓库配置步骤，并演示不同版本之间的兼容性和具体的安装命令流程。 适合人群：从事地理信息系统（GIS）工作的...
postgis 教程
06-25
postgis空间数据库权威教程 最完整的空间数据库手册
Postgis教程
12-18
PostgreSQL 是世界上技术最先进的开源数据库，其前身是1977年一个源于Berkeley名为Ingres的非关系型数据库，其项目领导人为Michael Stonebraker教授。1982年该教授商业化了Ingres；1985年，Michael Stonebraker教授回到Berkeley，开始对新的数据库设计进行研究，并于次年在美国防务高级研究项目局（DARPA）、陆军研究办公室 （ARO）、国家科学基金（NSF）以及ESL, Inc等机构的赞助下启动了Postgres（Post-Ingres）项目。
PostGIS教程（多套）
07-19
包含了postgis中的五套教程以及postgis的函数库
PostGIS源代码安装指南
04-15
### PostGIS源代码安装指南详解 #### 前言 本文档主要介绍如何从源代码安装PostGIS 2.5，并解决在使用`st_asmvt`函数返回MapBox格式矢量切片时遇到的`missinglibprotobuf-c`错误。通过详细步骤确保PostGIS能顺利...
PostGIS-入门实例
孙霸天的专栏
03-03 2474
PostgreSQL已经内置了很多空间特性：几何数据类型、几何类型函数与运算符、空间数据索引但对于现实世界的复杂需求仍力有不逮：有很多辅路支路的道路，需要用多条折线来表示；行政区域的飞地，需要用多个多边形的集合来表示;高效判断点是否在多边形内、两条折线是否有交点，两个区域相离、相邻、重叠还是包含。好在PostgreSQL强大 的扩展机制为解决问题留下了一个窗口，PostGIS 为此而生，它已经成为GIS行业的事实标准。
PostGIS入门与命令实战指南
Postgis学习笔记是一份详细的文档，旨在引导读者理解和掌握PostGIS，一个开源的地理信息系统扩展，用于在PostgreSQL数据库中存储和处理地理空间数据。以下是该文档中的关键知识点： 1. **PostGIS基础**: - 初步...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值