针对你遇到的OSM水体多边形质心落在外部的问题，这里有两种实用的解决方案，都可以基于GEOS库实现：

方法1：从边界最近点向内偏移（精确可控） #

这种方法先定位边界上离质心最近的点，再将其向内偏移得到有效内部点，适合需要精确控制位置的场景：

步骤代码示例： #

#include <geos/geom/Geometry.h>
#include <geos/geom/Point.h>
#include <geos/geom/GeometryFactory.h>
#include <cmath>
#include <limits>

using namespace geos::geom;

Coordinate getClosestInteriorPoint(Geometry* geometry, const Coordinate& centroidCoord, GeometryFactory* factory) {
    // 1. 创建质心点并检查是否在内部
    Point* centroidPoint = factory->createPoint(centroidCoord);
    if (geometry->contains(centroidPoint)) {
        Coordinate result = centroidCoord;
        delete centroidPoint;
        return result;
    }

    // 2. 获取几何体边界，找到离质心最近的边界点
    Geometry* boundary = geometry->getBoundary();
    Point* nearestBoundaryPoint = boundary->nearestPoint(centroidPoint);
    const Coordinate* nearestBoundaryCoord = nearestBoundaryPoint->getCoordinate();

    // 3. 计算向内偏移的向量（从质心指向边界点的方向，深入内部）
    double dx = nearestBoundaryCoord->x - centroidCoord.x;
    double dy = nearestBoundaryCoord->y - centroidCoord.y;
    double length = sqrt(dx*dx + dy*dy);
    if (length < 1e-9) { // 避免除以0
        delete centroidPoint;
        delete boundary;
        delete nearestBoundaryPoint;
        return *nearestBoundaryCoord;
    }
    dx /= length;
    dy /= length;

    // 4. 偏移微小距离（根据坐标单位调整，这里用0.00001度适配WGS84）
    const double offset = 0.00001;
    Coordinate interiorCoord(
        nearestBoundaryCoord->x + dx * offset,
        nearestBoundaryCoord->y + dy * offset
    );

    // 5. 验证偏移点是否在内部， fallback 到边界质心（如果偏移失败）
    Point* interiorPoint = factory->createPoint(interiorCoord);
    if (!geometry->contains(interiorPoint)) {
        Geometry* buffer = interiorPoint->buffer(offset);
        Geometry* intersection = geometry->intersection(buffer);
        if (!intersection->isEmpty()) {
            interiorCoord = intersection->getCentroid()->getCoordinate();
        } else {
            interiorCoord = *nearestBoundaryCoord;
        }
        delete buffer;
        delete intersection;
    }

    // 清理资源
    delete centroidPoint;
    delete boundary;
    delete nearestBoundaryPoint;
    delete interiorPoint;

    return interiorCoord;
}

方法2：多次随机内部点选最优（简单易实现） #

如果你不需要极致精确，只是需要一个接近质心的内部点，这种方法更简单：利用getInteriorPoint多次生成随机内部点，选择离质心最近的那个：

步骤代码示例： #

#include <geos/geom/Geometry.h>
#include <geos/geom/Point.h>
#include <cmath>
#include <limits>

using namespace geos::geom;

Coordinate getClosestRandomInteriorPoint(Geometry* geometry, const Coordinate& centroidCoord) {
    Coordinate bestCoord;
    double minDistance = std::numeric_limits<double>::max();
    const int maxAttempts = 10; // 可根据需求调整次数

    for (int i = 0; i < maxAttempts; ++i) {
        Point* interiorPoint = geometry->getInteriorPoint();
        if (!interiorPoint) continue;

        const Coordinate* coord = interiorPoint->getCoordinate();
        double distance = centroidCoord.distance(*coord);
        if (distance < minDistance) {
            minDistance = distance;
            bestCoord = *coord;
        }
        delete interiorPoint;
    }

    // 如果所有尝试都失败（理论上不会发生在有效多边形上）， fallback 到质心
    if (minDistance == std::numeric_limits<double>::max()) {
        bestCoord = centroidCoord;
    }

    return bestCoord;
}

注意事项 #

• 坐标单位：偏移量offset需要根据你的几何体坐标系统调整（比如WGS84用度，UTM用米），避免偏移过大或过小。
• 资源清理：GEOS的几何对象需要手动释放内存，避免内存泄漏。
• 复杂几何体：对于MultiPolygon等复杂类型，两种方法都能兼容，因为GEOS的getBoundary、contains等方法支持复杂几何体。

内容的提问来源于stack exchange，提问作者Walter