CN112200880A - Solid rocket engine combustion surface retreating image reconstruction method - Google Patents

Solid rocket engine combustion surface retreating image reconstruction method Download PDF

Info

Publication number
CN112200880A
CN112200880A CN202010646466.4A CN202010646466A CN112200880A CN 112200880 A CN112200880 A CN 112200880A CN 202010646466 A CN202010646466 A CN 202010646466A CN 112200880 A CN112200880 A CN 112200880A
Authority
CN
China
Prior art keywords
image
combustion surface
iteration
projection
microwave
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
Application number
CN202010646466.4A
Other languages
Chinese (zh)
Inventor
王世辉
陆明
乔文生
王欢欢
黄家骥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
INNER MONGOLIA AEROSPACE POWER MACHINERY TESTING INSTITUTE
Original Assignee
INNER MONGOLIA AEROSPACE POWER MACHINERY TESTING INSTITUTE
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by INNER MONGOLIA AEROSPACE POWER MACHINERY TESTING INSTITUTE filed Critical INNER MONGOLIA AEROSPACE POWER MACHINERY TESTING INSTITUTE
Priority to CN202010646466.4A priority Critical patent/CN112200880A/en
Publication of CN112200880A publication Critical patent/CN112200880A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention relates to a dynamic change process of a combustion surface in a solid rocket engine test, in particular to a method for reconstructing a combustion surface migration image of a solid rocket engine. The method comprises the steps of obtaining a magnetic line path of a microwave CT imaging field of the combustion surface of the engine, obtaining microwave projection data of the combustion surface of the engine, reconstructing a retreating image of the combustion surface of the engine, and judging the reconstructed quality. The image reconstruction method provided by the invention mainly determines a projection matrix according to a projection geometric relation, constructs a projection observation model according to a projection physical process, and realizes target reconstruction through optimization solution. The reconstruction method has low requirement on the completeness of projection data, can introduce prior constraint information according to specific imaging conditions, can inhibit noise and finally achieves high reconstruction quality.

Description

Solid rocket engine combustion surface retreating image reconstruction method
Technical Field
The invention relates to a dynamic change process of a combustion surface in a solid rocket engine test, in particular to a method for reconstructing a combustion surface migration image of a solid rocket engine.
Background
The solid rocket engine has become the preferred power device of the strategic and tactical missile weapon system at present due to the advantages of small volume, strong maneuverability, convenient maintenance and use and the like. With the continuous development of the technology of the solid missile weapon system in China, the solid rocket engine of the new generation missile weapon system adopts a propellant with higher fuel value and more sensitivity. However, in the process of developing the novel solid rocket engine, problems such as large combustion speed prediction deviation, erosion combustion, unstable combustion, difficult technical zero-setting positioning and the like occur from the root of combustion, and the important reason for the problems is that the dynamic change process of the combustion surface in the solid rocket engine test is not clear, the relevant key parameters and change rules of the solid rocket engine are difficult to accurately obtain, and the development process of the solid rocket engine in China is seriously hindered. So far, the change of the combustion surface of the explosive column in the solid rocket engine in China can not be detected or can not be accurately detected, and the development level, the performance improvement and the quality zeroing of the solid rocket engine in China are seriously influenced. The research on the combustion surface retrogradation testing technology of the solid rocket engine is urgently needed, and the research on the combustion surface retrogradation image reconstruction method of the solid rocket engine is one of the problems which must be overcome.
The analytic method can obtain better reconstruction quality for complete projection data, has high reconstruction speed, but has higher requirement on the completeness of the data. When projection data are rare and projection angles are missing, high-quality reconstruction effect cannot be obtained. The reconstruction by the analytic method requires complete projection information, and generally data with more than 180 projection angles are reconstructed, so that better reconstruction quality can be obtained for complete projection data.
Disclosure of Invention
Technical problem to be solved
The invention provides a method for reconstructing a combustion surface retreating image of a solid rocket engine, which is used for improving the image edge reconstruction quality; the number of magnetic lines of force of the microwave CT imaging field is increased to increase the imaging information quantity, so that the reconstruction that the error of the retrogradation images of the solid rocket engines with different combustion surface diameters is within 1mm is realized.
Technical scheme adopted for solving technical problem
The method for reconstructing the combustion surface retreating image of the solid rocket engine comprises the following steps:
1) magnetic line path for acquiring combustion surface microwave CT imaging field
According to the electromagnetic wave attenuation characteristic, battery magnetic wave tomography is simulated, and further through a transmitting device, 29 receiving devices work simultaneously, and firstly, the magnetic line path of a preliminary microwave CT imaging field is obtained. Then according to the propagation rule of the electromagnetic field, respectively adding 2 magnetic lines on two sides of each magnetic line of the preliminarily obtained magnetic line path of the microwave CT imaging field so as to increase the information quantity and obtain the final magnetic line path of the combustion surface microwave CT imaging field;
2) acquiring microwave projection data of combustion surface
Solving a projection matrix according to a magnetic line simulation path equation, wherein a in the projection matrix AijDefining the length of an intersection line of an i-ray and an j-pixel, discretizing a reconstruction region according to a projection matrix, uniformly arranging n transmitting/receiving devices on a circumference which takes the center of a measured object as a circle center, taking one of the n transmitting/receiving devices as a transmitter each time, switching a receiving state, and detecting scattered data together with all other receiving devices, wherein the propagation of microwaves in the measured medium meets the wave equation in a common form:
Figure BDA0002573207890000021
wherein
Figure BDA0002573207890000022
Representing the reflectivity of a medium
Figure BDA0002573207890000023
The scalar function of (a) is selected,
Figure BDA0002573207890000024
is a vector laplacian operator. Through a series of deformation operations, the formula can be converted into:
Figure BDA0002573207890000025
wherein
Figure BDA0002573207890000026
Called objective function,kbRepresents the wave number, mu, of a plane wave in a background mediumrAnd εrRepresenting the relative permeability and relative permittivity relative to the background medium due to the total field
Figure BDA0002573207890000027
Can be expressed as an incident field
Figure BDA0002573207890000031
And a scattered field
Figure BDA0002573207890000032
And incident wave of
Figure BDA0002573207890000033
Can be defined as:
Figure BDA0002573207890000034
thus, it is possible to obtain:
Figure BDA0002573207890000035
with the green's function, the illumination produced by all sources can be found as:
Figure BDA0002573207890000036
thus obtaining "projection" data a;
3) combustion surface migration image reconstruction
To fully exploit the gradient sparsity of an image, a minimized image can be obtained by a formula, i.e.
Figure BDA0002573207890000037
Constraint conditions are as follows: subject to y is Ax
Wherein | x | purpleTVIs the TV norm of the image x,
Figure BDA0002573207890000038
representing an image function, from N image pixels x1,x2,…,xNAnd (4) forming. y is the attenuation of the microwave through the combustion surface and A is the projection matrix.
The formula for the TV norm of image x is as follows,
Figure BDA0002573207890000039
wherein,
Figure BDA00025732078900000310
image reconstruction is carried out by taking minimization of TV norm as a target, and an iterative formula is as follows
Figure BDA00025732078900000311
Wherein p isiFor projecting the measured values,. mu.jFor the attenuation coefficient to be determined, λ is the relaxation factor, aijConstructing an image reconstruction algorithm based on minimization of a TV norm by performing convex set orthogonal projection and TV noise reduction twice correction processing on a calculation result in each step of iterative calculation process for the length of the i microwave intersection j pixels;
4) reconstruction quality discrimination
Comparing the obtained reconstructed values of the combustion surfaces with different diameters with the actual values, finishing iteration when the error between the two values is not more than 1mm, proving that the reconstruction quality of the combustion surfaces meets the requirement, otherwise, continuing the iteration until the error between the two values is not more than 1mm,
δ=max|fi-fi′|
where δ is the error between the reconstructed value and the actual value, fiIs the true numerical value of the fire face diameter when the image serial number is i, fi' is the reconstructed value of the combustion surface diameter when the image number is i.
Further, the specific algorithm steps of the combustion surface retrogradation image reconstruction are as follows:
1) initial value selection: let S represent a star undersampled matrix, an initial observation set
Figure BDA0002573207890000042
Fourier transform coefficients representing the original image taken only at the star sample line position (the Fourier-centric slice theorem indicates that sampling along the star sample line in the frequency domain is inherently consistent with the random sampling required by CS theory)
Figure BDA0002573207890000043
Representing the multiplication of corresponding elements of the matrix. Further, let the iteration initial value x0=DFT-1(y0);
2) And (3) iterative calculation: solving and calculating the following iterative formula by using a conjugate gradient method:
Figure BDA0002573207890000044
wherein λ isnDenotes the iteration step size, dnRepresents the direction of the conjugate gradient;
3) orthogonal projection: orthogonal projection is performed on y as Ax, and calculation is performed first
Figure BDA0002573207890000045
DFT of (1), preserving the frequency domain information y on the initial sample point0The frequency domain information of the non-sampling point is updated by using the current iteration result without changing, and then the inverse DFT of the frequency domain updating result is calculated, namely
Figure BDA0002573207890000046
4) Denoising TV: the frequency-domain star sampling causes the iteration result to contain randomly distributed noise, so that the iteration result of the above formula is matched
Figure BDA0002573207890000047
Carrying out noise reduction treatment by using the following formula to obtain the (n + 1) th step iterationEnd result x ofn+1I.e. by
Figure BDA0002573207890000048
Wherein for any given two-dimensional image x ═ x (x)ij) The TV noise reduction matrix is TV _ noise (x) ═ tij) And is and
Figure BDA0002573207890000051
5) and (3) updating the conjugation direction: search direction d of conjugate gradient methodnIs a combination of the direction of the negative gradient of the current iteration and the direction of the last search iteration, i.e.
dn+1=dTV(xn+1)+βdn
Wherein d isTV(xn+1) Express | | xn+1||TVIn the direction of the gradient of (a),
Figure RE-GDA0002788221550000052
is a proportional coefficient, | ·| non-woven phosphorFA Frobenius norm representing a matrix;
6) termination criteria: and stopping iteration when the iteration times reach a certain set value or the iteration result is converged, or returning to the step II to continue the process by iterative computation.
Advantageous effects
The image reconstruction method provided by the invention mainly determines a projection matrix according to the projection geometric relation, constructs a projection observation model according to the projection physical process, and realizes target reconstruction through optimization solution. The reconstruction method has low requirement on the completeness of projection data, can introduce prior constraint information according to specific imaging conditions, can inhibit noise and finally achieves high reconstruction quality.
Drawings
FIG. 1: a three-dimensional reconstruction algorithm flow chart of a combustion surface retreating image of a solid rocket engine;
FIG. 2: a strong electric field simulation result diagram observed at the radiation edge of the antenna;
FIG. 3: a primary magnetic line path simulation graph;
FIG. 4: finally, a magnetic line path simulation graph is obtained;
FIG. 5: microwave tomography projection matrix images;
FIG. 6: microwave projection data of different combustion surface diameters;
FIG. 7: a combustion surface retreating reconstruction algorithm flow chart;
FIG. 8: reconstructing result graphs of different combustion surface diameters;
FIG. 9: comparing the combustion surface diameter results;
FIG. 10: error plots.
Detailed Description
FIG. 1 is a flowchart of a method for reconstructing a combustion surface retrogradation image of a solid rocket engine. The main processes of the invention comprise the steps of obtaining the magnetic line path of the microwave CT imaging field of the combustion surface of the engine, obtaining the microwave projection data of the combustion surface of the engine, reconstructing the retreating image of the combustion surface of the engine and judging the reconstruction quality.
The invention is further described with reference to the accompanying drawings and the detailed description.
As shown in figure 1, the method for reconstructing the combustion surface retrogradation image of the solid rocket engine comprises four parts of obtaining a magnetic line path of a microwave CT imaging field of the combustion surface of the engine, obtaining microwave projection data of the combustion surface of the engine, reconstructing the combustion surface retrogradation image of the engine and judging the reconstruction quality. The method comprises the following specific steps:
1) magnetic line path for acquiring combustion surface microwave CT imaging field
According to the electromagnetic wave attenuation characteristic, battery magnetic wave tomography is simulated, and when the diameter of a grain is 200mm and the number of the loop antenna arrays is 30, the included angle between two adjacent receiving antennas at the left and right of a transmitting antenna (when the antenna is 10mm away from the arc surface of the grain) is about 168 degrees. The simulation result of the strong electric field observed at the radiation edge of the antenna is shown in fig. 2. As can be seen from fig. 2, microstrip antenna radiation ensures that each receiving antenna can receive signals. Further according to the mechanism of 1 transmitting and 29 receiving, the path of magnetic lines of force for primarily acquiring the microwave CT imaging field is shown in fig. 3. Then, according to the propagation rule of the electromagnetic field, 2 magnetic lines are added on two sides of each magnetic line in fig. 2 to increase the information amount, and the final magnetic line path of the combustion surface microwave CT imaging field is obtained, as shown in fig. 4.
2) Acquiring microwave projection data of combustion surface
Solving the projection matrix for the magnetic flux line simulation path equation of fig. 4 is shown in fig. 5. The reconstructed area was discretized using a circular geometry model as shown in fig. 5, and 29 transmitters/receivers were uniformly arranged on a circle centered on the center of the measured object. One of them is used as a transmitter each time, and then the receiving state is switched to detect the scattered data together with all other receivers. In the formula
Figure BDA0002573207890000061
On the basis of the obtained projection data, the microwave receiving values of 3 magnetic lines in the final magnetic line path of the combustion surface microwave CT imaging field are all set to be the same value, and microwave projection data of the combustion surface 30 × 90 is obtained, as shown in fig. 6.
3) Combustion surface migration image reconstruction
As shown in FIG. 7, the reconstruction algorithm includes six parts, namely initial value selection, iterative computation, orthogonal projection, TV noise reduction, conjugate direction update and termination criterion. According to the formula
Figure BDA0002573207890000071
Figure BDA0002573207890000072
Minimizing the obtained TV norm as a target according to a formula
Figure BDA0002573207890000073
The iterative formula of (2) is reconstructed, and the specific steps of the algorithm are as follows:
firstly, selecting an initial value: let S represent a star undersampled matrix, an initial observation set
Figure BDA0002573207890000074
Fourier transform coefficients representing the original image taken only at the star sample line position (the Fourier center slice theorem indicates that the star sample line is followed in the frequency domainThe sampling is inherently consistent with the random sampling required by Cs theory), where
Figure BDA0002573207890000075
Representing the multiplication of corresponding elements of the matrix. Further, let the iteration initial value x0=DFT-1(y0);
And (2) iterative calculation: solving and calculating the following iterative formula by using a conjugate gradient method:
Figure BDA0002573207890000076
wherein λ isnDenotes the iteration step size, dnThe conjugate gradient direction is indicated.
③ orthogonal projection: orthogonal projection is performed on y as Ax, and calculation is performed first
Figure BDA0002573207890000077
DFT of (1), preserving the frequency domain information y on the initial sample point0The frequency domain information of the non-sampling point is updated by using the current iteration result without changing, and then the inverse DFT of the frequency domain updating result is calculated, namely
Figure BDA0002573207890000078
Fourthly, noise reduction of the TV: the frequency-domain star sampling causes the iteration result to contain randomly distributed noise, so that the iteration result of the above formula is matched
Figure BDA0002573207890000079
Carrying out noise reduction treatment by using the following formula to obtain a final result x of the (n + 1) th iterationn+1I.e. by
Figure BDA00025732078900000710
Wherein for any given two-dimensional image x ═ x (x)ij) The TV noise reduction matrix is TV _ noise (x) ═ tij) And is and
Figure BDA0002573207890000081
conjugate direction updating: search direction d of conjugate gradient methodnIs a combination of the direction of the negative gradient of the current iteration and the direction of the last search iteration, i.e.
dn+1=dTV(xn+1)+βdn
Wherein d isTV(xn+1) Express | | xn+1||TVIn the direction of the gradient of (a),
Figure RE-GDA0002788221550000082
is a proportional coefficient, | ·| non-woven phosphorFThe Frobenius norm of the matrix is represented.
Sixthly, termination criterion: and stopping iteration when the iteration times reach a certain set value or the iteration result is converged, or returning to the step II to continue the process by iterative computation.
According to multiple times of calculation, the iteration is performed for 100 times, and the iteration result is converged. The reconstruction result is shown in FIG. 8.
4) Reconstruction quality discrimination
At this time, the reconstructed results of different combustion surface diameters are compared with the actual numerical values, fig. 9 shows that the error curve between the reconstructed numerical values and the actual numerical values of the combustion surfaces with different diameters can be obtained from the error curve, the error between the reconstructed numerical values and the actual numerical values is not more than 1mm, and the reconstructed quality of the combustion surfaces is proved to meet the requirements.

Claims (2)

1. The method for reconstructing the combustion surface retreating image of the solid rocket engine is characterized by comprising the following steps of:
1) magnetic line path for acquiring combustion surface microwave CT imaging field
According to the electromagnetic wave attenuation characteristic, battery magnetic wave tomography is simulated, 29 receiving devices further work simultaneously through one transmitting device, firstly, a magnetic line path of a preliminary microwave CT imaging field is obtained, then, according to the propagation rule of an electromagnetic field, 2 magnetic lines of force are added on two sides of each magnetic line of the magnetic line path of the preliminary obtained microwave CT imaging field respectively, so that the information quantity is increased, and the final magnetic line path of the combustion surface microwave CT imaging field is obtained;
2) acquiring microwave projection data of combustion surface
Solving a projection matrix according to a magnetic line simulation path equation, wherein a in the projection matrix AijDefining the length of an intersection line of i-ray and j-pixel, discretizing a reconstruction region according to a projection matrix, uniformly arranging n transmitting/receiving devices on a circumference taking the center of a measured object as a circle center, taking one of the n transmitting/receiving devices as a transmitter each time, switching a receiving state, and detecting scattered data together with all other receiving devices, wherein at the moment, the propagation of microwaves in a measured medium meets the wave equation in a common form:
Figure FDA0002573207880000011
wherein
Figure FDA0002573207880000012
Representing the reflectivity of a medium
Figure FDA0002573207880000013
The scalar function of (a) is selected,
Figure FDA0002573207880000014
for the vector laplacian, the formula can be converted into:
Figure FDA0002573207880000015
wherein
Figure FDA0002573207880000016
Called the objective function, kbRepresents the wave number, mu, of a plane wave in a background mediumrAnd εrIndicating relative permeability and relative to the background mediumDielectric constant due to total field
Figure FDA0002573207880000017
Can be expressed as an incident field
Figure FDA0002573207880000018
And a scattered field
Figure FDA0002573207880000019
And incident wave of
Figure FDA00025732078800000110
Can be defined as:
Figure FDA00025732078800000111
thus, it is possible to obtain:
Figure FDA00025732078800000112
with the green's function, the illumination produced by all sources can be found as:
Figure FDA0002573207880000021
thus obtaining "projection" data a;
3) combustion surface migration image reconstruction
To fully exploit the gradient sparsity of an image, a minimized image can be obtained by a formula, i.e.
Figure FDA0002573207880000022
Constraint conditions are as follows: subject to y is Ax
Wherein | x | purpleTVIs the TV norm of the image x,
Figure FDA0002573207880000023
representing an image function, from N image pixels x1,x2,L,xNThe composition, y is the attenuation value of the microwave passing through the combustion surface, A is a projection matrix,
the formula for the TV norm of image x is as follows,
Figure FDA0002573207880000024
wherein,
Figure FDA0002573207880000025
image reconstruction is carried out by taking minimization of TV norm as a target, and an iterative formula is as follows
Figure FDA0002573207880000026
Wherein p isiFor projecting the measured values,. mu.jFor the attenuation coefficient to be determined, λ is the relaxation factor, aijConstructing an image reconstruction algorithm based on TV norm minimization by performing convex set orthogonal projection and TV noise reduction twice correction processing on a calculation result in each step of iterative calculation process for the length of i microwave intersection j pixels;
4) reconstruction quality discrimination
Comparing the obtained reconstructed values of the combustion surfaces with different diameters with the actual values, finishing iteration when the error between the two values is not more than 1mm, proving that the reconstruction quality of the combustion surfaces meets the requirement, otherwise, continuing the iteration until the error between the two values is not more than 1mm,
δ=max|fi-f′i|
where δ is the error between the reconstructed value and the actual value, fiIs a true value of the combustion face diameter f 'of the image number i'iThe numerical value is the reconstructed value of the combustion surface diameter when the image serial number is i.
2. The method for reconstructing a combustion surface retrogradation image of a solid-rocket engine according to claim 1, wherein the algorithm comprises the following steps:
1) initial value selection: let S represent a star undersampled matrix, an initial observation set
Figure FDA0002573207880000031
Fourier transform coefficients representing the original image taken only at the star sample line position (the Fourier-centric slice theorem indicates that sampling along the star sample line in the frequency domain is inherently consistent with the random sampling required by CS theory)
Figure FDA0002573207880000032
Multiplying corresponding elements of the representation matrix, and further, making an initial value x of the iteration0=DFT-1(y0);
2) And (3) iterative calculation: solving and calculating the following iterative formula by using a conjugate gradient method:
Figure FDA0002573207880000033
wherein λ isnDenotes the iteration step size, dnRepresents the direction of the conjugate gradient;
3) orthogonal projection: orthogonal projection is performed on y as Ax, and calculation is performed first
Figure FDA0002573207880000034
DFT of (1), preserving the frequency domain information y on the initial sample point0The frequency domain information of the non-sampling point is updated by using the current iteration result without changing, and then the inverse DFT of the frequency domain updating result is calculated, namely
Figure FDA0002573207880000035
4) Denoising TV: the frequency-domain star sampling causes the iteration result to contain randomly distributed noise, so that the iteration result of the above formula is matched
Figure FDA0002573207880000036
Carrying out noise reduction treatment by using the following formula to obtain a final result x of the (n + 1) th iterationn+1I.e. by
Figure FDA0002573207880000037
Wherein for any given two-dimensional image x ═ x (x)ij) The TV noise reduction matrix is TV _ noise (x) ═ tij) And is and
Figure FDA0002573207880000038
5) and (3) updating the conjugation direction: search direction d of conjugate gradient methodnIs a combination of the direction of the negative gradient of the current iteration and the direction of the last search iteration, i.e.
dn+1=dTV(xn+1)+βdn
Wherein d isTV(xn+1) Express | | xn+1||TVIn the direction of the gradient of (a),
Figure FDA0002573207880000041
as a proportionality coefficient, | g | non-conducting phosphorFA Frobenius norm representing a matrix;
6) termination criteria: and (3) stopping iteration when the iteration times reach a certain set value or the iteration result is converged, or else, returning to 2) iterative computation to continue the process.
CN202010646466.4A 2020-07-07 2020-07-07 Solid rocket engine combustion surface retreating image reconstruction method Pending CN112200880A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010646466.4A CN112200880A (en) 2020-07-07 2020-07-07 Solid rocket engine combustion surface retreating image reconstruction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010646466.4A CN112200880A (en) 2020-07-07 2020-07-07 Solid rocket engine combustion surface retreating image reconstruction method

Publications (1)

Publication Number Publication Date
CN112200880A true CN112200880A (en) 2021-01-08

Family

ID=74006179

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010646466.4A Pending CN112200880A (en) 2020-07-07 2020-07-07 Solid rocket engine combustion surface retreating image reconstruction method

Country Status (1)

Country Link
CN (1) CN112200880A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115423890A (en) * 2022-09-15 2022-12-02 京心禾(北京)医疗科技有限公司 Tomographic image iterative reconstruction method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4662222A (en) * 1984-12-21 1987-05-05 Johnson Steven A Apparatus and method for acoustic imaging using inverse scattering techniques
CN105653754A (en) * 2014-12-05 2016-06-08 内蒙航天动力机械测试所 Solid rocket engine virtual test base system and testing method
CN208057259U (en) * 2018-02-08 2018-11-06 上海理工大学 Retire test system in a kind of Solid Rocket Engine Test three-dimensional combustion face
WO2019183584A1 (en) * 2018-03-23 2019-09-26 Memorial Sloan Kettering Cancer Center Deep encoder-decoder models for reconstructing biomedical images

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4662222A (en) * 1984-12-21 1987-05-05 Johnson Steven A Apparatus and method for acoustic imaging using inverse scattering techniques
CN105653754A (en) * 2014-12-05 2016-06-08 内蒙航天动力机械测试所 Solid rocket engine virtual test base system and testing method
CN208057259U (en) * 2018-02-08 2018-11-06 上海理工大学 Retire test system in a kind of Solid Rocket Engine Test three-dimensional combustion face
WO2019183584A1 (en) * 2018-03-23 2019-09-26 Memorial Sloan Kettering Cancer Center Deep encoder-decoder models for reconstructing biomedical images

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
刘玉等: "基于TTLS正则化方法的超声层析图像重建", 《电视技术》 *
史连艳等: "火箭弹MR/GPS制导系统应用仿真", 《系统仿真学报》 *
宫秀良等: "固体火箭发动机推力模拟装置研制", 《质量与可靠性》 *
柯丽等: "基于反投影的MIT动态图像重建方法研究", 《仪器仪表学报》 *
柯丽等: "基于滤波反投影的脑磁感应迭代重建算法研究", 《仪器仪表学报》 *
柳芳艳: "高质量光声图像重构方法研究", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115423890A (en) * 2022-09-15 2022-12-02 京心禾(北京)医疗科技有限公司 Tomographic image iterative reconstruction method
CN115423890B (en) * 2022-09-15 2023-09-19 京心禾(北京)医疗科技有限公司 Tomographic image iterative reconstruction method

Similar Documents

Publication Publication Date Title
CN108008385B (en) Interference environment ISAR high-resolution imaging method based on management loading
CN103713288B (en) Sparse Bayesian reconstruct linear array SAR formation method is minimized based on iteration
US7840053B2 (en) System and methods for tomography image reconstruction
Noghanian et al. Microwave tomography
CN105044722B (en) The full Bayesian Discriminating Features extracting method of synthetic aperture radar target
CN105388476B (en) A kind of chromatography SAR imaging methods based on joint sparse model
CN106918810B (en) A kind of microwave relevance imaging method when the amplitude phase error there are array element
CN115438571B (en) Ground wave radar wave field calculation method and device based on machine learning
CN116563096B (en) Method and device for determining deformation field for image registration and electronic equipment
CN111856362A (en) Magnetic resonance imaging method, device, system and storage medium
Shangguan et al. GF-3 polarimetric data quality assessment based on automatic extraction of distributed targets
Albert et al. Probabilistic direction-dependent ionospheric calibration for LOFAR-HBA
CN110554384A (en) imaging method based on microwave signal
Lee et al. Application of whale optimization algorithm to inverse scattering of an imperfect conductor with corners
CN112200880A (en) Solid rocket engine combustion surface retreating image reconstruction method
EP2817655A2 (en) Systems and methods for image sharpening
CN112773396A (en) Medical imaging method based on full waveform inversion, computer equipment and storage medium
CN110687535B (en) Rapid microwave imaging method
CN114690242B (en) Least square reverse time migration method with low noise
CN114137000B (en) Hybrid electromagnetic imaging algorithm based on Bayesian compressed sensing and Boen iteration
CN105608719A (en) Rapid CT image reconstruction method based on two-stage projection adjustment
Lei et al. Synthetic images augmentation for robust SAR target recognition
CN109633635B (en) Meter wave radar height measurement method based on structured recursive least squares
CN110146878A (en) The quantitative microwave imaging method that background media iteration based on the regularization of multiplying property updates
Hopp et al. Fat ray ultrasound transmission tomography: preliminary experimental results with simulated data

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
AD01 Patent right deemed abandoned

Effective date of abandoning: 20231017

AD01 Patent right deemed abandoned