气象卫星这个门类很多人第一反应是天气预报里的那张云图实际上它早就演化成一套分工明确的观测体系有的卫星守在固定位置盯着同一片区域每几分钟拍一次有的卫星沿着南北向轨道绕地球转、每天固定时刻从同一地点上方掠过。搞清这套体系的差别直接决定了你做遥感反演、灾害监测、农业估产还是资料同化时该选哪颗星、哪个通道、哪一级产品。下面这份整理是我这几年在数据处理和项目选型里反复用到的东西从轨道原理、载荷设计、国内外主要系列一直到数据落地和踩过的坑能直接拿去用。1. 气象卫星的轨道分类与体系设计逻辑1.1 静止轨道与极轨的物理约束差异气象卫星的轨道选择不是工程师拍脑袋定的而是被两个物理量死死卡住观测范围和时间分辨率。这一对指标天然互斥所以整个体系不得不拆成两条线。静止轨道的核心约束是轨道周期必须等于地球自转周期。按开普勒第三定律反推T 取 86164 秒恒星日半长轴 a (μT²/4π²)^(1/3)代入 μ 398600 km³/s²得到 a ≈ 42164 km减去地球半径 6371 km轨道高度约 35786 km。这个高度是算出来的不是选出来的——任何偏离都会让卫星相对地面漂移失去静止的意义。好处是一颗星能覆盖约三分之一地球表面全圆盘成像周期可以压到 10 分钟甚至更短代价是距离太远星下点像元几何尺寸做不细1 微弧度在 35786 km 处对应约 36 米一个 15 微弧度的探测器也只能做到 500 米级别红外通道普遍在 2 到 4 公里。极轨走的是另一条路。轨道高度控制在 800 到 900 公里配合约 98 度的倾角做成太阳同步轨道保证卫星每次经过同一纬度时的地方时基本一致——这一点极其关键否则不同日期的影像光照条件不同做时间序列对比就没法看了。按 836 公里高度算轨道周期a 6371 836 7207 kmT 2π√(a³/μ) ≈ 6090 秒也就是约 101.5 分钟一天绕地球约 14.2 圈。这个数字直接决定了单星的重访能力单颗极轨星对中低纬某地一天最多覆盖两次升轨、降轨各一次要缩短重访间隔只能靠多星组网或者扩大幅宽。我给新手解释这两者的差别时常用一个类比静止轨道卫星像小区门口固定机位的摄像头一直盯着一个方向能看到事情的完整过程但看不清细节极轨卫星像巡逻车能凑近看清楚但同一地点一天只能路过有限几次。做短时强对流、台风路径这种过程类监测必须靠静止星做地表分类、植被指数、精细地表温度只能靠极轨星。1.2 轨道参数怎么影响你的实际使用轨道参数不是仅供科普的数字它直接决定了你能拿到什么样的数据用的时候必须先想清楚三件事。地方时。太阳同步轨道有个降交点地方时的概念指卫星由北向南经过赤道时的当地时刻。这个值决定了观测的光照几何。下午星地方时 13:30 到 14:00 左右在午后对流发展旺盛时段过境对强对流、热力特征敏感上午星地方时 10:00 左右光照角度低、阴影长对地表纹理和云的立体结构更友好。晨昏轨道卫星则专门卡在日出日落附近过境太阳高度角接近 0可见光通道基本失效但对温度日变化的捕捉和夜间低云、雾的识别特别有价值还能给数值模式提供辐射收支的连续观测。选数据之前先查清降交点地方时比盲目试算有效得多。幅宽与重叠。极轨星单条扫描带的幅宽从 2000 多公里到 2900 公里不等相邻轨道在高纬度会大幅重叠在低纬度可能留缝。做区域拼接时如果不检查相邻轨道的重叠率很容易在低纬出现条带缺口。经验上纬度超过 50 度后重叠会非常充分30 度以下就要仔细核对轨道间隔。定标与几何精度。静止星长时间盯一个点几何定位相对稳定但会有卫星平台微振动、姿态漂移带来的像元抖动极轨星每轨都要重新算星历和姿态定位误差在山区和海岸线附近更容易被看出来。同一颗星的 L1 产品通常已经做过几何精校正但如果你要拿它和其他传感器做像素级配准还是得自己做一次控制点匹配。提示轨道类型决定了时间分辨率和空间分辨率这对矛盾如何分配做方案设计时先明确业务需要的是看得勤还是看得清不要在选型阶段就试图两者兼得。2. 载荷体系拆解一颗气象卫星到底在看什么2.1 成像通道的波长分配逻辑多通道成像仪是气象卫星绝对的主力载荷通道数量从早期的个位数发展到现在的十几到二十几个每一个通道都不是凑数的背后对应一类明确的地物或大气目标。可见光到近红外波段主要解决云和地表长什么样。0.47 微米附近的蓝光通道用于真彩色合成和海洋水色观测同时能反映气溶胶0.51 微米绿光配合红光做真彩色0.64 微米红光是云检测和植被指数的核心通道云在这个波段反射率高和海表、植被对比明显0.86 微米近红外利用植被叶片细胞结构的强反射特性是植被长势监测的经典通道同时用于水陆边界识别1.6 微米短波红外对云相态特别敏感冰云在这个波段吸收强、反射低水云则反射高是区分冰云和水云的关键2.1 微米通道则对云粒子有效半径敏感。红外波段解决温度多少和水汽在哪。3.7 到 3.9 微米中波红外对高温目标极其敏感火点、工业热源、火山口在这个波段亮温能比周围高出几十开尔文是火情监测的首选但同时它白天受反射太阳光污染只能夜间或做差分使用。6.2 到 7.3 微米是水汽吸收带探测的是对流层中上层的水汽分布图像上能直接看出急流位置、干湿区边界和卷云羽8.6 微米对薄卷云和火山灰敏感10.8 微米是窗区通道云顶温度、地表温度反演的主力12.0 微米与 10.8 微米做差分可以判断云粒子大小和薄卷云13.3 微米二氧化碳吸收带用于云顶高度估计和大气温度廓线反演。这些通道怎么组合是关键。雾的识别用 0.86 微米反射率和 10.8 微米亮温的差值因为雾滴在近红外反射率低、温度又接近地表夜间低云检测靠 3.9 微米和 10.8 微米的亮温差低云在这个差分上通常呈现负值或接近零而中高云接近正沙尘和云的区别要看 10.8 微米与 12.0 微米差分的符号和幅度沙尘的差分往往呈现与云相反的特征。这些通道组合知识比背参数重要得多实际做产品时靠的就是这些差分组合。2.2 垂直探测与高光谱的介入成像仪看的是平面垂直探测仪看的是剖面。早期用十几个到二十个红外通道做大气温度、湿度廓线反演垂直分辨率只有几公里误差也大。后来高光谱技术进入气象卫星一次观测能拿到上千个通道垂直分辨率提升到 1 公里左右这就是干涉式或者光栅式的高光谱红外探测仪。这类载荷的价值在于给数值天气预报提供初始场。模式对初值极其敏感尤其是湿度场和温度场的垂直结构如果只有成像仪资料模式只能拿到云顶温度和整层水汽总量垂直方向的信息严重不足。高光谱探测仪能以很高的光谱分辨率切开二氧化碳吸收带和水的吸收带反演出分层温湿度这对 3 到 5 天预报的改进是实打实的。静止轨道上的高光谱探测是近年的一大进步。因为卫星能长时间停留可以做到每小时甚至更快一次的温湿度廓线等于在固定位置放了一台大气 CT 机。这类产品对强对流触发前的环境场诊断特别有用——你能看到对流发生前几个小时低层暖湿、中层干冷这种不稳定层结是如何一步步建立的。注意高光谱产品数据量大、通道多直接全通道使用非常吃算力。实际业务里通常先做通道筛选把对目标敏感、彼此相关性低的通道挑出来再做反演或者同化通道数能压到原来的十分之一。2.3 微波、散射计与主动雷达的补位可见光和红外有个致命短板穿不透云。台风最强的对流云区、暴雨的降水核心在红外图像上只是一片白看不到内部结构。微波载荷就是来解决这个的。微波辐射计工作在 10 到 183 GHz 频段能穿透非降水云和薄云。低频通道10 到 37 GHz对地表和降水敏感高频通道89 GHz 以上对云中冰相粒子敏感。降水的散射信号在 89 GHz 上表现为亮温下降海洋上的暖云降水则表现为低频亮温升高两个方向结合就能估算降水强度。微波的温度权重函数分布在不同高度多个通道联立还能反演大气温度和水汽廓线而且不受云污染影响。散射计是另一个思路主动发射 C 波段或 Ku 波段脉冲测量海面粗糙度进而反演海面风场。海面风速越大毛细波越粗糙后向散射截面越大。这个原理让散射计成为海上大风、台风外围风场监测的核心手段缺点是空间分辨率较粗十几到几十公里近岸和冰区容易受污染。星载主动雷达目前数量不多但地位特殊。双频降水测量雷达能直接测得降水的三维反射率结构配合辐射计可以显著提升降水反演精度尤其是弱降水和固态降水的估计。还有一些小卫星搭载了云廓线雷达专门看云的垂直层次。载荷类型主要观测目标典型频段/通道分辨率量级典型局限多通道成像仪云、地表、气溶胶0.47-13.3 微米十余通道0.5-4 公里受云遮挡红外仅反映云顶高光谱红外探测仪温湿度垂直廓线3.7-15 微米千级通道十几公里足迹有云时反演困难微波辐射计降水、温湿廓线、海冰10-183 GHz5-50 公里陆地上地表发射率复杂散射计海面风场C/Ku 波段12-25 公里近岸、冰区污染星载雷达降水三维结构Ku/Ka 双频5 公里级幅宽窄覆盖有限掩星接收机大气折射率、温湿L1/L2 导航信号数百米垂直分辨率水平分辨率低事件稀疏3. 国内业务星座风云系列的定位与分工3.1 极轨系列的技术演进与选型要点国内的极轨气象卫星从早期的单通道扫描仪起步一路迭代到现在的多载荷综合观测平台。这个系列的演进逻辑很清楚从能看见到看得细再到看得准。第一代解决的是有没有的问题搭载的扫描辐射计通道很少分辨率在公里级主要提供云图和粗略的地表信息。第二代开始引入十通道量级的扫描辐射计增加了短波红外和水汽通道分组方式也做了优化可以为区域应用提供更灵活的数据下传。到第三代载荷数量明显增加形成了一套完整的组合中分辨率光谱成像仪达到二十多个通道、星下点 250 米分辨率、幅宽接近 3000 公里是主力产品高光谱红外探测仪提供上千通道的垂直探测微波温度计、微波湿度计、微波成像仪负责穿云观测地球辐射收支仪负责辐射平衡臭氧探测仪负责大气化学成分还有空间环境监测和掩星探测载荷。第三代还引入了晨昏轨道星专门填补清晨和傍晚的观测空档。用第三代极轨数据时我比较关注的几个点幅宽覆盖单轨幅宽接近 3000 公里中纬度相邻轨道重叠充分但低纬度拼图仍要核对轨道间隔尤其是想拼半个月的合成产品时。通道一致性同一载荷的可见光通道和红外通道分辨率不同250 米与 1 公里做彩色合成时要先重采样别直接把 1 公里通道插值到 250 米就当成高分辨率用那样只会让图像更糊。定标状态极轨星有太阳定标和星上黑体定标但可见光通道会随时间的衰减用长序列做趋势分析前先查定标说明必要时做交叉定标。晨昏星的用法晨昏轨道光照角度极低可见光基本没法用但红外通道在日出前后的温度对比对识别低云、雾和逆温层特别有效做晨雾预报时非常值得调出来看。3.2 静止系列的扫描能力与快速成像静止系列的第一代解决了在固定位置连续看的问题五通道配置一个可见光加四个红外在那个年代已经能支撑逐小时云图。第二代提升了通道数和定量化水平主要指标在可见光 1.25 公里、红外 5 公里左右扫描方式采用线扫描靠卫星自旋配合扫描镜完成二维成像。这套机制很经典但也有物理限制自旋稳定的卫星转速固定扫描节奏受限一次全圆盘成像要半小时左右。第三代换成了三轴稳定平台扫描镜可以自由摆动成像周期和灵活性大幅提升。核心成像仪通道数达到十几个可见光分辨率做到 1 公里以内红外 4 公里左右全圆盘成像周期缩短到 15 分钟上下区域扫描可以做到 5 分钟一次。同时搭载了干涉式大气垂直探测仪提供逐小时的温湿度廓线闪电成像仪可以在白天黑夜连续监测闪电事件对强对流初生和台风内部的对流活动判断很有帮助还有空间环境监测载荷。后续发射的改进型进一步增加了快速成像仪区域扫描可以做到 1 分钟、250 米分辨率这对中小尺度对流系统的追踪是质变——过去 15 分钟一次只能看到结果1 分钟一次能看到过程对流云的快速膨胀、云顶抬升速率都能算出来这些参数是判断是否即将出现强天气的关键。静止星有个容易忽视的技术点扫描与积分时间的权衡。想让图像信噪比高就要延长每个像元的积分时间但积分时间长了成像周期就慢而且地球自转和卫星微振动会带来拖影。所以快速扫描模式下的图像噪声通常比常规模式大做定量反演时要留意产品说明里的噪声指标别拿快速扫描的数据去做微弱信号检测。实操心得用静止星做短时预报时别只看单帧云图把连续多帧做成动画再看云顶亮温的降温速率、云团面积的扩张速度这两个指标对判断对流是否要爆发比单帧形态更可靠。3.3 降水测量与专用载荷的补充降水观测一直是个难点。地面雨量站分布不均海上基本空白雷达覆盖有限卫星就成了最有效的补位手段。国内在这块做了一条专门的降水测量卫星搭载双频降水测量雷达和微波成像仪雷达工作在 Ku 和 Ka 两个频段。双频搭配的道理在于Ku 频段对中等强度降水敏感、衰减小Ka 频段对弱降水敏感但强降水会饱和两者结合能扩大降水反演的量程同时能反演雨滴谱信息区分小雨和暴雨的微物理特征。这个频段选择和国际上主流的双频降水雷达思路一致目的是形成互补的全球降水观测网。除了降水雷达还有专门的低倾角轨道设计来提升中低纬降水观测的采样频率。低倾角轨道的思路是让卫星更频繁地扫过热带和副热带地区那里是全球降水最集中的区域但常规太阳同步轨道的覆盖频次不够。用降水雷达产品时要注意几个现实问题雷达幅宽通常只有几百公里一天覆盖范围有限单靠它做区域降水统计样本量不足得配合微波成像仪做宽幅补充雷达反射率到降水强度的转换依赖滴谱假设不同气候区、不同降水类型的偏差可以很大做定量应用前最好用地面雷达或者雨量站做一次本地化校验。4. 国外主要气象卫星体系横向对比4.1 北美业务体系从极轨到静止的完整链条北美地区的极轨体系经历了几代演进。早期业务星搭载的扫描辐射计是十通道量级、公里级分辨率服役时间长为长时间序列研究提供了宝贵资料。后续的联合极轨系统做了大幅升级主力载荷是一台二十多通道的成像仪星下点分辨率做到 375 米、幅宽超过 3000 公里配套的高光谱红外探测仪通道数在千级垂直探测精度高微波探测仪二十多通道覆盖 23 到 183 GHz温湿度廓线反演的核心还有掩星接收机和辐射收支载荷。这个体系的特点是成像加探测加微波三件套齐全且长期连续对做气候趋势分析很友好。有一点值得注意这套体系的成像仪有一个跨轨扫描加像元聚合的设计部分通道在沿轨和跨轨方向的分辨率不同直接投影时会出现像元拉长。处理的时候要么按最粗方向重采样要么按扫描几何自己重算坐标不能简单当成规则格点用。静止体系在近十年完成了换代新一代静止星的成像仪有 16 个通道可见光 0.5 公里、红外 2 公里全圆盘 10 分钟一次还支持 5 分钟一次的美国本土区域扫描和 1 分钟级别的中尺度扇区扫描。这个灵活扫描能力是很大的进步平时做全圆盘监测台风逼近时切换到中尺度扇区1 分钟一次盯着风暴中心能清楚看到眼壁替换、云顶暖心结构的演变。体系上采用东西两颗星分别负责不同经度加上在轨备份保证业务不中断。真彩色合成的体验值得一提16 通道配置里有 0.47 微米的蓝光通道这是做真彩色图像的必要条件。没有蓝光通道的卫星只能做伪彩色植被和云的区分度明显差一截。4.2 欧洲体系静止与极轨的双线布局欧洲的静止系列从早期的三通道载荷起步逐步升级到十二通道、3 公里分辨率的水平全圆盘 15 分钟一次还能做 5 分钟一次的欧洲区域快速扫描。第三代静止星改用了全新平台成像仪通道数增加到 16 个分辨率提升到 0.5 到 2 公里全圆盘 10 分钟一次、欧洲区域 2.5 分钟一次同时携带了红外高光谱探测仪和紫外可见光光谱仪后者对空气质量、臭氧、二氧化氮等痕量气体的监测是强项。这套体系的一个特色是成像加探测分星部署成像星和探测星搭配工作避免载荷之间争抢平台资源。极轨方面欧洲的极轨系列以载荷多样著称。成像仪是六通道、1 公里分辨率虽然通道数不多但定标质量好、序列长高光谱红外探测仪通道数超过 8000光谱分辨率极高是全球数值同化的主力资料之一微波湿度探测仪负责湿度廓线散射计提供海洋风场还有大气成分监测载荷。极轨轨道采用上午星设计地方时在 9:30 左右与下午星形成互补。用欧洲极轨数据做长时间序列时要特别注意跨代交叉定标。不同代际的同类载荷虽然设计指标接近但滤光片响应函数有细微差异直接拼接会出现台阶。正规做法是选取重叠运行期做交叉定标拟合出转换系数再拼接序列。4.3 东亚、南亚与极区体系的补充价值东亚地区的静止系列在国际上使用率很高其中一颗的核心指标是 16 通道、可见光 0.5 公里、红外 2 公里全圆盘 10 分钟一次还能做日本周边 2.5 分钟一次的快速扫描。另一个很强的点是真彩色图像的渲染——官方提供的真彩色合成产品做了瑞利散射校正海面颜色和云层的层次感做得非常自然做展示和科普极好。数据开放程度也高很多做云检测、气溶胶反演的项目组都拿它做过算法验证。南亚的静止系列有六通道成像仪和十九通道垂直探测仪的组合全圆盘半小时左右一次主要服务于季风区监测和热带气旋预警。韩国的静止系列有 16 通道成像仪和海洋水色成像仪兼顾陆海观测。俄罗斯的静止星通道数在十个左右重点覆盖高纬度地区极轨星搭载低分辨率多通道扫描仪和微波辐射计对极区海冰、积雪监测有独特价值。极区观测是个容易被忽略但很实用的方向。极轨卫星在高纬度轨道重叠密、观测频次高加上极区专用载荷对海冰边缘、冰面温度、积雪深度的敏感性对航运、气候研究都有直接用途。高纬度地区还有个额外好处极轨卫星在那里看得次数多一天能覆盖七八次甚至更多做逐小时序列完全够用。体系代表系列轨道类型成像通道数星下点最高分辨率全圆盘周期国内极轨风云第三代太阳同步20 以上250 米不适用国内静止风云第四代地球静止14-150.5-1 公里15 分钟北美极轨联合极轨系列太阳同步22成像仪375 米不适用北美静止GOES-R 系列地球静止160.5 公里10 分钟欧洲静止MTG 系列地球静止160.5 公里10 分钟欧洲极轨Metop 系列太阳同步6主成像仪1 公里不适用东亚静止Himawari 系列地球静止160.5 公里10 分钟南亚静止INSAT 系列地球静止61 公里30 分钟5. 数据落地实操从产品分级到读取处理5.1 产品分级体系与实际含义拿到气象卫星数据第一件事是搞清手里的文件是哪一级。这个分级在各家体系里命名一致但内容细节有差别。L0 是原始码流包含完整的仪器数据包和工程参数一般用户用不到除非要做定标重处理。L1 是加了辐射定标和几何定位的数据分 L1A未做几何校正和 L1B已做辐射和几何处理。真正拿来用的大部分是 L1 和 L2。L2 是反演产品云量、云顶温度、地表温度、海表温度、气溶胶光学厚度、降水率、大气廓线这些都在这一级每个产品文件里通常还带质量标记和置信度。L3 是投影到规则格网上的拼接产品分日、旬、月合成适合做区域分析和趋势研究。L4 是再分析或融合产品多源资料同化后的结果格点规整、时间连续但已经不是原始观测。选级的原则很简单做算法研究用 L1做应用分析用 L2做气候统计用 L3/L4。跳过 L1 直接拿 L3 去做反演改进基本没意义因为中间的重采样已经把原始信息平均掉了。反过来只想看某地一周的植被变化用 L3 的合成产品比自己去拼 L1 快得多。L2 产品有个坑不同产品的空间分辨率不同。云掩膜可能是 1 公里地表温度可能是 5 公里气溶胶可能是 10 公里把它们叠在一起做分析前一定要统一到同一格网而且要用最邻近或者面积加权不能用双线性插值把低分辨率产品糊到高分辨率上再当高分辨率用。5.2 分辨率、投影与重采样的取舍投影选择要看区域大小和纬度。赤道附近用等经纬度投影误差小中纬度大范围区域用兰勃特或者阿尔伯斯等积投影更合适极区必须用极射赤面投影用等经纬度会严重畸变。这个不是审美问题直接影响面积计算和梯度计算——用错投影算出来的面积能偏差百分之十几。重采样方法的选择也有讲究。连续场温度、湿度用双线性或者三次卷积能保持空间连续性分类场云掩膜、土地覆盖必须用最邻近否则会凭空造出不存在的类别辐射率场做重采样前最好先转成反射率或者亮温直接对原始计数值插值没有物理意义。计算一个具体例子极轨星某个载荷幅宽 2900 公里星下点 250 米扫描线上约有 11600 个可见光像元。如果全幅宽按 250 米格网存储一行就是 11600 个点。假设载荷有 25 个通道其中 5 个是 250 米、20 个是 1 公里量化位数按 12 到 14 比特算单轨数据量级在几个 GB 到十几 GB 之间。这个量级决定了你不能用单机内存硬啃得按扫描行分块读边读边处理。静止星的数据量算法不同。全圆盘按 2748×2748 像元、14 个通道、每通道 14 比特粗略估算单帧约 14 × 2748² × 14 ÷ 8 ≈ 1.85 亿字节接近 185 MB。10 分钟一帧一天 144 帧单星日数据量超过 25 GB。做长时间序列分析前先想清楚是全部下载还是按需裁剪——直接下全量会很快撑爆硬盘。5.3 读取、定标与可视化的完整流程下面这套流程是我自己常用的从文件读取到出图通用性比较强改改参数就能套到不同卫星上。第一步是环境准备。Python 生态比较顺手核心依赖是这几个pip install h5py netCDF4 numpy pyresample cartopy matplotlib scipy第二步是探测文件结构。HDF5 格式的数据集往往嵌套很深先看清维度名和属性再决定读哪些import h5py with h5py.File(satellite_file.h5, r) as f: def walk(name, obj): if isinstance(obj, h5py.Dataset): print(name, obj.shape, obj.dtype, dict(obj.attrs)) f.visititems(walk)第三步是辐射定标。可见光通道的计数值转反射率的通用形式是import numpy as np def dn_to_reflectance(dn, gain, offset, solar_zenith_deg, esd1.0): radiance dn * gain offset cos_sza np.cos(np.deg2rad(solar_zenith_deg)) cos_sza np.clip(cos_sza, 0.05, 1.0) # 避免低太阳高度角导致反射率爆炸 return radiance * np.pi * esd**2 / (cos_sza * 1361.0)这里有几个细节太阳高度角必须裁剪否则晨昏线附近会出现极大的假值日地距离修正esd在一年中变化约 3.4%做定量反演时不能忽略反射率超过 1 的情况在高反射云和雪面上是正常的不要粗暴截断。红外通道要转亮温用 Planck 反函数def radiance_to_bt(radiance, wl_um, c11.191042e-5, c21.4387752): # radiance 单位需与普朗克常数配套这里示例为常用简化式 return c2 / (wl_um * np.log(c1 / (wl_um**5 * radiance) 1.0))第四步是几何处理与重投影。极轨数据自带经纬度场可以用 pyresample 做面积定义和重采样静止数据用的是扫描角加卫星位置转换公式相对固定直接用官方给的查找表更快。from pyresample import geometry, kd_tree src geometry.SwathDefinition(lonslons, latslats) dst geometry.AreaDefinition( region, target area, region_id, {proj: lcc, lat_0: 30.0, lon_0: 110.0, lat_1: 25.0, lat_2: 45.0}, 2000, 2000, (-2000000.0, -2000000.0, 2000000.0, 2000000.0), ) result kd_tree.resample_nearest(src, data, dst, radius_of_influence3000)第五步是渲染。真彩色合成建议用 0.47、0.51、0.64 三个通道分别映射到蓝、绿、红先做瑞利散射的粗略校正简单的做法是按大气光学厚度经验扣除蓝光通道的偏移再做 Gamma 增强。红外图建议用亮温表示色标换成冷色代表低温的递减色带比彩虹色更容易读。水汽通道用 6.2 微米色带用蓝到白的渐变能直观看出干区下沉和湿区上升的气流结构。实操心得做动画比做单帧有价值得多。把连续时刻的亮温按固定色标渲染成序列用 ffmpeg 合成云顶快速冷却的区域会自动跳出来这个直观性远超任何数值指标。6. 常见问题与排查技巧实录6.1 条带噪声与定标异常怎么处理极轨和静止数据都会出现条带成因不同处理方式也不一样。探测器响应不一致造成的条带最常见尤其在扫描式成像仪上。现象是同一条通道的图像里出现规律性的竖条纹间隔对应探测器数量。排查方法把图像沿扫描方向做列平均画出曲线如果看到周期性的锯齿就基本确认。处理方式有两种一是用星上定标数据做逐探测器增益校正二是用场景统计做归一化——取一段均匀场景比如大片洋面或均匀云层计算每个探测器的均值把所有探测器拉伸到同一均值。扫描镜速度波动造成的条带间隔不规律通常和卫星姿态或扫描机构有关。这种条带用场景统计很难完全消除只能靠低通滤波压制代价是损失细节。如果业务对细纹理要求高建议换用经过地面系统重处理的 L1 产品。红外通道的非线性在高温端或者低温端会表现出来。现象是亮温反演出来的值在冷云顶附近偏差大或者水陆边界出现异常跳变。处理方法是用星上提供的查找表做非线性校正别自己拟合二次多项式因为不同批次的探测器响应不一样。可见光通道的衰减是长期问题。随着时间推移光学元件老化、紫外照射会让响应下降表现为同样的地表目标反射率逐年降低。做多年趋势分析时如果发现某个区域的植被指数整体单调下降先别高兴很可能不是植被变了是定标漂了。这时候要做交叉定标找一颗在轨时间重叠的新卫星在重叠期做同步观测比对拟合出校正系数。6.2 时间、坐标与几何定位的坑这类问题最隐蔽因为图像看起来没毛病但数据是错的。时间戳时区混用。有些产品用 UTC有些用本地时有些文件名里是观测开始时间、文件属性里是结束时间。做多源融合时我曾经把两个相差 8 小时的数据当成同一时刻拼在一起结果云图完全不匹配查了两天才发现是时区问题。建议在读取阶段统一转成 UTC并在代码里显式标注。轨道号与日期不对应。极轨卫星的轨道号是连续递增的跨日时轨道号会回绕用轨道号做时间索引容易出错。更稳妥的做法是用文件里的观测起始时间做索引。地理定位偏移。静止星在轨运行久了姿态会有微小漂移L1 产品的定位精度会下降。判断方法找一个海岸线明显、地形对比强的区域把图像上的海岸线和真实矢量线叠一下偏移超过一个像元就要考虑用地面控制点做二次配准。极轨星的定位误差在山区更明显因为地形起伏会影响定位精度。高纬度投影畸变。用等经纬度投影处理 60 度以上区域时东西方向的像元尺寸会被严重拉伸做面积统计会出错。解决方法改用极射赤面或者兰勃特投影或者在等经纬度格网上加一个面积权重因子乘以纬度的余弦。扫描角与天顶角混淆。静止星的数据里有两个角卫星天顶角和太阳天顶角。辐射定标要用太阳天顶角大气校正要用卫星天顶角两者混用会让反演结果系统性偏差。产品说明里通常写得很清楚但字段名缩写很像容易看错。6.3 常见问题速查表现象可能原因排查动作处理方式图像出现等间隔竖条纹探测器响应不一致列平均曲线看周期性逐探测器增益归一化亮温在冷云顶偏差大红外非线性未校正与星上查找表比对应用非线性校正表多年植被指数单调下降可见光通道定标漂移与同期新星交叉比对拟合交叉定标系数多源数据时间对不齐时区或时间字段含义不同逐个检查文件时间属性统一转换到 UTC海岸线错位一个以上像元姿态漂移或几何校正不足叠加矢量岸线目视检查控制点二次配准高纬区域面积统计偏大等经纬度投影畸变计算像元实际面积换投影或加余弦权重海上出现异常高反射率斑点太阳镜面反射耀斑检查太阳高度角与相对方位标记剔除或做耀斑掩膜夜间低云识别不出来仅用单一红外通道查看 3.9 与 10.8 微米差分改用双通道差分法降水产品在小雨区偏高滴谱假设不适配与地面雨量站比对做本地化系数订正拼图出现条带缺口相邻轨道重叠不足核对轨道间隔与幅宽扩大时间窗口或多星融合7. 场景选星与组合使用经验7.1 不同业务场景该拿哪颗星的数据选星这件事本质上是在时间分辨率、空间分辨率、穿透能力、覆盖范围四个维度里做取舍。我按几类常见场景整理了一下思路。短时强对流监测。核心需求是快10 分钟一次是底线1 分钟一次最好。首选静止星看云顶亮温的降温速率、云团面积的扩张速度、上冲云顶的出现。配合闪电成像仪的数据如果闪电密度在十几分钟内陡增同时云顶快速冷却强天气基本要来了。极轨星在这里只能做补充因为过境间隔太长。台风路径与强度。路径靠环境场和引导气流静止星的水汽通道能看清引导气流和干空气入侵强度靠暖心结构和眼壁特征微波成像仪能穿透云层看到暖心是判断快速增强的关键。实际操作中我一般这么组合静止星做逐 10 分钟的连续监测极轨微波载荷每天几次提供内部结构散射计提供海上风场实况三者叠加判断。地表温度与城市热岛。必须用极轨高分辨率数据250 米到 1 公里级别静止星的 4 公里红外太粗。选择晴空、高太阳高度角的过境时刻做分裂窗反演。注意城市地表比热和发射率与郊区差异大反演时要用地表发射率产品做校正否则白天城市温度会被系统性低估。农业长势与干旱。用极轨的可见光近红外通道算植被指数配合热红外算温度植被干旱指数。时间上按旬或者半月合成避开云污染。这里有个技巧用最大值合成法取一段时间内植被指数的最大值能有效抑制残余云的影响因为云会让植被指数偏低取最大值就自动把云筛掉了。海面风场与海浪。散射计是主力但它对近岸和冰区不友好。实际用的时候要把 25 公里分辨率的散射计风场和静止星的云导风、微波辐射计的海温结合起来做一致性检查。台风外围大风区经常出现散射计饱和这时候要看亮温异常或者用被动微波的风速反演做补充。大气环境与气溶胶。静止星的可见光多通道可以做到 10 分钟一次的气溶胶光学厚度对污染过程的追踪很有用极轨的高光谱和偏振载荷能提供更精细的粒子分类。做沙尘识别时10.8 与 12.0 微米差分的符号是最直接的判据沙尘区通常表现为负值或者明显小于云区。7.2 多源数据融合的几个实操原则融合不是简单叠加几个原则踩过坑之后总结出来的。先统一时空基准再做任何计算。所有数据先转 UTC、转同一投影、重采样到同一格网。这一步看似机械但出问题最多。我的习惯是每次融合前先做一次空跑把所有源数据的空间范围、时间戳、有效值范围打印出来人工核对一遍再往下走。保留质量标记别急着填洞。各家的 L2 产品都带质量控制字段云掩膜、置信度、有效性标记。做统计时先用质量标记筛一遍再决定是否填补。盲目用插值填云洞会让数据看起来完整实际上引入大量虚假信号做趋势分析时后果很严重。分辨率对齐取最粗而不是最细。融合不同分辨率的数据时统一到最细的格网看起来很美好但低分辨率产品被插值放大后会给人高分辨率的错觉实际上信息量没有增加。合理做法是统一到最粗的公共分辨率或者用降尺度方法依赖高分辨率辅助变量而不是简单插值。交叉验证是必须的。两套独立观测在同一时空点上的差异往往能暴露定标、几何或者反演算法的问题。我一般会挑几个典型天气过程做对比晴空地表、厚云区、海岸线、高山积雪这几类目标对定标误差和几何误差最敏感。长期序列要注意代际断点。卫星有寿命新老交替时载荷指标会变。做十年以上的趋势分析必须处理断点找出重叠运行期做交叉定标把老数据调整到新基准上。这一步做不好趋势分析的结论可能是完全反的。7.3 我踩过的几个典型坑第一个坑是把快速扫描数据当常规数据用。快速扫描模式下积分时间短、信噪比低拿去做微弱温度差异的分析噪声会盖过信号。后来养成习惯凡是定量反演先看产品说明里的噪声指标快速扫描数据只用于定性判断和动态追踪。第二个坑是忽视太阳高度角对可见光通道的影响。早期做真彩色合成没做太阳天顶角归一化结果同一区域不同季节的图像亮度差一大截看起来像是地表变了。后来在预处理里固定加上反射率归一化问题消失。第三个坑是用单帧判断对流强度。有一回看红外云图某块云顶亮温很低判断是强对流结果后续发展很弱。后来看动画才发现那块冷云顶是高层卷云砧下面的对流早就衰亡了。判断对流要看云顶亮温的变化率而不只是绝对值还要结合水汽通道看是否还有水汽输送。第四个坑是把投影误差当成数据错误。有次拼接区域产品发现北部边缘有一片空白以为是数据缺失查了半天才明白是投影转换后目标格网超出了源数据覆盖范围。这类问题现在我会在重采样前先算一次源数据的实际覆盖多边形和目标区域求交提前知道哪里会缺。第五个坑是忽视星下点与边缘的差异。静止星的边缘视场因为观测路径长、大气路径倾斜同样的地表目标在图像边缘和中心的亮温能差好几开尔文。做区域对比分析时如果不按卫星天顶角做归一化边缘区域的系统性偏差会被误读为真实的空间差异。提示所有卫星数据都有适用范围超出适用范围使用不一定报错但结果会悄悄偏掉。用之前花十分钟读一遍产品说明比事后查三天问题划算得多。最后分享一个自己常用的做法建一个卫星档案表把常用卫星的轨道类型、过境地方时、通道配置、分辨率、产品级别、数据延迟、已知问题都记在一张表里。项目一来先查表定方案比每次从头调研快得多也能避免重复踩同一个坑。这张表我自己维护了几年现在遇到任何新的观测需求基本能在半小时内定出用哪几颗星、拿哪一级产品、做什么预处理。
