如何从灰度卫星影像中提取矢量线段？技术问询

一、直接从ridge-detection的检测结果提取坐标 #

你当前使用的ridge-detection工具中，detectLines()返回的result对象已经存储了每条线段的原始坐标信息（而非仅生成可视化图像）。先打印result查看结构，一般包含线段的起点、终点或完整点序列：

# 先查看result的结构，确认坐标字段
print(type(result))
print(result[0])  # 打印第一条线段的内容

假设result是线段点列表的集合，可直接提取并转换为GeoJSON格式：

import json

vector_segments = []
for line_points in result:
    # 转换为GeoJSON要求的坐标格式
    coords = [[point[0], point[1]] for point in line_points]
    if len(coords) >= 2:
        vector_segments.append({
            "type": "Feature",
            "geometry": {"type": "LineString", "coordinates": coords},
            "properties": {}
        })

# 保存为GeoJSON文件
with open("satellite_edges.geojson", "w") as f:
    json.dump({
        "type": "FeatureCollection",
        "features": vector_segments
    }, f)

二、对输出灰度图做矢量化后处理 #

如果无法直接获取result的坐标，可对生成的img_only_lines二值图进行轮廓提取：

步骤1：用OpenCV提取线条轮廓 #

import cv2
import numpy as np
import json

# 将PIL图像转为OpenCV格式的二值图
img_np = np.array(img_only_lines)
if img_np.ndim == 3:
    img_np = cv2.cvtColor(img_np, cv2.COLOR_RGB2GRAY)
# 反转黑白，确保线条为前景区域
_, binary_img = cv2.threshold(img_np, 127, 255, cv2.THRESH_BINARY_INV)

# 提取轮廓
contours, _ = cv2.findContours(binary_img, cv2.RETR_EXTERNAL, cv2.CHAIN_APPROX_SIMPLE)

步骤2：转换为矢量格式并保存 #

vector_features = []
for idx, cnt in enumerate(contours):
    # 简化轮廓点，减少冗余坐标
    epsilon = 0.001 * cv2.arcLength(cnt, True)
    approx_points = cv2.approxPolyDP(cnt, epsilon, True)
    # 转换为坐标列表
    coords = [point[0].tolist() for point in approx_points]
    if len(coords) >= 2:
        vector_features.append({
            "type": "Feature",
            "geometry": {"type": "LineString", "coordinates": coords},
            "properties": {"segment_id": idx}
        })

# 保存为GeoJSON
with open("edges_from_image.geojson", "w") as f:
    json.dump({"type": "FeatureCollection", "features": vector_features}, f)

# 转换为SHP格式（需提前安装geopandas）
import geopandas as gpd
gdf = gpd.GeoDataFrame.from_features(vector_features)
gdf.to_file("satellite_edges.shp")

三、专业遥感工具辅助矢量化 #

如果需要结合卫星影像的地理投影信息，可采用以下方式：

• 将处理后的灰度图导入QGIS，使用提取线工具（Processing Toolbox → 向量提取 → 提取线）完成矢量化，同时可直接赋予图像的地理投影参数
• 使用GDAL命令行工具gdal_polygonize.py将线条区域转为面要素，再提取边界得到线矢量

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