【干货】免费12.5米精度ALOS DEM:从批量下载到镶嵌重投影全流程自动化

pip install requests rasterio numpy

import osimport requests# 1. 你的下载链接列表 (从第一阶段获取的链接放这里)download_urls = [    "https://datapool.asf.alaska.edu/RTC/ALOS/AP_14756_FBS_F0500_RT1.zip",    "https://datapool.asf.alaska.edu/RTC/ALOS/AP_14121_FBS_F3090_RT1.zip",    # ... 粘贴你导出的所有链接]# 2. 配置 Cookie (在ASF官网按F12 -> Network -> 刷新 -> 找个请求看Request Headers)# 或者直接使用 Earthdata 的账号密码登录接口（此处演示Cookie法，最简单）headers = {    'Cookie': '你的_urs_user_login_session_cookie',    'User-Agent': 'Mozilla/5.0 ...'}save_path = r"D:\download"def download_file(url):    local_filename = url.split('/')[-1]    path = os.path.join(save_path, local_filename)    # 流式下载，内存不爆炸    with requests.get(url, headers=headers, stream=True) as r:        r.raise_for_status()        with open(path, 'wb') as f:            for chunk in r.iter_content(chunk_size=8192):                f.write(chunk)    print(f"[√] {local_filename} 下载完毕")for url in download_urls:    download_file(url)

import osimport zipfileTARGET_DIR = r"E:\库区消落带\ALOS_DEM"ZIP_FILES = [r"D:\download\AP_14756_FBS_F0500_RT1.zip", ...] # 你的压缩包路径def full_unzip():    if not os.path.exists(TARGET_DIR):        os.makedirs(TARGET_DIR)    for zip_path in ZIP_FILES:        with zipfile.ZipFile(zip_path, 'r') as zf:            # 这里我们暴力全解压，后续脚本会自动识别DEM文件            zf.extractall(TARGET_DIR)        print(f"[√] 解压完成: {os.path.basename(zip_path)}")if __name__ == "__main__":    full_unzip()

import osimport globimport rasteriofrom rasterio.merge import mergefrom rasterio.vrt import WarpedVRTfrom rasterio.enums import Resampling# 输入文件夹 (解压后的目录)INPUT_FOLDERS = [    r"E:\库区消落带\ALOS_DEM\AP_14756_FBS_F0500_RT1",    # ... 你的所有文件夹路径]OUTPUT_DIR = r"E:\库区消落带\ALOS_DEM"OUTPUT_FILENAME = "ALOS_Mosaic_WGS84_All.tif"def mosaic_final():    print("--- 正在搜寻 DEM 文件 ---")    tif_files = []    for folder in INPUT_FOLDERS:        # 自动寻找 .dem.tif 文件        found = glob.glob(os.path.join(folder, "*.dem.tif"))        if found: tif_files.append(found[0])    print(f"--- 开始处理 {len(tif_files)} 个文件 ---")    src_files = []    try:        for fp in tif_files:            src = rasterio.open(fp)            # 【核心】建立虚拟投影，强制转为 WGS84            vrt = WarpedVRT(src, crs='EPSG:4326', resampling=Resampling.bilinear)            src_files.append(vrt)        # 执行内存镶嵌        mosaic, out_trans = merge(src_files)        # 更新元数据并输出        out_meta = src_files[0].meta.copy()        out_meta.update({            "driver": "GTiff",            "height": mosaic.shape[1],            "width": mosaic.shape[2],            "transform": out_trans,            "crs": "EPSG:4326",            "compress": "lzw"        })        with rasterio.open(os.path.join(OUTPUT_DIR, OUTPUT_FILENAME), "w", **out_meta) as dest:            dest.write(mosaic)        print("[成功] 完美拼接！")    finally:        for src in src_files: src.close()if __name__ == "__main__":    mosaic_final()