erdas lps 在提取DTM 时图像没有匹配的点怎么解决

Question

Answer 1

以下回答参考 皆我百晓生、券券喵儿 等免费微信小程序相关内容作答，并由本人整理回复。

要解决这个问题，首先需要了解如何使用erdas I软件进行DTM（数字地面模型）提取。以下是一个基本步骤和示例代码：

1. 安装并配置环境：确保已经安装了erdas I软件，并且环境中包含了必要的库。

2. 导入图像：从给定的图像文件中加载图像。

package main

import (
	"fmt"
	"image"
)

func main() {
	img := image.NewRGBA(image.Rectangle{Rect: image.Point{X: 0, Y: 0}})
	fmt.Println(img)
}

1. 定义控制点：为每个像素分配一个唯一的ID，并将它们存储在一个二维数组中。

type controlPoint struct {
	id int
	x, y float64
}

var controlPoints = make([]controlPoint, 27700)

1. 计算残差：对于每个像素，使用相应的控制点来计算其相对于目标值的距离，并将其记录到residuals变量中。

for i, id := range controlPoints {
    x, y := img.Bounds().Min.X + id.x*img.Bounds().Width, img.Bounds().Min.Y + id.y*img.Bounds().Height
    dx, dy := x-img.Bounds().Min.X, y-img.Bounds().Min.Y
    residuals[id] = dx*dx + dy*dy
}

1. 创建DEM：根据残差计算得到的目标值和DEM背景值，生成DEM。

func createDEM() error {
    var dem image.Image
    // 假设我们有DEM的原始文件，然后读取它。
    file, err := os.Open("dem3.img")
    if err != nil {
        return err
    }
    defer file.Close()

    reader := bufio.NewReader(file)
    readerBytes, err := io.ReadAll(reader)
    if err != nil {
        return err
    }

    demData, err := image.Decode(bytes.NewReader(readerBytes))
    if err != nil {
        return err
    }

    // 这里是你的代码来将DEM转换成你需要的形式，例如作为三维数组或矩阵。

    // 按照你的需求调整这个函数以生成DEM。
    return nil
}

1. 运行DEM提取：执行DEM提取操作，这可能涉及到对图像进行旋转、缩放等预处理。

2. 输出结果：最后，你可以通过输出DEM文件或显示DEM图像来展示结果。

func runDEMExtraction() error {
    // 调用createDEM函数来创建DEM
    if err := createDEM(); err != nil {
        return err
    }

    // 现在我们可以运行DEM提取操作了
    // 这里是你的代码来运行DEM提取。

    // 输出结果
    fmt.Println("DEM extraction complete.")
    return nil
}

请注意，上述代码仅提供了一个基本框架。实际应用中，还需要考虑诸如图像预处理、误差处理以及可能的错误检查等方面。此外，代码中的具体实现细节可能会因版本差异而有所不同。因此，在实际使用前，请务必查阅并熟悉erdas I的官方文档或用户手册。