GeoDataFrame提取Geometry的XY值变为整数的异常问题求助

可能原因 #

• 几何数据导入时被隐式截断：从SQL读取过程中，数据库驱动或读取工具（如read_postgis）可能错误地将浮点坐标转换为整数类型，或截断了小数部分
• GeoPandas x/y 属性的类型推断问题：当几何对象的坐标存储为整数类型时，直接提取x/y会返回整数；即使WKT显示有小数，底层Shapely对象可能已丢失精度
• Pandas自动类型转换：提取后的数值被Pandas自动推断为整数类型，覆盖了原本的浮点精度

排查步骤 #

1. 验证几何对象的实际精度：

# 查看第一条记录的原始坐标值（非显示值）
print(df_kunde.geometry.iloc[0].x)
# 查看GeoSeries的底层数据类型
print(df_kunde.geometry.dtype)

2. 检查SQL导入环节：确认数据库中存储的坐标是浮点类型，且读取时未设置截断精度的参数（如某些驱动的integer_only类选项）

解决方案 #

1. 强制浮点类型提取坐标：
   避免直接使用geometry.x/y属性，改用apply显式转换为浮点类型，绕过自动类型推断：

df_kunde['utmx'] = df_kunde.geometry.apply(lambda geom: float(geom.x))
df_kunde['utmy'] = df_kunde.geometry.apply(lambda geom: float(geom.y))

2. 修正WKT转换的精度控制：
   to_wkt()默认trim=True会修剪小数末尾的零和小数点，导致浮点值看起来像整数。如果需要保留完整精度，使用trim=False；若需强制转为整数，可在提取后显式转换：

# 保留完整小数精度的WKT
df_kunde['wkt_full'] = df_kunde.geometry.to_wkt(trim=False)
# 强制转为整数的坐标列
df_kunde['utmx_int'] = df_kunde['utmx'].astype(int)

3. 重新导入数据确保精度：
   使用geopandas.read_postgis时，指定geom_col参数并确保数据库连接驱动正确解析浮点类型，例如：

import geopandas as gpd
from sqlalchemy import create_engine

engine = create_engine("postgresql://user:pass@host/db")
df_kunde = gpd.read_postgis("SELECT * FROM kunde_table", engine, geom_col="geometry")

相关反馈 #

该问题在GeoPandas 0.12.x和Pandas 1.4.x版本中已有用户反馈，核心是底层Shapely与Pandas类型交互的细节问题。即使升级新版本，若数据导入环节存在精度丢失，问题仍会复现。

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