Python Shapely/PyProj面积计算与Google Earth引擎或Geojson.io不匹配

计算多边形的面积时，使用shapely/pyproj的结果与google earth engine或geojson.io的结果截然不同。我认为这与投影有关，当我研究它时，我发现了一些例子，比如Alders等面积投影，但在PyProj的2.1.3版中，这些似乎不是一个选项。我还想知道google earth engine和geojson.io是否有一种根据多边形区域自动选择适合区域的区域投影的方法，但我不确定他们会如何选择这些

例如：

coords=[-97.59238135821987，43.474565304017， -97.59244690469288, 43.47962399877412, -97.59191951546768, 43.47962728271748, -97.59185396090983, 43.47456565304017, -97.59238135821987, 43.47456565304017] projection=partialpyproj.transform，pyproj.Projinit='epsg:4326'，pyproj.Projinit='epsg:3857' shapely.ops.transformprojection、shapely.geometry.Polygonsample\u coords.area 返回45573.993884405005 m^2的面积

---

Google Earth引擎返回23944.14737277293和Geojson.io返回23997.77

我发现了这个问题，并与shapely提出了相应的github问题

简而言之，shapely使用的是一种不同的方法，其来源不清楚，其他人使用的是JPL论文中描述的方法