NetLogo中点到线段垂直距离计算问题及代码修正咨询

问题分析 #

原代码的核心问题在于上下文混淆：在ask links块中使用myself时，myself指代的是当前执行的link而非call海龟，导致传递给perpendicular-distance的点参数错误。此外，原函数未处理线段重合的边界情况，且计算的是到直线而非线段的距离，不符合实际需求。

修正方案 #

1. 重构垂直距离计算函数 #

改用坐标参数替代海龟对象，彻底避免上下文依赖，同时处理线段重合、投影超出线段范围的边界情况：

to-report perpendicular-distance [point-x point-y node1-x node1-y node2-x node2-y]
  ; 计算点(x,y)到线段(node1-x,node1-y)-(node2-x,node2-y)的垂直距离
  let dx1 node2-x - node1-x
  let dy1 node2-y - node1-y
  let dx2 point-x - node1-x
  let dy2 point-y - node1-y
  let line-length-sq (dx1 * dx1) + (dy1 * dy1)
  
  ; 处理线段两端点重合的情况：直接返回点到该点的距离
  if line-length-sq = 0 [
    report sqrt ((dx2 * dx2) + (dy2 * dy2))
  ]
  
  ; 计算投影参数t，将t限制在[0,1]内确保投影在线段上
  let t ((dx2 * dx1) + (dy2 * dy1)) / line-length-sq
  let t-clamped clamp t 0 1
  
  ; 计算投影点坐标
  let proj-x node1-x + t-clamped * dx1
  let proj-y node1-y + t-clamped * dy1
  
  ; 返回点到投影点的距离
  report sqrt ((point-x - proj-x) * (point-x - proj-x) + (point-y - proj-y) * (point-y - proj-y))
end

2. 修正连接逻辑 #

从call海龟的上下文遍历所有links，直接获取坐标计算距离，避免上下文混淆，并实现合理的连接逻辑：

to connect-calls
  ask calls [
    let closest-edge nobody
    let min-distance 99999
    let my-x xcor
    let my-y ycor
    
    ; 遍历所有links，计算当前call到线段的垂直距离
    foreach links [
      let current-link ?
      let n1 end1 current-link
      let n2 end2 current-link
      let dist perpendicular-distance my-x my-y [xcor] of n1 [ycor] of n1 [xcor] of n2 [ycor] of n2
      
      if dist < min-distance [
        set min-distance dist
        set closest-edge current-link
      ]
    ]
    
    if closest-edge != nobody [
      ; 方案1：创建垂点node并连接call（模拟连接到线段）
      let n1 end1 closest-edge
      let n2 end2 closest-edge
      let dx1 [xcor] of n2 - [xcor] of n1
      let dy1 [ycor] of n2 - [ycor] of n1
      let dx2 my-x - [xcor] of n1
      let dy2 my-y - [ycor] of n1
      let line-length-sq (dx1 * dx1) + (dy1 * dy1)
      let t ((dx2 * dx1) + (dy2 * dy1)) / line-length-sq
      let t-clamped clamp t 0 1
      let proj-x [xcor] of n1 + t-clamped * dx1
      let proj-y [ycor] of n1 + t-clamped * dy1
      
      create-nodes 1 [
        setxy proj-x proj-y
        set color red
        create-link-with myself [set color blue]
      ]

      ; 方案2：连接到线段的其中一个端点（根据距离选择更近的）
      ; let dist-n1 distance end1 closest-edge
      ; let dist-n2 distance end2 closest-edge
      ; if dist-n1 <= dist-n2 [
      ;   create-link-with end1 closest-edge [set color blue]
      ; ] else [
      ;   create-link-with end2 closest-edge [set color blue]
      ; ]
    ]
  ]
end

关键调整说明 #

• 上下文规避：使用foreach links替代ask links，从call海龟的上下文直接访问所有links，彻底避免myself指代错误。
• 边界处理：新增线段重合时的距离计算，避免除以0错误；通过clamp函数将投影限制在线段范围内，确保计算的是到线段而非无限直线的距离。
• 合理连接：原代码中重复创建node间的link无意义，提供两种符合需求的连接方案：创建垂点node连接，或连接到线段更近的端点。

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