为何GeoDataFrame几何列是POINT Z(long, lat)而非POINT(long, lat)？

两者的核心区别 #

• POINT：二维平面点，仅包含水平坐标（如经度、纬度），无垂直维度信息。
• POINT Z：三维空间点，在水平坐标基础上额外存储了Z值（通常代表高程、海拔等垂直数据）。你看到的POINT Z (long, lat)只是显示简化，实际该点的Z值是存在的（可能为默认0或未在输出中展示）。

转换可行性 #

二者完全可以互相转换：

• POINT Z → POINT：丢弃Z值，保留二维坐标即可。
• POINT → POINT Z：为每个点添加指定的Z值（如固定0或自定义高程数据）。

Python转换解决方案 #

你修改EPSG无效是因为坐标系（EPSG）和几何维度是独立属性，修改坐标系不会移除Z值。以下是针对geopandas的处理方案：

方法1：自定义函数移除Z值 #

import geopandas as gpd
from shapely.geometry import Point

# 读取目标shapefile
geo_points = gpd.read_file('../point_data/groundtruth.shp')

# 定义移除Z值的处理函数
def strip_z(geom):
    if geom.has_z:
        return Point(geom.x, geom.y)
    return geom

# 应用函数转换几何列
geo_points['geometry'] = geo_points['geometry'].apply(strip_z)

# 验证转换结果
print(geo_points['geometry'].geom_type.unique())  # 输出应为['Point']

方法2：使用shapely内置工具（shapely 1.8+版本支持） #

import geopandas as gpd
from shapely.force_2d import force_2d

geo_points = gpd.read_file('../point_data/groundtruth.shp')
# 直接调用force_2d移除Z/M值
geo_points['geometry'] = geo_points['geometry'].apply(force_2d)

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