基於百度地圖的區域重疊判斷

原文連接: Fyerl's Blog

前些日的一個小需求：
用戶在後臺劃不規則區域，區域之間不能重疊，如圖

判斷分兩步：
一、判斷多變形是否有相交線段，無則進行第二步判斷（公式）
二、判斷多變形之間是否存在頂點與多邊形的包含關係（BMapLib.GeoUtils.isPointInPolygon）

代碼以下：

/**
 * 點: { lat: xxx, lng: xxx }
 * 線: [{ lat: xxx, lng: xxx }, { lat: xxx, lng: xxx }]
 * 面: [{ lat: xxx, lng: xxx }, { lat: xxx, lng: xxx }, { lat: xxx, lng: xxx }...]
 * */
const { BMap, BMapLib } = window;

/**
 * 線段是否相交
 * seg: [{ lat: xxx, lng: xxx }, { lat: xxx, lng: xxx }]
 * */
function isSegmentsIntersectant(segA, segB) {
  const abc = (segA[0].lat - segB[0].lat) * (segA[1].lng - segB[0].lng) - (segA[0].lng - segB[0].lng) * (segA[1].lat - segB[0].lat);
  const abd = (segA[0].lat - segB[1].lat) * (segA[1].lng - segB[1].lng) - (segA[0].lng - segB[1].lng) * (segA[1].lat - segB[1].lat);
  if (abc * abd >= 0) {
    return false;
  }

  const cda = (segB[0].lat - segA[0].lat) * (segB[1].lng - segA[0].lng) - (segB[0].lng - segA[0].lng) * (segB[1].lat - segA[0].lat);
  const cdb = cda + abc - abd;
  return !(cda * cdb >= 0);
}

/**
 * 判斷兩多邊形邊界是否相交
 */
function isPolygonsIntersectant(plyA, plyB) {
  for (let i = 0, il = plyA.length; i < il; i++) {
    for (let j = 0, jl = plyB.length; j < jl; j++) {
      const segA = [plyA[i], plyA[i === il - 1 ? 0 : i + 1]];
      const segB = [plyB[j], plyB[j === jl - 1 ? 0 : j + 1]];
      if (isSegmentsIntersectant(segA, segB)) {
        return true;
      }
    }
  }
  return false;
}

/**
 * 判斷兩多變形是否存在點與區域的包含關係(A的點在B的區域內或B的點在A的區域內)
 */
function isPointInPolygonBidirectional(plyA, plyB) {
  const [pA, pB] = [[], []];
  plyA.forEach((item) => {
    pA.push(new BMap.Point(item.lng, item.lat));
  });

  plyB.forEach((item) => {
    pB.push(new BMap.Point(item.lng, item.lat));
  });


  let [a, b] = [false, false];
  a = pA.some(item => BMapLib.GeoUtils.isPointInPolygon(item, new BMap.Polygon(pB)));
  if (!a) {
    b = pB.some(item => BMapLib.GeoUtils.isPointInPolygon(item, new BMap.Polygon(pA)));
  }

  return a || b;
}


/**
 * 判斷多邊形是否重疊
 * */
export function isPolygonsOverlap(plyA, plyB) {
  return isPolygonsIntersectant(plyA, plyB) || isPointInPolygonBidirectional(plyA, plyB);
}