能否参照示例缩放SQL geography多边形？需实现内外x米缓冲，可借助DotSpatial/R

当然可以！下面我会分别介绍几种满足你需求的实现方式，从优先推荐的SQL原生方案，到DotSpatial和R工具的替代方案：

SQL Server的geography类型自带STBuffer()方法，直接就能实现向内/向外的缓冲操作，而且是基于椭圆体计算的，精度有保障。

基础用法 #

-- 向外缓冲100米（正数代表向外扩展）
SELECT original_geo.STBuffer(100) AS buffered_geo_out
FROM your_table;

-- 向内缓冲50米（负数代表向内收缩）
SELECT original_geo.STBuffer(-50) AS buffered_geo_in
FROM your_table;

关键注意点 #

• 确保你的geography对象的空间参考标识符（SRID）是正确的（比如常用的WGS84，SRID=4326），这样STBuffer()的距离参数才会以米为单位计算。
• 如果向内缓冲的距离过大，可能会导致多边形变成无效几何甚至消失，建议先通过STIsValid()检查原始几何的有效性，必要时用MakeValid()修复：

  SELECT original_geo.MakeValid().STBuffer(-50) AS buffered_geo_in
  FROM your_table
  WHERE original_geo.STIsValid() = 0;
  

如果SQL方案受限于权限或版本，DotSpatial作为.NET生态的空间处理库，也能轻松实现缓冲操作。

核心步骤与示例代码 #

using DotSpatial.Topology;
using DotSpatial.Projections;

// 1. 从SQL读取地理数据（这里假设是WKT格式）
string originalWkt = "POLYGON((-122.34900 47.65100, -122.34900 47.65200, ...))";
IGeometry originalGeom = Geometry.FromWkt(originalWkt);

// 2. 转换到投影坐标系（确保米单位的准确性，比如UTM带区）
ProjectionInfo sourceCrs = ProjectionInfo.FromEpsgCode(4326); // WGS84
ProjectionInfo targetCrs = ProjectionInfo.FromEpsgCode(32610); // 示例UTM 10N
originalGeom = originalGeom.Project(sourceCrs, targetCrs);

// 3. 执行缓冲操作
IGeometry bufferedOut = originalGeom.Buffer(100); // 向外100米
IGeometry bufferedIn = originalGeom.Buffer(-50);  // 向内50米

// 4. 转换回原始地理坐标系
bufferedOut = bufferedOut.Project(targetCrs, sourceCrs);
bufferedIn = bufferedIn.Project(targetCrs, sourceCrs);

小提示 #

如果对精度要求不高，也可以直接在地理坐标系下执行缓冲，但投影到平面坐标系能获得更准确的米级计算结果。

R的sf包是处理空间数据的神器，语法简洁且支持地理坐标系的椭圆体缓冲计算。

核心步骤与示例代码 #

library(sf)
library(odbc)

# 1. 连接SQL Server并读取空间数据
conn <- dbConnect(
  odbc(),
  Driver = "SQL Server",
  Server = "你的服务器地址",
  Database = "你的数据库名"
)

spatial_df <- dbGetQuery(conn, "SELECT original_geo FROM your_table") %>%
  st_as_sf(wkt = "original_geo", crs = 4326) # 转换为sf对象，指定坐标系

# 2. 执行缓冲操作
buffered_out <- st_buffer(spatial_df, dist = 100) # 向外100米
buffered_in <- st_buffer(spatial_df, dist = -50)  # 向内50米

# 3. 关闭数据库连接
dbDisconnect(conn)

优势说明 #

sf包的st_buffer()在地理坐标系（如4326）下会自动使用椭圆体计算距离，不需要手动投影，非常省心。

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