### 批量读取TXT文件并下载ICESat-2 ATL03数据 要实现通过Python批量读取TXT文件中的信息来下载ICESat-2的ATL03数据，可以按照以下方法操作： #### 1. 准备工作 首先需要安装必要的库。以下是主要依赖项： - `pandas`：用于处理CSV/TXT文件。 - `requests` 或 `urllib`：用于发送HTTP请求以下载数据。 - `h5py` 和 `numpy`：用于后续处理HDF5格式的数据。 可以通过以下命令安装所需库： ```bash pip install pandas requests h5py numpy ``` --- #### 2. 批量读取TXT文件 假设TXT文件存储的是每行一条记录，其中包含轨道编号或其他可用于查询ICESat-2数据的关键字段。下面是一段代码示例，展示如何批量读取TXT文件的内容： ```python import os import pandas as pd def read_txt_files(folder_path): """ 批量读取指定目录下的所有 TXT 文件，并返回一个 DataFrame。 :param folder_path: 存储 TXT 文件的文件夹路径 :return: 包含所有 TXT 文件内容的 Pandas DataFrame """ all_data = [] txt_files = [f for f in os.listdir(folder_path) if f.endswith('.txt')] for file_name in txt_files: file_path = os.path.join(folder_path, file_name) try: df = pd.read_csv(file_path, header=None, sep='\s+', comment='#') df['source_file'] = file_name # 添加来源文件名作为新列 all_data.append(df) except Exception as e: print(f"Error reading {file_name}: {str(e)}") combined_df = pd.concat(all_data, ignore_index=True) return combined_df ``` 此函数会遍历给定文件夹内的所有TXT文件，将其内容加载到Pandas DataFrame中，并附加源文件名称以便追踪[^4]。 --- #### 3. 下载ICESat-2 ATL03数据 NASA Earthdata 提供了一个API接口，允许开发者通过脚本访问和下载数据。为了使用该服务，需先注册Earthdata账号并获取认证令牌。 ##### (a) 获取认证令牌 登录[NASA Earthdata](https://urs.earthdata.nasa.gov/)网站后，创建账户并通过以下方式生成认证令牌： - 使用用户名密码进行基本身份验证； - 或者配置OAuth2授权流程。 ##### (b) 编写下载逻辑 以下代码展示了如何根据轨道号（RGT）和其他参数自动下载对应的ATL03 HDF5文件： ```python import requests from urllib.parse import quote_plus def download_atl03(rgt_list, output_folder, username, password): """ 根据 RGT 列表下载 ICESat-2 的 ATL03 数据。 :param rgt_list: 轨道编号列表 :param output_folder: 输出文件夹路径 :param username: NASA Earthdata 用户名 :param password: NASA Earthdata 密码 """ base_url = "https://n5eil01u.ecs.nsidc.org/egi/request" session = requests.Session() session.auth = (username, password) for rgt in rgt_list: query_params = { "short_name": "ATL03", "version": "006", "temporal": "", # 时间范围可自定义 "bbox": "-180,-90,180,90", # 地理边界框 "rgt": str(rgt), "processing_level": "L2", } encoded_query = "&".join([f"{k}={quote_plus(str(v))}" for k, v in query_params.items()]) request_url = f"{base_url}?{encoded_query}" response = session.get(request_url) if response.status_code != 200: print(f"Failed to fetch data for RGT {rgt}. Status code: {response.status_code}") continue links = [ link.split(",")[1].strip('"') for link in response.text.strip().split("\n") if link.startswith("https:") ] for url in links: filename = url.split("/")[-1] save_path = os.path.join(output_folder, filename) with open(save_path, 'wb') as f: res = session.get(url, stream=True) if res.status_code == 200: for chunk in res.iter_content(chunk_size=1024 * 1024): # 按 MB 分块下载 if chunk: f.write(chunk) print(f"Downloaded {filename}") else: print(f"Failed to download {url}, status code: {res.status_code}") ``` 在此代码片段中，`download_atl03` 函数接受一组轨道编号以及目标输出文件夹位置，依次尝试从NSIDC服务器拉取匹配的ATL03数据文件[^5]。 --- #### 4. 结合两部分功能 最后一步是将上述两个模块结合起来，形成完整的解决方案： ```python if __name__ == "__main__": input_folder = "/path/to/txt/files" # 替换为实际存放TXT文件的位置 output_folder = "/path/to/output/folder" # 替换为目标存储空间 earthdata_username = "your_earthdata_username" earthdata_password = "your_earthdata_password" # Step 1: 读取TXT文件 track_info = read_txt_files(input_folder) unique_rghts = track_info.iloc[:, 0].unique() # 假设第一列为RGTS # Step 2: 开始下载 download_atl03(unique_rghts.tolist(), output_folder, earthdata_username, earthdata_password) ``` 以上程序实现了从本地TXT文件解析需求至远程ICESat-2数据集检索的过程[^6]。 --- ###