CMEMS全球海洋产品数据–超级详细下载教程【python】

copernicusmarine login# 按提示输入用户名和密码

# 搜索海面高度相关数据集copernicusmarine describe --contains "sea level" | grep dataset_id

copernicusmarine subset \  --dataset-id cmems_obs-sl_glo_phy-ssh_nrt_allsat-l4-duacs-0.125deg_P1D \  --variable sla \  --start-datetime 2026-04-15 \  --end-datetime 2026-04-15 \  --minimum-longitude 110 \  --maximum-longitude 120 \  --minimum-latitude 10 \  --maximum-latitude 20 \  --output-directory ./test_data

#!/usr/bin/env python3"""CMEMS 数据下载测试脚本（自动登录版）内置账号密码，一键运行"""import osimport sysimport subprocessfrom datetime import datetime, timedelta# ==================== 账号配置（直接写入）====================USERNAME = "183@qq.com"PASSWORD = "."   # 密码# ==================== 下载参数 ====================DATASET_ID = "cmems_obs-sl_glo_phy-ssh_nrt_allsat-l4-duacs-0.125deg_P1D"# 下载日期（取3天前，确保数据已发布）TARGET_DATE = (datetime.utcnow() - timedelta(days=3)).strftime("%Y-%m-%d")# 若上述日期无数据，可手动指定：TARGET_DATE = "2026-04-15"LON_MIN, LON_MAX = 110, 120LAT_MIN, LAT_MAX = 10, 20VARIABLES = ["sla"]OUTPUT_DIR = "./cmems_test_output"OUTPUT_FILENAME = f"sla_test_{TARGET_DATE}.nc"# ==================== 函数定义 ====================def auto_login():    """自动登录 CMEMS"""    print("🔐 正在自动登录 CMEMS...")    login_cmd = [        "copernicusmarine", "login",        "--username", USERNAME,        "--password", PASSWORD,        "--force-overwrite"    ]    try:        subprocess.run(login_cmd, check=True, capture_output=True, text=True)        print("✅ 登录成功")    except subprocess.CalledProcessError as e:        print(f"❌ 登录失败: {e.stderr}")        sys.exit(1)def download_test_data():    """执行数据下载"""    print(f"\n📥 开始下载测试数据...")    print(f"   数据集: {DATASET_ID}")    print(f"   日期: {TARGET_DATE}")    print(f"   区域: 经度 {LON_MIN}~{LON_MAX}, 纬度 {LAT_MIN}~{LAT_MAX}")    print(f"   变量: {VARIABLES}")    print(f"   输出目录: {OUTPUT_DIR}")    os.makedirs(OUTPUT_DIR, exist_ok=True)    try:        from copernicusmarine import subset        subset(            dataset_id=DATASET_ID,            start_datetime=TARGET_DATE,            end_datetime=TARGET_DATE,            minimum_longitude=LON_MIN,            maximum_longitude=LON_MAX,            minimum_latitude=LAT_MIN,            maximum_latitude=LAT_MAX,            variables=VARIABLES,            output_directory=OUTPUT_DIR,            output_filename=OUTPUT_FILENAME,        )        return True    except Exception as e:        print(f"\n❌ 下载失败: {e}")        return Falsedef verify_file():    """验证下载的文件"""    file_path = os.path.join(OUTPUT_DIR, OUTPUT_FILENAME)    if not os.path.exists(file_path):        print("❌ 文件未生成")        return False    file_size = os.path.getsize(file_path) / 1024    print(f"\n✅ 文件已保存: {file_path}")    print(f"   文件大小: {file_size:.2f} KB")    try:        import xarray as xr        ds = xr.open_dataset(file_path)        print(f"   变量: {list(ds.data_vars.keys())}")        print(f"   时间维度: {ds.time.values}")        ds.close()        print("✅ 文件格式验证通过")    except ImportError:        print("⚠️ 未安装 xarray，跳过格式验证")    except Exception as e:        print(f"⚠️ 读取文件时出错: {e}")    return True# ==================== 主程序 ====================def main():    print("="*60)    print("   CMEMS 数据下载测试脚本（自动登录版）")    print("="*60)    auto_login()    if download_test_data():        verify_file()        print("\n🎉 测试完成！环境配置正确，可以正常下载数据。")    else:        print("\n⚠️ 测试失败，请检查网络或数据集ID")if __name__ == "__main__":    main()

python download_cmems_test.py