|
|
|
|
|
/**
|
|
|
|
|
|
文件描述：空间分析方法
|
|
|
|
|
|
创建时间：2024/4/15 9:52
|
|
|
|
|
|
创建人：Zhaipeixiu
|
|
|
|
|
|
*/
|
|
|
|
|
|
|
|
|
|
|
|
import {getDistance, getElevation} from "@/utils/map/geocomputation.ts";
|
|
|
|
|
|
import {Cartesian3, Viewer} from "cesium";
|
|
|
|
|
|
type ProfileResult = {
|
|
|
|
|
|
distanceArray:number[],
|
|
|
|
|
|
elevationArray:number[],
|
|
|
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
|
|
|
* 两点间地形剖面分析数值计算
|
|
|
|
|
|
* @param viewer 地图Viewer
|
|
|
|
|
|
* @param start 起点
|
|
|
|
|
|
* @param end 终点
|
|
|
|
|
|
* @param interval 线段采样间隔 m为单位
|
|
|
|
|
|
* @return 从起点至终点剖面线的海拔高度数组
|
|
|
|
|
|
*/
|
|
|
|
|
|
export function elevationProfile(viewer: Viewer, start:Cartesian3, end:Cartesian3, interval: number)
|
|
|
|
|
|
{
|
|
|
|
|
|
let breakPointsHeight:number[] = [] //断点的高程m
|
|
|
|
|
|
let distanceFromStart:number[] = [] //断点至起点的距离m
|
|
|
|
|
|
// 计算首尾点距离 m
|
|
|
|
|
|
let totalLen = getDistance(start, end) * 1000
|
|
|
|
|
|
// 获取起点高度
|
|
|
|
|
|
breakPointsHeight.push(getElevation(viewer, start))
|
|
|
|
|
|
distanceFromStart.push(0)
|
|
|
|
|
|
//断点数量
|
|
|
|
|
|
let breakNum = Math.floor(totalLen/interval)
|
|
|
|
|
|
// 如果采样间隔小于首尾点距离，则获取每个断点的坐标 并获取其高度
|
|
|
|
|
|
if(breakNum>=1){
|
|
|
|
|
|
for (let i = 1; i <= breakNum; i++) {
|
|
|
|
|
|
let breakP = Cartesian3.lerp(start, end, i/breakNum, new Cartesian3())
|
|
|
|
|
|
breakPointsHeight.push(getElevation(viewer, breakP))
|
|
|
|
|
|
distanceFromStart.push(getDistance(start,breakP)*1000)
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
|
|
|
|
|
// 获取终点高度
|
|
|
|
|
|
breakPointsHeight.push(getElevation(viewer, end))
|
|
|
|
|
|
distanceFromStart.push(totalLen)
|
|
|
|
|
|
return { distanceArray:distanceFromStart, elevationArray:breakPointsHeight }
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/**
|
|
|
|
|
|
* 折线段地形剖面分析数值计算
|
|
|
|
|
|
* @param viewer 地图Viewer
|
|
|
|
|
|
* @param polyline 折线点数组
|
|
|
|
|
|
* @param interval 线段采样间隔 m为单位
|
|
|
|
|
|
* @return 从折线起点至终点剖面线的海拔高度数组
|
|
|
|
|
|
*/
|
|
|
|
|
|
export function profileAnalyse(viewer: Viewer, polyline:Cartesian3[],interval: number){
|
|
|
|
|
|
let result:ProfileResult = { distanceArray:[], elevationArray:[] }
|
|
|
|
|
|
for (let i = 0; i < polyline.length - 2; i++) {
|
|
|
|
|
|
let temp = elevationProfile(viewer,polyline[i],polyline[i+1],interval)
|
|
|
|
|
|
result.elevationArray = result.elevationArray.concat(temp.elevationArray)
|
|
|
|
|
|
result.distanceArray = result.distanceArray.concat(temp.distanceArray)
|
|
|
|
|
|
}
|
|
|
|
|
|
return result
|
|
|
|
|
|
}
|
