CN114046888B - 基于卷积神经网络的波束合成推扫辐射计定标方法 - Google Patents
基于卷积神经网络的波束合成推扫辐射计定标方法 Download PDFInfo
- Publication number
- CN114046888B CN114046888B CN202111274687.4A CN202111274687A CN114046888B CN 114046888 B CN114046888 B CN 114046888B CN 202111274687 A CN202111274687 A CN 202111274687A CN 114046888 B CN114046888 B CN 114046888B
- Authority
- CN
- China
- Prior art keywords
- receiving
- radiometer
- broom
- push
- neural network
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000013527 convolutional neural network Methods 0.000 title claims abstract description 47
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 41
- 238000003786 synthesis reaction Methods 0.000 title claims abstract description 41
- 244000007853 Sarothamnus scoparius Species 0.000 title claims description 20
- 230000005855 radiation Effects 0.000 claims abstract description 22
- 238000012549 training Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 21
- 230000001427 coherent effect Effects 0.000 claims description 11
- 230000002068 genetic effect Effects 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 6
- 230000001537 neural effect Effects 0.000 abstract description 2
- 238000011176 pooling Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 210000002569 neuron Anatomy 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000013528 artificial neural network Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/52—Radiation pyrometry, e.g. infrared or optical thermometry using comparison with reference sources, e.g. disappearing-filament pyrometer
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明涉及基于卷积神经网络的波束合成推扫辐射计定标方法,该方法基于卷积神经网络的全链路定标方法,结合其他步骤得到的先验信息,再根据推扫辐射计周期性的对已知微波辐射亮温信息的定标场进行观测,利用定标场的亮温信息以及推扫辐射计系统输出的功率信号,对卷积神经网络的模型参数进行循环、迭代,不断训练,直至搜索得到满足误差阈值时的最优模型参数。推扫辐射计工作时,通过对微波辐射亮温已知的定标场的观测,获取定标数据样本,可以周期性的进行卷积神经网格模型参数训练。
Description
技术领域
本发明涉及基于卷积神经网络的波束合成推扫辐射计定标方法,属于空间微波遥感技术领域。
背景技术
不同于实孔径辐射计和综合孔径辐射计体制,波束合成推扫辐射计系统在数字域进行波束合成,以此实现天线电性能的超高波束效率,避免机械扫描与大口径天线之间的矛盾,能够用于弥补近岸高精度数据的空白。波束合成推扫辐射计系统的星载应用尚属空白,且定标技术也无公开文献提出。波束合成过程中的加权系数的获取直接关系到波束合成后的推扫辐射计系统性能。初步得到的幅相加权系数是通过测试得到的馈源次级天线方向图数据、耦合相干噪声得到接收通道的幅度和相位不一致性,再利用遗传算法和序列二次规划算法进行最优值搜索得到的,然而系统在星载或机载应用中,馈源次级天线方向图和接收通道的幅相特性会发生变化,如果沿用初始的幅相加权系数值,会导致波束合成后的波束天线方向图性能(主要包括主波束效率、旁瓣性能)不能满足指标要求。
目前,对于多个接收通道之间的幅度和相位不一致性定标方法主要有三种:(1)在辐射计系统外部远场条件下放置噪声源,对所有接收链路的输出信号进行两两进行复相关运算输出复相关矩阵,结合外部噪声源相对于辐射计系统的位置,对相关矩阵进行最小二乘求解,得到所有接收通道之间的幅相不一致性。(2)通过在接收通道输入端注入相同的相干噪声信号,同样是对两两接收通道进行复相关运算,并相关矩阵进行最小二乘求解得到幅相不一致性。(3)对接收通道注入相干噪声信号,通过自适应滤波器结合LMS算法,以均方误差最小为判定准则,不断循环迭代,对波束合成的幅相加权系数进行更新。不足之处:方法(1)~(3)对相干噪声信号的功率要求高,所有接收通道两两复相关对数据处理资源、运算能力提出了极高的要求。为了获取高主波束效率以及大幅宽观测,波束合成推扫辐射计系统的接收通道数量庞大,通常需要上千路接收通道进行实现,而辐射计系统的功耗、运算资源、存储空间等高要求都面临巨大的挑战,系统应用受到限制。
发明内容
本发明的技术解决问题是:克服上述现有技术的不足,提供基于卷积神经网络的波束合成推扫辐射计定标方法,提升波束合成推扫辐射计系统的探测精度,为高精度的亮温获取提供了方法。
本发明的技术解决方案为:基于卷积神经网络的波束合成推扫辐射计定标方法,该方法步骤如下:
S1、测量得到推扫辐射计的环焦反射面结合密集馈源阵列条件下的所有接收链路的馈源次级天线方向图Ffeed(i),i=1~M,M为推扫辐射计接收链路数目;
S2、以期望得到的各波束合成后天线方向图作为优化目标,基于推扫辐射计各接收链路的馈源次级天线方向图Ffeed(i),i=1~M,得到最优波束合成幅相加权系数矩阵C0(K×N),K为合成波束数目,N为每个合成波束中所选择接收链路数目;
S3、采用耦合器对推扫辐射计所有接收链路馈入幅度相位相同的相干噪声信号,得到推扫辐射计各波束合成选择的接收链路中接收通道的幅相不一致性Δg′(K×N);
S4、采用推扫辐射计各波束合成选择的接收链路中接收通道的幅相不一致性Δg′(K×N),更新波束合成的幅相加权系数C1(K×N),得到定标后的波束合成天线方向图Fb'eam(k),k=1~K;
S5、通过推扫辐射计系统观测微波辐射亮温已知的定标场,以推扫辐射计观测微波辐射亮温已知的定标场输出的功率信号PModel作为输入层、定标场的微波辐射亮温信息TModel作为输出层,步骤S4所得到的定标后的波束合成天线方向图Fb'eam(k),k=1~K为卷积神经网络模型卷积层特征图K个特征图的初始值,对卷积神经网络模型进行训练,确定卷积神经网络模型参数,得到与推扫辐射计等效的卷积神经网络模型,实现波束合成推扫辐射计的全链路误差的定标。
优选地,所述步骤S2先采用遗传算法初步搜索得到波束合成幅相加权系数矩阵的全局最优解;再利用序列二次规划算法加强局部搜索,最终得到最优波束合成幅相加权系数矩阵C0(K×N)。
优选地,所述步骤S3中接收链路中接收通道的幅相不一致性包括接收链路中接收通道之间的相位不一致性和幅度不一致性。
优选地,所述接收链路中接收通道之间的相位不一致性通过如下方法得到:
以某一个接收链路作为参考链路,将所有接收链路中接收通道输出的电压信号与参考链路中接收通道输出的电压信号作复相关,得到的相关系数的相位即为接收链路中接收通道之间的相位不一致性。
优选地,所述接收链路中接收通道之间的幅度不一致性通过如下方法得到:
以某一个接收链路作为参考链路,所有接收链路中接收通道输出的电压信号作自相关得到接收链路中接收通道的输出功率,各接收链路中接收通道输出功率与参考链路中接收通道的输出功率相除,则得到接收链路中接收通道之间的幅度不一致性。
优选地,所述步骤S4中,更新后的幅相加权系数矩阵C1(K×N)为:
C1(K×N)=C0(K×N)·*Δg′(K×N)
式中,“·*”表示两个矩阵的对应元素相乘,Δg′(K×N)的每一行元素对应每一个波束合成选择的N个接收链路中接收通道的幅相不一致性。
优选地,所述步骤S6采用梯度优化方法对卷积神经网络模型参数进行优化,直至卷积神经网络的误差小于预设阈值。
本发明与现有技术相比的优点在于:
(1)、本发明建立了基于卷积神经网络的波束合成推扫辐射计定标方法,通过地面测量得到的馈源次级天线方向图数据,周期性的耦合相干噪声粗略得到接收链路中接收通道的幅相不一致性。再通过推扫辐射计系统观测微波辐射亮温已知的定标场,利用卷积神经网络算法的多层监督、深度学习特性,建立推扫辐射计系统的反向模型,完成波束合成推扫辐射计系统的全链路误差定标;
(2)、本发明将幅度相位相同的相干噪声信号利用耦合的方式注入到所有接收链路中的接收通道,耦合方式对噪声功率要求小,再利用所有接收链路中接收通道输出的电压信号之间作复相关和接收通道自身信号作自相关,从相关系数和功率信号中得到接收链路中接收通道的相位和幅度不一致性,能够减小波束合成的幅相加权系数的不确定性。
(3)、本发明所有接收链路中接收通道输出信号仅需与参考通道(例如:选择第一个接收通道)作一次复相关运算和接收通道自身自相关运算,就能够从辐射计系统输出的相关系数和功率信号中得到所有接收链路中接收通道的幅度和相位不一致性,这样,利用极少的运算次数初步完成接收链路中接收通道之间的幅相不一致性定标,能够初步对后续定标方法的定标结果进行约束,避免后续的卷积神经网络模型搜索过程中落入局部最优值,有利于卷积神经网络最优参数的搜索。
(4)、本发明提出了基于卷积神经网络的方法对波束合成的加权系数进行定标,弥补了传统多个接收链路之间仅对接收通道部分存在的幅相不一致性进行定标的不足,能够实现对推扫辐射计系统的全链路误差定标。
(5)、本发明是针对多路接收链路之间的幅相不一致性导致推扫辐射计系统直接输出性能下降的问题,提出的同时对多个接收链路误差进行定标的方法,可以应用于星载、机载及地面的多波束多通道接收的辐射计及雷达系统,提高系统应用性能。
附图说明
图1为本发明实施例推扫辐射计系统组成示意图;
图2为本发明实施例波束合成链路组成图;
图3(a)为本发明实施例测试得到的天线方向图示意图;
图3(b)为本发明实施例期望得到的天线方向图示意图;
图4为本发明实施例方法流程图;
图5为本发明实施例卷积神经网络流程图;
图6为本发明实施例卷积神经网络模型训练过程。
具体实施方式
下面结合附图对本发明的具体实施方式进行进一步的详细描述。
区别于传统实孔径辐射计系统,波束合成推扫辐射计无需进行圆锥扫描,一次观测就能获取高分辨率、宽观测视场目标场景的微波辐射亮温,克服了高分辨率对大口径天线机械扫描的工程实现难题,成为微波遥感领域的研究热点。但该辐射计系统在应用时由于环境、系统等不理想性导致系统的接收链路(主要包括天线、接收通道)存在误差,此时通过系统输出的功率信号与辐射计系统在地面获取的馈源次级天线方向图等信息无法正确求解出观测目标场景的微波辐射亮温信息。为了获取高精度的亮温信息,必须对推扫辐射计系统的误差进行定标。而传统的辐射计系统没有波束合成的环节,故传统的定标方法无法适用于波束合成推扫辐射计系统的误差定标。
推扫辐射计通过星载的手段,对地球进行观测,能够实现对大气、海洋、极区等地球目标场景的微波辐射亮温的测量,本方法发明了基于卷积神经网络的波束合成推扫辐射计定标方法,能够提升系统的测量精度,高精度测量结果将为气象预报、海洋环境监测、灾害预报等提供重要的信息,对提高短期极端天气以及长期气候预测的准确性有着非常重要的作用,对我国的国防建设和国民生活有着极其重要的意义。
如图1所示,波束合成推扫辐射计系统包括:环焦反射面、密集馈源阵列、M个馈源、M个垂直极化接收通道、M个水平极化接收通道、数据处理器。以下以M个接收链路为对象进行描述,接收链路由馈源和接收通道组成,所述接收通道为垂直极化接收通道或者水平极化接收通道。
对于每个接收链路而言,环焦反射面将观测目标场景的微波辐射亮温信号反射、聚焦至馈源阵列的馈源,馈源接收微波辐射亮温信号转换成射频信号,送至接收通道对信号放大、下变频、滤波,变成中频信号;所有接收链路中接收通道输出的电压信号送至数据处理器,由数据处理器对所有接收链路中接收通道输出的电压信号进行复相关或自相关运算。
如图2所示,推扫辐射计系统中共有M个接收链路,为了形成满足需求的波束,每个波束由N个接收链路中接收通道的输出信号加权求和得到,而M个馈源共能形成K个波束,K个波束间的接收链路存在交叉复用。数据处理器按照每个波束选择的N个接收链路中接收通道输出信号完成波束合成并作自相关,得到输出每个波束的功率信号,根据每个波束的功率信号,采用两点定标方法,得到目标场景的微波辐射亮温信息。
本发明提供了基于卷积神经网络的波束合成推扫辐射计定标方法,具体步骤如下:
S1、测量得到推扫辐射计的环焦反射面结合密集馈源阵列条件下的所有接收链路的馈源次级天线方向图Ffeed(i),i=1~M,M为推扫辐射计接收链路数目;
优选方案具体如下:
每个馈源次级天线方向图,具体为:在球面近场将推扫辐射计系统的环焦反射面和密集馈源阵列以馈源阵列中心单元几何位置中心放置于球面近场机械臂规定的中心靶点上(由测量场地提供);通过控制电机扫描的方式旋转和移动球面近场的机械臂,使球面近场射频发射信号(由测量场地提供)位于馈源阵列坐标系下的不同位置,覆盖馈源的全空间立体角,即可获得所有馈源在全空间立体角的射频信号;将所有馈源的射频信号与球面近场射频发射信号(由测量场地提供)的幅度相除、相位相减,获得每个馈源的次级天线方向图Ffeed。
S2、以期望得到的各波束合成后天线方向图作为优化目标,基于推扫辐射计各接收链路的馈源次级天线方向图Ffeed(i),i=1~M,得到最优波束合成幅相加权系数矩阵C0(K×N),K为合成波束数目,N为每个合成波束中所选择接收链路数目;
优选方案具体如下:
如图3(a)和图3(b)所示,单个馈源的次级天线方向图性能很难进一步提升,为了获得波束窄、低旁瓣的天线方向图,采用多个馈源合成一个波束的形式,得到波束合成后的天线方向图。馈源编号1、2……M的次级天线方向图分别为Ffeed(1)、Ffeed(2)……Ffeed(M),合成波束编号1、2……K对应的天线方向图为Fbeam1、Fbeam2……FbeamK。波束合成天线方向图为期望得到的目标(如图3(b)所示波束合成后天线方向图),馈源次级天线方向图由步骤S1得到,波束合成天线方向图之间的关系有:
式中,Ffeed(i)表示第i个接收链路的馈源次级天线方向图,上式可简写为矩阵相乘:
Fbeam(K×1)=C0(K×N)·*F′feed(K×N)
式中,“·*”表示两个矩阵的对应元素相乘,F′feed(K×N)矩阵中的每一行表示每一个波束所选择的N个接收链路中对应的N个馈源次级天线方向图。
求解上式的波束合成幅相加权系数矩阵C0,未知数个数远大于方程个数,故直接用矩阵求逆的方法获取加权系数矩阵是一个病态的数学过程,所得的解不稳定,存在极大的误差。
本步骤使用的遗传算法,是一种借鉴生物界自然选择和进化机制发展起来的高度并行、随机、自适应搜索算法。使用群体搜索技术,将种群代表一组问题解,通过对当前种群施加如上述波束合成天线方向图之间的关系式中表达的约束、选择、交叉和变异等一系列遗传操作,产生新一代的种群,并逐步使种群进化到包含近似最优解的状态,快速找到全局最优解;
再利用序列二次规划算法,在每个迭代点构造一个二次规划子问题,以子问题的解作为迭代搜索方向并沿该方向进行一维搜索,通过反复迭代最终逼近最优解C0(K×N)。
通过遗传算法和序列二次规划算法搜索得到的波束加权系数矩阵C0(K×N),结合馈源天线方向图数据,就可以得到每个波束合成后实际的天线方向图Fbeam。
波束合成的幅相加权系数初始值C0是基于地面测量得到的馈源天线方向图所得到的,而没有考虑接收链路中馈源、接收通道的幅相特性发生变化时带来的影响。由于工程实现的不理想性和工作环境温度的变化,会引起接收通道的幅相不一致性发生变化。同时,地面研制的天线在轨后由于展开、环境温度等的变化,馈源次级天线方向图会发生变化,建立星载应用时实际波束合成天线方向图模型为:
式中,Δg(i)表示第i个接收链路的幅相特性误差量,ΔF(i)表示第i个接收链路馈源的次级天线方向图误差。可见,在轨后波束合成天线方向图与地面测试获得的波束合成天线方向图存在较大差异。
S3、采用耦合器对推扫辐射计所有接收链路馈入幅度相位相同的相干噪声信号,得到推扫辐射计各波束合成选择的接收链路中接收通道的幅相不一致性Δg′(K×N);
接收链路中接收通道的幅相不一致性包括接收链路中接收通道之间的相位不一致性和幅度不一致性。
本发明通过耦合器对所有接收链路中的接收通道馈入幅度相位相同的相干噪声信号,以某一个接收链路作为参考链路(选择第一个接收通道),将所有接收链路中接收通道输出的电压信号与参考链路中接收通道输出的电压信号作复相关,得到相关系数的相位即为接收通道之间的相位不一致性;
以某一个接收链路作为参考链路(选择第一个接收通道),将所有接收链路中接收通道输出的电压信号作自相关,输出的功率信号与参考链路中接收通道相除,则得到接收链路中接收通道之间的幅度不一致性,至此,得到所有接收链路中接收通道的幅相不一致性Δg,优选方案具体如下:
从接收通道的工程实现可知,接收通道的幅相不一致性特征比较大,为了获取全局最优值,需要通过耦合器对所有接收通道馈入幅度相位相同的相干噪声信号。将所有接收通道输出电压信号与参考通道(选择第一个接收通道)输出的电压信号相乘,得到复相关系数的实部Vreal;将所有接收通道输出电压信号进行90°移相之后再与参考通道输出的电压信号相乘,得到互相关系数的虚部Vimag;接收通道之间的相位不一致性可以从复相关系数中求解:
式中,atan表示三角函数的反正切函数。
所有接收链路中接收通道输出的电压信号与自身信号相乘,得到接收通道的自相关值,即功率信号Power。各接收链路中接收通道的功率信号Poweri,i=1~M与参考接收通道(选择第一个接收通道)功率信号Power1的比值,即为接收链路中接收通道的幅度不一致性:
式中,pi为第i个接收链路中接收通道的幅度不一致性。
接收通道的幅相不一致性Δg为:
推扫辐射计各波束合成选择的接收链路中接收通道的幅相不一致性矩阵Δg′(K×N)表达如下式,矩阵中的每一行表示每一个波束所选择的N个接收链路中接收通道的幅相不一致性。
S4、采用推扫辐射计各波束合成选择的接收链路中接收通道的幅相不一致性Δg′(K×N),更新波束合成的幅相加权系数C1(K×N),得到定标后的波束合成天线方向图Fb'eam(k),k=1~K;
优选方案具体如下:
接收链路中接收通道之间的幅相不一致性定标后波束合成天线方向图为:
更新后的幅相加权系数矩阵C1:
C1(K×N)=C0(K×N)·*Δg′(K×N)
式中,“·*”表示两个矩阵的对应元素相乘,Δg′(K×N)中的每一行表示每一个波束所选择的N个接收链路中接收通道的幅相不一致性。获取接收通道的幅相不一致性特性之后,在轨波束合成的实际天线方向图与步骤S4定标后的波束合成方向图之间仍然存在误差,残留的误差为在轨后馈源天线方向图的误差以及接收通道的幅相不一致性定标后的残差。步骤S3对波束合成天线方向图进行了初步定标,避免了步骤S5采用卷积神经网络搜索过程中落入局部最优值,有利于该方法获取波束合成推扫辐射计系统全链路优化的结果。
S5、通过推扫辐射计系统观测微波辐射亮温已知的定标场,以推扫辐射计观测微波辐射亮温已知的定标场输出的功率信号PModel作为输入层、定标场的微波辐射亮温信息TModel作为输出层,步骤S4所得到的定标后的波束合成天线方向图Fb'eam(k),k=1~K为卷积神经网络模型卷积层特征图K个特征图的初始值,对卷积神经网络模型进行训练,确定卷积神经网络模型参数,得到与推扫辐射计等效的卷积神经网络模型,实现波束合成推扫辐射计的全链路误差的定标。
优选方案具体如下:
如图5所示,卷积神经网络分为:输入层、卷积层、池化层、全连接层以及输出层。输入层:推扫辐射计系统观测目标场景输出的功率信号,即波束合成后的自相关值。卷积层用于特征提取,是卷积神经网络实现的核心,不同的卷积核用于不同特征的提取,卷积核越多,能够对输入数据提取的特征就越多。每一个特征提取层的神经元的输入与上一层的局部相连,通过神经元得到该局部区域的特征。池化层的作用是在保证信息完整的基础上,减少卷积层数据量,提升卷积神经网络的运算速度。全连接层实际上就是神经网络中隐藏层的部分,全连接层的神经元与上一层的池化层的神经元上的节点相互连接,但是同一全连接层中的各个神经元是不会相互连接的。输出层的即为观测目标场景的微波辐射亮温信息。
卷积层的计算:
式中,为第l层的第j'个特征图,i'表示选择的卷积核函数行数,Kil'j'为第l层的卷积核函数,f(·)为激励函数,/>为偏置参数,Mj'为选择输入特征图的集合。结合步骤S4得到的初步定标后的波束合成天线方向图,获得卷积层第一层的共K个特征图的初始值卷积层l-1之后连接到下一个卷积层l的梯度为:
式中,为卷积层中l层第j'个特征图与卷积层l-1层的误差信号,up表示升力采样操作,u、v表示矩阵的每个元素的位置坐标。
池化层的计算,原理就是每个输出特征图的大小就是输入特征图的缩小版,如下式所示:
式中,down(·)为下采样函数,β为乘性偏置参数,b为加性偏置参数。
池化层梯度计算:
全连接层使各层的神经元之间保持全连接,用以模拟卷积神经网络模型误差结果oj':
式中,k'表示第k'层卷积层,为卷积层中第k'层每个特征图j'的误差信号。当oj'满足期望值时(自定义误差阈值:误差阈值越小,得到的模型越准确,定标效果越好,但是运算量增加导致运算时间长),卷积层与池化层停止梯度更新,完成全局最优值搜索,返回的/>即为适用于波束合成推扫辐射计系统的卷积神经网络模型。
如图6所示,通过推扫辐射计系统观测微波辐射亮温已知的定标场(可选择冷空、海面定标场、雨林或者沙漠等定标场),定标场的微波辐射亮温信息为TModel时,推扫辐射计系统输出功率信号为PModel,PModel作为输入层,TModel作为输出层,步骤S4定标后的波束合成天线方向图F′beam作为卷积神经网络模型初始值。通过不断的循环、迭代,训练模型,直至搜索得到满足误差阈值时的最优模型参数,包括卷积层和池化层的参数选取。
综上所述,本发明基于卷积神经网络的全链路定标方法,结合其他步骤得到的先验信息,再根据推扫辐射计周期性的对已知微波辐射亮温信息的定标场进行观测,利用定标场的亮温信息以及推扫辐射计系统输出的功率信号,对卷积神经网络的模型参数进行循环、迭代,不断训练,直至搜索得到满足误差阈值时的最优模型参数。推扫辐射计工作时,通过对微波辐射亮温已知的定标场的观测,获取定标数据样本,可以周期性的进行卷积神经网格模型参数训练。
基于定标后获得的最优卷积神经网络模型参数,应用于波束合成推扫辐射计观测其他微波辐射计亮温未知的目标场景时,结合推扫辐射计系统输出的功率信号,可以获得高探测精度的微波辐射亮温信息,以此实现波束合成推扫辐射计的全链路误差高精度定标。
本发明说明书中未作详细描述的内容属于本领域技术人员的公知技术。
Claims (7)
1.基于卷积神经网络的波束合成推扫辐射计定标方法,其特征在于步骤如下:
S1、测量得到推扫辐射计的环焦反射面结合密集馈源阵列条件下的所有接收链路的馈源次级天线方向图Ffeed(i),i=1~M,M为推扫辐射计接收链路数目;
S2、以期望得到的各波束合成后天线方向图作为优化目标,基于推扫辐射计各接收链路的馈源次级天线方向图Ffeed(i),i=1~M,得到最优波束合成幅相加权系数矩阵C0(K×N),K为合成波束数目,N为每个合成波束中所选择接收链路数目;
S3、采用耦合器对推扫辐射计所有接收链路馈入幅度相位相同的相干噪声信号,得到推扫辐射计各波束合成选择的接收链路中接收通道的幅相不一致性Δg′(K×N);
S4、采用推扫辐射计各波束合成选择的接收链路中接收通道的幅相不一致性Δg′(K×N),更新波束合成的幅相加权系数C1(K×N),得到定标后的波束合成天线方向图F′beam(k),k=1~K;
S5、通过推扫辐射计系统观测微波辐射亮温已知的定标场,以推扫辐射计观测微波辐射亮温已知的定标场输出的功率信号PModel作为输入层、定标场的微波辐射亮温信息TModel作为输出层,步骤S4所得到的定标后的波束合成天线方向图F′beam(k),k=1~K为卷积神经网络模型卷积层特征图K个特征图的初始值,对卷积神经网络模型进行训练,确定卷积神经网络模型参数,得到与推扫辐射计等效的卷积神经网络模型,实现波束合成推扫辐射计的全链路误差的定标。
2.根据权利要求1所述的基于卷积神经网络的波束合成推扫辐射计定标方法,其特征在于所述步骤S2先采用遗传算法初步搜索得到波束合成幅相加权系数矩阵的全局最优解;再利用序列二次规划算法加强局部搜索,最终得到最优波束合成幅相加权系数矩阵C0(K×N)。
3.根据权利要求1所述的基于卷积神经网络的波束合成推扫辐射计定标方法,其特征在于所述步骤S3中接收链路中接收通道的幅相不一致性包括接收链路中接收通道之间的相位不一致性和幅度不一致性。
4.根据权利要求3所述的基于卷积神经网络的波束合成推扫辐射计定标方法,其特征在于所述接收链路中接收通道之间的相位不一致性通过如下方法得到:
以某一个接收链路作为参考链路,将所有接收链路中接收通道输出的电压信号与参考链路中接收通道输出的电压信号作复相关,得到的相关系数的相位即为接收链路中接收通道之间的相位不一致性。
5.根据权利要求3所述的基于卷积神经网络的波束合成推扫辐射计定标方法,其特征在于所述接收链路中接收通道之间的幅度不一致性通过如下方法得到:
以某一个接收链路作为参考链路,所有接收链路中接收通道输出的电压信号作自相关得到接收链路中接收通道的输出功率,各接收链路中接收通道输出功率与参考链路中接收通道的输出功率相除,则得到接收链路中接收通道之间的幅度不一致性。
6.根据权利要求1所述的基于卷积神经网络的波束合成推扫辐射计定标方法,其特征在于所述步骤S4中,更新后的幅相加权系数矩阵C1(K×N)为:
C1(K×N)=C0(K×N).*Δg′(K×N)
式中,“.*”表示两个矩阵的对应元素相乘,Δg′(K×N)的每一行元素对应每一个波束合成选择的N个接收链路中接收通道的幅相不一致性。
7.根据权利要求1所述的基于卷积神经网络的波束合成推扫辐射计定标方法,其特征在于所述步骤S6采用梯度优化方法对卷积神经网络模型参数进行优化,直至卷积神经网络的误差小于预设阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111274687.4A CN114046888B (zh) | 2021-10-29 | 2021-10-29 | 基于卷积神经网络的波束合成推扫辐射计定标方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111274687.4A CN114046888B (zh) | 2021-10-29 | 2021-10-29 | 基于卷积神经网络的波束合成推扫辐射计定标方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114046888A CN114046888A (zh) | 2022-02-15 |
CN114046888B true CN114046888B (zh) | 2024-02-20 |
Family
ID=80206732
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111274687.4A Active CN114046888B (zh) | 2021-10-29 | 2021-10-29 | 基于卷积神经网络的波束合成推扫辐射计定标方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114046888B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117196418B (zh) * | 2023-11-08 | 2024-02-02 | 江西师范大学 | 一种基于人工智能的阅读教学质量评估方法及系统 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004028731A (ja) * | 2002-06-25 | 2004-01-29 | Natl Space Development Agency Of Japan | 人工衛星用トータルパワー方式マイクロ波放射計の校正方法 |
CN101936778A (zh) * | 2010-08-05 | 2011-01-05 | 中国兵器工业第二○六研究所 | 一种辐射计定标方法 |
CN101975971A (zh) * | 2010-09-30 | 2011-02-16 | 中国科学院国家天文台 | 探月卫星微波探测仪的星上定标方法 |
CN105988146A (zh) * | 2015-01-29 | 2016-10-05 | 中国科学院空间科学与应用研究中心 | 一种星载微波辐射计的应用数据处理方法 |
CN109212496A (zh) * | 2018-09-26 | 2019-01-15 | 西安空间无线电技术研究所 | 一种星载微波辐射计天线误差校正方法 |
CN110617889A (zh) * | 2019-08-28 | 2019-12-27 | 西安空间无线电技术研究所 | 一种应用于综合孔径微波辐射计的高稳定性测试方法 |
CN110716185A (zh) * | 2019-10-22 | 2020-01-21 | 国家卫星气象中心 | 一种用于星载微波辐射计的在轨辐射校正方法 |
CN111538000A (zh) * | 2020-03-30 | 2020-08-14 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 均匀圆阵列综合孔径辐射计亮温反演成像方法 |
CN111929651A (zh) * | 2020-07-30 | 2020-11-13 | 西安空间无线电技术研究所 | 一种综合孔径辐射计系统相位残差标定方法 |
CN112285659A (zh) * | 2020-07-30 | 2021-01-29 | 西安空间无线电技术研究所 | 一种基于综合孔径辐射计在轨更新亮温重构矩阵的方法 |
DE112019003384T5 (de) * | 2018-07-02 | 2021-03-25 | Elecs Industry Co., Ltd. | Verfahren zum kalibrieren eines mikrowellen-radiometers |
CN113218510A (zh) * | 2021-05-13 | 2021-08-06 | 上海航天测控通信研究所 | 一种星载一维综合孔径辐射计的馈源口面外定标装置 |
CN113484917A (zh) * | 2021-06-30 | 2021-10-08 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种微波辐射计自动化定标的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI434057B (zh) * | 2011-11-24 | 2014-04-11 | Ind Tech Res Inst | 輻射計之校正裝置、校正系統與校正方法 |
-
2021
- 2021-10-29 CN CN202111274687.4A patent/CN114046888B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004028731A (ja) * | 2002-06-25 | 2004-01-29 | Natl Space Development Agency Of Japan | 人工衛星用トータルパワー方式マイクロ波放射計の校正方法 |
CN101936778A (zh) * | 2010-08-05 | 2011-01-05 | 中国兵器工业第二○六研究所 | 一种辐射计定标方法 |
CN101975971A (zh) * | 2010-09-30 | 2011-02-16 | 中国科学院国家天文台 | 探月卫星微波探测仪的星上定标方法 |
CN105988146A (zh) * | 2015-01-29 | 2016-10-05 | 中国科学院空间科学与应用研究中心 | 一种星载微波辐射计的应用数据处理方法 |
DE112019003384T5 (de) * | 2018-07-02 | 2021-03-25 | Elecs Industry Co., Ltd. | Verfahren zum kalibrieren eines mikrowellen-radiometers |
CN109212496A (zh) * | 2018-09-26 | 2019-01-15 | 西安空间无线电技术研究所 | 一种星载微波辐射计天线误差校正方法 |
CN110617889A (zh) * | 2019-08-28 | 2019-12-27 | 西安空间无线电技术研究所 | 一种应用于综合孔径微波辐射计的高稳定性测试方法 |
CN110716185A (zh) * | 2019-10-22 | 2020-01-21 | 国家卫星气象中心 | 一种用于星载微波辐射计的在轨辐射校正方法 |
CN111538000A (zh) * | 2020-03-30 | 2020-08-14 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 均匀圆阵列综合孔径辐射计亮温反演成像方法 |
CN111929651A (zh) * | 2020-07-30 | 2020-11-13 | 西安空间无线电技术研究所 | 一种综合孔径辐射计系统相位残差标定方法 |
CN112285659A (zh) * | 2020-07-30 | 2021-01-29 | 西安空间无线电技术研究所 | 一种基于综合孔径辐射计在轨更新亮温重构矩阵的方法 |
CN113218510A (zh) * | 2021-05-13 | 2021-08-06 | 上海航天测控通信研究所 | 一种星载一维综合孔径辐射计的馈源口面外定标装置 |
CN113484917A (zh) * | 2021-06-30 | 2021-10-08 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种微波辐射计自动化定标的方法 |
Non-Patent Citations (4)
Title |
---|
In-Orbit External Calibration Method for Synthetic Aperture Radiometer;Jiakun Wang 等;《 IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing》;第14卷;3763 - 3773 * |
一种一维综合孔径微波辐射计的定标方法;吴琼, 吴季;遥感技术与应用(03);53-56 * |
基于神经网络控制的星载微波辐射计;赵晓雪;张晓辉;;微计算机信息(13) * |
旋转圆环阵列综合孔径微波辐射计系统研究;韩东浩;《中国博士学位论文全文数据库 (信息科技辑)》(第09期);I140-62 * |
Also Published As
Publication number | Publication date |
---|---|
CN114046888A (zh) | 2022-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Stoffelen et al. | Ambiguity removal and assimilation of scatterometer data | |
US5345539A (en) | Radar apparatus using neural network for azimuth and elevation detection | |
CN109725317B (zh) | 一种基于一维综合孔径微波辐射计的海面亮温成像仿真方法 | |
CN111538000B (zh) | 均匀圆阵列综合孔径辐射计亮温反演成像方法 | |
CN114046888B (zh) | 基于卷积神经网络的波束合成推扫辐射计定标方法 | |
CN114894819A (zh) | 一种星载gnss-r土壤湿度反演方法 | |
CN115079109A (zh) | 一种合成孔径雷达成像质量指标评价系统 | |
CN115184916A (zh) | 一种海面风速联合反演方法、装置、介质及计算设备 | |
Wu et al. | Performance analysis of circular antenna array for microwave interferometric radiometers | |
Jia et al. | Rice biomass estimation using radar backscattering data at S-band | |
CN108344968A (zh) | 一种基于多峰量子布谷鸟搜索机制的正交传播算子测向方法 | |
CN114187533A (zh) | 一种基于随机森林时序分类的GB-InSAR大气改正方法 | |
Kerr et al. | New radiometers: SMOS-a dual pol L-band 2D aperture synthesis radiometer | |
Dodson et al. | The science case for simultaneous mm-wavelength receivers in radio astronomy | |
Liu et al. | A combined L-band synthetic aperture radiometer and fan-beam scatterometer for soil moisture and ocean salinity measurement | |
Zhang et al. | Global sea surface height measurement from CYGNSS based on machine learning | |
CN113687446B (zh) | 一种基于星载综合孔径微波辐射计的近场外定标方法 | |
Tournadre et al. | High-resolution imaging of the ocean surface backscatter by inversion of altimeter waveforms | |
Radhakrishnan et al. | Rainfall Estimation from Tempest-D Cubesat Observations | |
CN116184343A (zh) | 基于相控阵雷达的三维空间蜂群目标检测及信息估计方法 | |
Vasudevan et al. | Backpropagation neural-network-based retrieval of atmospheric water vapor and cloud liquid water from IRS-P4 MSMR | |
Alam et al. | Microwave radiometer calibration using deep learning with reduced reference information and 2-d spectral features | |
Wang et al. | Coastal Sea Surface Temperature Inversion from Microwave Radiometer using Radial Basis Function Neural Network | |
CN114325625B (zh) | 海洋-电离层混沌动力学特征提取方法、装置、介质、应用 | |
Zhang et al. | Results Analysis of Coastal Regions Sea Surface Salinity Retrieval from Aquarius Mission Using Deep Neural Network |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |