feat:基于最大斜率的两点间通视算法实现

src/components/ProfileAnalysis.vue

@@ -4,7 +4,7 @@
   创建人：Zhaipeixiu
 -->
 <script setup lang="ts">
-import {onMounted, ref} from 'vue'
+import {ref} from 'vue'
 import * as echarts from 'echarts'
 const emit = defineEmits(['confirmDia'])
 let showDia = ref<boolean>(false)
@@ -22,12 +22,12 @@ const confirmDia = (): void => {
 }
 
 /**
- *  绘制折线图
+ *  绘制折线图（航线碰撞检测）
  * @param xData x数组
  * @param yData y数组，以面积线绘制
  * @param yData2 y数组，以折线绘制
  */
-const drawChart = (xData:number[],yData:number[],yData2:number[]) => {
+const drawChart_AirlineDetect = (xData:number[],yData:number[],yData2:number[]) => {
   let myChart = echarts.init(document.getElementById('profileChart'))
   // 绘制图表
   myChart.setOption({
@@ -119,13 +119,106 @@ const drawChart = (xData:number[],yData:number[],yData2:number[]) => {
     ]
   });
 }
+/**
+ *  绘制折线图（地形剖面）
+ * @param xData x数组
+ * @param yData y数组，以面积线绘制
+ */
+const drawChart_TerrainProfile = (xData:number[],yData:number[]) => {
+  let myChart = echarts.init(document.getElementById('profileChart'))
+  // 绘制图表
+  myChart.setOption({
+    legend:{
+      show: true,
+      type: 'plain',
+      top: '7%',
+      data:[
+        {
+          name: 'groundLine',
+          itemStyle: 'inherit',
+          lineStyle: 'inherit',
+        },
+        {
+          name: 'airLine',
+          itemStyle: 'inherit',
+          lineStyle: 'inherit',
+        }]
+    },
+    tooltip: {
+      show: true,
+      trigger: 'axis',
+      axisPointer: {
+        type: 'cross'
+      },
+      formatter:'地表高度: {c0}'
+    },
+    xAxis: {
+      data: xData,
+      name: '距离',
+      nameTextStyle: {
+        fontWeight:'bolder',
+        fontSize: 14
+      },
+      nameLocation: 'end',
+      axisLine:{
+        onZero: false,
+        show: true,    // 是否显示坐标轴轴线
+        symbol: ['none', 'arrow'],
+        symbolSize: [7, 10]
+      },
+      axisLabel: {
+        formatter: '{value} m',
+        margin: 5,
+      },
+      axisTick: {
+        show: true,    // 是否显示坐标轴刻度
+        inside: true,     // 坐标轴刻度是否朝内，默认朝外
+        alignWithLabel: true,
+        lineStyle: {
+          color: '#000000',     //刻度线的颜色
+          type: 'solid',    //坐标轴线线的类型（solid实线类型；dashed虚线类型；dotted点状类型）
+        },
+      }
+    },
+    yAxis: {
+      type: 'value',
+      name: '高度',
+      nameTextStyle: {
+        fontWeight:'bolder',
+        fontSize: 14
+      },
+      nameLocation: 'end',
+      position: 'left',
+      axisLabel: {
+        formatter: '{value} m'
+      },
+      axisLine: {
+        show: true,
+        symbol: ['none', 'arrow'],
+        symbolSize: [7, 10]
+      }
+    },
+    series: [
+      {
+        name:'groundLine',
+        type: 'line',
+        data: yData,
+        areaStyle: {
+          color: '#37a5fb',
+          opacity: 0.5
+        }
+      }
+    ]
+  });
+}
 
 // vue3中规定，使用了 <script setup> 的组件是默认私有的：
 // 一个父组件无法访问到一个使用了 <script setup> 的子组件中的任何东西，除非子组件在其中通过 defineExpose 宏显式暴露
 defineExpose({
   openDia,
   closeDia,
-  drawChart,
+  drawChart_AirlineDetect,
+  drawChart_TerrainProfile
 })
 
 

src/utils/map/SpatialAnalysis.ts

@@ -59,3 +59,44 @@ export function profileAnalyse(viewer: Viewer, polyline:Cartesian3[],interval: n
  }
  return result
 }
+
+/**
+ * 两点间通视分析（基于最大斜率的算法） <br>
+ * 返回值：通视为true，不通视为false，异常返回undefined
+ * @param viewer 地图
+ * @param viewpoint 视点
+ * @param target 目标点
+ * @param h 视点地面抬高，默认为1m
+ * @param breakNum 采样间隔，默认为100个断点
+ */
+export function visibilityAnalyse(viewer: Viewer, viewpoint:Cartesian3, target:Cartesian3,
+                                  h= 1, breakNum = 100) {
+ // 获取视点高度和目标点高度
+ let viewpointH = getElevation(viewer, viewpoint) + h
+ let targetH = getElevation(viewer, target)
+ if (viewpointH === -9999.2024 + h ) {
+  console.log("无法获取视点海拔高度！")
+  return undefined
+ }
+ if (targetH === -9999.2024 ) {
+  console.log("无法获取目标点海拔高度！")
+  return undefined
+ }
+ // 计算首尾点距离 m
+ let totalLen = getDistance(viewpoint, target) * 1000
+ // 计算首尾点斜率
+ let visibleK = (targetH - viewpointH)/totalLen
+
+ // 逐个计算断点与视点的斜率
+ for (let i = 1; i <= breakNum; i++) {
+  let breakP = Cartesian3.lerp(viewpoint, target, i/breakNum, new Cartesian3())
+  let breakPH = getElevation(viewer, breakP)  //断点海拔
+  let breakPDis = getDistance(viewpoint,breakP)*1000  //断点与视点的距离
+  // 计算断点与视点的斜率
+  let breakPK = (breakPH-viewpointH)/breakPDis
+  if(breakPK>visibleK){
+   return false
+  }
+ }
+ return true
+}

src/utils/map/geocomputation.ts

@@ -142,7 +142,7 @@ function getPolygonArea(vertexs: Cartesian3[]) {
 }
 
 /**
- * 获取坐标点的海拔高度
+ * 获取坐标点的海拔高度 m为单位。 -9999.2024为错误值。
  * @param viewer 地图Viewer
  * @param pos 坐标点 Cartographic|Cartesian3|[lon,lat]
  */
@@ -150,14 +150,14 @@ export function getElevation(viewer: Viewer, pos: Cartographic|Cartesian3|number
   let cartographic = undefined
   if(pos instanceof Array){
      cartographic = Cartographic.fromDegrees(Angle.degree2rad(pos[0]), Angle.degree2rad(pos[1]))
-    return viewer.scene.globe.getHeight(cartographic)??-99.2024
+    return viewer.scene.globe.getHeight(cartographic)??-9999.2024
   }
   else if(pos instanceof Cartesian3){
     cartographic = Cartographic.fromCartesian(pos)
-    return viewer.scene.globe.getHeight(cartographic)??-99.2024
+    return viewer.scene.globe.getHeight(cartographic)??-9999.2024
   }
   else{
-    return viewer.scene.globe.getHeight(pos)??-99.2024
+    return viewer.scene.globe.getHeight(pos)??-9999.2024
   }
 }