如何从KML作业区多边形计算机器人导航用中心线？

需求背景 #

手里有一份拖拉机生成的KML文件，里面包含多个固定宽度为3米的农田作业区多边形，需要计算这些多边形的中心线，作为1.5米宽玉米地除草机器人的导航航点。

过往尝试的局限 #

之前用PHP写代码，通过计算多边形对侧坐标的中点来生成中心线，但这个方法有明显缺陷：仅适用于对侧坐标数量为偶数的情况，碰到高度弯曲的多边形时完全失效，生成的中心线偏差极大。

当前实现的初步方案 #

开发了一套PHP脚本，核心思路如下：

1. 先检测多边形的长边（也就是作业区的长度方向边）
2. 针对每个坐标点，计算相对于前一个坐标点90°方位、距离1.5米的点
3. 将这些计算出的点串联，形成作业区的中心线

目前这套脚本生成的结果准确性尚可，但还有优化空间。

可优化方向 #

• 长边检测逻辑优化：复杂弯曲多边形中，可结合边的长度阈值+相邻边的角度变化来精准识别作业区的长度方向边，减少短边误判
• 中心线平滑处理：在生成的航点间加入插值（比如线性插值、贝塞尔曲线），让中心线更贴合弯曲作业区的轮廓，提升机器人导航的流畅度
• 坐标精度校正：针对KML的WGS84坐标，短距离下地球曲率影响可忽略，但高精度场景可加入校正逻辑，避免累积误差
• 闭合多边形适配：对于闭合的作业区多边形，要确保中心线首尾衔接，避免出现断点或重叠

内容的提问来源于stack exchange，提问作者Bart De Reu