从GEE下载sentinel-1土壤湿度双极化数据python代码

在Google Earth Engine (GEE) 中，你可以使用Python API来下载Sentinel-1 SAR（ Synthetic Aperture Radar）数据，包括土壤湿度信息，特别是对于双极化数据。以下是一个简单的示例代码片段，展示了如何设置查询、过滤条件和下载数据：

# 导入必要的库
import ee
from google.colab import drive
import os

# 登录到GEE
ee.Authenticate()

# 链接到Google Drive以便存储数据
drive.mount('/content/drive')

# 设置工作区
os.environ['GEE_DATA'] = '/content/drive/MyDrive/gee_data'

# 初始化GEE
ee.Initialize()

# 查询Sentinel-1数据
def get_sentinel1_data(start_date, end_date, area):
    # 创建 Sentinel-1 模型
    s1 = ee.ImageCollection('COPERNICUS/S1_GRD')
    
    # 过滤只包含VV和 VH极化的数据
    polarizations = ['VV', 'VH']
    s1 = s1.filter(ee.Filter.listContains('transmitterReceiverPolarization', polarizations))
    
    # 设置时间范围
    s1 = s1.filterDate(start_date, end_date)
    
    # 定义兴趣区域
    geometry = ee.Geometry.Rectangle(area)
    
    return s1.filterBounds(geometry).sort('CLOUDY_PIXEL_PERCENTAGE').limit(1)

# 使用函数并替换实际日期和地理位置
start_date = '2022-01-01'
end_date = '2022-01-31'
area = ee.Geometry.Polygon([[-122.5, 47], [-121.5, 47], [-121.5, 48], [-122.5, 48]])

# 下载数据
download_url = get_sentinel1_data(start_date, end_date, area).first().getDownloadURL({
    'scale': 10,  # 数据分辨率
    'crs': 'EPSG:4326',  # 空间参考系统
    'filePerBand': False  # 合并所有波段到一个文件
})

# 执行下载
drive.mount("/content/drive", force_remount=True)  # 可能需要再次挂载
data_path = f'/content/drive/MyDrive/{download_url.split("/")[-1]}'
ee.data.downloadToDrive(download_url=download_url, description=data_path)