python写用tif立体像对提取DEM的代码

以下是一个简单的 Python 代码示例，用于从立体影像对中提取数字高程模型（DEM）：

import numpy as np
import rasterio as rio

def extract_dem(left_image_path, right_image_path):
    # Load the left and right images as numpy arrays
    with rio.open(left_image_path) as left:
        left_image = np.transpose(left.read(), (1, 2, 0))
    with rio.open(right_image_path) as right:
        right_image = np.transpose(right.read(), (1, 2, 0))

    # Compute the disparity map using a disparity algorithm
    disparity_map = compute_disparity(left_image, right_image)

    # Compute the DEM from the disparity map and camera parameters
    dem = compute_dem(disparity_map, camera_params)

    return dem

其中，`left_image_path` 和 `right_image_path` 是左右立体影像对的文件路径，`compute_disparity()` 是用于计算视差图的算法函数，`compute_dem()` 是用于从视差图和相机参数计算 DEM 的函数，`camera_params` 是相机的内部和外部参数。请注意，这只是一个简单的示例，您需要根据自己的数据和要求进行修改。