“京师一号”极地小卫星宽幅影像辐射校正方法研究_张卓宇

《“京师一号”极地小卫星宽幅影像辐射校正方法研究_张卓宇》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

北京师范大学学报（自然科学版）2023-02104Journal of Beijing Normal University（Natural Science）59（1）“京师一号”极地小卫星宽幅影像辐射校正方法研究*张卓宇†　章　影（北京师范大学全球变化与地球系统科学研究院，遥感科学国家重点实验室，100875，北京）摘要　“京师一号”小卫星是中国自主研制的首颗民用极地光学卫星，搭载有中高分辨率的多光谱宽幅相机，兼顾重访周期和空间分辨率，具有高覆盖、中高分辨率的优势，对极地变化监测研究具有重大意义；但其0级数据的条带噪声掩盖图像的地物特征，无法反映地物的真实辐射特性，严重影响了数据应用.为此，本文提出了一种基于统计学原理，针对宽幅影像去除条带的方法.实验结果表明，相较于原始数据，该方法较能稳健地去除数据中的条带噪声，还原地表真实辐射信息，提升了“京师一号”宽幅相机多光谱数据的应用价值.关键词　京师一号；南极；条带噪声；辐射校正中图分类号　TP752.1DOI：10.12202/j.0476-0301.20212092019年9月12日上午11时25分，中国首颗专用生产应用的关键所在.于极地气候与环境监测的遥感小卫星“京师一号”经过30多年的发展，已经形成了几十种不同的（研制代码：BNU-1）在中国太原卫星发射中心成功发去除条带影像的方法.现有去除条带噪声的处理方射并入轨.它由3个有效载荷组成：多光谱宽幅相机法可分为3大类[3].在可见光和近红外光谱范围内工作，每像素12 bit辐1）定标法.分为实验室定标、星上定标、场地定射分辨率，可用于监测极区变化；另外2个传感器是标法.实验室定标是使用积分球[4]和黑体进行定标；窄视场高分辨率遥感相机和船舶AIS接收机，能够实星上定标分为主动内置定标和被动向阳定标，还可引现遥感同步观测和船舶定位[1].宽幅相机包括4个光入漫射板进行定标[5]；场地定标作为弥补在轨传感器谱波段的CMOS面阵传感器，采用推扫式成像方式，发生变化的替代方法，选择大范围地势平坦、能见度分为PMS1和PMS2左右2个焦平面.每波段开窗选高、反射率均一的场景代替标准光源[6].取成像质量最好的一排用于采集信息，每排有5 0562）基于数学模型的方法.常见的空间域统计方法个电子元件，分景的0级产品每列有5 000个像素.有直方图匹配、直方图均衡化[7]、矩匹配等，这类方法宽幅相机具有高覆盖、中高分辨率的优点，可以假设不同探测器件所探测的地物应具有相同的辐射获得高质量的光学影像.但由于光学镜头造成的不分布特征.通过计算条带与正常影像灰度统计信息均匀性、国产卫星镀膜技术受限、CMOS各传感器的的相关关系来去除条带噪声[8]，计算效率高，但易出响应不一致、CMOS的暗噪声及热环境的不均匀等原现条带噪声去除不够彻底或处理过度等问题[9]；基于因，导致小卫星各焦平面的各波段0级产品有垂直于变分的条带去除是利用正则化技术对条带特征进行航迹的条带存在，从而影响了对影像的视觉解释；同刻画，并采用整体求解的方式解算无噪声的影像[10]；时，由于缺少星上定标设备，冰雪的反照率高而小卫滤波法分为频率域滤波和空间域滤波2种，频率域滤星辐射分辨率低，条带噪声表现出复杂的非线性特征[2]，波法是利用傅里叶变换、小波变换[11]将空间域的遥以及卫星的无控几何定位精度远大于1个像元给影感影像变换到频率域，通过不同滤波器滤除噪声对应像的条带处理带来了挑战.的频率，最后将滤波结果变换到空间域后得到去除条卫星发射后，往往需要对遥感器进行在轨测试，带后的影像.以评估轨运行效果，之后可以投入正式使用.其中，3）综合法.科学家将多种校正方法和算法进行综对遥感器的辐射性能定量评估是后续遥感产品投入合运用，例如形状检测算子结合空间插值算法对影像*国家自然科学基金资助项目（41925027）† 通信作者：张卓宇（1997—），女，硕士，助理工程师.研究方向：极地遥感. E-mail：201921490047@mail.bnu.edu.cn收稿日期：2021-08-27

1第1期张卓宇等：“京师一号”极地小卫星宽幅影像辐射校正方法研究105的坏线进行补充[12]，又如综合运用实验室定标和均匀反差极小，导致影像的噪声给图面质量带来不良影场景统计的校正方法.响.目前专门针对南极地区影像条带噪声去除的文作为新发射的小型光学卫星，“京师一号”极地卫献相对较少，且年代久远[14]，选择南极地区部分影像星具有不同太阳高度角的成像条件，宽幅相机的各波作为研究重点有助于填补这一方面的研究空白.段间纵向条带无明显周期性，而不同焦平面各波段间条带具有相似的空间分布的特点，同时因极地小卫星2    数据和方法镜头制造、镀膜及实验室相对校正定标不准确造成运用“京师一号”卫星影像，采用统计和均匀场景了卫星影像相对辐射精度不高，在影像中间出现偏色法对南极地貌的0级真彩色影像进行辐射校正，并利和条纹，影像边缘出现彩色等具体情况.因此“京师用同时期相同波段的MODIS影像对辐射校正效果进一号”的辐射校正问题亟待解决，故本文选取了覆盖行精度验证.“京师一号”的多光谱宽幅相机，幅宽南极典型地貌的区域作为研究区，旨在探究去除小744 km，星下点分辨率约74 m，总计有4个波段，地面卫星影像彩色条带的方法，提升“京师一号”的应用分辨率均为80 m，该数据来自北京师范大学全球变化与价值.地球系统科学研究院参与研制的地面引接系统（http://earth.bnu.edu.cn/index）.此外，本文使用了HEG工具对自1    研究区域https://ladsweb.modaps.eosdis.nasa.gov/search/网站获取南极洲作为地球的最南端，包括南极大陆的冰的MODIS影像进行处理，进行波段合成，显示出地物盖、冰架及周边岛屿.四周被大西洋、印度洋、太平的真实色彩. MODIS为搭载在Terra和Aqua这2颗洋包围，南极大陆边缘有隶属南太平洋的别林斯高晋卫星上的中分辨率成像光谱仪，是美国地球观测系统海、罗斯海、阿蒙森海和属于南大西洋的威德尔海（EOS）计划中用于观测全球生物和物理过程的重要等.作为当今淡水的重要储存地，南极洲蕴藏全球仪器.它拥有36个中等分辨率（0.25~1.00 μm）的光谱72%的淡水资源，景观以冰雪地貌为主.南极大陆气波段，每1~2 d可以对地球表面观测1次.本文仅使候严寒干燥，是世界上平均海拔最高、最冷、风力最用了其中红绿蓝3个波段的数据，详细记录可参考大的陆地[13].冰雪在卫星影像上亮度极大，对比微弱，表1，辐射校正前的示例影像记录见表2和图1.表 1 卫星参数对照MODISBNU-1波段波段序号波长范围/nm地面分辨率/m波段序号波长范围/nm地面分辨率/m红1620~6702503638~676/642~677绿4545~5655002560~599/562~60280蓝3459~4795001471~506/473~508全色0527 ~769/532~764表 2 本文使用的示例数据序号BNU-1MODIS特征地物MOD02QKM.A2019281.0405.061.20192811317511BNU1_PMS2_20191008044851_200214500_102_0001_001云、海岸MOD02HKM.A2019281.0405.061.2019281131751MOD02QKM.A2019321.2125.061.20193221743252BNU1_PMS1_20191117201830_200001045_102_0003_001云、海岸MOD02HKM.A2019321.2125.061.2019322174325MOD02QKM.A2019321.1810.061.20193221741333BNU1_PMS2_20191117201940_200018633_101_0001_001冰盖MOD02HKM.A2019321.1810.061.2019322174133MOD02QKM.A2019323.1935.061.20193251713404BNU1_PMS1_20191119202826_200018637_102_0004_001海岸MOD02HKM.A2019323.1935.061.20193251713402.1    辐射校正原理对“京师一号”小卫星的传感器像模块：CMOS传感器和后处理电子链路.每个探测来讲，从入射光信号到输出数字灰度值要经过光学成器上测得的电压与探测器在积分时间内收集的光子

2106北京师范大学学报（自然科学版）第59卷62°64°66°E175°180°W175°170°165°160°155°150°ES79°00′S70°00′78°00′70°30′77°00′76°00′71°00′75°00′050100200km74°00′050100200km170°180°W170°160°150°140°130°E180°170°160°150°ES78°00′S84°00′76°00′82°00′050100200km050100200km74°00′图 1    辐射和几何校正前的真彩色合成影像数成正比，每个探测器测得的电压通过模数转换器转式中：Sλ，i和bλ，i分别为第i个探测器在飞行前确定的换成DN（digital number）值，理想的CMOS传感器，在绝对辐射定标系数和偏置差.没有入射光时也应没有输出电平，但在实际使用时，适用于“京师一号”大量遥感数据的处理流程：由于热环境、探测器本身性能差异总会有或多或少1）自动筛选并统计大量0°太阳高度角影像，发现在全的输出信号[15].每波段的每1个探测器i的输出值为黑的条件下，影像有位置和大小相对固定的暗噪声，该像素入射辐射的响应转换比率与暗电流之和，初始DN值较大，不能忽略，可用于提取各传感器的偏置模型为值，其大小与探测器本身的性能和在轨环境相关，是条带与随机噪声产生的原因之一；2）通过对比大量有Si=Ri
Li+Ni;（1）效影像，发现同焦平面的同波段影像具有相近的明暗式中：Ri为传感器的灵敏度；Li为入射辐射值；Ni为暗空间分布趋势，且冰盖影像的真彩色合成图明显不呈电平；i为第i个传感器.经过放大、量化、模数转换白色，说明卫星发射后的在轨辐射情况和实验室定标等环节，将电信号转换为数字信号，表达式为结果存在差距，需在去除暗噪声后通过均匀地表影像Draw=Ds+Db;（2）重新计算增益值；3）经过目视筛选了天气良好、反射率分布均一的冰盖影像，经过掩膜、裁剪等操作，最式中：Draw为原始灰度值；Ds为输出电平转换成的灰终计算出各传感器的增益值.度值；Db为灰度偏置值；条带去除是通过改变、调节规范式中的增益系Ds=Gi
Si;（3）数Gi的大小和偏置值Db的大小，使CMOS传感器阵式中：Si为CMOS的输出电平；Gi为量化增益系数.列中各传感器具有新的增益系数和偏置值.当入射将式（2）、（3）整合后，得辐照度均匀一致时，得到相同的灰度值，地物的反射率得到真实的反映，对于推扫式影像，任意一行的第Draw=Gi
Si+Db;（4）i个像素的相对辐射校正模型为式中：Db为灰度偏置值.第i个探测器的原始DN值ci=(DiDb)
Ui;（6）最终会转换为第j个图像像素值和波长λ相关的大气层顶辐射值，这称为绝对辐射校正，模型为式中：ci为定标法校正后的DN值；Di为原始DN值；Db为暗噪声值即灰度偏置值；Ui为探元的归一化增D;i;j=G;i;j
S;i+b;i;（5）益值，均由实验获得.算法的具体处理流程如图2所

3第1期张卓宇等：“京师一号”极地小卫星宽幅影像辐射校正方法研究107示，流程主要包括：交集，目视筛选出没有过曝，少云、岩石少、没有海水1）统计全色波段直方图，自动筛选出直方图峰值的冰盖影像作为提取增益系数的参考影像.<5且横坐标最大值≠0的影像集；2次筛选出最大太3）将筛选后的冰盖影像每行减去暗噪声，采用单阳高度角为0°的RGB影像集，提取各焦平面各波段波段阈值法去除冰盖上的阴影、水体、裸岩等非冰雪的偏置值.像元，使用滑动窗口计算列中值后取其平均值，视作2）制作BNU的数据索引图，与极地冰盖边界取标准，提取增益系数.BNU-1L0数据1)筛选用于相对辐射2)小卫星相对辐射校正参考的BNU-1数据校正参数获取A：暗噪声参考影像筛选B：增益修正参考A：提取各传感器的B：提取各传感器的影像筛选偏置值增益系数绘制全色波段影像制作索引图与冰盖边统计列暗噪声将冰盖参考影像直方图界取交集影像众数减去偏置值目视筛选没有过曝,单波段阈值法去除水筛选直方图峰值＜5,存储偏置值云、岩石、水体少体、岩石和阴影并保Xmax＞0的影像的冰盖影像留条带噪声平滑计算列中值后筛选中心太阳高度角C：数据预处理求均值min影像裁剪有云、水体的计算各列中值与均值筛选暗噪声的参考影影像行的比值并存储像集3)小卫星相对辐射校正匹配影像和校正系数批量相对辐射校正使用定量评价指标评价图 2    适用于“京师一号”影像的自动相对辐射校正的流程2.2    辐射校正方法挑选θ最大值为0°的数据集，查找不同焦平面RGB2.2.1计算暗噪声波段的数据集.1）通过计算全色波段的直方图，首先筛选直方3）计算每景列均值后求各景列均值的平均值作图峰值≤5，且直方图横坐标最大值（全图灰度最大为暗噪声，公式为值）>0的影像，使用IDL语言的函数实现，公式为∑Nj1Hmax(I)≤5且Imax>0;（7）N;i=ND;i;j;（9）jj=1式中：I为图像的二维数组；H为直方图统计的函数，1∑nB;q=N;q;（10）可得到1个离散的一维数组；max为求最大值的函数，nq=1这一步筛选可以排除全图无效的情况.式中：Nλ，i为每景影像的列均值，λ为不同波段，i为每2）在不考虑天空光的影响下，根据太阳高度角个传感器的序列号；Nj为一景影像共有多少行，j为每和影像灰度值的经验式行的像素序列号；Bλ，q为待求的暗噪声，即前文的Oi，gi;j=fi;j
sin;（8）n为符合要求的影像数量；q为影像的序列号.需要注式中：gi，j为影像在太阳斜射时得到的DN值；fi，j为影意南北半球扫描成像方向相反，所以，南北半球的暗像在太阳直射时得到的DN值；θ为太阳高度角，2次噪声数据应分开计算.暗噪声影像及剖面见图3.

4108北京师范大学学报（自然科学版）第59卷a50b4030DN值2010010002000300040005000列号PMS2B2暗噪声剖面暗噪声影像图 3    暗噪声剖面（a）及影像（b）M;i=medC;i;j;（12）2.2.2计算增益系数1）利用ArcPy工具，根据“京师一号”的地理数1∑i+14R;i=M;q;（13）据，生成空间索引图，与南极的海岸线相交.进一步Nq=i14目视筛选出没有过度曝光、少云、岩石少、没有海水式中：Mλ，i代表列中值的结果；med代表中值函数；的冰盖影像作为提取增益系数的参考影像.计算参q代表滑动的窗口位置.考影像去除暗噪声的二维数组，公式为4）计算滑动后的特征向量的平均值，并将其视C;i;j=D;i;jB;q;（11）为标准，最终计算每个传感器的增益系数，公式为式中：Cλ，i，j为去除暗噪声的灰度值；Dλ，i，j为原灰∑i1度值.S=R;i;（14）ip=12）因小卫星的波段较少，且地理产品精度有限，故利用经验值，掩膜去除冰盖影像上为水体、裸岩、SU;i=;（15）阴影的像元，同时尽可能保留条带噪声，并手动裁掉M;i云像元，效果如图4所示.式中：Sλ代表特征向量的平均值的结果；Uλ，i代表每3）计算掩膜后的列中值，并使用大小为29的滑个传感器的增益系数.图5为PMS1焦面辐射校正效动窗口平滑计算列中值，公式为果对比，图6为PMS2焦面辐射校正对比.abca. PMS2焦面B1的原始影像；b. 去除原始影像中不均匀的少量冰面湖、阴影和岩石，提取增益系数的过程；c. 校正系数作用于原图后的整体效果.图 4    辐射校正过程2.3    辐射校正结果精度评价精度验证使用列均值较大变化，侧面反映出光学镜头的制造与打磨，以及分布图，并在对应同期MODIS影像中均匀取1 000个焦平面的设计方案对成像效果造成较大影响.由图8点对值，最终以线性回归方程和线性回归拟合优度可知，BNU-1和MODIS影像之间存在明显的线性关R2作为定量评价指标.系，证明辐射校正的效果的真实性与可靠性；结合表3由图7可知，“京师一号”辐射校正后的数值在影可知，序号1和4的拟合效果最佳，拟合优度R2均值像中部位置基本保持不变，随着向镜头边缘靠近，影分别为0.79和0.87，序号2和3的拟合效果较低，拟像的清晰程度受到较大影响，影像的列均值也产生了合优度R2均值分别为0.55和0.36，也反映出基于均

5第1期张卓宇等：“京师一号”极地小卫星宽幅影像辐射校正方法研究109a辐射校正前PMS1:Bl影像c辐射校正前PMS1:B2影像e辐射校正前PMS1:B3影像b辐射校正后PMS1:Bl影像d辐射校正后PMS1:B2影像f辐射校正后PMS1:B3影像图 5    PMS1焦面辐射校正效果对比a辐射校正前PMS2:Bl影像c辐射校正前PMS2:B2影像e辐射校正前PMS2:B3影像b辐射校正后PMS2:Bl影像d辐射校正后PMS2:B2影像f辐射校正后PMS2:B3影像图 6    PMS2焦面辐射校正效果对比匀场景的辐射校正方法在不同的地域，产生的效果也进，通过SIFT算子将“京师一号”影像与同期MODIS不尽相同.影像配准. MODIS数据在利用NASA研发的HEG工具完成数据格式转换后在ArcGIS中投影为极投影平3    结果和讨论面坐标系.“京师一号”卫星的几何定位精度仍需要人工改图9中左上角的红星表示试验区域，红框中为

6110北京师范大学学报（自然科学版）第59卷2400表 3    BNU-1和MODIS拟合方程和拟合优度（R2）R-L02200R-DestripeG-L0G-Destripe序号波段回归方程R22000B-L0B-Destripe红y=21.21x+847.5770.7918001绿y=20.57x+310.520.79平均DN1600蓝y=17.99x+561.040.79红y=7.45x+427.280.5314002绿y=7.53x+958.140.531200010002000300040005000蓝y=8.83x+1092.830.58列号红y=12.08x−620.9050.38注：来源于BNU1_PMS2_20191002041454_200220500_102_0001_001影像．3绿y=10.36x+730.670.34图 7    “京师一号”RGB波段列均值分布蓝y=10.03x+1039.520.37红y=12.92x−615.390.87红4绿y=13.43x+99.840.878000蓝y=15.32x+348.920.896000相对辐射校正系数，根据各纠正系数文件名和各焦平MODIS4000面的各波段文件名进行匹配，进行批量相对辐射校正.从结果来看，适用于国产极地小卫星影像的自动2000相对辐射纠正技术在以下几个方面具有明显优势：01502002503003504004505005501）该技术能够实现自动化筛选暗噪声参考影像，绿避免因为卫星数据量大，又要兼备影像空间范围、拍8000摄时间、数据质量等方面带来的人工干预.70002）该技术利用极地均匀雪面影像进行了相对辐6000射校正系数的提取，使得在多景非均匀雪面的影像上5000MODIS4000也能够得到有效的校正效果.从原始真彩色合成图3000像可以看出，0级产品的海冰分布红、黄绿色，而校正2000后的海冰已经没有了偏色的效果；从冰盖影像辐射校1000150200250300350400450500正前后对比效果来看，焦面PMS2的0级影像在真彩蓝色合成后，会在大范围雪面上形成红、绿背景，而校7000正后的影像已经没有了偏色效果，证明该算法的可行性.600050004    结论40001）针对“京师一号”小卫星的0级数据辐射性能MODIS3000问题导致的条纹现象提出了一种去噪声算法，首先进2000行暗噪声的处理，其次使用均匀场景掩膜法在高反射的纯净冰盖上提取出每个像元探测器的增益系数.150200250300350400450500BNU-12）辐射校正效果在南极得到了广泛应用，得到注：影像来源于BNU1_PMS2_20191117201940_了较好的目视辐射校正效果.列均值的图片表明，影200018633_101_0001_001．像在中间部位色差不大，而在边缘位置数值变化较图 8    “京师一号”和MODIS散点大；0级影像中的大范围泛红、泛绿和在海冰等地区“京师一号”经过辐射校正和几何校正后的效果，底发黄的问题，均得到了一定程度的改善，通过和图为MODIS同期真彩色合成影像.对比图1和图9MODIS同期影像对比，二者的拟合优度在海岸影像可知经过本文的辐射校正后，影像的辐射得到了有效较高（>0.53），而在冰盖上的拟合优度较低，最小仅提升.“京师一号”影像利用步骤2中的方案得出的为0.37，证明该算法在不同区域存在差异.

7第1期张卓宇等：“京师一号”极地小卫星宽幅影像辐射校正方法研究11162°64°66°E180°175°170°165°160°155°150°ES70°00′S78°00′70°30′71°00′76°00′02550100150200km02550100150200km71°30′74°00′160°170°180°W170°160°150°140°135°130°125°E175°180°W175°170°165°160°155°150°ES78°00′S84°00′76°00′82°00′74°00′02550100150200km02550100150200km图 9    辐射和几何校正后真彩色合成影像3）在南极的经典区域，冰盖、冰架、海岸等区域，images by histogram modification[J].Computer Graphics均得到了较好的辐射校正效果.但其中海岸影像的Image Processing，1979，10（1）：69边缘没有对齐，且存在黄色、红色、蓝色的窄边，且经［8］NICHOL J E，VOHORA V. Noise over water surfaces in过校正前的影像形状均不规则，而经过SIFT算法与Landsat TM images[J].International Journal of RemoteMODIS配准后，影像呈规则的方形，且边缘的彩色边Sensing，2004，25（11）：2087缘也已去除，证明影像的几何问题也会影响辐射校正［9］张兵，张浩，陈正超，等.一种基于图像统计量的相对辐射纠正算法[J].遥感学报，2006，10（5）：630的效果.［10］皮原征，储栋，管小彬，等.顾及水陆差异的高分五号影5　参考文献像条带去除[J].遥感学报，2020，24（4）：360［11］INFANTE，OMAR S. Wavelet analysis for the elimination［1］袁勤，张帝，何志鹏，等.微景星座首发星京师一号技术特of striping noise in satellite images[J].Optical点分析[J].国际太空，2019（10）：4Engineering，2001，40（7）：1309［2］鄂栋臣，沈强.南极冰貌卫星影像条带噪声的剔除[J].极［12］MARKHAM B L，BONCYK W C，HELDER D L，et al.地研究，2004，16（1）：1Landsat-7 Enhanced Thematic Mapper Plus：radiometric［3］段依妮，张立福，晏磊，等.遥感影像相对辐射校正方法及calibration and prelaunch performance[J].Canadian适用性研究[J].遥感学报，2014，18（3）：597Journal of Remote Sensing，1997，23（4）：318［4］王淑荣，邢进，李福田.利用积分球光源定标空间紫外遥感光谱辐射计[J].光学精密工程，2006，14（2）：185［13］陈立奇.南极和北极地区在全球变化中的作用研究[J].地学前缘，2002，9（2）：245［5］李德仁，王密.“资源三号”卫星在轨几何定标及精度评估[J].航天返回与遥感，2012，33（3）：1［14］甘信铮，孙家抦.南极卫星影像条带噪声的消除[J].武汉测绘科技大学学报，1994，19（4）：332［6］刘银年，孙德新，曹开钦，等.高分五号可见短波红外高光谱相机在轨辐射性能评估[J].遥感学报，2020，24（4）：352［15］高正清，杨志高，王险峰.相对辐射定标与相对辐射校正［7］HORN B K P，WOODHAM R J. Destriping Landsat MSS场[J].影像技术，2009，21（4）：48

8112北京师范大学学报（自然科学版）第59卷Radiometric correction of BNU-1 wide-swath imageryZHANG ZhuoyuZHANG Ying（ State Key Laboratory of Remote Sensing Science，College of Global Change and Earth System Science，Beijing Normal University，100875，Beijing，China）AbstractBNU-1，the first civil polar optical micro-sat developed by Chinese scientists，is equipped with amedium and high-resolution multi-spectral wide-swath camera，and takes into account both revisit cycle andresolution. BNU-1 is of great significance for the study of climate and environment change and monitoring in polarregions，at high coverage and medium resolution. However，the band noise of level 0 data masks the ground feature ofimage and does not reflect the real radiation characteristics of ground features，seriously affecting data application.Therefore，in this paper a method is proposed to remove the strip based on statistical principle for wide-swath image.The algorithm uses the dark noise extraction method corresponding to the image of 0° solar altitude angle and theuniform scene mask method to extract offset and gain values of each detector respectively，and realizes the stripremoval. Experimental data indicate that compared to original data，this algorithm robustly removes band noise，restores real surface radiation information，and gives full play to the application value of multispectral data of BNU-1.KeywordsBNU-1；the Antarctic；strip noise；radiometric correction【责任编辑：刘先勤】