持续更新。。。。。。。。。。。。。。。
【第三版】
MOD021KM
是 MODIS L1B 级数据产品的核心标识符,代表 Terra 卫星获取的 1 公里分辨率的校准辐射数据。以下是关于该数据的详细解析与应用指南:
一、MOD021KM 数据详解
属性 | 说明 |
---|---|
卫星平台 | Terra (上午星) |
产品级别 | L1B (Level 1B) |
分辨率 | 1 公里 (实际包含 250m/500m/1km 波段,以 1km 网格封装) |
核心内容 | 36 个波段的辐射亮度值 (Calibrated Radiance) 和 地理定位信息 |
文件格式 | HDF-EOS (.hdf) |
波段范围 | 可见光(1-7) ~ 热红外(20-36) |
典型用途 | 真彩色合成、云检测、火灾点提取、大气校正输入数据 |
二、关键数据获取方式
1. 官方下载源
LAADS DAAC (推荐)
https://ladsweb.modaps.eosdis.nasa.gov
路径:/allData/61/MOD021KM/
(61 为 Collection 6.1 版本)
优势:下载速度快,按轨道号组织Earthdata Search
https://search.earthdata.nasa.gov
优势:可通过地图框选区域搜索
2. 文件名解析示例
MOD021KM.A2024185.0230.061.2024185234323.hdf
A2024185
:2024年第185天获取0230
:UTC时间02:30061
:Collection 6.12024185234323
:处理时间戳
三、数据处理流程(Python示例)
步骤 1:读取辐射亮度值
import pyhdf.SD as hdf
# 打开HDF文件
file = hdf.SD('MOD021KM.A2024185.0230.061.hdf')
datasets = file.datasets()
# 获取500m分辨率波段8(近红外)数据
band_index = 'EV_500_Aggr1km_RefSB' # 500m波段组名称
data_band8 = file.select(band_index).get()[:,:,7] # 第8波段在数组中的索引
步骤 2:辐射定标 → 表观反射率
# 从属性获取缩放因子
attrs = file.select(band_index).attributes()
scale = attrs['reflectance_scales'][7]
offset = attrs['reflectance_offsets'][7]
# 计算表观反射率
reflectance_band8 = (data_band8 - offset) * scale
步骤 3:生成真彩色影像 (RGB:波段1-4-3)
import numpy as np
import matplotlib.pyplot as plt
# 读取波段1(红),4(绿),3(蓝)
red = file.select('EV_250_Aggr1km_RefSB').get()[:,:,0] * 0.0001
green = file.select('EV_500_Aggr1km_RefSB').get()[:,:,3] * 0.0001
blue = file.select('EV_500_Aggr1km_RefSB').get()[:,:,2] * 0.0001
# 合成并增强
rgb = np.dstack([red, green, blue])
rgb_normalized = np.clip(rgb * 3.5, 0, 1) # 对比度增强
# 显示
plt.imshow(rgb_normalized)
plt.savefig('MODIS_TrueColor.jpg', dpi=300)
四、专业工具推荐
MRT (MODIS Reprojection Tool)
- 官方工具:地理校正/格式转换
- 下载:https://lpdaac.usgs.gov/tools/modis_reprojection_tool
HEG (HDF-EOS to GeoTIFF)
- 批量转换HDF为GeoTIFF
- 下载:https://earthdata.nasa.gov/esdis/eso/standards-and-references/heg
QGIS + HDF插件
- 可视化查看波段结构
五、典型应用场景
火点检测
使用 波段21+22(4μm) 热辐射值:
if rad_band21 > 330 K and rad_band22/rad_21 > 0.8: 判定为火点
云掩膜生成
结合波段26(1.38μm)和31(11μm)区分云/晴空气溶胶反演
深蓝波段(3)与红波段(1)组合计算光学厚度
六、注意事项
- Bowtie效应
近边缘数据需用MRT
进行去重叠处理 - 数据量预警
单文件约 ~150 MB,全球每日数据量超1TB - 定标版本
Collection 6.1 比旧版改进了海洋耀斑校正
替代方案:若只需快速可视化,直接使用 NASA Worldview 在线平台可免去本地处理。
通过正确处理 MOD021KM
,您将获得科研级遥感分析的基石数据!
内容 | 地址 链接 |
---|---|
总览 | 【第三版】系统集成项目管理工程 |
十五至尊图 | 第三版 |
版本记录:
- 2025年6月27第一版