PostgreSQL 中的多边形区域与 QGIS 不同

Question

我下载了 OpenStreetMap 数据并使用 osm2pgsql 上传到本地 PostgreSQL 数据库。 由于 ST_Union function 只有一条记录，我创建了一个新表（目前）。

CREATE TABLE IF NOT EXISTS m_temp.budapest_merge
(
    id text COLLATE pg_catalog."default" NOT NULL,
    st_union geometry(Polygon,3857),
    ter double precision,
    CONSTRAINT budapest_merge_pkey PRIMARY KEY (id)
)

然后我检查了多边形的面积。

SELECT ST_Area(st_union) FROM m_temp.budapest_merge

结果是1151156216.758265 sqm ，这对我来说很奇怪，因为我期待大约520 sqkm平方米。 我在 QGIS 中打开了 PostgreSQL 图层，使用字段计算器工具的$area function 计算了面积。 结果是526031443.36087507 sqm ，这也符合我的期望和现实。

我的问题是这里有什么问题？ QGIS 项目位于同一投影中 (EPSG:3857)。

我也知道命名几何列 st_union 并不是很好，但这不是这里的主要问题。

Answer 1

在 Qgis 中， $area是（如果定义的话）椭圆体区域，即真实的地面区域，而area(geometry)是平面区域，即受几何投影影响很大。

在 Postgis 中， st_area总是使用几何投影，所以它与 QGIS area(geometry)等价。

话虽如此，当您远离赤道时， 3857会严重扭曲距离和面积。 使用此投影来计算距离或面积是错误的。