【GIS】算法原理：几何作图与计算

最新推荐文章于 2025-08-25 16:51:38 发布

原创最新推荐文章于 2025-08-25 16:51:38 发布 · 412 阅读

6 ·

CC 4.0 BY-SA版权

文章标签：

#算法

GIS 专栏收录该内容

17 篇文章

订阅专栏

在GIS（地理信息系统）中，几何作图与计算是空间分析的基础。本文将从线段延长、线段交点计算、线段与圆的交点计算、中心点计算、过点作垂线、过点作平行线、三点共圆、线段打断等方面，详细讲解GIS中的几何作图与计算算法，并提供C++代码实现。

几何作图与计算算法

线段延长

给定线段 AB 和点 P，输入延长距离 d，求延长线。
步骤：
1. 计算线段 AB 的长度 L。
2. 判断 P 点距离 A 或 B 更近。
3. 根据最近点计算延长点 D 的坐标。
4. 连接 D 点与最近点，得到延长线。

计算线段或直线与线段的交点

给定两条线段 L0 = P1P2 和 L1 = Q1Q2，计算它们的交点。
步骤：
1. 判断两条线段是否相交。
2. 将线段视为直线，联立直线方程求解交点。
3. 判断交点是否在线段范围内。

计算线段或直线与圆的交点

给定圆心 O、半径 r 和线段 L = P1P2，计算交点。
步骤：
1. 判断线段是否完全在圆内。
2. 计算直线与圆的交点。
3. 判断交点是否在线段范围内。

中心点的计算

多边形的中心点（质心）通过分割多边形为三角形，计算各三角形中心点的加权平均。
公式：
Cx = (Σ(Cx_i * Area_i)) / ΣArea_i
Cy = (Σ(Cy_i * Area_i)) / ΣArea_i

过点作垂线

给定点 C 和线段 AB，求过 C 点垂直于 AB 的垂线。
步骤：
1. 计算点 C 到直线 AB 的垂足 P。
2. 连接 C 和 P，得到垂线。

过点作平行线

给定点 C 和线段 AB，求过 C 点平行于 AB 的平行线。
步骤：
1. 计算 AB 的方向向量。
2. 根据方向向量和平移距离计算新线段的端点。

三点共圆

给定三点 A、B、C，求通过三点的圆。
步骤：
1. 计算两条边的垂直平分线，求交点作为圆心。
2. 计算圆心到任意一点的距离作为半径。

线段打断

给定线段 AB 和距离 d，在线段上插入点 C，将线段分为两部分。
步骤：
1. 计算线段 AB 的长度 L。
2. 根据比例计算点 C 的坐标。

C++代码实现

#include <iostream>
#include <cmath>
using namespace std;

// 二维点结构体
struct Point {
    double x, y;
    Point(double x = 0, double y = 0) : x(x), y(y) {}
};

// 计算两点距离
double distance(const Point& a, const Point& b) {
    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}

// 线段延长
Point extendSegment(const Point& A, const Point& B, const Point& P, double d) {
    double L = distance(A, B);
    double dx = B.x - A.x;
    double dy = B.y - A.y;
    if (distance(P, A) < distance(P, B)) {
        return {A.x - dx * d / L, A.y - dy * d / L};
    } else {
        return {B.x + dx * d / L, B.y + dy * d / L};
    }
}

// 计算线段交点
Point computeIntersection(const Point& P1, const Point& P2, const Point& Q1, const Point& Q2) {
    double A1 = P2.y - P1.y;
    double B1 = P1.x - P2.x;
    double C1 = A1 * P1.x + B1 * P1.y;

    double A2 = Q2.y - Q1.y;
    double B2 = Q1.x - Q2.x;
    double C2 = A2 * Q1.x + B2 * Q1.y;

    double det = A1 * B2 - A2 * B1;
    if (det == 0) {
        return {INFINITY, INFINITY}; // 平行或共线
    } else {
        double x = (B2 * C1 - B1 * C2) / det;
        double y = (A1 * C2 - A2 * C1) / det;
        return {x, y};
    }
}

// 计算线段与圆的交点
void computeCircleIntersection(const Point& O, double r, const Point& P1, const Point& P2) {
    // 判断线段是否在圆内（省略）
    // 计算直线与圆的交点（省略）
}

// 计算多边形中心点
Point computeCentroid(const vector<Point>& polygon) {
    double area = 0, Cx = 0, Cy = 0;
    int n = polygon.size();
    for (int i = 0; i < n; i++) {
        int j = (i + 1) % n;
        double temp = (polygon[i].x * polygon[j].y - polygon[j].x * polygon[i].y);
        area += temp;
        Cx += (polygon[i].x + polygon[j].x) * temp;
        Cy += (polygon[i].y + polygon[j].y) * temp;
    }
    area /= 2;
    Cx /= (6 * area);
    Cy /= (6 * area);
    return {Cx, Cy};
}

int main() {
    // 测试线段延长
    Point A(0, 0), B(4, 4), P(2, 2);
    Point D = extendSegment(A, B, P, 2);
    cout << "延长点坐标: (" << D.x << ", " << D.y << ")" << endl;

    // 测试线段交点
    Point Q1(0, 4), Q2(4, 0);
    Point intersection = computeIntersection(A, B, Q1, Q2);
    cout << "交点坐标: (" << intersection.x << ", " << intersection.y << ")" << endl;

    // 测试多边形中心点
    vector<Point> polygon = {{0, 0}, {4, 0}, {4, 4}, {0, 4}};
    Point centroid = computeCentroid(polygon);
    cout << "多边形中心点: (" << centroid.x << ", " << centroid.y << ")" << endl;

    return 0;
}