CN114663530A - 一种遥感影像偏航90度的相对辐射定标方法和装置 - Google Patents
一种遥感影像偏航90度的相对辐射定标方法和装置 Download PDFInfo
- Publication number
- CN114663530A CN114663530A CN202210301439.2A CN202210301439A CN114663530A CN 114663530 A CN114663530 A CN 114663530A CN 202210301439 A CN202210301439 A CN 202210301439A CN 114663530 A CN114663530 A CN 114663530A
- Authority
- CN
- China
- Prior art keywords
- yaw data
- ccd
- probe
- sub
- original
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 239000000523 sample Substances 0.000 claims abstract description 81
- 238000001914 filtration Methods 0.000 claims abstract description 27
- 238000012545 processing Methods 0.000 claims abstract description 20
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 16
- 230000005855 radiation Effects 0.000 claims description 18
- 238000003491 array Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 6
- 201000009482 yaws Diseases 0.000 abstract description 2
- 238000004891 communication Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
-
- G06T5/70—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/90—Determination of colour characteristics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Abstract
本发明提供了一种遥感影像偏航90度的相对辐射定标方法和装置,涉及图像处理的技术领域,包括:获取原始偏航数据,并将原始偏航数据旋转45度,得到目标偏航数据,其中,原始偏航数据为CCD线阵90度偏航时CCD线阵中的CCD探元采集到的地物遥感影像;对目标偏航数据进行Savitzky‑Golay滤波处理,得到CCD探元的理论灰度值;基于最小二乘算法、CCD探元的理论灰度值和目标偏航数据,计算CCD探元的探元参数,其中,探元参数包括:增益参数和偏置参数;基于探元参数,对原始偏航数据进行相对辐射定标,得到所述原始偏航数据的相对辐射定标结果,解决了现有的遥感偏航影像的相对辐射定标方法的准确性较低的技术问题。
Description
技术领域
本发明涉及图像处理的技术领域,尤其是涉及一种遥感影像偏航90度的相对辐射定标方法和装置。
背景技术
现有的遥感影像定标方式如下:
1.在轨星上定标。在轨星上定标精度高但操作易受限制,且难以校正传感器因发射震动和在轨响应状态衰减等原因产生的辐射响应差异,无法满足卫星整个运行周期的高精度辐射定标。
2.实验室积分球定标。在遥感器发射之前对其进行的波长位置、辐射精度、空间定位等的定标,将仪器的输出值转换为辐射值。有的仪器内有内定定标系统。但是在仪器运行之后,还需要定期定标,以监测仪器性能的变化,相应调整定标参数。
3.均匀场地地面定标人工地面定标场需要大量人力财力维护,易受天气、轨道等环境的条件影像。基于自然地物的沙漠、水体等均匀场地的在轨场地定标和在轨统计定标均存在数据符合数学统计特性的假设,且需要海量样本或均匀场数据,在卫星入轨初期难以满足高频次的定标需求,在统计数据量不足时,会因地物差异大而产生校正过度或校正不完全的情况,均匀性校正效果通常不理想。
遥感卫星在轨控制技术的提高,使得基于90°偏航数据的相对辐射定标技术得到广泛应用,如珠海一号、Landsat 8、QuickBird等。常规的偏航数据相对辐射定标方法主要是将各个探元校正到所有探元平均响应水平,但影像中的亮暗条纹在有时候未必服从正态分布,因此计算定标参数会不准确。
针对上述问题,还未提出有效的解决方案。
发明内容
有鉴于此,本发明的目的在于提供一种遥感影像偏航90度的相对辐射定标方法和装置,以缓解了现有的遥感影像相对辐射定标方法的准确性较低的技术问题。
第一方面,本发明实施例提供了一种遥感影像偏航90度的相对辐射定标方法,包括:获取原始偏航数据,并将所述原始偏航数据旋转45度,得到目标偏航数据,其中,所述原始偏航数据为CCD线阵90度偏航时所述CCD线阵中的CCD探元采集到的地物遥感影像;对所述目标偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的理论灰度值;基于最小二乘算法、所述CCD探元的理论灰度值和所述目标偏航数据,计算所述CCD探元的探元参数,其中,所述探元参数包括:增益参数和偏置参数;基于所述探元参数,对所述原始偏航数据进行相对辐射定标,得到所述原始偏航数据的相对辐射定标结果。
进一步地,所述原始偏航数据包含不同波段的子偏航数据。
进一步地,对所述目标偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的理论灰度值,包括:对所述子偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的子理论灰度值;基于不同波段的子偏航数据对应的CCD探元的子理论灰度值,得到所述CCD探元的理论灰度值。
进一步地,基于最小二乘算法、所述CCD探元的理论灰度值和所述目标偏航数据,计算所述CCD探元的探元参数,包括:基于所述最小二乘算法、所述子理论灰度值和所述目标偏航数据,计算所述CCD探元的子探元参数;基于不同波段对应的子探元参数,得到所述CCD探元的探元参数。
进一步地,基于所述CCD探元的探元参数,对所述原始偏航数据进行相对辐射定标,得到所述原始偏航数据的相对辐射定标结果,包括:对所述CCD探元的子探元参数,对所述子探元参数对应的子偏航数据进行相对辐射标定,得到子相对辐射标定结果;基于所述子相对辐射标定结果,得到所述原始偏航数据的相对辐射定标结果。
第二方面,本发明实施例还提供了一种遥感影像偏航90度的相对辐射定标装置,包括:获取单元,滤波单元,计算单元和标定单元,其中,所述获取单元,用于获取原始偏航数据,并将所述原始偏航数据旋转45度,得到目标偏航数据,其中,所述原始偏航数据为CCD线阵90度偏航时所述CCD线阵中的CCD探元采集到的地物遥感影像;所述滤波单元,用于对所述目标偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的理论灰度值;所述计算单元,用于基于最小二乘算法、所述CCD探元的理论灰度值和所述目标偏航数据,计算所述CCD探元的探元参数,其中,所述探元参数包括:增益参数和偏置参数;所述标定单元,用于基于所述探元参数,对所述原始偏航数据进行相对辐射定标,得到所述原始偏航数据的相对辐射定标结果。
进一步地,所述原始偏航数据包含不同波段的子偏航数据。
进一步地,所述滤波单元,用于:对所述子偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的子理论灰度值;基于不同波段的子偏航数据对应的CCD探元的子理论灰度值,得到所述CCD探元的理论灰度值。
第三方面,本发明实施例还提供了一种电子设备,包括存储器以及处理器,所述存储器用于存储支持处理器执行上述第一方面中所述方法的程序,所述处理器被配置为用于执行所述存储器中存储的程序。
第四方面,本发明实施例还提供了一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序。
在本发明实施例中,通过获取原始偏航数据,并将原始偏航数据旋转45度,得到目标偏航数据,其中,原始偏航数据为CCD线阵90度偏航时CCD线阵中的CCD探元采集到的地物遥感影像;对目标偏航数据进行Savitzky-Golay滤波处理,得到CCD探元的理论灰度值;基于最小二乘算法、CCD探元的理论灰度值和目标偏航数据,计算CCD探元的探元参数,其中,探元参数包括:增益参数和偏置参数;基于探元参数,对原始偏航数据进行相对辐射定标,得到所述原始偏航数据的相对辐射定标结果,达到了不依赖地面均匀地物或均匀定标场对遥感偏航影像进行相对定标的目的,进而解决了现有的遥感偏航影像的相对辐射定标方法的准确性较低的技术问题,从而实现了提高遥感偏航影像的相对辐射定标方法的准确性的技术效果。
本发明的其他特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的种遥感影像偏航90度的相对辐射定标方法的流程图;
图2为本发明实施例提供的种遥感影像偏航90度的相对辐射定标装置的示意图;
图3为本发明实施例提供的一种电子设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一:
根据本发明实施例,提供了一种遥感影像偏航90度的相对辐射定标方法的实施例,需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
图1是根据本发明实施例的一种遥感影像偏航90度的相对辐射定标方法的流程图,如图1所示,该方法包括如下步骤:
步骤S102,获取原始偏航数据,并将所述原始偏航数据旋转45度,得到目标偏航数据,其中,所述原始偏航数据为CCD线阵90度偏航时所述CCD线阵中的CCD探元采集到的地物遥感影像;
需要说明的是,所述原始偏航数据包含不同波段的子偏航数据。
步骤S104,对所述目标偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的理论灰度值;
步骤S106,基于最小二乘算法、所述CCD探元的理论灰度值和所述目标偏航数据,计算所述CCD探元的探元参数,其中,所述探元参数包括:增益参数和偏置参数;
步骤S108,基于所述探元参数,对所述原始偏航数据进行相对辐射定标,得到所述原始偏航数据的相对辐射定标结果。
在本发明实施例中,通过获取原始偏航数据,并将原始偏航数据旋转45度,得到目标偏航数据,其中,原始偏航数据为CCD线阵90度偏航时CCD线阵中的CCD探元采集到的地物遥感影像;对目标偏航数据进行Savitzky-Golay滤波处理,得到CCD探元的理论灰度值;基于最小二乘算法、CCD探元的理论灰度值和目标偏航数据,计算CCD探元的探元参数,其中,探元参数包括:增益参数和偏置参数;基于探元参数,对原始偏航数据进行相对辐射定标,得到所述原始偏航数据的相对辐射定标结果,达到了不依赖地面均匀地物或均匀定标场对遥感偏航影像进行相对定标的目的,进而解决了现有的遥感偏航影像的相对辐射定标方法的准确性较低的技术问题,从而实现了提高遥感偏航影像的相对辐射定标方法的准确性的技术效果。
在本发明实施例中,步骤S104包括如下步骤:
步骤S11,对所述子偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的子理论灰度值;
步骤S12,基于不同波段的子偏航数据对应的CCD探元的子理论灰度值,得到所述CCD探元的理论灰度值。
在本发明实施例中,获取CCD探元的子理论灰度值的方式如下:
设每个波段子偏航数据的行方向像素值为xi={x1,x2,…,xn}(n为列数),若考虑一组以n为中心的2M+1个数据,可以用k阶多项式表示:
式中,k表示多项式的阶次,ak表示多项式系数,p(n)表示拟合后的值。
最小二乘拟合的残差为:
当n<M时,参与计算的点只有x0,x1,…,xM,当n>M时需要计算x0,…,x(-M),同样,n>n-M时,需要计算n+M个数据点的值,计算公式如下:
xn+m=xn-m。
下面说一下如何获得ak。
由基本的微积分知识可知,若要ε对各个参数的偏导数都应为0,也就是:
化简后可得:
i=0,1,…N。因此,Savitzky-Golay滤波卷积平滑法是移动平滑算法的改进:
需要说明的是,Savitzky-Golay滤波需要根据滤波的程度输入多项式拟合阶次和窗口的点数,实验表明多项式拟合阶次小于2,点数大于100的效果较好。
在本发明施例中,步骤S106包括如下步骤:
步骤S21,基于所述最小二乘算法、所述子理论灰度值和所述目标偏航数据,计算所述CCD探元的子探元参数;
步骤S22,基于不同波段对应的子探元参数,得到所述CCD探元的探元参数。
在本发明实施例中,计算CCD探元的子探元参数的过程如下:
最小二乘法原理是给定的二维数据点Pi(xi,yi),根据公式计算拟合线y=ax+b与给定数据点的拟合偏差并使偏差的平方达到最小,拟合线在点Pi处的偏差δj=axj+b-yi,其中j=1,2,3,…,m。依据以下公式,计算拟合直线y=ax+b中斜率和截距参数。
对理论灰度值与目标偏航数据采用最小二乘法拟合计算每个CCD探元的增益参数和偏置参数。
在本发明实施例中,步骤S108包括如下步骤:
步骤S31,对所述CCD探元的子探元参数,对所述子探元参数对应的子偏航数据进行相对辐射标定,得到子相对辐射标定结果;
步骤S32,基于所述子相对辐射标定结果,得到所述原始偏航数据的相对辐射定标结果。
在本发明实施例中,确定子相对辐射标定结果如下:
根据计算所得各个CCD探元的增益和偏置两项辐射定标参数依据以下公式对影像进行偏航相对辐射定标:
y=ax+b。
其中x表示原始偏航数据的像元值,y为定标后的像元值。
由于不同波段的定标系数不一样,按照上述方法对原始偏航数据中所有波段进行逐一定标参数的获取,再进行各个波段辐射定标处理,获得原始偏航数据中所有波段相对辐射定标的结果。
实施例二:
本发明实施例还提供了一种遥感影像偏航90度的相对辐射定标装置,该遥感影像偏航90度的相对辐射定标装置用于执行本发明实施例上述内容所提供的遥感影像偏航90度的相对辐射定标方法,以下是本发明实施例提供的遥感影像偏航90度的相对辐射定标装置的具体介绍。
如图2所示,图2为上述遥感影像偏航90度的相对辐射定标装置的示意图,该遥感影像偏航90度的相对辐射定标装置包括:获取单元10,滤波单元20,计算单元30和标定单元40。
所述获取单元10,用于获取原始偏航数据,并将所述原始偏航数据旋转45度,得到目标偏航数据,其中,所述原始偏航数据为CCD线阵90度偏航时所述CCD线阵中的CCD探元采集到的地物遥感影像;
所述滤波单元20,用于对所述目标偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的理论灰度值;
所述计算单元30,用于基于最小二乘算法、所述CCD探元的理论灰度值和所述目标偏航数据,计算所述CCD探元的探元参数,其中,所述探元参数包括:增益参数和偏置参数;
所述标定单元40,用于基于所述探元参数,对所述原始偏航数据进行相对辐射定标,得到所述原始偏航数据的相对辐射定标结果。
在本发明实施例中,通过获取原始偏航数据,并将原始偏航数据旋转45度,得到目标偏航数据,其中,原始偏航数据为CCD线阵90度偏航时CCD线阵中的CCD探元采集到的地物遥感影像;对目标偏航数据进行Savitzky-Golay滤波处理,得到CCD探元的理论灰度值;基于最小二乘算法、CCD探元的理论灰度值和目标偏航数据,计算CCD探元的探元参数,其中,探元参数包括:增益参数和偏置参数;基于探元参数,对原始偏航数据进行相对辐射定标,得到所述原始偏航数据的相对辐射定标结果,达到了不依赖地面均匀地物或均匀定标场对遥感偏航影像进行相对定标的目的,进而解决了现有的遥感偏航影像的相对辐射定标方法的准确性较低的技术问题,从而实现了提高遥感偏航影像的相对辐射定标方法的准确性的技术效果。
实施例三:
本发明实施例还提供了一种电子设备,包括存储器以及处理器,所述存储器用于存储支持处理器执行上述实施例一中所述方法的程序,所述处理器被配置为用于执行所述存储器中存储的程序。
参见图3,本发明实施例还提供一种电子设备100,包括:处理器50,存储器51,总线52和通信接口53,所述处理器50、通信接口53和存储器51通过总线52连接;处理器50用于执行存储器51中存储的可执行模块,例如计算机程序。
其中,存储器51可能包含高速随机存取存储器(RAM,Random Access Memory),也可能还包括非不稳定的存储器(non-volatile memory),例如至少一个磁盘存储器。通过至少一个通信接口53(可以是有线或者无线)实现该系统网元与至少一个其他网元之间的通信连接,可以使用互联网,广域网,本地网,城域网等。
总线52可以是ISA总线、PCI总线或EISA总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示,图3中仅用一个双向箭头表示,但并不表示仅有一根总线或一种类型的总线。
其中,存储器51用于存储程序,所述处理器50在接收到执行指令后,执行所述程序,前述本发明实施例任一实施例揭示的流过程定义的装置所执行的方法可以应用于处理器50中,或者由处理器50实现。
处理器50可能是一种集成电路芯片,具有信号的处理能力。在实现过程中,上述方法的各步骤可以通过处理器50中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器50可以是通用处理器,包括中央处理器(Central Processing Unit,简称CPU)、网络处理器(Network Processor,简称NP)等;还可以是数字信号处理器(Digital SignalProcessing,简称DSP)、专用集成电路(Application Specific Integrated Circuit,简称ASIC)、现成可编程门阵列(Field-Programmable Gate Array,简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成,或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器,闪存、只读存储器,可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器51,处理器50读取存储器51中的信息,结合其硬件完成上述方法的步骤。
实施例四:
本发明实施例还提供了一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,计算机程序被处理器运行时执行上述实施例一中所述方法的步骤。
另外,在本发明实施例的描述中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本申请所提供的几个实施例中,应该理解到,所揭露的系统、装置和方法,可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,又例如,多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种遥感影像偏航90度的相对辐射定标方法,其特征在于,包括:
获取原始偏航数据,并将所述原始偏航数据旋转45度,得到目标偏航数据,其中,所述原始偏航数据为CCD线阵90度偏航时所述CCD线阵中的CCD探元采集到的地物遥感影像;
对所述目标偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的理论灰度值;
基于最小二乘算法、所述CCD探元的理论灰度值和所述目标偏航数据,计算所述CCD探元的探元参数,其中,所述探元参数包括:增益参数和偏置参数;
基于所述探元参数,对所述原始偏航数据进行相对辐射定标,得到所述原始偏航数据的相对辐射定标结果。
2.根据权利要求1所述的方法,其特征在于,所述原始偏航数据包含不同波段的子偏航数据。
3.根据权利要求2所述的方法,其特征在于,对所述目标偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的理论灰度值,包括:
对所述子偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的子理论灰度值;
基于不同波段的子偏航数据对应的CCD探元的子理论灰度值,得到所述CCD探元的理论灰度值。
4.根据权利要求3所述的方法,其特征在于,基于最小二乘算法、所述CCD探元的理论灰度值和所述目标偏航数据,计算所述CCD探元的探元参数,包括:
基于所述最小二乘算法、所述子理论灰度值和所述目标偏航数据,计算所述CCD探元的子探元参数;
基于不同波段对应的子探元参数,得到所述CCD探元的探元参数。
5.根据权利要求4所述的方法,其特征在于,基于所述CCD探元的探元参数,对所述原始偏航数据进行相对辐射定标,得到所述原始偏航数据的相对辐射定标结果,包括:
对所述CCD探元的子探元参数,对所述子探元参数对应的子偏航数据进行相对辐射标定,得到子相对辐射标定结果;
基于所述子相对辐射标定结果,得到所述原始偏航数据的相对辐射定标结果。
6.一种遥感影像偏航90度的相对辐射定标装置,其特征在于,包括:获取单元,滤波单元,计算单元和标定单元,其中,
所述获取单元,用于获取原始偏航数据,并将所述原始偏航数据旋转45度,得到目标偏航数据,其中,所述原始偏航数据为CCD线阵90度偏航时所述CCD线阵中的CCD探元采集到的地物遥感影像;
所述滤波单元,用于对所述目标偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的理论灰度值;
所述计算单元,用于基于最小二乘算法、所述CCD探元的理论灰度值和所述目标偏航数据,计算所述CCD探元的探元参数,其中,所述探元参数包括:增益参数和偏置参数;
所述标定单元,用于基于所述探元参数,对所述原始偏航数据进行相对辐射定标,得到所述原始偏航数据的相对辐射定标结果。
7.根据权利要求6所述的装置,其特征在于,所述原始偏航数据包含不同波段的子偏航数据。
8.根据权利要求7所述的装置,其特征在于,所述滤波单元,用于:
对所述子偏航数据进行Savitzky-Golay滤波处理,得到所述CCD探元的子理论灰度值;
基于不同波段的子偏航数据对应的CCD探元的子理论灰度值,得到所述CCD探元的理论灰度值。
9.一种电子设备,其特征在于,包括存储器以及处理器,所述存储器用于存储支持处理器执行权利要求1至5任一项所述方法的程序,所述处理器被配置为用于执行所述存储器中存储的程序。
10.一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,其特征在于,计算机程序被处理器运行时执行上述权利要求1至5任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210301439.2A CN114663530A (zh) | 2022-03-24 | 2022-03-24 | 一种遥感影像偏航90度的相对辐射定标方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210301439.2A CN114663530A (zh) | 2022-03-24 | 2022-03-24 | 一种遥感影像偏航90度的相对辐射定标方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114663530A true CN114663530A (zh) | 2022-06-24 |
Family
ID=82030706
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210301439.2A Pending CN114663530A (zh) | 2022-03-24 | 2022-03-24 | 一种遥感影像偏航90度的相对辐射定标方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114663530A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115830146A (zh) * | 2023-02-10 | 2023-03-21 | 武汉玄景科技有限公司 | 一种航天光学遥感相机的在轨相对辐射定标与校正方法 |
-
2022
- 2022-03-24 CN CN202210301439.2A patent/CN114663530A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115830146A (zh) * | 2023-02-10 | 2023-03-21 | 武汉玄景科技有限公司 | 一种航天光学遥感相机的在轨相对辐射定标与校正方法 |
CN115830146B (zh) * | 2023-02-10 | 2023-05-09 | 武汉玄景科技有限公司 | 一种航天光学遥感相机的在轨相对辐射定标与校正方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109522516B (zh) | 基于随机森林回归算法的土壤湿度检测方法、装置及电子设备 | |
CN106887016B (zh) | 一种gf-4卫星序列图像自动相对配准方法 | |
CN103226819A (zh) | 一种基于分段统计的相对辐射校正方法 | |
CN112927306B (zh) | 拍摄装置的标定方法、装置及终端设备 | |
CN113450353B (zh) | 一种叶面积指数精度的优化方法和装置 | |
CN116046692B (zh) | 一种基于高光谱的土壤重金属污染监测方法和装置 | |
CN111508011A (zh) | 一种飞行时间相机的深度数据校准方法 | |
CN114663530A (zh) | 一种遥感影像偏航90度的相对辐射定标方法和装置 | |
CN114241064B (zh) | 一种遥感卫星内外方位元素实时几何定标方法 | |
CN108983172B (zh) | 一种随机进化约束的sar辐射定标方法 | |
CN112345090A (zh) | 点目标红外辐射测量的方法、设备及可读存储介质 | |
CN113077384B (zh) | 一种数据空间分辨率提高方法、装置、介质及终端设备 | |
Lindsley et al. | Adapting the SIR algorithm to ASCAT | |
CN107220947B (zh) | 一种遥感图像相对辐射校正方法及系统 | |
CN113553549B (zh) | 一种植被覆盖度反演方法、装置、电子设备及存储介质 | |
CN116430069A (zh) | 机器视觉流体流速测量方法、装置、计算机设备及存储介质 | |
CN114820376A (zh) | 条带噪声的融合校正方法、装置、电子设备及存储介质 | |
CN116379935A (zh) | 基于多光谱影像的岛礁浅海水深预测方法、装置以及设备 | |
CN111246190B (zh) | 基于核辐射定位仪的图像处理方法、装置及电子设备 | |
CN113409366B (zh) | 一种遥感卫星影像的辐射校正方法及装置 | |
CN115235625A (zh) | 一种光谱共焦回波位置标定方法、装置、介质及计算机设备 | |
CN110930344B (zh) | 一种目标的质量确定方法、装置、系统及电子设备 | |
CN111751006A (zh) | 一种数据标定方法、目标辐射特性分析方法和装置 | |
CN111751286A (zh) | 基于变化检测算法的土壤水分提取方法 | |
Tilton et al. | On-orbit measurement of the focal length of the SNPP VIIRS instrument |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |