无人机照片贴地实现

来源：cnblogs　　作者：当时明月在曾照彩云归　　时间：2024/3/27 8:39:03　　对本文有异议

1. 引言

将无人机拍摄拍摄的影像作为底图可以实现快速甚至实时的更新，这对于应急指挥（如，森林防火）有着很大的意义

常规的做法是使用无人机拍摄一组照片，然后将这一组照片放入生产软件（如，ContextCapture）生产正射影像，然后地图切片发布

这里记录的是另一种做法，直接将无人机照片贴地，这种方式速度很快，几乎无需后台处理，照片细节完整保留，当然这种方式存在着误差，尤其地势起伏较大的地方可能并不适用

2. 实现过程

这里先有以下设定：

无人机拍摄的影像都是竖直向下的类似正射的照片（不是类正射的照片贴地视角会很奇怪）
无人机的朝向是正北方向
地面是接近于平面

无人机照片的位置关系如下图所示，可以根据这些数据快速简单地算出照片的贴地范围

笔者这里使用的无人机是大疆无人机型号H20T，参考大疆官网给出的H20T的数据：技术参数_禅思 H20 系列无人机负载_DJI大疆行业应用

这里的DFOV是82.9°（广角相机）

上图中的width和height是图片的宽高（单位：px），用以计算与斜边的角度

上图中的relative altitude是相对高度，表示地面与无人机（相机）的距离

width、height和relative altitude都可以直接从照片中读取

要进行计算还有一点需要做的，就是把经纬度坐标转换为投影坐标系下的坐标，计算后再将投影坐标数据转换为经纬度，这样就完成了计算

详细的代码实现如下：

/**
 * 根据照片的数据信息，将照片粘贴到地图上（最好是正射照片，无人机指向正北方）
 * @param imgWidth 图片宽度（像素）
 * @param imgHeight 图片高度（像素）
 * @param distance 相机与地面的距离（米）
 * @param url 图片的URL
 * @param latitude 照片中心点的纬度
 * @param longitude 照片中心点的经度
 * @param dfov 相机的视场角（默认为82.9）
 */
function pasteImageByHeight(imgWidth, imgHeight, distance, url, latitude, longitude, dfov = 82.9) {
  console.log(imgWidth, imgHeight, distance, url, latitude, longitude, dfov)
  const thisDistance = parseFloat(distance)
  if (isNaN(thisDistance) || thisDistance <= 0) {
    globalAlert("照片与地面的距离信息有误")
    return
  }
  const thisLatitude = parseFloat(latitude)
  const thisLongitude = parseFloat(longitude)
  if (isNaN(thisLatitude) || isNaN(thisLongitude)) {
    globalAlert("照片的经纬度信息有误")
    return
  }
  const hypotenuse = distance * Math.tan(dfov / 2 * Math.PI / 180)
  const angle = Math.atan(imgHeight / imgWidth)
  const halfWidth = hypotenuse * Math.sin(angle)
  const halfHeight = hypotenuse * Math.cos(angle)
  const projCrs = "+proj=utm +zone=51 +datum=WGS84 +units=m +no_defs"
  const projectedPos = Proj4("EPSG:4326", projCrs, [thisLongitude, thisLatitude])
  const [x, y] = projectedPos
  const leftTop = Proj4(projCrs, "EPSG:4326", [x + halfWidth, y + halfHeight])
  const rightBottom = Proj4(projCrs, "EPSG:4326", [x - halfWidth, y - halfHeight])
  console.log(leftTop, rightBottom)
  console.log({ xmin: rightBottom[0], xmax: leftTop[0], ymin: rightBottom[1], ymax: leftTop[1] })
  const imageLayer = new mars3d.layer.ImageLayer({
    name: "无人机航拍影像",
    url: url,
    rectangle: { xmin: rightBottom[0], xmax: leftTop[0], ymin: rightBottom[1], ymax: leftTop[1] },
    zIndex: 20,
    opacity: 0.6
  })
  map.addLayer(imageLayer)
  
  map.flyToExtent({
    xmin: longitude - 0.001,
    xmax: longitude + 0.001,
    ymin: latitude - 0.001,
    ymax: latitude + 0.001
  }, {
    duration: 1
  })
  return imageLayer.id as string
}