Mirror synthetic aperture radiometer imaging method and system based on weight vector optimization

Abstract

The present invention discloses a mirror-image synthetic aperture radiometer imaging method and system based on weight vector optimization, and relates to the field of microwave remote sensing and detection technology. The present invention includes: receiving antenna array arrangement data, generating a cosine transform matrix according to the antenna array arrangement data, performing inverse cosine transform on the cosine transform matrix to generate an initial weight vector; calculating the product of the initial weight vector and the cosine transform matrix to obtain an array factor, optimizing the array factor to suppress the influence of the side lobes or interference areas on the imaging direction; solving the transfer equation output by the antenna array to obtain cosine visibility; using the optimized weight vector to multiply the cosine visibility to obtain a reconstructed image. The present invention can regulate the influence of the non-imaging direction on the imaging direction in the array factor by optimizing the weight vector in the image reconstruction process, thereby reducing its side lobes, and in particular, can suppress the side lobes in the interference direction, thereby improving the quality of the reconstructed image.

Description

Technical Field

[0001] The invention relates to the field of microwave remote sensing and detection technology, and in particular to a mirror synthetic aperture radiometer imaging method and system based on weight vector optimization. Background Art

[0002] Spatial resolution is a key indicator of microwave radiometers. Existing microwave radiometers mainly use real aperture or synthetic aperture systems, which have difficulties in improving spatial resolution. The main difficulty of the real aperture system is the large-diameter antenna: large-diameter antennas are engineering challenges in terms of development, transportation and installation, and their surface deformation control, scanning system design and accuracy assurance are also difficult. The main difficulty of the synthetic aperture system is the high system complexity: the high spatial resolution indicator places high demands on the number of antennas of the synthetic aperture system radiometer. The number of subsequent receivers will increase accordingly, and the number of related devices is approximately the square of the number of antennas. Complex systems also face reliability assurance issues. To solve the above problems, researchers have proposed the concept of mirror synthetic aperture, and have preliminarily verified its principle. The spatial resolution can be improved by placing a reflector next to the antenna array, which is an option to solve the problem of improving spatial resolution in microwave radiometers.

[0003] However, the existing image reconstruction methods of mirror-synthetic aperture do not consider the relationship between the pixel value of the reconstructed image and the scene brightness temperature distribution. In fact, the pixel value of the mirror-synthetic aperture reconstructed image is related to the scene brightness temperature in all directions, and this relationship can be characterized by the array factor. The array factor of the existing reconstruction method is related to the pixel position, that is, the relationship between the pixel value in different brightness temperature inversion directions and the scene brightness temperature distribution is different, which manifests as spatial variation. The spatial variation of the array factor causes the pixel value of the reconstructed image to be disturbed by the brightness temperature in other directions in the scene, and this interference lacks an effective way to suppress it, causing the image quality to deteriorate. For this reason, we propose a mirror-synthetic aperture radiometer imaging method and system based on weight vector optimization. Summary of the invention

[0004] The object of the present invention is to provide a mirror aperture synthesis radiometer imaging method and system based on weight vector optimization to solve the problems raised in the above background technology.

[0005] According to a first aspect of the present invention, in order to achieve the above-mentioned purpose, the present invention provides the following technical solution: a mirror synthetic aperture radiometer imaging method based on weight vector optimization, comprising the following steps:

[0006] receiving antenna array arrangement data, generating a cosine transform matrix according to the antenna array arrangement data, and performing an inverse cosine transform on the cosine transform matrix to generate an initial weight vector;

[0007] Calculate the product of the initial weight vector and the cosine transform matrix to obtain the array factor, optimize the array factor, and reduce the impact of the side lobes or interference areas on the imaging direction;

[0008] Solving the transfer equation of the antenna array output yields the cosine visibility;

[0009] The optimized weight vector is multiplied by the cosine visibility to obtain the reconstructed image.

[0010] Furthermore, a cosine transform matrix is ​​generated according to the antenna array arrangement data, and an inverse cosine transform is performed on the cosine transform matrix to generate an initial weight vector, as follows:

[0011] (21) The antenna array is arranged as follows: [x1, x2, ...x N ], where x i represents the antenna position and is a positive integer multiple of 0.5, and N is the number of antennas;

[0012] (22) The cosine transform matrix CM is a matrix that describes the relationship from the original scene brightness temperature distribution vector to the cosine visibility vector. The following relationship is established based on CM:

[0013] CV=T·CM

[0014] Where CV is the cosine visibility matrix, T is the original scene brightness temperature distribution vector;

[0015] (23) Assume that the original scene brightness temperature distribution is discretized into a 1×P vector, and the number of cosine visibility obtained by the antenna array arrangement is Q, where Q is the number of elements in the following set:

[0016] {|x i -x j |}∪{|x i +x j |}

[0017] Where x i and x j are any two antennas in the antenna array;

[0018] (24) Then the cosine transform matrix CM can be expressed as:

[0019]

[0020] where ξ i , i = 1, 2, ... P is the direction cosine of the imaging field of view, P is the number of field divisions, u i For any two antennas x in the antenna array i and x j The baseline formed is equal to |x i -x j|,=1,2,..,N, j=1,2,...,N;

[0021] (25) For a certain direction ξ within the imaging field of view i , whose inverse cosine transform is defined as:

[0022]

[0023] Specifically:

[0024]

[0025] The initial weight vector Win0 obtained is:

[0026]

[0027] In the formula, ξ i , i = 1, 2, ... P is the direction cosine of the imaging field of view, u i For any two antennas x in the antenna array i and x j The baseline formed is equal to |x i -x j |,=1,2,...,N,j=1,2,...,N.

[0028] Furthermore, the product of the weight vector and the cosine transform matrix, i.e., the array factor, is optimized to reduce the influence of the side lobes or interference areas on the imaging direction, as follows:

[0029] (31)ξ i The relationship between the reconstructed brightness temperature of the direction and the brightness temperature distribution of the original scene can be expressed as:

[0030]

[0031] That is, the reconstructed brightness temperature in this direction is the product of the original scene brightness temperature distribution and the array factor AF;

[0032] (32) Establish an optimization model, optimize win0 in the above formula, and obtain the optimized weight vector Winop:

[0033]

[0034] Where ξ0 is the imaging direction, ξ j It is necessary to control the imaging interference direction of the array factor side lobe to be less than ε, ξ k These are other imaging directions in which the array factor side lobes need to be controlled to be smaller than μ.

[0035] Furthermore, the transfer equation of the antenna array output is solved to obtain the cosine visibility, which is as follows:

[0036] (41) For any antenna i and antenna j in the antenna array, the correlation value R ij It is expressed as:

[0037] R ij =CV(d i -d j )±CV(d i +d j )

[0038] The ± sign on the right side of the equation is determined by the polarization of the antenna in the antenna array, and CV is the cosine visibility, which is related to the scene brightness temperature distribution as follows:

[0039] Where ξ=sinθ is the direction cosine, for

[0040] Corrected scene brightness temperature distribution function per unit solid angle, u is the wavelength-normalized spatial frequency;

[0041] (42) The correlation values ​​of all antennas in the antenna array are combined to obtain the correlation value matrix R, and their corresponding cosine visibility values ​​are combined to obtain the cosine visibility matrix CV. The relationship between the correlation value matrix and the cosine visibility matrix can be described by a linear equation system, namely, the transfer equation:

[0042] R = P·CV, where P is called the transfer matrix, which is composed of the antenna array

[0043] Column decision;

[0044] The correlation value matrix R is obtained by mirroring the synthetic aperture radiometer, which is substituted into the transfer equation and solved to obtain the cosine visibility matrix CV.

[0045] Furthermore, the optimized weight vector is multiplied by the cosine visibility to obtain the reconstructed image, as follows:

[0046] Get ξ i The reconstructed brightness temperature in the direction is:

[0047]

[0048] By traversing all directions of the scene brightness temperature distribution, a reconstructed image of the scene brightness temperature distribution can be obtained.

[0049] According to a second aspect of the present invention, the present invention provides an imaging system of a one-dimensional mirror aperture synthesis radiometer based on weight vector optimization, which is used to implement the above-mentioned imaging method of the mirror aperture synthesis radiometer based on weight vector optimization, comprising:

[0050] An initial weight vector generation module is used to receive antenna array arrangement data, generate a cosine transform matrix according to the antenna array arrangement data, and perform an inverse cosine transform on the cosine transform matrix to generate an initial weight vector;

[0051] The first calculation module is used to calculate the product of the initial weight vector and the cosine transform matrix to obtain the array factor, optimize the array factor, and reduce the influence of the side lobe or interference area on the imaging direction;

[0052] The second calculation module is used to solve the transfer equation output by the antenna array to obtain the cosine visibility;

[0053] The third calculation module is used to multiply the optimized weight vector by the cosine visibility to obtain a reconstructed image.

[0054] According to a third aspect of the present invention, the present invention provides a terminal device, comprising a memory, a processor, and a computer program stored in the memory and capable of running on the processor, wherein the memory stores a computer program capable of running on the processor, and when the processor loads and executes the computer program, the above-mentioned mirror synthetic aperture radiometer imaging method based on weight vector optimization is adopted.

[0055] According to a fourth aspect of the present invention, the present invention provides a storage medium comprising computer executable instructions, which, when executed by a computer processor, are used to perform the above-mentioned mirror synthetic aperture radiometer imaging method based on weight vector optimization.

[0056] The present invention has at least the following beneficial effects:

[0057] The present invention can regulate the influence of the non-imaging direction on the imaging direction in the array factor by optimizing the weight vector in the image reconstruction process, thereby reducing its side lobe. In particular, the side lobe in the interference direction can be suppressed, thereby improving the quality of the reconstructed image.

[0058] Of course, any product implementing the present invention does not necessarily need to achieve all of the advantages described above at the same time. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] Figure 1 A schematic diagram of a process of the method according to an embodiment of the present invention;

[0060] Figure 2 A schematic diagram of the structure of a one-dimensional mirror aperture synthesis microwave radiometer provided by an embodiment of the present invention;

[0061] Figure 3 A schematic diagram of an antenna array receiving a signal in an embodiment of the present invention;

[0062] Figure 4 A schematic diagram of the antenna array arrangement in an embodiment of the present invention;

[0063] Figure 5 Schematic diagram comparing the conventional array factor and the optimized array factor described in the embodiments of the present invention. DETAILED DESCRIPTION

[0064] The following will be combined with the drawings in the embodiments of the present disclosure to clearly and completely describe the technical solutions in the embodiments of the present disclosure. Obviously, the described embodiments are only part of the embodiments of the present disclosure, rather than all the embodiments. Based on the embodiments in the present disclosure, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present disclosure.

[0065] Figure 2 shows a schematic diagram of the structure of a one-dimensional mirror aperture synthesis radiometer. Figure 3 shows a schematic diagram of a one-dimensional mirror-synthetic aperture antenna array receiving signal. Figure 1 A flow chart of a mirror-type synthetic aperture radiometer imaging method based on weight vector optimization is shown. For ease of explanation, only the part related to the present invention is shown, which is described in detail as follows:

[0066] like Figure 2 As shown, the one-dimensional mirror aperture synthesis radiometer includes an antenna array, a reflector, a receiving channel array, an A / D array and a correlator, wherein the antenna array simultaneously receives the radiation signal from the scene and the reflected signal via the reflector, the output high-frequency signal of the antenna array is amplified and down-converted by the receiving channel array to obtain an intermediate frequency signal, which is then processed by the A / D array to obtain a digital signal, and finally the correlator performs correlation operation on the digital signal to obtain a correlation value matrix, and finally a reconstructed image is obtained from the correlation value matrix through a computer program.

[0067] like Figure 3 As shown, for any antenna i and antenna j in the antenna array, the correlation value Rij is expressed as:

[0068] Rij=CV(di-dj)±CV(di+dj)

[0069] The ± sign on the right side of the equation is determined by the polarization of the antenna in the antenna array, and CV is the cosine visibility, which is related to the scene brightness temperature distribution as follows:

[0070]

[0071] Among them, ξ=sinθ is the direction cosine, is the scene brightness temperature distribution function per unit solid angle, and u is the wavelength-normalized spatial frequency.

[0072] The correlation values ​​of all antennas in the antenna array are combined to obtain the correlation value matrix R, and their corresponding cosine visibility is combined to obtain the cosine visibility matrix CV. Due to the linear relationship between a single correlation value and two cosine visibility, the relationship between the correlation value matrix and the cosine visibility matrix can be described by a linear equation system, namely the transfer equation:

[0073] R=P·CV

[0074] Among them, P is called the transfer matrix, which is determined by the antenna array.

[0075] The correlation value matrix R is obtained by mirroring the synthetic aperture radiometer, and it is substituted into the transfer equation to solve the cosine visibility matrix CV. The discrete inverse cosine transform is performed on them to obtain the reconstructed image.

[0076]

[0077] Where Win0 is the discrete inverse cosine transform matrix.

[0078] See also Figure 1 , the imaging method of mirror aperture radiometer based on weight vector optimization includes the following steps:

[0079] S1: Generate a cosine transform matrix based on the antenna array arrangement of the one-dimensional mirror aperture radiometer, as follows:

[0080] (21) The antenna array is arranged as follows: [x1, x2, ...x N ], where x i represents the antenna position and is a positive integer multiple of 0.5, and N is the number of antennas;

[0081] (22) The cosine transform matrix CM is a matrix that describes the relationship from the original scene brightness temperature distribution vector to the cosine visibility vector. Based on CM, the following relationship can be established:

[0082] CV=T·CM

[0083] Where CV is the cosine visibility matrix, T is the original scene brightness temperature distribution vector;

[0084] Assume that the original scene brightness temperature distribution is discretized into a 1×P vector, and the number of cosine visibility obtained by the antenna array arrangement is Q, where Q is the number of elements in the following set:

[0085] {|x i -x j |}∪{|x i +x j |}

[0086] Where x i and xj are any two antennas in the antenna array;

[0087] Then the discrete cosine transform matrix CM can be expressed as:

[0088]

[0089] Due to the cosine transform relationship between the original scene brightness temperature distribution and the cosine visibility, the inverse cosine transform can be used to reconstruct the brightness temperature distribution from the cosine visibility. Specifically, after solving the transfer moment equation to obtain the cosine visibility, if we want to reconstruct ξ i The brightness temperature in the direction is calculated by constructing the initial weight vector Win0 based on the inverse cosine transform. Then the reconstructed brightness temperature in this direction can be calculated as:

[0090]

[0091] Where Win0 is:

[0092]

[0093] In the formula, ξ i , i = 1, 2, ... P is the direction cosine of the imaging field of view, u i For any two antennas x in the antenna array i and x j The baseline formed is equal to |x i -x j |,=1,2,..,N, j=1,2,...,N;

[0094] S2: Optimize the product of the weight vector and the cosine transform matrix, that is, the array factor, to reduce the impact of the side lobes or certain interference areas on the imaging direction. From this, we can get, ξ i The relationship between the reconstructed brightness temperature of the direction and the brightness temperature distribution of the original scene is expressed as:

[0095]

[0096] That is, the reconstructed brightness temperature in this direction is the product of the brightness temperature distribution of the original scene and the array factor AF, that is, the reconstructed brightness temperature in a certain direction is related to the brightness temperatures in all directions of the imaging scene, and only one of these directions is the reconstructed brightness temperature direction, and the brightness temperatures in the remaining directions can be regarded as interference;

[0097] In order to make the reconstructed brightness temperature in a certain direction reflect the actual brightness temperature in the original scene as realistically as possible, win0 in the above formula can be optimized so that the value of this direction in AF is constant and the values ​​of other directions are as small as possible. Specifically, the following optimization model is established:

[0098]

[0099] Where ξ0 is the imaging direction, ξ j It is necessary to control the imaging interference direction of the array factor side lobe to be less than ε, ξ k These are other imaging directions in which the array factor side lobes need to be controlled to be smaller than μ.

[0100] It should be noted that based on the above optimization model, intelligent algorithms such as simulated annealing algorithm, particle swarm algorithm, genetic algorithm, etc. can be used for optimization, or ready-made tools of mathematical software can be used directly for optimization, and finally the optimized weight vector Winop can be obtained;

[0101] S3: Solve the transfer equation of the antenna array output to obtain the cosine visibility;

[0102] CV=A·R

[0103] Where A is the pseudo-inverse of P, and the MP generalized inverse method, Tikhonov regularization method, truncated singular value method and other algorithms can be used. The relevant output of the mirror comprehensive aperture belongs to the first kind of Fredholm integral equation. The first kind of integral equation with a fully continuous operator is unstable even if it has a solution. In addition, the actual measurement data inevitably contains systematic errors and noise, which will cause the approximate solution to deviate seriously from the true solution. Due to the underdetermined nature of the matrix, the matrix inversion based on the MP generalized inverse is very ill-posed. In order to obtain a stable approximate solution to the ill-posed problem, a regularization method is needed. The basic idea of ​​the regularization method is to redefine the concept of the solution to the ill-posed problem using some additional information of the specific problem, and then introduce a positive definite functional to give a stable method to approximate the solution of the original problem, that is, to find a suitable inverse operator to make a stable approximate solution. Therefore, considering the systematic error and noise, the truncated singular value decomposition and Tikhonov regularization are used to obtain a stable approximate solution of the cosine visibility function;

[0104] S4: multiply the optimized weight vector by the cosine visibility to obtain the reconstructed image;

[0105] Based on the optimized weight vector Winop obtained in step S2 and the cosine visibility CV obtained in step S3, we get ξ i The reconstructed brightness temperature in the direction is:

[0106]

[0107] By traversing all directions of the scene brightness temperature distribution, a reconstructed image of the scene brightness temperature distribution can be obtained.

[0108] In order to further illustrate the imaging method of the mirror aperture radiometer based on weight vector optimization provided by the embodiment of the present invention, it is now described in detail as follows in conjunction with a specific example:

[0109] This embodiment uses the array factor sidelobe reduction of the mirror synthesis aperture of a 12-element uniform array. In this embodiment, Figure 4 For the array used in the simulation, a 12-unit uniform array is located on the right side of the reflector. Figure 5 For comparison of the two array factors, the array factor without weight vector optimization is drawn with a red line, and the array factor with weight vector optimization is drawn with a blue line.

[0110] The specific steps are as follows:

[0111] (1) For a 12-element uniform array, the distance between the reflector and the nearest antenna in the antenna array is set to 1 times the minimum antenna spacing, and the corresponding transfer matrix, cosine visibility number, and discrete cosine transform matrix are obtained accordingly;

[0112] (2) Starting from the first direction in the field of view, take it as the imaging direction, calculate the initial weight vector, multiply the initial weight vector by the cosine transform matrix, and obtain the array factor optimized by the unused weight vector in this direction;

[0113] (3) Starting from the initial weight vector, the influence of other directions in the array factor on the imaging direction is optimized to obtain the optimized weight vector, and further obtain the optimized array factor;

[0114] (4) For (2) and (3), calculate the array factors without weight vector optimization and with weight vector optimization for all imaging directions within the field of view;

[0115] From the simulation results, it can be seen that Figure 5 As shown in the figure, at the direction of 0.05, the conventional array factor reaches the maximum sidelobe of -10.8621, while the optimized array factor is -16.5409, a decrease of about 6dB; at the interference direction of 0.105714, the sidelobe of the conventional array factor reaches -22.0674, while the optimized array factor is -53.4878, a decrease of more than 30dB.

[0116] In summary, the present invention can optimize the weight vector to control the influence of the non-imaging direction on the imaging direction in the array factor, thereby reducing its side lobe, and can suppress the side lobe in the interference direction, thereby improving the quality of the reconstructed image.

[0117] It should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device.

[0118] For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to the specific circumstances. When an element is referred to as being "assembled on", "installed on", "fixed on" or "set on" another element, it can be directly on the other element or there can also be a centered element. When an element is considered to be "connected" to another element, it can be directly connected to the other element or there may be a centered element at the same time. The terms "vertical", "horizontal", "up", "down", "left", "right" and similar expressions used herein are only for illustrative purposes and are not intended to be the only implementation method.

[0119] Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the appended claims and their equivalents.

[0120] In the description of this specification, the description with reference to the terms "one embodiment", "example", "specific example", etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present disclosure. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

Claims

1. A mirror aperture radiometer imaging method based on weight vector optimization, characterized in that: The following steps are involved: receiving antenna array arrangement data, generating a cosine transform matrix according to the antenna array arrangement data, and performing an inverse cosine transform on the cosine transform matrix to generate an initial weight vector; Calculate the product of the initial weight vector and the cosine transform matrix to obtain the array factor, optimize the array factor, and reduce the impact of the side lobes or interference areas on the imaging direction; Solving the transfer equation of the antenna array output yields the cosine visibility; The optimized weight vector is multiplied by the cosine visibility to obtain the reconstructed image.

2. The image synthesis aperture radiometer imaging method based on weight vector optimization according to claim 1, characterized in that: Generate a cosine transform matrix based on the antenna array arrangement data, and perform an inverse cosine transform on the cosine transform matrix to generate an initial weight vector, as follows: (21) The antenna array is arranged as follows: [x1, x2, ...x N ], where x i represents the antenna position and is a positive integer multiple of 0.5, and N is the number of antennas; (22) The cosine transform matrix CM is a matrix that describes the relationship from the original scene brightness temperature distribution vector to the cosine visibility vector. The following relationship is established based on CM: CV=T·CM Where CV is the cosine visibility matrix, T is the original scene brightness temperature distribution vector; (23) Assume that the original scene brightness temperature distribution is discretized into a 1×P vector, and the number of cosine visibility obtained by the antenna array arrangement is Q, where Q is the number of elements in the following set: {|x i -x j |}∪{|x i +x j |} Where x i and x j are any two antennas in the antenna array; (24) Then the cosine transform matrix CM can be expressed as: where ξ i , i = 1, 2, ... P is the direction cosine of the imaging field of view, P is the number of field divisions, u i For any two antennas x in the antenna array i and x j The baseline formed is equal to |x i -x j |,=1,2,..,N, j=1,2,...,N; (24) For a certain direction ξ within the imaging field of view i , whose inverse cosine transform is defined as: Specifically: The initial weight vector Win0 obtained is: In the formula, ξ i , i = 1, 2, ... P is the direction cosine of the imaging field of view, u i For any two antennas x in the antenna array i and x j The baseline formed is equal to |x i -x j |,=1,2,...,N,j=1,2,...,N.

3. The image synthesis aperture radiometer imaging method based on weight vector optimization according to claim 2, characterized in that: The product of the weight vector and the cosine transform matrix, i.e. the array factor, is optimized to reduce the impact of the side lobes or interference areas on the imaging direction, as follows: (31)ξ i The relationship between the reconstructed brightness temperature of the direction and the brightness temperature distribution of the original scene can be expressed as: That is, the reconstructed brightness temperature in this direction is the product of the original scene brightness temperature distribution and the array factor AF; (32) Establish an optimization model, optimize win0 in the above formula, and obtain the optimized weight vector Winop: Where ξ0 is the imaging direction, ξ j It is necessary to control the imaging interference direction of the array factor side lobe to be less than ε, ξ k These are other imaging directions in which the array factor side lobes need to be controlled to be smaller than μ.

4. The imaging method of mirror aperture radiometer based on weight vector optimization according to claim 1, characterized in that: Solving the transfer equation for the antenna array output yields the cosine visibility as follows: (41) For any antenna i and antenna j in the antenna array, the correlation value R ij It is expressed as: R ij =CV(d i -d j )±CV(d i +d j ) The ± sign on the right side of the equation is determined by the polarization of the antenna in the antenna array, and CV is the cosine visibility, which is related to the scene brightness temperature distribution as follows: Where ξ=sinθ is the direction cosine, is the scene brightness temperature distribution function for the modified unit solid angle, and u is the wavelength-normalized spatial frequency; (42) The correlation values ​​of all antennas in the antenna array are combined to obtain the correlation value matrix R, and their corresponding cosine visibility values ​​are combined to obtain the cosine visibility matrix CV. The relationship between the correlation value matrix and the cosine visibility matrix can be described by a linear equation system, namely, the transfer equation: R = P·CV, where P is called the transfer matrix, which is determined by the antenna array; The correlation value matrix R is obtained by mirroring the synthetic aperture radiometer, which is substituted into the transfer equation and solved to obtain the cosine visibility matrix CV.

5. The imaging method of mirror aperture radiometer based on weight vector optimization according to claim 1, characterized in that: The optimized weight vector is multiplied by the cosine visibility to obtain the reconstructed image, as follows: Get ξ i The reconstructed brightness temperature in the direction is: By traversing all directions of the scene brightness temperature distribution, a reconstructed image of the scene brightness temperature distribution can be obtained.

6. An imaging system of a one-dimensional mirror aperture synthesis radiometer based on weight vector optimization, used to implement the imaging method of a mirror aperture synthesis radiometer based on weight vector optimization according to any one of claims 1 to 5, characterized in that: include: An initial weight vector generation module is used to receive antenna array arrangement data, generate a cosine transform matrix according to the antenna array arrangement data, and perform an inverse cosine transform on the cosine transform matrix to generate an initial weight vector; The first calculation module is used to calculate the product of the initial weight vector and the cosine transform matrix to obtain the array factor, optimize the array factor, and reduce the influence of the side lobe or interference area on the imaging direction; The second calculation module is used to solve the transfer equation output by the antenna array to obtain the cosine visibility; The third calculation module is used to multiply the optimized weight vector by the cosine visibility to obtain a reconstructed image.

7. A terminal device, comprising a memory, a processor, and a computer program stored in the memory and capable of running on the processor, characterized in that: The memory stores a computer program that can be run on the processor. When the processor loads and executes the computer program, the image synthesis aperture radiometer imaging method based on weight vector optimization described in any one of claims 1 to 5 is adopted.

8. A storage medium containing computer executable instructions, characterized in that: When the computer executable instructions are executed by a computer processor, they are used to perform the image synthesis aperture radiometer imaging method based on weight vector optimization according to any one of claims 1 to 5.

Images