[cesium] 卫星扫描 + 地面雷达 + 实时通信过境 效果

学习cesium时看到github有一个可视化卫星轨道的项目，闲着没事模拟了一下可是效果。做的比较简陋。

更多在线效果：在线地址 gis用户

效果图

示例代码：

**
 * 卫星通过效果
 */
const SatellitePassEffect = class {
  constructor(opt) {
    // this._state = CV.INITIALIZED
    this.viewer = opt.viewer
    this.clock = opt.viewer.clock;
    this.handler = new Cesium.ScreenSpaceEventHandler(this.viewer.canvas)
    this.layer = new Cesium.CustomDataSource('cv-sate')
    this.viewer && this.viewer.dataSources.add(this.layer)
    this._init();
  }

  _init() {
    //时间速率
    this.STSpeed = 10;
    // 卫星实体
    this.satelliteAll = [];
    //通信
    this.Transits = [];
    //格式数据存储
    //查看方式区分开 以免冲突
    this.passTime = [];
    this.passTime_old = [];
    //刷新的时间
    this.renderTime = null;
    //拾取的对象
    this.pickedFeature = null;
    //标识符
    this.iden = true;
    //实体
    this.EntityObj = [];
    //线实体
    this.lineEntity = [];

    this._bindEvent()
    this._loadCZML();
  }

  _bindEvent() {
    this.handler.setInputAction(e => {
      try {
        if (!e.position) {
          return false;
        }
        this.pickedFeature = this.viewer.scene.pick(e.position);
        if (!Cesium.defined(this.pickedFeature) && this.pickedFeature == undefined) {
          return false;
        };
        if (this.pickedFeature.id.description == undefined) { //自己创建的
          return false;
        };
        let f_name = this.pickedFeature.id.name, f_position, position;
        if (this.pickedFeature.id.idea == "radar") {
          position = this.pickedFeature.id.position.getValue(this.clock.currentTime);
          this.pickedFeature.type = "radar";
        } else {
          position = this.pickedFeature.id.position.getValue(this.clock.currentTime);
          this.pickedFeature.type = "satellite";
        }
        f_position =  Cesium.Ellipsoid.WGS84.cartesianToCartographic(position)
        this.pickedFeature.id.description = this.infoTable_2(f_name, f_position);
        console.log(f_name, f_position);
        //viewer.selectedEntity = selectedEntity;
        this.pickedFeature.id.name = f_name;
        this.pickedFeature.temp = this.pickedFeature.id.description + "<h2> Passes <h2>";
        this.iden = false; //点击事件改变标识符
      } catch (error) {
        console.log(error);
      }
    }, Cesium.ScreenSpaceEventType.LEFT_CLICK)
  }

  _loadCZML() { //load数据
    this.viewer.dataSources.add(Cesium.CzmlDataSource.load('data/test.czml')).then((dataSource) => {
      this.dataSource = dataSource;
      this.radar = dataSource.entities.getById("Facility/AGI");
      let satellite1 = dataSource.entities.getById("Satellite/ISS");
      let satellite2 = dataSource.entities.getById("Satellite/Geoeye1");
      let transit1 = dataSource.entities.getById("Facility/AGI-to-Satellite/ISS");
      let transit2 = dataSource.entities.getById("Facility/AGI-to-Satellite/Geoeye1");
      this.satelliteAll = [satellite1, satellite2];
      this.Transits = [transit1, transit2];
      this.radar.idea = "rader";
      try {
        this.scan();
        this.satelliteInfo();
        //_self.createEntity([_self.radar]);
        this.communication();
      } catch (e) {
        console.log(e);
      }
    })
  }
  /**
   * 卫星过境通信效果
   */
  communication() {
    this.p_line = []; //创建线的点

    if (this.radar == null) { //雷达
      return false;
    }

    let r_position = this.radar.position.getValue(this.clock.currentTime);
    let r_point = Cesium.Ellipsoid.WGS84.cartesianToCartographic(r_position)
    this.rr_point = [parseInt(r_point.longitude / Math.PI * 180), parseInt(r_point.latitude / Math.PI * 180)];
    if (this.satelliteAll.length == 0) {
      return false;
    }
    //遍历卫星
    for (let i in this.satelliteAll) {
      this.p_line[i] = [];
      let sate = this.satelliteAll[i];
      let s_position = sate.position.getValue(this.clock.currentTime);
      this.s_point =  Cesium.Ellipsoid.WGS84.cartesianToCartographic(s_position)
      this.p_line[i].push(r_position.clone());
      this.p_line[i].push(s_position.clone());
      let lineObj = CV.E.createDynamicPolyline({ positions: this.p_line[i], width: 1 })

      this.lineEntity.push(this.layer.entities.add(lineObj));
    }

  }
  /** 
  *  判断一个点是否在圆的内部 
  *  @param point  测试点坐标 
  *  @param circle 圆心坐标 
  *  @param r 圆半径 
  *  返回true为真，false为假 
  *  */
  pointInsideCircle(point, circle, r) {
    if (r === 0) return false
    var dx = circle[0] - point[0]
    var dy = circle[1] - point[1]
    return dx * dx + dy * dy <= r * r
  }
  /**
   * 卫星信息
   */
  satelliteInfo() {
    if (this.Transits.length == 0) {
      return false;
    };
    this.formatTransit();
    this.selectedEntity = new Cesium.Entity();
    this.selectedEntity.name = "PASS";
    this.selectedEntity.description = this.infoTable_1(dayjs(Cesium.JulianDate.addHours(this.clock.currentTime, -8, new Cesium.JulianDate())).format("YYYY-MM-DD HH:mm:ss"));
    this.viewer.selectedEntity = this.selectedEntity;

    this.clock.onTick.addEventListener(clock => {
      if (!clock.shouldAnimate) return;
      if (this.iden) this.selectedEntity.description = this.infoTable_1(this.pass(clock));//标识符  进来展示所有卫星信息
      if (!this.iden) {
        if (this.pickedFeature == null) return;
        if ("radar" == this.pickedFeature.type) {
          this.pickedFeature.id.description = this.pickedFeature.temp + this.infoTable_1(this.pass(clock));
        } else {
          let position = this.pickedFeature.id.position.getValue(clock.currentTime);
          let f_position = Cesium.Ellipsoid.WGS84.cartesianToCartographic(position)
          this.pickedFeature.id.description = this.infoTable_2(this.pickedFeature.id.name, f_position) + ' <h2> Passes </h2>' + this.infoTable_3(this.pass(clock), this.pickedFeature.id.name);
          this._addRadarScen(position)
        }
      }
      //判断是否在地面雷达通信范围
      if (this.satelliteAll.length == 0) {
        return false;
      }
      for (let i in this.satelliteAll) { //遍历卫星
        let sate = this.satelliteAll[i];
        let position = sate.position.getValue(clock.currentTime);
        let s_point =  Cesium.Ellipsoid.WGS84.cartesianToCartographic(position)
        this.p_line[i].pop();
        this.p_line[i].push(position.clone());
        this.ss_point = [parseInt(s_point.longitude / Math.PI * 180), parseInt(s_point.latitude / Math.PI * 180)]
        let flag = this.pointInsideCircle(this.ss_point, this.rr_point, 30);
        if (flag) {
          this.lineEntity[i].show = true;
        } else {
          this.lineEntity[i].show = false;
        }
      }
    })
    //时间轴结束
    this.clock.onStop.addEventListener(clock => {
      if (this.Transits.length == 0) {
        return false;
      };
      //格式化卫星数据  
      this.formatTransit();
    });
  }

   //新增雷达探照
   _addRadarScen(position) {

      //console.log(CV.U.courseAngle(this.radar.position.getValue(viewer.clock.currentTime),position))

  }
  /**
   * 格式化通信数据
   */
  formatTransit() {
    if (this.Transits.length == 0) {
      return false;
    };
    this.passTime = [], this.passTime_old = []; //查看方式区分开 以免冲突
    for (let i in this.Transits) {
      let transit = this.Transits[i]
      let intervals = [], intervals_old = [];
      let n_interval = transit.availability._intervals;
      for (let ii in n_interval) {
        let interval = n_interval[ii]
        let start = dayjs(Cesium.JulianDate.addHours(interval.start, -8, new Cesium.JulianDate())).format("YYYY-MM-DD HH:mm:ss");
        let stop = dayjs(Cesium.JulianDate.addHours(interval.stop, -8, new Cesium.JulianDate())).format("YYYY-MM-DD HH:mm:ss");
        intervals.push({ name: transit.name, "startTime": start, "stopTime": stop, "interval": dayjs(stop).diff(dayjs(start), 'millisecond') });
        intervals_old.push({ name: transit.name, "startTime": start, "stopTime": stop, "interval": dayjs(stop).diff(dayjs(start), 'millisecond') });
      }
      this.passTime.push(intervals);
      this.passTime_old.push(intervals_old);
    }
  }
  /**
   * 创建雷达实体
   */
  createEntity(radars) {
    try {
      if (radars.length == 0) {
        return false;
      }
      for (let i in radars) {
        let radar = radars[i], l, r;
        let positions = radar.position.getValue(this.clock.currentTime);
        this.r_point = positions;
        if (positions.length == 0) {
          return false;
        };
        let cartographic =  Cesium.Ellipsoid.WGS84.cartesianToCartographic(positions)
        let lat = Cesium.Math.toDegrees(cartographic.latitude), lon = Cesium.Math.toDegrees(cartographic.longitude), height = cartographic.height;
        //radarscan
        r = new Cesium.HeadingPitchRoll(Cesium.Math.toRadians(90),
          Cesium.Math.toRadians(0), Cesium.Math.toRadians(0));
        l = Cesium.Cartesian3.fromDegrees(lon, lat, height);
        this.EntityObj.push(this.layer.entities.add(CV.E.getCustomRadar(l, r)));
      }
    } catch (e) {
      console.log(e);
    }
  }
  /**
   * 添加扫描物
   */
  scan() {
    if (this.satelliteAll.length == 0) {
      return false;
    };
    for (let i in this.satelliteAll) {
      let entity = this.satelliteAll[i]
      let cartesian = entity.position.getValue(this.clock.currentTime);
      let position = CV.T.transformCartesianToWSG84(cartesian)
      //let positions = _self.mouseManager.worldToLonlat(cartesian);
      this.bindScan(position, entity);
    }
  }

  pass(clock) { //当前时间
    let currentTime = dayjs(Cesium.JulianDate.addHours(clock.currentTime, -8, new Cesium.JulianDate())).format("YYYY-MM-DD HH:mm:ss");
    return currentTime;
  }
  index_rm(n) {
    if (this.passTime_old.length == 0) {
      return false;
    };
    this.passTime_old[n].splice(0, 1); //删除第一个
  }
  /**
  * 绑定扫描物
  */
  bindScan(positions, entityObj) {
    //let modelMatrix = this.Primitives.countModelMatrix(positions);
    //this.scanEntity = new entityFactory({type:"createScan",data:{modelMatrix:modelMatrix,positions:positions,v:this.CoreV}});
    this.scanEntity = CV.E.createDynamicCylinder({ positions: positions, entity: entityObj, viewer: this.viewer, cylinder: { legnth: 600000, slices: 4, bottomRadius: 600000 / 2 } })
    this.EntityObj.push(this.layer.entities.add(this.scanEntity));
  }
  /*
 * table 1
 */
  infoTable_1(currentTime) {
    try {
      let _self = this, renderTime = _self.renderTime;
      if (_self.passTime_old.length == 0) {
        return false;
      };
      var tr = "", table = `<table class="cesium-infoBox-defaultTable"><thead><tr><th>卫星</th><th>倒计时(ms)</th><th>通信开始(date)</th><th>通信结束(date)</th><th>通信时长(ms)</th></tr></thead><tbody>`;
      for (var n in _self.passTime_old) {
        if (_self.passTime_old[n].length == 0) continue;
        var interval_pass = _self.passTime_old[n][0]; //始终取第一个 
        renderTime = dayjs(interval_pass.startTime).diff(dayjs(currentTime));
        if (renderTime <= 0) {
          if (renderTime <= -(interval_pass.interval)) {
            _self.index_rm(n);
          } else {
            renderTime = "PASS";
          }
        }
        tr += `<tr><td>${interval_pass.name}</td><td>${renderTime}</td><td>${interval_pass.startTime}</td><td> ${interval_pass.stopTime}</td><td> ${interval_pass.interval}</td></tr>`;
      }
      return table + tr + `</tbody></table>`;
    } catch (e) {
      console.log(e);
    }
  }
  /*
  * table 2
  */
  infoTable_2(f_name, cartesian) {
    if (f_name == undefined && cartesian == undefined) {
      return false;
    };
    let tr = "", table = `<h2> Position </h2><table class="cesium-infoBox-defaultTable"><thead><tr><th>Name</th><th>Latitude</th><th>Longitude</th><th>Elevation</th></tr></thead><tbody>`;
    let f_point = [parseInt(cartesian.longitude / Math.PI * 180), parseInt(cartesian.latitude / Math.PI * 180)];
    tr = `<tr><td>${f_name}</td><td>${f_point[0]}°</td><td>${f_point[1]}°</td><td> ${parseInt(cartesian.height)}</td></tr>`;
    return table + tr + `</tbody></table>`;
  }
  /*
  * table 3
  */
  infoTable_3(currentTime, featureName) {
    let _self = this, renderTime = _self.renderTime;
    if (_self.passTime.length == 0 && featureName == undefined) {
      return false;
    };
    let t_interval = function () {
      for (var i in _self.passTime) { if (_self.passTime[i][0].name.indexOf(featureName) != -1) return _self.passTime[i]; }
    }
    let intervals = t_interval();
    var tr = "", table = `<table class="cesium-infoBox-defaultTable"><thead><tr><th>卫星</th><th>倒计时(ms)</th><th>通信开始(date)</th><th>通信结束(date)</th><th>通信时长(ms)</th></tr></thead><tbody>`;
    for (let i in intervals) {
      let interval = intervals[i]
      renderTime = dayjs(interval.startTime).diff(dayjs(currentTime));
      if (renderTime <= 0) renderTime = 0;
      tr += `<tr><td>${interval.name}</td><td>${renderTime}</td><td>${interval.startTime}</td><td> ${interval.stopTime}</td><td> ${interval.interval}</td></tr>`;
    }
    return table + tr + `</tbody></table>`;
  }
  /**
  * 向后飞行
  */
  back() {
    this.STSpeed = -50;
    this.clock.multiplier = this.STSpeed;
  }
  /**
   * 向前飞行
   */
  forward() {
    this.STSpeed = 50;
    this.clock.multiplier = this.STSpeed;
  }
  /**
   * 飞行速度
   */
  speed(speedNum) {
    console.log(speedNum);
    this.clock.multiplier = speedNum;
  }
  /**
   * 开始飞行
   */
  start() {
    this.clock.shouldAnimate = true;
  }
  /**
   * 暂停飞行
   */
  pause() {
    this.clock.shouldAnimate = false;
  }
  /**
   * 清除飞行
   */
  remove() {
    if (this.EntityObj.length == 0) {
      return false;
    }

    this.layer.entities.removeAll();
    this.clock.dataSources.remove(this.dataSource);
    this.clock.clock.shouldAnimate = false;
    this.handlerAction.destroy();
    this.handlerAction = null;
    this.EntityObj = [];
  }
}