如何在EPSG:3857的OpenLayers中绘制200公顷的圆形？

我需要在采用EPSG:3857投影的OpenLayers地图上绘制一个覆盖200公顷面积的圆形Polygon要素。找到了一段参考代码，但看不懂其中的数学逻辑，也不知道怎么适配到自己的代码里。

参考代码 #

new ol.geom.Circle(centerLongitudeLatitude,200*1000/ol.proj.getPointResolution('EPSG:3857', 1, centerLongitudeLatitude))

当前代码 #

const geometry = myFeature.getGeometry();
const circle = new Circle(geometry.getCoordinates(), myFeature.get("SIZE_HECTARES"));
myFeature.setGeometry(fromCircle(circle));

核心逻辑拆解 #

1. EPSG:3857投影的变形问题：Web墨卡托投影（EPSG:3857）会随纬度产生距离变形——越靠近两极，地图上1个单位对应的实际地面距离越短。直接用米作为半径画圆，实际面积会出现偏差。
2. 参考代码的作用：ol.proj.getPointResolution('EPSG:3857', 1, center)返回中心点位置，1个地图单位对应的实际地面距离（米）。用**实际需要的地面距离（米）**除以这个值，就能得到适配投影变形的地图单位半径，保证实际面积准确。
3. 面积转半径的计算：200公顷等于200×10000=2000000平方米，根据圆的面积公式 S=πr²，反推实际地面半径 r = √(S/π)。

修改后的代码 #

const geometry = myFeature.getGeometry();
const center = geometry.getCoordinates();

// 1. 从面积计算实际地面半径（米）
const areaHectares = myFeature.get("SIZE_HECTARES"); // 或直接写200
const areaSquareMeters = areaHectares * 10000;
const radiusMeters = Math.sqrt(areaSquareMeters / Math.PI);

// 2. 根据投影变形修正半径，得到地图单位的半径
const projection = 'EPSG:3857';
const pointResolution = ol.proj.getPointResolution(projection, 1, center);
const mapUnitRadius = radiusMeters / pointResolution;

// 3. 创建圆形并转为Polygon
const circle = new ol.geom.Circle(center, mapUnitRadius);
myFeature.setGeometry(ol.geom.Polygon.fromCircle(circle));

关键说明 #

• 如果fromCircle是ol.geom.Polygon.fromCircle的简写，直接用全称可避免作用域问题。
• 若SIZE_HECTARES是动态传入的值，直接保留myFeature.get("SIZE_HECTARES")即可适配不同面积需求。

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