如何解决缓冲后质心点坐标精度不够的问题？

Question

buffer生成多边形（此处为正方形）时，生成的几何点坐标与.centroid方法在多边形after他们这一代。

这里举个例子，只有一点。

from shapely.ops import transform
import geopandas as gpd
import shapely.wkt
import pyproj
from math import sqrt

def edge_size(area): return sqrt(area)*1e3

point = "POINT (4379065.583907348 2872272.254645019)"
point = shapely.wkt.loads(point)

center = gpd.GeoSeries(point)

project = pyproj.Transformer.from_proj(
    pyproj.Proj('epsg:3395'), 
    pyproj.Proj('epsg:4326'),
    always_xy=True) 

center = center.apply(lambda p: transform(project.transform, p))

print(center.iloc[0])

square = point.buffer(
    edge_size(3), cap_style=3) #distance of 3km2

square = gpd.GeoSeries(square)
square = square.apply(lambda p: transform(project.transform, p))
square = square.apply(lambda p: p.centroid)

print(square.iloc[0])

#POINT (39.33781544185747 25.11929860805248)
#POINT (39.33781544185747 25.11929777802279)

这会导致之后的处理错误。
首先，这正常吗？以及如何解决这个问题？

我也报告了我的问题here。感谢您的关注。

Answer 1

从 GitHub 复制我的答案以供后代使用。

这不是bug，而是对坐标变换的误解。你必须记住，一个投影中的正方形在另一个投影中不是正方形。

如果坚持使用相同的 CRS，则缓冲区质心的输出等于初始点。但是重新投影的多边形的质心略微偏离，特别是因为您进行的重新投影使几何体在一个方向上倾斜。

How to overcome this problem?

在一个 CRS 中完成所有操作，并在完成后重新投影。

如何解决缓冲后质心点坐标精度不够的问题？

How to solve the lack of precision of the coordinates of the centroid points after a buffer?

geopandas

shapely

pyproj

fiona