Disclosure of Invention
The technical problem to be solved by the invention is as follows: the method overcomes the defects of the prior art, provides the beam forming push-broom radiometer calibration method based on the convolutional neural network, improves the detection precision of the beam forming push-broom radiometer system, and provides a method for obtaining high-precision brightness and temperature.
The technical solution of the invention is as follows: a beam synthesis push-broom radiometer calibration method based on a convolutional neural network comprises the following steps:
s1, measuring to obtain a feed secondary antenna directional pattern F of the ring focal reflection surface of the push-broom radiometer combined with all receiving links under the condition of dense feed arrayfeed(i) I is 1-M, and M is the number of receiving links of the push-broom radiometer;
s2, taking each antenna directional diagram after beam synthesis expected to be obtained as an optimization target, and taking a feed source secondary antenna directional diagram F of each receiving link based on push-broom radiometerfeed(i) I is 1 to M, and an optimal beam forming amplitude-phase weighting coefficient matrix C is obtained0(K×N)K is the number of synthesized beams and N is the number of selected receive links in each synthesized beam;
s3, feeding coherent noise signals with the same amplitude and phase to all receiving links of the push-broom radiometer by adopting the coupler to obtain the amplitude and phase inconsistency delta g 'of the receiving channels in the receiving links selected by the wave beam synthesis of the push-broom radiometer'(K×N);
S4, adopting push-scan radiometer beamsSynthesizing amplitude phase disparity Δ g 'for receive channels in selected receive links'(K×N)Updating the amplitude-phase weighting coefficient C of the beam forming1(K×N)Obtaining a calibrated beam forming antenna directional pattern Fb'eam(k),k=1~K;
S5, observing the calibration field with known microwave radiation brightness and temperature through the push-scan radiometer system, and observing the power signal P output by the calibration field with known microwave radiation brightness and temperature through the push-scan radiometer systemModelMicrowave radiation brightness and temperature information T as input layer and scaling fieldModelAs an output layer, the scaled beam-forming antenna pattern F obtained in step S4b'eam(k) And K is 1-K and is the initial value of K characteristic graphs of the convolutional neural network model convolutional layer characteristic graph, the convolutional neural network model is trained, the parameters of the convolutional neural network model are determined, the convolutional neural network model equivalent to the push-broom radiometer is obtained, and the calibration of the full link error of the beam synthesis push-broom radiometer is realized.
Preferably, in step S2, a genetic algorithm is first used to initially search for a global optimal solution of the beam-forming amplitude-phase weighting coefficient matrix; then, local search is enhanced by utilizing a sequential quadratic programming algorithm, and finally, an optimal beam forming amplitude-phase weighting coefficient matrix C is obtained0(K×N)。
Preferably, the amplitude-phase inconsistency of the receiving channels in the receiving chain in step S3 includes phase inconsistency and amplitude inconsistency between the receiving channels in the receiving chain.
Preferably, the phase inconsistency between the receiving channels in the receiving chain is obtained by the following method:
and taking one receiving link as a reference link, and performing complex correlation on the voltage signals output by the receiving channels in all the receiving links and the voltage signals output by the receiving channels in the reference link to obtain the phase of the correlation coefficient, namely the phase inconsistency among the receiving channels in the receiving links.
Preferably, the amplitude inconsistency between the receiving channels in the receiving chain is obtained by the following method:
and taking one receiving link as a reference link, and performing autocorrelation on voltage signals output by the receiving channels in all the receiving links to obtain the output power of the receiving channels in the receiving links, wherein the amplitude inconsistency between the receiving channels in the receiving links is obtained by dividing the output power of the receiving channels in each receiving link by the output power of the receiving channels in the reference link.
Preferably, in the step S4, the updated amplitude-phase weighting coefficient matrix C1(K×N)Comprises the following steps:
C1(K×N)=C0(K×N)·*Δg′(K×N)
where "· denotes the multiplication of corresponding elements of the two matrices, Δ g'(K×N)Each row of elements of (a) corresponds to the amplitude disparity of the receive channels in the N receive chains selected by each beamforming.
Preferably, in the step S6, the convolutional neural network model parameters are optimized by using a gradient optimization method until an error of the convolutional neural network is smaller than a preset threshold.
Compared with the prior art, the invention has the advantages that:
(1) the invention establishes a beam synthesis push-broom radiometer calibration method based on a convolutional neural network, and the amplitude-phase inconsistency of a receiving channel in a receiving link is roughly obtained through periodic coupling coherent noise through the directional diagram data of a secondary antenna of a feed source obtained through ground measurement. Then, observing a calibration field with known microwave radiation brightness and temperature through the push-broom radiometer system, and establishing a reverse model of the push-broom radiometer system by utilizing the multilayer supervision and deep learning characteristics of the convolutional neural network algorithm to finish the whole link error calibration of the beam synthesis push-broom radiometer system;
(2) the coherent noise signals with the same amplitude and phase are injected into the receiving channels of all receiving chains in a coupling mode, the coupling mode has low requirement on noise power, complex correlation is carried out between voltage signals output by the receiving channels in all the receiving chains, self-correlation is carried out on the signals of the receiving channels, the phase and amplitude inconsistency of the receiving channels in the receiving chains is obtained from the correlation coefficient and the power signals, and the uncertainty of the amplitude-phase weighting coefficient of beam forming can be reduced.
(3) The amplitude and phase inconsistency of the receiving channels in all receiving links can be obtained from the correlation coefficient and the power signal output by the radiometer system only by performing one-time complex correlation operation on the output signal of the receiving channel in all receiving links and a reference channel (for example, selecting the first receiving channel) and performing self-correlation operation on the receiving channel, so that the amplitude and phase inconsistency calibration between the receiving channels in the receiving links is preliminarily completed by using few operation times, the calibration result of a subsequent calibration method can be preliminarily restricted, the phenomenon that the subsequent convolutional neural network model search process falls into a local optimal value is avoided, and the search of the optimal parameter of the convolutional neural network is facilitated.
(4) The invention provides a convolutional neural network-based method for calibrating the weighting coefficients of beam forming, overcomes the defect that the traditional multiple receiving links only calibrate the amplitude-phase inconsistency existing in the receiving channel part, and can realize the calibration of the full link error of the push-broom radiometer system.
(5) The invention provides a method for calibrating errors of a plurality of receiving links simultaneously aiming at the problem that the direct output performance of a push-broom radiometer system is reduced due to amplitude-phase inconsistency among a plurality of receiving links, and the method can be applied to spaceborne, airborne and ground multi-beam multi-channel receiving radiometers and radar systems and improves the application performance of the system.
Detailed Description
The following describes embodiments of the present invention in further detail with reference to the accompanying drawings.
Different from the traditional real-aperture radiometer system, the beam synthesis push-broom radiometer can acquire the microwave radiation brightness temperature of a target scene with high resolution and wide observation field of view without cone scanning, overcomes the engineering realization problem of mechanical scanning of a large-aperture antenna by high resolution, and becomes a research hotspot in the field of microwave remote sensing. However, when the radiometer system is applied, due to the non-idealities of the environment, the system and the like, an error exists in a receiving link (mainly comprising an antenna and a receiving channel) of the system, and at the moment, microwave radiation brightness and temperature information of an observation target scene cannot be correctly solved through a power signal output by the system and information such as a feed source secondary antenna directional diagram and the like acquired by the radiometer system on the ground. In order to obtain high precision brightness temperature information, the error of the push-scan radiometer system must be calibrated. However, the conventional radiometer system has no beam forming step, so the conventional calibration method cannot be applied to error calibration of the beam forming push-broom radiometer system.
The push-broom radiometer observes the earth by means of satellite-borne, and can measure the microwave radiation brightness and temperature of earth target scenes such as atmosphere, ocean, polar regions and the like.
As shown in fig. 1, the beam-forming push-scan radiometer system includes: the device comprises a ring focus reflecting surface, a dense feed source array, M feed sources, M vertical polarization receiving channels, M horizontal polarization receiving channels and a data processor. The following description is given by taking M receiving links as objects, where the receiving link is composed of a feed source and a receiving channel, and the receiving channel is a vertical polarization receiving channel or a horizontal polarization receiving channel.
For each receiving link, the ring focus reflecting surface reflects and focuses microwave radiation brightness and temperature signals of an observation target scene to a feed source of a feed source array, the feed source receives the microwave radiation brightness and temperature signals, converts the microwave radiation brightness and temperature signals into radio frequency signals, and sends the radio frequency signals to a receiving channel to amplify, down-convert and filter the signals into intermediate frequency signals; the voltage signals output by the receiving channels in all the receiving links are sent to the data processor, and the data processor carries out complex correlation or autocorrelation operation on the voltage signals output by the receiving channels in all the receiving links.
As shown in fig. 2, the push-scan radiometer system has M receiving chains, and in order to form beams meeting requirements, each beam is obtained by weighted summation of output signals of receiving channels in the N receiving chains, and the M feed sources can form K beams, and the receiving chains among the K beams have cross multiplexing. The data processor completes beam synthesis and autocorrelation according to the output signals of the receiving channels in the N receiving links selected by each beam to obtain a power signal for outputting each beam, and microwave radiation brightness and temperature information of a target scene is obtained by adopting a two-point calibration method according to the power signal of each beam.
The invention provides a beam forming push-broom radiometer calibration method based on a convolutional neural network, which comprises the following specific steps of:
s1, measuring to obtain a feed secondary antenna directional pattern F of the ring focal reflection surface of the push-broom radiometer combined with all receiving links under the condition of dense feed arrayfeed(i) I is 1-M, and M is the number of receiving links of the push-broom radiometer;
the preferred scheme is as follows:
each feed source secondary antenna directional diagram specifically comprises: placing a ring focal reflecting surface and a dense feed source array of the push-broom radiometer system on a central target point (provided by a measuring field) specified by the spherical near field mechanical arm according to the geometric position center of a feed source array central unit in the spherical near field; the mechanical arm of the spherical near field is rotated and moved by controlling the scanning mode of the motor, so that the radio frequency emission signals (provided by a measurement field) of the spherical near field are positioned at different positions under the array coordinate system of the feed source, and the full space solid angle of the feed source is covered, thus obtaining all the feed sourcesRadio frequency signals at a full spatial solid angle; dividing the radio frequency signals of all the feed sources by the amplitude and phase of the spherical near-field radio frequency transmitting signal (provided by the measuring field) to obtain a secondary antenna directional diagram F of each feed sourcefeed。
S2, taking each antenna directional diagram after beam synthesis expected to be obtained as an optimization target, and taking a feed source secondary antenna directional diagram F of each receiving link based on push-broom radiometerfeed(i) I is 1 to M, and an optimal beam forming amplitude-phase weighting coefficient matrix C is obtained0(K×N)K is the number of synthesized beams and N is the number of selected receive links in each synthesized beam;
the preferred scheme is as follows:
as shown in fig. 3(a) and fig. 3(b), the performance of the secondary antenna pattern of a single feed is difficult to further improve, and in order to obtain an antenna pattern with a narrow beam and low side lobes, a beam is synthesized by a plurality of feeds to obtain an antenna pattern after beam synthesis. The secondary antenna patterns of the feed source numbers 1 and 2 … … M are respectively Ffeed(1)、Ffeed(2)……Ffeed(M) the antenna pattern F corresponding to the synthesized beam numbers 1 and 2 … … Kbeam1、Fbeam2……FbeamK. The beamformed antenna pattern is the desired target (e.g., the beamformed antenna pattern shown in fig. 3 (b)), the feed secondary antenna pattern is obtained in step S1, and the beamformed antenna pattern has the following relationship:
in the formula, Ffeed(i) The feed secondary antenna pattern representing the ith receive chain can be abbreviated as matrix multiplication:
Fbeam(K×1)=C0(K×N)·*F′feed(K×N)
in the formula (I), the compound is shown in the specification,". denotes the multiplication of corresponding elements of the two matrices, F'feed(K×N)Each row in the matrix represents a corresponding N feed secondary antenna patterns in the N receive chains selected for each beam.
Solving the beam forming amplitude-phase weighting coefficient matrix C of the above formula0The unknown number is far greater than the equation number, so that the method of directly using matrix inversion to obtain the weighting coefficient matrix is a pathological mathematical process, the obtained solution is unstable, and great errors exist.
The genetic algorithm used in the step is a highly parallel, random and self-adaptive search algorithm developed by taking the natural selection and evolution mechanism of the biology as a reference. Using a group search technology, representing the group as a group of problem solutions, generating a new generation of group by applying a series of genetic operations such as constraint, selection, intersection, variation and the like expressed in the relational expression between the beam-forming antenna directional diagrams on the current group, gradually evolving the group to a state containing an approximately optimal solution, and quickly finding out an overall optimal solution;
and then constructing a quadratic programming subproblem at each iteration point by using a sequential quadratic programming algorithm, taking the solution of the subproblem as an iteration search direction, performing one-dimensional search along the direction, and finally approaching an optimal solution C through repeated iteration0(K×N)。
Wave beam weighting coefficient matrix C obtained by searching through genetic algorithm and sequence quadratic programming algorithm0(K×N)Combining the antenna directional pattern data of the feed source, the actual antenna directional pattern F after each wave beam is synthesized can be obtainedbeam。
Initial value C of amplitude-phase weighting coefficient of beam forming
0The antenna directional diagram is obtained based on a feed source antenna directional diagram obtained through ground measurement, and influence brought by changes of amplitude-phase characteristics of a feed source and a receiving channel in a receiving link is not considered. Due to non-idealities in engineering implementation and variations in operating environment temperature, variations in amplitude-phase non-uniformity of the receive channel may result. Meanwhile, the directional diagram of the secondary antenna of the feed source can be changed due to the change of the expansion, the environmental temperature and the like of the ground-developed antenna after the antenna is in orbit, and the actual beam combination is established during the satellite-borne applicationAntenna direction diagram model
Comprises the following steps:
in the formula, Δ g (i) represents the magnitude-phase characteristic error of the ith reception link, and Δ f (i) represents the secondary antenna pattern error of the ith reception link feed. Therefore, the on-orbit beam synthesis antenna directional pattern has a larger difference from the beam synthesis antenna directional pattern obtained by the ground test.
S3, feeding coherent noise signals with the same amplitude and phase to all receiving links of the push-broom radiometer by adopting the coupler to obtain the amplitude and phase inconsistency delta g 'of the receiving channels in the receiving links selected by the wave beam synthesis of the push-broom radiometer'(K×N);
The amplitude-phase disparity of the receive channels in the receive chain includes phase disparity and amplitude disparity between the receive channels in the receive chain.
The method comprises the steps of feeding coherent noise signals with the same amplitude and phase into receiving channels in all receiving links through a coupler, taking one receiving link as a reference link (selecting a first receiving channel), performing complex correlation on voltage signals output by the receiving channels in all the receiving links and voltage signals output by the receiving channels in the reference link, and obtaining the phase of a correlation coefficient, namely the phase inconsistency between the receiving channels;
using a certain receiving link as a reference link (selecting a first receiving channel), performing autocorrelation on voltage signals output by receiving channels in all receiving links, and dividing an output power signal by the receiving channels in the reference link to obtain amplitude inconsistency among the receiving channels in the receiving links, and thus obtaining the amplitude inconsistency Δ g of the receiving channels in all receiving links, wherein the preferable scheme is as follows:
as known from the engineering realization of the receiving channel, the amplitude-phase inconsistency characteristic of the receiving channel is large, and in order to obtain the global optimumIt is necessary to feed coherent noise signals of the same amplitude and phase to all reception channels via couplers. Multiplying the voltage signals output by all the receiving channels with the voltage signal output by the reference channel (selecting the first receiving channel) to obtain the real part V of the complex correlation coefficientreal(ii) a All the receiving channel output voltage signals are subjected to 90-degree phase shift and then multiplied by the voltage signals output by the reference channel to obtain an imaginary part V of the cross-correlation coefficientimag(ii) a The phase inconsistency between the receive channels can be solved from the complex correlation coefficients:
in the formula, atan represents an arctangent function of a trigonometric function.
The voltage signals output by the receiving channels in all receiving links are multiplied by the signals of the receiving channels to obtain the autocorrelation value of the receiving channels, namely the Power signal Power. Power signal Power of receiving channel in each receiving chainiI 1-M and a reference receiving channel (selecting the first receiving channel) Power signal Power1The ratio of (a) to (b) is the amplitude inconsistency of the receiving channels in the receiving link:
in the formula, piThe amplitude disparity of the receive channels in the ith receive chain.
The amplitude-phase inconsistency Δ g of the receive channel is:
amplitude-phase inconsistency matrix delta g 'of receive channels in receive chain selected by push-scan radiometer beam synthesis'(K×N)Expressing the following formula, each row in the matrix represents the amplitude phase inconsistency of the receiving channels in the N receiving chains selected by each beam.
S4 amplitude phase inconsistency deltag 'of receiving channels in receiving chain selected by adopting beam synthesis of push-scan radiometer'(K×N)Updating the amplitude-phase weighting coefficient C of the beam forming1(K×N)Obtaining a calibrated beam forming antenna directional pattern Fb'eam(k),k=1~K;
The preferred scheme is as follows:
the beam combination antenna directional diagram after the amplitude-phase inconsistency among the receiving channels in the receiving link is calibrated is as follows:
updated amplitude-phase weighting coefficient matrix C1:
C1(K×N)=C0(K×N)·*Δg′(K×N)
Where "· denotes the multiplication of corresponding elements of the two matrices, Δ g'(K×N)Each row in (a) represents the amplitude phase disparity of the receive channels in the N receive chains selected for each beam. After the amplitude-phase inconsistency of the receiving channel is obtained, an error still exists between the actual antenna directional pattern of the on-track beam synthesis and the beam synthesis directional pattern calibrated in step S4, and the residual error is an error of the on-track feed source antenna directional pattern and a residual error after the amplitude-phase inconsistency of the receiving channel is calibrated. Step S3 carries out preliminary calibration to the beam forming antenna directional diagram, avoids the situation that the partial optimal value falls in the convolution neural network searching process adopted in step S5, and is beneficial to obtaining the result of the full link optimization of the beam forming push-broom radiometer system by the method.
S5, observing the calibration field with known microwave radiation brightness and temperature through the push-scan radiometer system, and observing the power signal P output by the calibration field with known microwave radiation brightness and temperature through the push-scan radiometer systemModelMicrowave radiation brightness and temperature information T as input layer and scaling fieldModelAs an output layer, the scaled beam-forming antenna pattern F obtained in step S4b'eam(k) And K is 1-K and is the initial value of K characteristic graphs of the convolutional neural network model convolutional layer characteristic graph, the convolutional neural network model is trained, the parameters of the convolutional neural network model are determined, the convolutional neural network model equivalent to the push-broom radiometer is obtained, and the calibration of the full link error of the beam synthesis push-broom radiometer is realized.
The preferred scheme is as follows:
as shown in fig. 5, the convolutional neural network is divided into: an input layer, a convolutional layer, a pooling layer, a full-link layer, and an output layer. An input layer: the push-broom radiometer system observes a power signal output by a target scene, namely an autocorrelation value after beam synthesis. The convolution layer is used for feature extraction and is a core for realizing the convolutional neural network, different convolution kernels are used for extracting different features, and the more convolution kernels are, the more features can be extracted from input data. The input of the neuron of each feature extraction layer is connected with the local part of the previous layer, and the feature of the local area is obtained through the neuron. The pooling layer has the functions of reducing the data volume of the convolution layer and improving the operation speed of the convolution neural network on the basis of ensuring the integrity of information. The fully-connected layer is actually part of a hidden layer in the neural network, and the neurons of the fully-connected layer are connected with nodes on the neurons of the pooling layer of the previous layer, but the neurons in the same fully-connected layer are not connected with each other. The output layer is the microwave radiation brightness and temperature information of the observation target scene.
Calculation of convolutional layer:
in the formula (I), the compound is shown in the specification,
for the j 'th feature map of the l layer, i' represents the number of selected convolution kernel rows, K
il
'j'Is the convolution kernel of the l-th layer, f (-) is the excitation function,
as a bias parameter, M
j'To select a set of input feature maps. Combining the preliminarily scaled beam-forming antenna directional pattern obtained in step S4 to obtain the initial values of K characteristic patterns of the first layer of the convolutional layer
The gradient for the convolutional layer l-1 followed by the connection to the next convolutional layer l is:
in the formula (I), the compound is shown in the specification,
for the jth feature map of layer l of the convolutional layer and the error signal of layer l-1 of the convolutional layer, up represents the lift sampling operation, and u, v represent the position coordinates of each element of the matrix.
The principle of the pooling layer calculation is that the size of each output feature map is a reduced version of the input feature map, as shown in the following formula:
where down (-) is a downsampling function, β is a multiplicative bias parameter, and b is an additive bias parameter.
Gradient calculation of the pooling layer:
the full-connection layer keeps full connection between neurons of each layer to simulate a convolutional neural network modelType error result oj':
Wherein k 'represents a k' th convolutional layer,
error signal for each profile j 'of the k' th layer in the convolutional layer. When o is
j'When the expected value is met (the user-defined error threshold value is smaller, the obtained model is more accurate, the calibration effect is better, but the operation time is long due to the increase of the operation amount), the gradient updating of the convolution layer and the pooling layer is stopped, the global optimum value searching is completed, and the returned value
The convolution neural network model is suitable for the beam-forming push-broom radiometer system.
As shown in fig. 6, a calibration field (optionally, a cold air calibration field, a sea surface calibration field, a rainforest or a desert calibration field) with known microwave radiation brightness temperature is observed through the push-broom radiometer system, and the microwave radiation brightness temperature information of the calibration field is TModelWhen the output power signal of the push-broom radiometer system is PModel,PModelAs input layer, TModelAs the output layer, beam-synthesized antenna pattern F 'scaled in step S4'beamAs the initial value of the convolutional neural network model. And training the model through continuous circulation and iteration until the optimal model parameters meeting the error threshold are obtained by searching, wherein the optimal model parameters comprise parameter selection of the convolutional layer and the pooling layer.
In summary, the full-link calibration method based on the convolutional neural network of the present invention combines the prior information obtained in other steps, periodically observes the calibration field of the known microwave radiation brightness and temperature information according to the push-broom radiometer, and performs loop and iteration on the model parameters of the convolutional neural network by using the brightness and temperature information of the calibration field and the power signal output by the push-broom radiometer system, and continuously trains until the optimal model parameters satisfying the error threshold are obtained by searching. When the push-broom radiometer works, calibration data samples are obtained through observation of a calibration field with known microwave radiation brightness and temperature, and parameter training of a convolution neural grid model can be periodically carried out.
Based on the optimal convolutional neural network model parameters obtained after calibration, when the optimal convolutional neural network model parameters are applied to a beam forming push-broom radiometer to observe a target scene with unknown brightness temperature of other microwave radiometers, the microwave radiation brightness temperature information with high detection precision can be obtained by combining power signals output by a push-broom radiometer system, and therefore high-precision calibration of full link errors of the beam forming push-broom radiometer is achieved.
Those skilled in the art will appreciate that the details of the invention not described in detail in this specification are well within the skill of those in the art.