python镶嵌dem

根据提供的引用内容，以下是一个用Python进行DEM数据批量镶嵌的示例：

import os
import glob
from osgeo import gdal

# 设置输入DEM数据文件夹路径和输出镶嵌后的数据路径
input_folder = 'path/to/input/folder'
output_file = 'path/to/output/file.tif'

# 获取输入文件夹中的所有DEM数据文件
dem_files = glob.glob(os.path.join(input_folder, '*.tif'))

# 读取第一幅DEM数据，获取其空间参考和像素大小等信息
first_dem = gdal.Open(dem_files[0])
projection = first_dem.GetProjection()
geotransform = first_dem.GetGeoTransform()
pixel_size = geotransform[1]

# 计算输出镶嵌后的图像的大小
x_size = int((first_dem.RasterXSize * pixel_size) / abs(geotransform[1]))
y_size = int((first_dem.RasterYSize * pixel_size) / abs(geotransform[5]))

# 创建输出镶嵌后的图像
driver = gdal.GetDriverByName('GTiff')
output_dem = driver.Create(output_file, x_size, y_size, 1, gdal.GDT_Float32)
output_dem.SetProjection(projection)
output_dem.SetGeoTransform(geotransform)

# 遍历所有DEM数据文件，将其读取并写入输出镶嵌后的图像
for dem_file in dem_files:
    dem = gdal.Open(dem_file)
    dem_data = dem.GetRasterBand(1).ReadAsArray()
    output_dem.GetRasterBand(1).WriteArray(dem_data, 0, 0)

# 关闭文件
output_dem = None

print("DEM数据镶嵌完成！")

请注意，这只是一个示例脚本，具体的DEM数据镶嵌过程可能因数据格式和需求而有所不同。你需要根据实际情况进行适当的修改和调整。