如何在Geopandas中实现类似Scikit-image的Swirl变换及相关问题

1. 是否可以对GeoPandas DataFrame执行类似Scikit-image的Swirl变换？ #

可以实现。核心逻辑是用Scikit-image的swirl函数对几何图形的坐标点做变换，再把变换后的坐标重新构建成GeoPandas支持的几何对象（如Polygon、Point等）。

2. 自定义GeoDataFrame变换后出现GEOMETRYCOLLECTION EMPTY的原因 #

常见诱因有这几点：

• 坐标范围与变换中心不匹配：Swirl默认以类似图像中心（常为[0,0]或图像尺寸中点）为原点，若你的GeoDataFrame坐标远离这个中心，变换后坐标可能超出合理范围，导致几何失效。
• 几何类型不兼容：如果你的swirl_poly函数只支持特定几何类型（比如仅处理Polygon），但自定义数据包含MultiPolygon、Point等其他类型，处理时会返回空几何。
• 变换函数逻辑错误：提取坐标、重构几何时出错——比如没把变换后的坐标转为可识别的几何格式，或者重构时x/y坐标顺序搞反。
• 坐标参考系（CRS）问题：若GeoDataFrame用的是地理坐标系（如WGS84经纬度），Swirl基于笛卡尔坐标系的变换会导致坐标过度扭曲，生成无效几何。

3. 如何对整个GeoDataFrame执行Swirl变换 #

无需逐个处理单个几何，直接对整个geometry列应用变换函数即可，示例代码如下：

import geopandas as gpd
from skimage.transform import swirl
from shapely.geometry import Polygon, Point, MultiPolygon

def swirl_geometry(geom, center=(0,0), strength=1, radius=100):
    # 按几何类型分别处理
    if isinstance(geom, Polygon):
        coords = list(geom.exterior.coords)
        xs, ys = zip(*coords)
        # skimage的swirl输入是(y, x)格式数组，注意顺序
        transformed = swirl([ys, xs], center=center, strength=strength, radius=radius)
        new_coords = list(zip(transformed[1], transformed[0]))
        return Polygon(new_coords)
    elif isinstance(geom, Point):
        x, y = geom.x, geom.y
        transformed = swirl([[y], [x]], center=center, strength=strength, radius=radius)
        return Point(transformed[1][0], transformed[0][0])
    elif isinstance(geom, MultiPolygon):
        transformed_polys = [swirl_geometry(poly, center, strength, radius) for poly in geom.geoms]
        return MultiPolygon(transformed_polys)
    else:
        return geom

# 加载自定义GeoDataFrame
gdf = gpd.read_file("your_data.geojson")

# 地理坐标系转笛卡尔投影坐标系（比如UTM，根据数据区域选对应EPSG）
gdf = gdf.to_crs(epsg=32633)

# 用数据中心作为变换原点，避免默认[0,0]导致的异常
center_x, center_y = gdf.geometry.centroid.x.mean(), gdf.geometry.centroid.y.mean()

# 对整个geometry列应用变换，生成新列
gdf['swirled_geometry'] = gdf.geometry.map(lambda g: swirl_geometry(g, center=(center_y, center_x), strength=2, radius=5000))

# 可视化结果
gdf.plot(column='swirled_geometry')

关键注意事项 #

• 坐标转换：必须先把经纬度的地理坐标系转为投影坐标系（如UTM），Swirl是平面变换，直接处理球面坐标会出问题。
• 变换中心：用数据自身的中心坐标（比如所有几何质心的平均值）替代默认原点，变换效果更合理。
• 几何类型全覆盖：函数要兼容数据中所有存在的几何类型，避免遗漏导致空几何。

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