基于气象预测数据优化卫星成像参数

基于气象预测数据优化卫星成像参数 Optimize satellite imaging parameters using meteorological forecast data 1 first-author 白　杨（1983-），女，辽宁鞍山人，博士，研究员，2006年于东北师范大学获得学士学位， 2011年于中国科学院长春光学精密机械与物理研究所获得博士学位，主要从事卫星遥感数据处理等方面的研究。E-mail： baiy776@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225299&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225308&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225303&type= 22.01333618 32.00399780 白　杨 （1983-），女，辽宁鞍山人，博士，研究员，2006年于东北师范大学获得学士学位， 2011年于中国科学院长春光学精密机械与物理研究所获得博士学位，主要从事卫星遥感数据处理等方面的研究。E-mail： baiy776@163.com baiy776@163.com 1 corresp E-mail： wll920707@sina.cn E-mail： wll920707@sina.cn 王灵丽（1989-），女，河南柘城人，工程师，2016年于大连理工大学获得硕士学位，主要从事卫星遥感影像数据辐射校正处理、图像去噪与增强算法、大气校正与定量遥感等研究。E-mail： wll920707@sina.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225311&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225319&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225316&type= 22.01333618 32.00401306 王灵丽 （1989-），女，河南柘城人，工程师，2016年于大连理工大学获得硕士学位，主要从事卫星遥感影像数据辐射校正处理、图像去噪与增强算法、大气校正与定量遥感等研究。E-mail： wll920707@sina.cn wll920707@sina.cn 1 1 1 1 2 2 长光卫星技术有限公司，吉林 长春 130051 长光卫星技术有限公司，吉林 长春 130051 Chang Guang Satellite Technology co.， Ltd.， Changchun 130051， China Chang Guang Satellite Technology co.， Ltd.， Changchun 130051， China 中国气象局数值预报中心， 北京 100081 中国气象局数值预报中心， 北京 100081 Numerical Prediction Center， China Meteorological Administration， Beijing 100081， China Numerical Prediction Center， China Meteorological Administration， Beijing 100081， China 卫星传感器成像参数包括增益、积分级数、行积分时间等，需要在传感器成像前通过地面站指令上注，合理的成像参数设置是获取高质量遥感影像数据的重要前提。本文提出了一种基于气象预测数据的卫星传感器成像参数优化与计算方法，首先，根据成像链路模型，结合气象预测数据中能够较好反映地表状态变化的地表反照率数据，解算目标地物成像时刻天顶辐亮度；然后，利用传感器绝对辐射定标系数解算不同成像参数下对应的图像灰度；最后，根据图像灰度适宜的策略选取合理的成像参数。通过吉林一号视频04星推扫成像模式验证本文提出的算法对全色和多光谱谱段图像信息熵分别提升4.37%~23.71%，采用本文提出的优化成像参数可以有效改善卫星遥感影像的有效动态范围，满足高质量遥感影像数据获取需求。 Satellite sensor imaging parameters include gain， integration level， integration time， etc.， which need to be input by command through a ground station before imaging. Reasonable imaging parameter settings are important prerequisites for obtaining high-quality remote sensing image data. In this paper， a satellite sensor imaging parameter optimization method that employs meteorological prediction data is proposed. First， according to an imaging link model combined with surface albedo data， which can better reflect surface state changes in meteorological prediction data， the radiance at the top of the atmosphere of a target during the imaging time is determined. Then， the corresponding grayscale image under different imaging parameters is obtained using the absolute radiation calibration coefficient of the sensor. Finally， reasonable imaging parameters are selected based on the appropriate strategy of obtaining grayscale images. The algorithm proposed in this paper is verified using the line-push broom imaging mode of the Jilin SP04 satellite. The image comentropy is improved by 4.37%-23.71% for panchromatic and multispectral bands. The results show that the imaging parameter optimization strategy proposed in this paper can improve the effective dynamic range of images and meet the demand for high-quality remote sensing image data acquisition. 卫星成像参数 气象预测数据 成像链路模型 地表反照率 绝对辐射定标 satellite imaging parameters meteorological forecast data image-link model surface albedo data absolute radiation calibration 1 引　言 1 光学卫星在对地观测获取遥感影像前，需要通过上注指令设置传感器成像参数，用于获取高动态范围的高质量遥感影像数据，如视频相机一般需要设置增益、帧频、单帧图像曝光时间，推扫相机需要设置积分级数、行积分时间、增益等必要的成像参数 ［ 1 1 - 2 2 1-2 ］ 。成像参数设置不合理会极大的降低数据质量，当成像参数设置过低时，获取的影像数据有效动态范围小，降低了影像对地物的解析和判读能力，当成像参数设置过高时，会导致图像的最大灰度超出传感器有效范围，而产生图像整体或部分区域灰度值饱和现象，使得数据缺失饱和区域地物的纹理及细节信息 ［ 3 3 - 5 5 3-5 ］ 。因此合理设置卫星成像参数是卫星获得高质量影像数据的重要保障和前提。 目前常规的光学卫星传感器成像参数设置方式是依据目标拍摄时刻的太阳高度角建立卫星成像参数查找表，这种设置方式简便易行，然而相同太阳高度角下拍摄目标点的地表反照率、大气环境以及卫星观测角度均会影响传感器的灰度响应，在不同地物拍摄过程中，难以保证图像均获得较好的拍摄效果。Skybox采用MODIS 大气光学厚度（Atmospheric Optical Thichness， AOT）历史平均数据，利用大气辐射传输模型6SV（Second Simulation of a Satellite Signal in the Solar Spectrum-Vector），预测曝光参数 ［ 2 2 ］ 。北京一号依据拍摄目标点太阳高度角，基于历史数据底图，最终根据不同的拍摄目标，如沙漠、海洋、城镇等，利用6S辐射传输模型预测曝光参数 ［ 6 6 ］ 。上述这两种方式均能在一定程度上改善图像的拍摄效果，但是基于历史数据预测的方法，对于地表随时间变化较大，或是难以根据历史数据预估的天气变化（如降雪等）的区域则无法获得理想的预测效果。如何将具有时变特性的地表反照率，以及大气变化产生的影响加入到光学卫星成像参数设置中，成为了光学卫星对地观测高效率获取高质量遥感影像数据的难题。 为了解决这一难题，本文提出了一种基于气象预测数据的卫星成像参数设置方法，方法基于气象预测数据数值预测模型获取未来10天不同时刻的地表反照率、总云量和雪层厚度数据，再根据拍摄目标点地表反照率和云量，结合MODTRAN（Moderate Resolution Atmospheric Transmission）大气辐射传输模型建立从地表到传感器的辐射传输模型，模拟不同成像参数下图像的平均灰度响应，依据图像灰度动态范围合理的策略选取合适的成像参数。为了验证本文提出的成像参数设置方法，利用吉林一号视频04星测试成像参数预测效果，测试结果表明，本文提出的方法预测精度能够较好地保证成像参数设置的合理有效。 2 影像检测系统的构成及其工作原理 1 在光学卫星对地面观测成像前，需通过任务规划设置卫星的成像参数，成像参数主要包括传感器的积分级数、增益、行积分时间等，并最终经地面站指令上注上传至卫星，用于保障卫星正常成像任务。卫星传感器的成像方式主要分为线阵传感器推扫成像（ 图1 图1 （a））和面阵传感器凝视成像（ 图1 图1 （b）），不同的传感器成像方式具有不同的成像参数及设置方式，本文主要研究线阵推扫成像方式传感器的成像参数获取方法。 1;2; 线阵推扫成像TDI CCD（Time Delay Integration Charge Coupled Device）传感器是将入瞳处辐亮度通过传感器光电转换获取图像灰度响应DN（Digital Number） ［ 7 7 ］ ，其定义如 式（1） 式（1） 所示： <math display="block" id="M1"><mi>D</mi><mi>N</mi><mo>=</mo><mfrac><mrow><mi mathvariant="normal">π</mi><msub><mrow><mi>A</mi></mrow><mrow><mi>d</mi></mrow></msub><mi>M</mi><msub><mrow><mi>t</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">t</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">4</mn><msup><mrow><mi>F</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>h</mi><mi>v</mi></mrow></mfrac><msub><mrow><mi>L</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>η</mi><mi>G</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225358&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225353&type= 37.59199905 7.78933334 ， 式中： <math id="M2"><mi>M</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225368&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225363&type= 3.38666677 2.87866688 为积分级数， <math id="M3"><msub><mrow><mi>t</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">t</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225378&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225373&type= 3.38666677 3.72533321 为传感器行积分时间， <math id="M4"><mi>G</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225386&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225380&type= 2.79399991 2.79399991 为转换增益，单位为 <math id="M5"><mi>D</mi><mi>N</mi><mo>/</mo><msup><mrow><mi>e</mi></mrow><mrow><mo>-</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225396&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225392&type= 11.09133339 3.47133350 ，可以在成像时进行调整； <math id="M6"><msub><mrow><mi>L</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225408&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225401&type= 3.55599999 3.72533321 为入瞳处谱段辐亮度，与拍摄地物相关； <math id="M7"><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225416&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225412&type= 2.70933342 3.72533321 为光学系统透过率，受大气参数影响； <math id="M8"><msub><mrow><mi>A</mi></mrow><mrow><mi>d</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225427&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225422&type= 4.14866638 3.72533321 为探测像元面积， <math id="M9"><mi>F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225436&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225431&type= 2.53999996 2.79399991 为光学系统相对孔径数， <math id="M10"><mi>η</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225446&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225441&type= 1.86266661 3.64066648 为传感器的量子效率，这些参数与传感器光学系统、电子学性能相关，在传感器设计完成后不再改变； <math id="M11"><mi>h</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225461&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225456&type= 1.94733334 2.79399991 为普朗克常量， <math id="M12"><mi>v</mi><mo>=</mo><mi>c</mi><mo>/</mo><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225470&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225466&type= 11.68400002 3.30199981 ， <math id="M13"><mi>c</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225481&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225476&type= 1.52399993 2.79399991 为光速。 根据 公式（1） 公式（1） ，在推扫成像中可调节的卫星成像参数为积分级数、增益和行积分时间。为了获取清晰的图像，成像参数行时间需要同像元运动速度相匹配 ［ 8 8 ］ 。由于卫星飞行过程中纬度的变化，该参数需要实时进行调整，因此卫星在成像期间会自主对该参数进行实时运算，无需通过成像前任务规划进行设置。 卫星规划推扫成像任务时，需要设置的卫星成像参数为积分级数和增益，由于不同传感器的设计差异或基于图像质量的考虑，也可能只调整积分级数或增益。若成像参数设置过高，会导致局部或整体图像饱和，丢失被观测地物的纹理细节信息；若成像参数设置过低，使得图像偏暗，从而降低数据的地物识别和解译能力。因此合理的设置卫星成像参数是充分利用传感器辐射分辨率，获取灰度动态范围适宜的图像，保障遥感影像数据质量的重要前提。 3 成像参数推演 1 3.1　气象预测数据 2 成像参数推演的目的是为了确定规划拍摄目标的积分级数和增益，保证图像灰度响应范围适宜。根据成像参数与图像灰度响应的关系，拍摄目标点辐亮度直接影响图像的灰度响应，而辐亮度主要受到地表地物反照率的影响。目前虽然能够基于历史地物反照率建立全球地物反照率底图用于卫星成像参数设置，然而地物的反照率并非固定不变，它会随着地表覆盖地物类型以及季节发生变化，且不同年份也会有差异。因此本文提出基于国家气象局提供的气象预测数据中的反照率数据进行成像参数优化的方法。 气象预测数据是中国气象局联合国内多家单位，在科技部“十五”国家重点科技攻关项目《中国气象数值预报系统创新研究》支持下，自主研究建立了新一代多尺度通用资料同化与数值预报系统GRAPES ［ 9 9 ］ 。从2016年6月开始，GRAPES全球同化预报系统正式业务化运行 ［ 10 10 ］ 。目前中国气象局数值预报中心实时业务系统提供的全球气象预测数据，能够预测未来10天的全球地表反照率、低云量、总云量和表层雪深，经过产品后处理后的全球数据栅格大小为1 440×2 880 pixels，网格水平分辨率为0.125°，预测数据时间分辨率为1 h。 图2 图2 为2017年3月2日19时的全球气象预测数据。 全球地表反照率气象预测数据能够较好地反映世界范围内不同区域的地表反照率，对于冰雪、沙漠、绿植都能较好地区分。同时，全球地表反照率预测数据对于随时间变化较大的地区同样具有较好的预测效果， 图3 图3 为长春地区2017年全年上午11时地表反照率变化曲线。 图3 图3 中横坐标为年积日，即当日在一年中的日序，纵坐标为地表反照率。根据 图3 图3 所示在积日为346，即2017年12月12日，长春出现降雪，地表反照率发生阶跃式地突变。可见，地表反照率预测数据能够很好地反映地表状态的改变，尤其对于降雪有较好的预测效果。 3.2　成像参数推演 2 建立遥感图像成像链路模型是进行成像参数推演的重要基础。在光学成像系统中，传感器获取的能量来源于地表反射和大气散射，太阳光从大气层外入射到地表，地物反射辐射太阳光，透过大气层进入卫星传感器，传感器将入瞳处辐亮度转化为图像灰度响应。若忽略遥感成像过程中由大气引起的光子吸收与散射，则其物理过程可如 图4 图4 表示。 根据遥感图像成像物理过程，结合气象预测数据中的地物反照率，进行成像参数推演。基于气象预测数据的成像参数推演流程如 图5 图5 所示。成像参数推演的具体流程为： （1）根据卫星任务规划输入规划拍摄的目标点经纬度，以及该目标点成像时刻； （2）根据成像时刻、目标点经纬度读取气象预测数据的地表反照率fal，并计算大气辐射传输模型需要的参数，包括太阳高度角和方位角、卫星观测高度角和方位角等。模型输入中，大气模型由拍摄点经纬度、拍摄季节确定 ［ 11 11 ］ 。由于不同地物有明显不同的反照率，故可参考地表反照率预报来粗略设置气溶胶类型，对于反照率小于0.08的地区，其气溶胶模型设置为海洋型；对于反照率在0.08~0.28区间的地区，气溶胶类型设置为乡村型；反照率在0.28~0.45区间的地区气溶胶类型设置为沙漠型，反照率大于0.45的地区设置为城市型。由于缺少能见度的预测数据及足够完整的全球能见度历史数据集，同时不同地区气溶胶分布具有明显局地性 ［ 12 12 ］ ，故参考MODIS的历史MOD08_M3大气产品中各年每月平均550 nm气溶胶光学厚度值，通过输入拍摄时刻前一年相同月份的平均550 nm气溶胶光学厚度来粗略代替表示能见度的特性 ［ 13 13 ］ 。最后将前述参数输入MODTRAN4大气辐射传输模型计算天顶光谱辐亮度； （3）根据各谱段的光谱响应曲线利用 公式（2） 公式（2） 计算对应谱段的辐亮度，并计算不同成像参数组合对应谱段辐亮度的灰度响应，其中图像灰度响应根据传感器绝对辐射定标系数图像灰度与辐亮度线性系数进行计算，计算公式如 式（3） 式（3） 所示，所计算出的灰度响应是成像目标区域的平均灰度响应： <math display="block" id="M14"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi></mrow></msub><mo>=</mo><mfrac><mrow><mo>∫</mo><msub><mrow><mi>R</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo><mi>L</mi><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>λ</mi></mrow><mrow><mo>∫</mo><msub><mrow><mi>R</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>λ</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225553&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225548&type= 39.28533173 14.30866623 ， 其中： <math id="M15"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225563&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225558&type= 6.43466663 3.72533321 为谱段辐亮度（band=1，2，…， n ）， L （ λ ）为利用大气辐射传输模型计算得到的天顶光谱辐亮度，单位均为 <math id="M16"><mi>W</mi><mo>/</mo><mo stretchy="false">(</mo><msup><mrow><mi mathvariant="normal">m</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>⋅</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">r</mi><mo>⋅</mo><mi mathvariant="normal">μ</mi><mi mathvariant="normal">m</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225572&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225568&type= 25.48466682 4.14866638 ， <math id="M17"><msub><mrow><mi>R</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225582&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225579&type= 12.19199944 4.48733330 为谱段 <math id="M18"><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225592&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225587&type= 7.19666672 2.79399991 的相对光谱响应曲线。 <math display="block" id="M19"><mi>D</mi><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>=</mo><mo stretchy="false">(</mo><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi></mrow></msub><mo>-</mo><mi>o</mi><mi>f</mi><mi>f</mi><mi>s</mi><mi>e</mi><msub><mrow><mi>t</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi><mo>,</mo><mi>i</mi></mrow></msub><mo stretchy="false">)</mo><mo>/</mo><mi mathvariant="normal">g</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">i</mi><msub><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi><mo>,</mo><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225602&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225598&type= 57.48866653 4.82599974 ， 其中： <math id="M20"><mi>D</mi><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi><mo>,</mo><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225612&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225607&type= 11.17599964 3.81000018 为对应谱段 <math id="M21"><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225622&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225617&type= 7.19666672 2.79399991 的辐亮度 <math id="M22"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225632&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225626&type= 6.43466663 3.72533321 在成像参数组合 <math id="M23"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225643&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225638&type= 1.10066664 2.79399991 （i=1，2，…， m ）时对应的灰度响应， <math id="M24"><mi mathvariant="normal">g</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">i</mi><msub><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>,</mo><mi>o</mi><mi>f</mi><mi>f</mi><mi>s</mi><mi>e</mi><msub><mrow><mi>t</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">d</mi><mo>,</mo><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225650&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225647&type= 27.09333420 4.57200003 为对应谱段和成像参数组合传感器绝对辐射定标线性系数的斜率和截距。 （4）设置图像灰度响应阈值，从不同参数组合对应的灰度响应中，查找与阈值最接近的成像组合参数作为该谱段最优成像组合参数，成像参数选取的原则是在图像不过曝的前提下尽量保证图像动态范围适宜，即避免获取的遥感影像数据在低地表反照率时过暗、在高地表反射率时过亮，故采用分段线性函数，并通过大量测试验证的效果对比设定阈值 thre ，其计算方法如 公式（4） 公式（4） 所示： <math display="block" id="M25"><mi>t</mi><mi>h</mi><mi>r</mi><mi>e</mi><mo>=</mo><mfenced open="{" close="" separators="|"><mrow><mtable><mtr><mtd><mo stretchy="false">(</mo><mn mathvariant="normal">0.75</mn><mi>f</mi><mi>a</mi><mi>l</mi><mo>+</mo><mn mathvariant="normal">0.1</mn><mo stretchy="false">)</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi><mo>_</mo><mi>D</mi><mi>N</mi><mo>,</mo><mtable><mtr><mtd/><mtd/></mtr></mtable><mi>f</mi><mi>a</mi><mi>l</mi><mo>≤</mo><mn mathvariant="normal">0.4</mn></mtd></mtr><mtr><mtd><mi>f</mi><mi>a</mi><mi>l</mi><mo>⋅</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi><mo>_</mo><mi>D</mi><mi>N</mi><mo>,</mo><mtable><mtr><mtd/><mtd/></mtr></mtable><mn mathvariant="normal">0.4</mn><mo><</mo><mi>f</mi><mi>a</mi><mi>l</mi><mo>≤</mo><mn mathvariant="normal">0.6</mn></mtd></mtr><mtr><mtd><mo stretchy="false">(</mo><mn mathvariant="normal">0.5</mn><mi>f</mi><mi>a</mi><mi>l</mi><mo>+</mo><mn mathvariant="normal">0.3</mn><mo stretchy="false">)</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi><mo>_</mo><mi>D</mi><mi>N</mi><mo>,</mo><mtable><mtr><mtd/><mtd/></mtr></mtable><mn mathvariant="normal">0.6</mn><mo><</mo><mi>f</mi><mi>a</mi><mi>l</mi></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225658&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225654&type= 71.79733276 16.67933464 ， 其中： <math id="M26"><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi><mo>_</mo><mi>D</mi><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225666&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225662&type= 13.88533401 3.30199981 为图像的最大灰度值， <math id="M27"><mi>f</mi><mi>a</mi><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225673&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225670&type= 4.14866638 3.64066648 为气象预测数据中获取的地表反照率。 图像灰度响应反映了通过成像链路模型推演的平均图像灰度，本文提出的图像参数选取策略是根据地表反照率设置图像灰度响应阈值，即在地表反照率偏低时，调高阈值，避免图像过暗，在地表反照率偏高时，降低阈值，避免高地表反照率引起图像饱和。在图像成像参数组合中选取灰度响应最接近阈值的一组作为成像参数来实现高动态范围的遥感影像数据获取。 为了定量化分析成像参数优化效果，利用信息熵的指标进行评估。信息熵主要反映图像包含地物信息的详细程度，一般用熵来表达，熵是从信息论角度反映影像信息丰富程度的一种度量方式，通过对图像信息熵的比较可以对比出影像的细节表现能力，信息熵Comentropy的计算如 公式（5） 公式（5） 所示： <math display="block" id="M28"><mi>C</mi><mi>o</mi><mi>m</mi><mi>e</mi><mi>n</mi><mi>t</mi><mi>r</mi><mi>o</mi><mi>p</mi><mi>y</mi><mo>=</mo><mo>-</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>k</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi><mo>_</mo><mi>D</mi><mi>N</mi></mrow></munderover></mstyle><mrow><mi>p</mi><mo stretchy="false">(</mo><msub><mrow><mi>b</mi></mrow><mrow><mi>k</mi></mrow></msub><mo stretchy="false">)</mo></mrow><mo>⋅</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><msub><mrow><mi mathvariant="normal">g</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mi>p</mi><mo stretchy="false">(</mo><msub><mrow><mi>b</mi></mrow><mrow><mi>k</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225679&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225676&type= 60.11333084 8.63599968 ， 其中： k 为灰度级数， <math id="M29"><msub><mrow><mi>b</mi></mrow><mrow><mi>k</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225685&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225682&type= 2.96333337 3.72533321 为第k级的灰度值， <math id="M30"><mi>p</mi><mo stretchy="false">(</mo><msub><mrow><mi>b</mi></mrow><mrow><mi>k</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225692&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225688&type= 7.87400007 4.65666676 为 <math id="M31"><msub><mrow><mi>b</mi></mrow><mrow><mi>k</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225698&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10225695&type= 2.96333337 3.72533321 出现的概率。 4 实验结果分析 1 为了验证本文提出的基于气象预测数据优化卫星成像参数算法的有效性，使用吉林一号视频04星在轨进行推扫成像获取的标准影像产品，并对获取的影像数据进行统计分析。 吉林一号视频04星（JL104B）是长光卫星技术有限公司于2017年11月成功发射入轨的高分辨率光学卫星，卫星通过搭载2台0.92 m高分辨率相机构成了双相机结构，每台相机均具备推扫和视频成像能力，本文使用其推扫模式成像验证所提出的算法，推扫模式下成像主要参数如 表1 表1 所示。推扫图像包含全色和多光谱谱段，其中多光谱所包含的5个谱段分别为Blue，Green，Red，Red Edge和NIR，双相机幅宽20 km。 成像参数 参数值 地面像元分辨率 0.92 m PAN, 3.68m MSS 光谱响应范围 PAN：450~800 nm Blue: 460~520 nm Green: 510~580 nm Red: 630~690 nm Red Edge: 700~750 nm NIR: 770~890 nm 量化位数 12 bits 幅宽 10.6 km(单相机),20 km(双相机) 吉林一号视频04星单景影像幅宽与气象预测数据像元分辨率0.125°基本相当，因此能够利用气象预测数据获取对应拍摄目标点的地表反照率。 根据文中提出的基于气象预测数据优化卫星成像参数的方法对比成像参数优化前后的传感器成像效果，实现优化参数的定量评价，其中 图6 图6 为长春区域雪景拍摄对比图， 表2 表2 为成像参数优化前后图像灰度与信息熵统计结果。 表2 表2 中饱和率（Saturation factor）为图像饱和像元的比率，Max（1%）和Min（1%）分别为对应图像灰度直方 图1 图1 %的灰度级，灰度响应范围Range为Max（1%）-Min（1%）。 1;2; 谱段 项目 优化前 优化后 比例 6 PAN Mean 3 797 1 679 Max(1%) 4 086 2 497 Min(1%) 2 201 1 011 Range 1 885 1 487 -21.13% Comentropy 8.31 10.27 23.71% Saturation factor 77.80% 0 图6 图6 （a）和（b）分别为长春区域（Lat：44.53，Lon： 125.60）2018年12月30日和2019年01月03日雪景全色图像，2018年12月30日成像根据目标点成像时刻太阳高度角19.3°设置增益为8.0x，积分级数为10级，2019年01月03日成像采用本文提出的成像参数优化策略，结合地表反照率预测数据，成像增益设置为4.0x，积分级数为10级。 图6 图6 （a）中出现了雪景局部过曝的现象， 表2 表2 列出图像饱和率为77.80%，图像部分细节信息丢失，而 图6 图6 （b）为采用成像参数优化后的成像效果，较好得保留了雪表面的特征，消除了图像局部饱和的现象，相比于成像参数优化前信息熵提升了23.71%。可见，本文提出的方法对于改善图像局部过曝有较好的效果。 同样地，对低照度区域采用本文提出的方法进行在轨成像测试， 图7 图7 为亚松森拍摄成像参数优化前后对比图，并列出各个谱段的灰度分布直方图， 表3 表3 为亚松森成像参数优化前后图像灰度与信息熵对比结果。 1;2;3;4; 谱段 项目 优化前 优化后 比例 3 Blue Mean 114 166 Range 132 127 -3.79% Comentropy 6.05 6.39 5.54% 3 Green Mean 165 233 Range 230 232 0.87% Comentropy 6.90 7.37 6.79% 3 Red Mean 103 147 Range 252 263 4.37% Comentropy 7.09 7.39 4.24% 3 Red Edge Mean 100 171 Range 154 259 68.18% Comentropy 6.58 7.81 18.63% 3 NIR Mean 291 483 Range 403 810 100.99% Comentropy 8.31 9.58 15.32% 图7 图7 （a）和（c）分别为亚松森（Lat： - 24.84，Lon： - 57.38）2018年12月5日和2019年3月21日植被假彩色图像（NIR-Red-Green）， 图7 图7 中（b）和（d）两次成像对应的各个谱段灰度直方图，根据灰度直方图，图9（c）图像灰度动态范围更大，而 图7 图7 （a）中Red谱段的平均灰度仅为103，整体灰度值偏低， 图7 图7 （c）成像效果优于 图7 图7 （a）。该目标点成像根据太阳高度角（65.8°和63.9°）设置成像增益为1.0x，积分级数为4级，而采用本文提出的成像参数优化策略，由于该目标点地表反照率偏低，2019年3月21日（图9（c））将成像增益调整为2.0x，积分级数4级，使得成像效果得到显著提升。根据 表3 表3 中统计结果，在成像参数优化后各个谱段的平均灰度响应增加，信息熵分别提高4.37%~18.63%。可见，本文提出的方法对于低亮度图像灰度有明显地改善，使得低亮度区域细节信息更为清晰（彩图见期刊电子版）。 根据视频04星推扫成像模式的实际在轨成像效果，采用本文提出的成像参数优化算法能够有效保证图像在不过曝的前提下获取适宜的图像灰度动态范围。 5 结　论 1 光学遥感卫星是获取对地观测遥感数据的重要手段，为了保障在轨运行的卫星对地成像能够获得合理的动态范围的高质量遥感数据，提升卫星获取数据的有效率及数据解译能力，本文提出了将气象卫星数据用于陆地观测卫星规划中，通过结合气象预测数据和卫星传感器绝对辐射定标结果利用辐射传输模型核算在轨拍摄图像灰度响应，为选取卫星合理的成像参数提供有效依据，从而在有效避免图像饱和的同时极大的提升了低亮度区域的图像灰度，从根本上提升了卫星获取有效数据的能力及所获取数据的解译能力。通过吉林一号视频04星在轨成像进行验证，实验证明本文提出的算法对全色和多光谱谱段图像信息熵分别提升4.37%~23.71%，对比成像参数优化策略使用前后获取的影像数据，充分证明了本文提出的方法对于图像过曝与图像灰度过暗具有良好的改善效果。进一步的研究计划将关注提升图像灰度预测精度，根据卫星传感器长期在轨成像特性跟踪精细设置灰度动态选取范围。 参考文献 MURTHY K ， SHEARN M ， SMILEY B D ， et al .. SkySat-1： very high-resolution imagery from a small satellite ［C］. Sensors ， Systems ， and Next-Generation Satellites XVIII. International Society for Optics and Photonics ， 2014 ， 9241 ： 92411E . SkySat-1： very high-resolution imagery from a small satellite 袁航飞 . 航天相机星上自适应调光的研究 ［D］. 长春 ： 中国科学院研究生院 ， 2016 . 航天相机星上自适应调光的研究 YUAN H F . Research on Automatic Adaptive Exposure for Space-borne Remote Sensing Camera ［D］. Changchun ： Graduate school of Chinese academy of sciences ， 2016 . （in Chinese） Research on Automatic Adaptive Exposure for Space-borne Remote Sensing Camera 曹旗磊 ， 董丽丽 ， 赵明 等 . 空间相机自适应曝光 ［J］. 光子学报 ， 2016 ， 45 （ 8 ）： 0822003 - 0822008 . 空间相机自适应曝光 CAO Q L ， DONG L L ， ZHAO M ， et al . . Adaptive exposure of space camera ［J］. Acta Photonica Sinca ， 2016 ， 45 （ 8 ）： 0822003 - 0822008 . （in Chinese） .Adaptive exposure of space camera 曹旗磊 ， 董丽丽 ， 赵明 等 . 低分辨率测光指导高分辨成像自适应调光 ［J］. 激光与光电子学进展 ， 2017 ， 54 （ 9 ）： 374 - 380 . 低分辨率测光指导高分辨成像自适应调光 CAO Q L ， DONG L L ， ZHAO M ， et al . . Low resolution metering guide for high resolution imaging adaptive exposure ［J］. Advances in laser and optoelectronics ， 2017 ， 54 （ 9 ）： 374 - 380 . （in Chinese） . Low resolution metering guide for high resolution imaging adaptive exposure 何春良 ， 李斌桥 ， 刘振旺 等 . TDI CMOS图像传感器曝光时间优化方法研究 ［J］. 光学学报 ， 2015 ， 35 （ 2 ）： 34 - 40 . TDI CMOS图像传感器曝光时间优化方法研究 HE CH L ， LI B Q ， LIU ZH W ， et al . . Integration time optimization for TDI CMOS image senor ［J］. Acta Optica Sinca ， 2015 ， 35 （ 2 ）： 34 - 40 . （in Chinese） . Integration time optimization for TDI CMOS image senor CHEN S T ， CHENG J H ， GAO W . Exposure adjustment of satellite cameras ［C］. Proceedings of 2009 IEEE International Geoscience and Remote Sensing Symposium ， July 12-17 ， 2009 ， Cape Town . New York ： IEEE ， 2009 ， 978 ： 857 - 860 . Exposure adjustment of satellite cameras 韩昌元 . 航天相机成像质量评价 ［J］. 光电工程 ， 2002 ， （ S1 ）： 11 - 24 . 航天相机成像质量评价 HAN CH Y . Imaging quality evaluation of space camera ［J］. Opto-Electronic Engineering ， 2002 ， （ S1 ）： 11 - 24 . （in Chinese） Imaging quality evaluation of space camera 曹彬才 ， 朱述龙 ， 邱振戈 ， 等 . TDI行积分时间跳变对视差的影响及纠正方法 ［J］. 测绘科学技术学报 ， 2015 ， 32 （ 6 ）： 63 - 67 . TDI行积分时间跳变对视差的影响及纠正方法 CAO B C ， ZHU SH L ， QIU ZH G ， et al . . Analysis of integration time adjustment’s impact on parallax and it’s correction ［J］. Journal of Geomatics Science and Technology ， 2015 ， 32 （ 6 ）： 63 - 67 . （in Chinese） . Analysis of integration time adjustment’s impact on parallax and it’s correction 王金成 ， 陆慧娟 ， 韩威 ， 等 . GRAPES 全球三维变分同化业务系统性能 ［J］. 应用气象学报 ， 2017 ， 28 （ 1 ）： 11 - 24 . GRAPES 全球三维变分同化业务系统性能 WANG J CH ， LU H J ， HAN W ， et al . . Improvements and performances of the operational GRAPES_GFS 3DVar system ［J］. Journal of Applied Meteorological Science ， 2017 ， 28 （ 1 ）： 11 - 24 . （in Chinese） . Improvements and performances of the operational GRAPES_GFS 3DVar system ​HAN W ， BORMANN N . Constrained adaptive bias correction for satellite radiances assimilation in the ECMWF 4D-Var ［C］. EGU General Assembly Conference Abstracts ， 2016 . Constrained adaptive bias correction for satellite radiances assimilation in the ECMWF 4D-Var BERK A ， ANDERSON G P ， ACHARYA P K ， et al .. MODTRAN 4 User ’ s Manual ［M］. Washington ： Air Force Research Laboratory ， 1999 . MODTRAN4User’s Manual 刘泽阳 ， 李学彬 ， 孙刚 ， 等 . 中国部分地区的春季气溶胶光学特性 ［J］. 光学 精密工程 ， 2020 ， 28 （ 5 ）： 1021 - 1028 . 中国部分地区的春季气溶胶光学特性 LIU Z Y ， LI X B ， SUN G ， et al . . Optical properties of spring aerosol in typical areas of China ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 5 ）： 1021 - 1028 . （in Chinese） . Optical properties of spring aerosol in typical areas of China KING M D ， PLATNICK S ， MENZEL W P ， et al . . Spatial and temporal distribution of clouds observed by MODIS onboard the terra and aqua satellites ［J］. IEEE Transactions on Geoscience & Remote Sensing ， 2013 ， 51 （ 7 ）： 3826 - 3852 . .Spatial and temporal distribution of clouds observed by MODIS onboard the terra and aqua satellites 10.37188/OPE.20212903.0628 V321.2+53 A 1004-924X（2021）03-0628-09 白杨,王灵丽,武红宇等.基于气象预测数据优化卫星成像参数[J].光学精密工程,2021,29(03):628-636. 白杨,王灵丽,武红宇等.基于气象预测数据优化卫星成像参数[J].光学精密工程,2021,29(03):628-636. BAI Yang,WANG Ling-li,WU Hong-yu,et al.Optimize satellite imaging parameters using meteorological forecast data[J].Optics and Precision Engineering,2021,29(03):628-636. BAI Yang,WANG Ling-li,WU Hong-yu,et al.Optimize satellite imaging parameters using meteorological forecast data[J].Optics and Precision Engineering,2021,29(03):628-636. 吉林省科技发展计划资助项目（No.20180201108GX）；海南省重大科技计划资助项目（No.ZDKJ2019007） 2020-11-20 2021-01-06 2021-03-15 2021-04-08 光学精密工程 3 29