pnpoly函数用(verty[i] > testy) != (verty[j] > testy)判断边是否跨越水平射线y=testy,自动排除水平边;交点横坐标比较testx
直接用
pnpoly函数最省事,但得小心顶点顺序、浮点边界和水平边——它不自动处理点在线段上、也不做包围盒预筛,硬套容易在地图或图形引擎里翻车。为什么
pnpoly的判断条件是(verty[i] > testy) != (verty[j] > testy)这个表达式本质是“边是否跨越水平射线 y = testy”,只对一上一下的边计数,跳过水平边(
verty[i] == verty[j]时两边比较结果相同,整个条件为 false)。它比y0 >= testy && y1 更鲁棒,因为避免了把顶点重复计入两次(比如两个相邻边共用一个顶点,都把它当“上端点”)。注意:它隐含假设多边形顶点按**环绕顺序**给出(顺/逆时针皆可),且首尾自动闭合;如果顶点乱序或未闭合,结果完全不可信。
nvert是顶点总数,不是边数;传入 4 就会自动连vert[3]→vert[0]- 该条件不处理
testy == verty[i]的情况,靠后续交点横坐标比较来规避重复- 若你用的是整数坐标且无精度担忧,可直接用;浮点坐标建议加 epsilon 容差,比如改用
std::abs(verty[i] - testy) > eps配合符号判断
pnpoly里交点横坐标计算为什么写成testx这是从参数方程推出来的交点 x 坐标:边从
(vertx[i], verty[i])到(vertx[j], verty[j]),代入 y = testy 解出 t,再算 x。关键点在于:
- 分母
verty[j]-verty[i]已由前一条件保证非零(否则(verty[i] > testy) != (verty[j] > testy)不成立)- 不等式用
而不是,确保当点恰好落在多边形左边界(x 最小处)时不被误判为内部- 这个比较天然排除了射线起点(即测试点自身)被当作交点的情况
- 若你希望包含左边界,可改为
,但必须同步处理点在线段上的逻辑,否则会漏掉边界点实际调用前必须做的三件事
很多人抄了
pnpoly就跑,结果在凹多边形交点附近、或高缩放地图坐标下频繁误判。真正落地时,这三步不能跳:
- 先用包围盒粗筛:
if (testx max_x || testy max_y) return false;—— 对批量判断能砍掉 70%+ 的无效计算- 若需支持“点在边上也算内部”,不能只改
pnpoly,得单独调用onSegment检查:对每条边,用叉积 + 点积判断测试点是否在线段范围内- 对带孔多边形(如岛屿+湖泊),要分别对外环跑一次
pnpoly得 true,再对每个内环跑一次——任一内环返回 true,就说明点在孔里,最终结果为 false射线法真正的坑不在主循环,而在你怎么准备输入数据和怎么解释输出——顶点顺序错一位、漏掉包围盒、或者把“点在线上”当成“点在内”,都会让结果在生产环境里静默出错。
C++免费学习笔记(深入):立即使用
在学习笔记中,你将探索 C++ 的入门与实战技巧!












