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。关键点在于：

C++14

C++14 对 C++11 的修正与增强版本，适合旧系统维护和较老工具链兼容。

下载

立即学习“C++免费学习笔记（深入）”；

• 分母 verty[j]-verty[i] 已由前一条件保证非零（否则 (verty[i] > testy) != (verty[j] > testy) 不成立）
• 不等式用 而不是 <code>，确保当点恰好落在多边形左边界（x 最小处）时不被误判为内部
• 这个比较天然排除了射线起点（即测试点自身）被当作交点的情况
• 若你希望包含左边界，可改为 ，但必须同步处理点在线段上的逻辑，否则会漏掉边界点

实际调用前必须做的三件事

很多人抄了 pnpoly 就跑，结果在凹多边形交点附近、或高缩放地图坐标下频繁误判。真正落地时，这三步不能跳：

• 先用包围盒粗筛：if (testx max_x || testy max_y) return false; —— 对批量判断能砍掉 70%+ 的无效计算
• 若需支持“点在边上也算内部”，不能只改 pnpoly，得单独调用 onSegment 检查：对每条边，用叉积 + 点积判断测试点是否在线段范围内
• 对带孔多边形（如岛屿+湖泊），要分别对外环跑一次 pnpoly 得 true，再对每个内环跑一次——任一内环返回 true，就说明点在孔里，最终结果为 false

射线法真正的坑不在主循环，而在你怎么准备输入数据和怎么解释输出——顶点顺序错一位、漏掉包围盒、或者把“点在线上”当成“点在内”，都会让结果在生产环境里静默出错。