阅读说明：本技术 一种组合扫描cl成像方法 (Combined scanning CL imaging method ) 是由 刘丰林 席雅睿 周鹏伍 于 2021-07-12 设计创作，主要内容包括：本发明涉及一种组合扫描CL成像方法,属于扫描成像技术领域。该方法包括：S1：构建直线-圆周组合扫描CL成像系统,并获得几何成像模型；S2：直线扫描：确定直线扫描等效投影角度θ,并计算出X射线源和平板探测器的步进距离,以及系统放大比K；S3：圆周扫描：确定圆周扫描偏角并计算出X射线源的移动距离和平板探测器的移动距离；S4：分别采集一组直线扫描和圆周扫描过程中的一组投影数据,并利用SIRT迭代算法对采集的投影数据进行图像重建。本发明解决了传统扫描CL成像方法数据缺失、纵向空间分辨力缺失等问题,提高投影数据的成像质量,同时使得成像设备体积足够小。(The invention relates to a combined scanning CL imaging method, belonging to the technical field of scanning imaging. The method comprises the following steps: s1: constructing a linear-circumferential combined scanning CL imaging system and obtaining a geometric imaging model; s2: linear scanning: determining a linear scanning equivalent projection angle theta, and calculating the step distance between an X-ray source and a flat panel detector and the system amplification ratio K; s3: circumferential scanning: determining a circular scan declination Calculating the moving distance of the X-ray source and the moving distance of the flat panel detector; s4: and respectively acquiring a group of projection data in the processes of a group of linear scanning and a group of circular scanning, and carrying out image reconstruction on the acquired projection data by using a SIRT iterative algorithm. The invention solves the problems of data loss, longitudinal spatial resolution loss and the like of the traditional scanning CL imaging method, improves the imaging quality of projection data and simultaneously ensures that the volume of imaging equipment is small enough.)

1. A combined scanning CL imaging method is characterized by comprising the following steps:

s1: constructing a linear-circumferential combined scanning CL imaging system and obtaining a geometric imaging model;

s2: linear scanning: determining a linear scanning equivalent projection angle theta, and calculating the step distance between an X-ray source and a flat panel detector and the system amplification ratio K;

s3: circumferential scanning: determining a circular scan declinationCalculating the moving distance of the X-ray source and the moving distance of the flat panel detector;

s4: and respectively acquiring a group of projection data in the processes of a group of linear scanning and a group of circular scanning, and carrying out image reconstruction on the acquired projection data by using a SIRT iterative algorithm.

2. The combined scan CL imaging method of claim 1, wherein in step S1, the constructed linear-circumferential combined scan CL imaging system comprises: the X-ray source, the flat panel detector and the detection object; the X-ray source is positioned below the detection object and emits X-rays; the detection object is positioned above the X-ray source; the flat panel detector is positioned above the detection object and receives the X-rays attenuated after passing through the detection object.

3. The combined-scan CL imaging method as claimed in claim 1, wherein in step S1, obtaining a geometric imaging model specifically includes: the position of the X-ray source target point at any moment is recorded as y_SObtaining a geometric imaging model as follows:

R＝S_O·tanγ

wherein, omega is the included angle between any projection ray and the x-y plane, gamma is the included angle between any two rays with the same omega angle, alpha is the included angle between any two rays with the same omega angle and the projection central ray on the x-y plane, and S_OIn order to detect the distance from the object to the track of the X-ray source, beta is the included angle between the ray and the ray in the projection center, t is the distance from the projection center to the ray, and t is the element [ -R, R]And R is the field radius.

4. The combined-scan CL imaging method according to claim 3, characterized in that in step S2, the linear scan specifically includes: the X-ray source and the detector perform relatively parallel linear scanning on a plate-shaped detection object along a certain linear direction; determining the equivalent projection angle theta of linear scanning, wherein the theta belongs to [0, pi ],the total scanning stroke L of the X-ray source_SComprises the following steps:

meanwhile, the total movement stroke L of the flat panel detector_DComprises the following steps:

and setting the number p of required sampling points before scanning, wherein the sampling unit angle of the linear equiangular scanning is as follows: Δ θ ═ θ/p; then the X-ray source step distance ay_SiComprises the following steps:

step distance delta y of flat panel detector_DiComprises the following steps:

wherein i is the number of sampling points scanned by the X-ray source or the flat panel detector;

the system amplification ratio K is:

5. the combined-scan CL imaging method according to claim 3, characterized in that in step S3, the circular scan includes: the X-ray source and the flat panel detector move to the designated positions, and 360-degree projection data are acquired through the indexing movement of the detection object; determining a circular scan declinationThe moving distance of the X-ray sourceIs far from L_S1Comprises the following steps:

similarly, the moving distance L of the flat panel detector_D1Comprises the following steps:

wherein S is_DThe distance from the track of the X-ray source to the track of the flat panel detector;

setting the number of sampling points of circumferential scanning as p, and then detecting the rotation indexing angle delta a of the object as 2 pi/p during sampling;

the system is adjusted by S_O、S_DThe FOV is adjusted.

6. The combined scan CL imaging method of claim 1, wherein in step S4, the acquired projection data is a linear matrix equation AX ═ b, whereA projection measurement matrix is adopted, and M is the total number of rays;for reconstructing the object, N is the total number of voxel points; a ═ a_mn) A system measurement matrix, M1., M, N1., N;

the method for reconstructing the image of the acquired projection data by using the SIRT iterative algorithm specifically comprises the following steps:

s41: calculating a correction term of an equation corresponding to the first ray for each pixel point, and storing the correction term in an array; calculating a correction term of an equation corresponding to the second ray for each pixel point, and adding the correction term into the array; until the correction term of the equation corresponding to the last ray to each voxel point is calculated and added into the array, one-time iteration updating utilizes all projection equations to complete the updating processing of the iteration solution under all projection angles;

s42: step S41 is applied to the projection data acquired in steps S2 and S3 until the reconstructed image satisfies the criterion requirement.

7. The combined scan CL imaging method of claim 6, characterized in that in step S41, the update processing step of the iterative solution at one projection angle is:

the iterative formula isWherein λ^kIs a relaxation factor for suppressing overcorrection, k is the number of iterations, i is 1. j 1.. N, N is the total number of voxel points, p_iIs the projection value of the ith ray,is an estimate of the ith ray, a_ijIs a projection coefficient reflecting the contribution of the jth voxel point to the ith ray; the iteration process specifically comprises the following steps:

input projection data p_iAnd an initial value is given to the user,whereinRepresenting an initial value of a jth voxel;

calculating estimated projection values of all rays:

thirdly, calculating the error of the image to be processed,

and (iv) repeating the steps S2 and S3 until all projection angles are finishedAfter error calculation for all equations, cumulative summation is performed:

calculating the correction value of the jth unknown quantity:

sixthly, carrying out the correction,performing one-time correction on all voxel points of the reconstructed image and then completing one-time iteration;

and seventhly, taking the result of the iteration as a temporary solution, wherein k is k +1, and repeating the steps from the step two to the step six until the criterion requirement is met.

8. The combined scanning CL imaging method according to one of the claims 1 to 7, characterized in that an imaging system suitable for the method comprises: mechanical systems, computer control systems, X-ray sources, and flat panel detectors.

Technical Field

The invention belongs to the technical field of scanning imaging, and relates to a combined scanning CL imaging method.

Background

A Computed Tomography (CT) technique is a nondestructive testing technique capable of effectively testing three-dimensional structural information inside an object, and is widely applied to various fields, such as the petroleum and natural gas field, the intelligent manufacturing field, the medical field, and the like, because of its advantages of high visualization degree, high resolution, no damage, and the like. When detecting the objectWhen the body is a circuit board, a chip, an airplane wing plate and other plate-shaped components (the length and the width are far larger than the thickness), the CT technology has certain limitations due to the limitation of the geometric structure: on one hand, for a structure with a special shape, the conventional rotational scanning method of CT has a distance S from the X-ray source to the rotation center of the object to be detected_OMay be large, causing a reduction in the spatial resolution of the image; on the other hand, the beam is severely attenuated in the direction in which the physical size of the object to be detected is large, and the projection data has a poor signal-to-noise ratio. For this reason, a computed tomography (CL) technique is studied and developed for the radiographic detection of the plate-shaped member.

In 1916, a classical layered imaging (CL) method was proposed by Andre Bocage, a french dermatologist, which enables film imaging of sections at specific depths; in 1932, the first layered imaging experiments were performed by D planes, which demonstrated that theoretically it was possible to obtain a multilayer radiation cross-section of an object by a series of X-ray irradiations. Early CL imaging typically incorporated a camera and image intensifier, and had slow imaging speed, low signal-to-noise ratio, and only one layer per scan. With the advent of digital detectors in the last 80's, CL was rapidly developed by applying appropriate offsets before superimposing images to obtain multi-layer sharp digital images (i.e., tomosynthesis techniques). Currently, CL technologies are mainly classified into linear type, circular type, swing type, C-arm type, and the like according to a system scanning manner.

In 2002, based on a micro-focus X-ray imaging system, lie politics and the like propose an oscillating CL scanning mode aiming at thin-layer structure objects such as multilayer large-area composite materials and the like, the scanning speed is high, and the real-time performance is high; in 2010, MAISL et al also mentioned wiggle-type CL of the same scanning structure for non-destructive inspection of large plate-like members. In a swing-type CL scanning configuration, projection data is acquired by swinging the object back and forth over a range of angles (less than 180 °), which is achievable with standard CT scanners, which are essentially a finite angle CT. The swing CL has the advantages that the detector and the X-ray source are kept relatively static in the scanning process, and the mechanical manufacturing difficulty is reduced; the limitation is that the image magnification varies during the scan, and thus the image cannot be reconstructed using the tomosynthesis method.

In 2010, Fu et al propose a large-field CL imaging method with an asymmetric rotational scanning structure, and the rotation offset angle is adoptedThe object to be detected or the rotating X-ray source and the offset detector to acquire projection data, and it is proved that the method can expand the imaging area and improve the imaging spatial resolution. The German IZFP company develops another scanning mode of rotary CL for the nondestructive detection of large heavy objects, a ray source and a detector rotate relative to a central point, and the system has the advantages that only an object rotates in the scanning process, the X ray source and the detector synchronously realize offset through a swinging frame, the space requirement is small, the illumination intensity of the detector is always the same as that of a plane type detector, the calibration of the detector at any sample position is not needed, and the system can detect the object which weighs 300kg and has the width of 140m at maximum. In 2015, Liu and the like, a Chinese academy of sciences, developed an industrial computer tomography (ICL) system, in which a detector is mounted on a C-shaped arm and can rotate 360 degrees around a z-axis, a detection object is placed on a turntable with three degrees of freedom, and the system can flexibly adjust spatial resolution, imaging field of view size and region of interest size, and has multiple imaging modes. The above method or system needs to tilt the detected object (or detector, X-ray source) to a certain angle with the z-axis, and the space structure of the whole system is complex.

In 2018, wangsu et al propose a simple and fast relative parallel linear scanning CL (PTCL) system, where the detection object is placed on the middle plane, and the detector and the X-ray source move in anti-parallel to obtain projection data; comparing an FDK three-dimensional analytic reconstruction algorithm with a combined algebraic reconstruction technology (SART + TV) based on image total variation minimization by taking a chip as a detection object, proving that the system can realize high-quality image reconstruction based on truncated projection data; on the basis, an orthogonal linear translation scanning (OTCL) system is proposed by ran et al, which realizes orthogonal scanning by rotating a detection object, and solves the problem of finite angle projection artifact existing in the PTCL to a certain extent; in 2020, field faithful construction and the like are combined with simulation and actual experiments, and a Simultaneous Iterative Reconstruction Technology (SIRT) is adopted to compare a PTCL with an OTCL, so that the higher image resolution and fewer artifacts of the OTCL are proved under the same scanning angle no matter global reconstruction or local reconstruction is adopted, but the image reconstruction quality is still to be improved. The planar linear CL system has the advantages of high scanning speed, simple system structure and certain limitation.

Although the above CL imaging methods have been used in various applications, the CL imaging quality still needs to be improved.

Disclosure of Invention

In view of the above, the present invention provides a combined scanning CL imaging method, which adopts a linear-circumferential combined scanning imaging manner, solves the problems of data loss, longitudinal spatial resolution loss, and the like in the conventional scanning CL imaging method, improves the imaging quality of projection data, and makes the volume of an imaging device sufficiently small.

In order to achieve the purpose, the invention provides the following technical scheme:

a combined scanning CL imaging method specifically comprises the following steps:

s1: constructing a linear-circumferential combined scanning CL imaging system and obtaining a geometric imaging model;

s2: linear scanning: determining a linear scanning equivalent projection angle theta, and calculating the step distance of an X-ray source and a flat panel detector and the system amplification ratio K;

s3: circumferential scanning: determining a circular scan declinationCalculating the moving distance of the X-ray source and the moving distance of the flat panel detector;

s4: and respectively acquiring a group of projection data in the processes of a group of linear scanning and a group of circular scanning, and carrying out image reconstruction on the acquired projection data by using a SIRT iterative algorithm.

Further, in step S1, the linear-circumferential combined scanning CL imaging system is constructed, including: the X-ray source, the flat panel detector and the detection object; the X-ray source is positioned below the detection object and emits X-rays; the detection object is positioned above the X-ray source; the flat panel detector is positioned above the detection object and receives the X-rays attenuated after passing through the detection object.

Further, in step S1, obtaining a geometric imaging model specifically includes: the position of the X-ray source target point at any moment is recorded as y_SAnd the distance between the flat panel detector and the center of the view field in the y-axis direction is recorded as y_D(ii) a Omega is the included angle between any projection ray and the x-y plane, gamma is the included angle between any two rays with the same omega angle, alpha is the included angle between any two rays with the same omega angle and the projection central ray on the x-y plane, and S_OIn order to detect the distance from the object to the track of the X-ray source, beta is the included angle between the ray and the ray in the projection center, t is the distance from the projection center to the ray, and t is the element [ -R, R]R is the radius of the field of view; the geometric imaging model is obtained as follows:

R＝S_O·tanγ

further, in step S2, the linear scanning specifically includes: the X-ray source and the detector perform relatively parallel linear scanning on a plate-shaped detection object along a certain linear direction (taking the y direction as an example in the figure); determining the equivalent projection angle theta of linear scanning, theta belongs to [0, pi), and then the total scanning stroke L of the X-ray source_SComprises the following steps:

meanwhile, the total movement stroke L of the flat panel detector_DComprises the following steps:

the line scan may be divided into an equidistant scan and an equiangular scan, and the equiangular scan is used here to ensure uniformity of the projected data. Setting the number p of sampling points before scanning, linear equiangular scanningThe unit angle of the sampling is as follows: Δ θ ═ θ/p; the X-ray source step distance ay_SiComprises the following steps:

step distance delta y of flat panel detector_DiComprises the following steps:

wherein i is the number of sampling points scanned by the X-ray source or the flat panel detector;

the system amplification ratio K is:

further, in step S3, the circular scanning includes: the X-ray source and the flat panel detector move to the designated positions, and 360-degree projection data are acquired through the indexing movement of the detection object; determining a circular scan declinationThe moving distance L of the X-ray source_S1Comprises the following steps:

similarly, the moving distance L of the flat panel detector_D1Comprises the following steps:

wherein S is_DThe distance from the track of the X-ray source to the track of the flat panel detector;

setting the number of sampling points of circumferential scanning as p, and then detecting the rotation indexing angle delta a of the object as 2 pi/p during sampling;

the system is adjusted by S_O、S_DThe FOV is adjusted.

Further, in step S4, the acquired projection data is the linear matrix equation AX ═ b, whereFor the projection measurement matrix, M is the total amount of data (or total number of rays);for reconstructing the object, N is the total number of voxel points; a ═ a_mn) A system measurement matrix, M1., M, N1., N;

the method for reconstructing the image of the acquired projection data by using the SIRT iterative algorithm specifically comprises the following steps:

s41: calculating a correction term of an equation corresponding to the first ray for each pixel point, and storing the correction term in an array; calculating a correction term of an equation corresponding to the second ray for each pixel point, and adding the correction term into the array; until the correction term of the equation corresponding to the last ray to each element point is calculated and added into the array, one-time iteration updating utilizes all projection equations to complete the updating processing of the iteration solution under all projection angles;

s42: step S41 is applied to the projection data acquired in steps S2 and S3 until the reconstructed image meets certain criteria requirements.

Further, in step S41, the update processing step of the iterative solution under one projection angle is:

the iterative formula isWherein λ^kIs a relaxation factor for suppressing overcorrection, k is the number of iterations, i is 1. j 1.. N, N is the total number of voxel points, p_iIs the projection value of the ith ray,is an estimate of the ith ray, a_ijIs a projection coefficient reflecting the jth voxel point to the ith rayA contribution; the iterative process specifically comprises:

input projection data p_iAnd an initial value is given to the user,whereinRepresenting an initial value of a jth voxel;

calculating estimated projection values of all rays:

thirdly, calculating the error of the image to be processed,

and (4) repeating the steps S2 and S3 until the error calculation of all equations under all projection angles is completed, and performing cumulative summation:

calculating the correction value of the jth unknown quantity:

sixthly, carrying out the correction,performing one correction on all voxel points of the reconstructed image to complete one iteration;

and seventhly, taking the result of the iteration as a temporary solution, wherein k is k +1, and repeating the steps from the step two to the step six until the criterion requirement is met.

Further, an imaging system suitable for use in the method includes: mechanical system, computer control system, X-ray source and flat panel detector.

The invention has the beneficial effects that: the method adopts a linear-circumferential combined scanning imaging mode, solves the problems of data loss, longitudinal spatial resolution loss and the like of the traditional scanning CL imaging method, improves the imaging quality of projection data, and simultaneously ensures that the imaging equipment has small volume.

Additional advantages, objects, and features of the invention will be set forth in part in the description which follows and in part will become apparent to those having ordinary skill in the art upon examination of the following or may be learned from practice of the invention. The objectives and other advantages of the invention may be realized and attained by the means of the instrumentalities and combinations particularly pointed out hereinafter.

Drawings

For purposes of promoting a better understanding of the objects, aspects and advantages of the invention, reference will now be made to the following detailed description taken in conjunction with the accompanying drawings in which:

FIG. 1 is a block diagram of a combined scanning CL imaging system of the invention;

FIG. 2 is a schematic diagram of a line-circle combination scan CL system;

FIG. 3 is a combined scan geometry model of the system of the present invention at any time;

FIG. 4 is a rectilinear, circular scan geometric model;

FIG. 5 is a 32 nd layer reconstructed image of a variety of CL scan simulations;

FIG. 6 is a diagram of a layer 32 vertical center waveform;

FIG. 7 is a schematic diagram of the motion of the combined scanning CL imaging system of the invention;

FIG. 8 is a block diagram of a combined scanning CL imaging system according to the invention;

reference numerals: the system comprises a system framework 1, a flat panel detector z-direction movement mechanism 2, a flat panel detector y-direction movement mechanism 3, a flat panel detector 4, a detection object 5, an objective table 6, a wheel disc 7, an X-ray source y-direction movement mechanism 8, an X-ray source 9 and an X-ray source z-direction movement mechanism 10.

Detailed Description

The embodiments of the present invention are described below with reference to specific embodiments, and other advantages and effects of the present invention will be easily understood by those skilled in the art from the disclosure of the present specification. The invention is capable of other and different embodiments and of being practiced or of being carried out in various ways, and its several details are capable of modification in various respects, all without departing from the spirit and scope of the present invention. It should be noted that the drawings provided in the following embodiments are only for illustrating the basic idea of the present invention in a schematic way, and the features in the following embodiments and examples may be combined with each other without conflict.

Wherein the showings are for the purpose of illustrating the invention only and not for the purpose of limiting the same, and in which there is shown by way of illustration only and not in the drawings in which there is no intention to limit the invention thereto; to better illustrate the embodiments of the present invention, some parts of the drawings may be omitted, enlarged or reduced, and do not represent the size of an actual product; it will be understood by those skilled in the art that certain well-known structures in the drawings and descriptions thereof may be omitted.

The same or similar reference numerals in the drawings of the embodiments of the present invention correspond to the same or similar components; in the description of the present invention, it should be understood that if there is an orientation or positional relationship indicated by terms such as "upper", "lower", "left", "right", "front", "rear", etc., based on the orientation or positional relationship shown in the drawings, it is only for convenience of description and simplification of description, but it is not an indication or suggestion that the referred device or element must have a specific orientation, be constructed in a specific orientation, and be operated, and therefore, the terms describing the positional relationship in the drawings are only used for illustrative purposes, and are not to be construed as limiting the present invention, and the specific meaning of the terms may be understood by those skilled in the art according to specific situations.

Referring to fig. 1-8, the present invention discloses a linear-circumferential combined scanning CL imaging system and method, as shown in fig. 1-2, the data acquisition mode of the imaging system is based on the X-ray source/detector relative parallel linear movement scanning and the object rotation around the z-axis circumferential scanning. As shown in FIG. 2, the X-ray source is located at the lowermost end; the object stage is positioned above the X-ray source and used for placing a plate-shaped detection object; the flat panel detector is positioned right above the object stage and used for receiving the attenuated X-rays.

1. The scanning geometric model of the linear-circumferential combined scanning CL imaging system is shown in figures 3-4, and the system consists of an X-ray source S, a flat panel detector D and a detection object. The X-ray source is positioned below the detection object and emits X-rays; the detection object is positioned above the X-ray source; the flat panel detector is positioned above the detection object and receives the X-rays attenuated after passing through the detection object.

As shown in FIG. 3, the position of the target point of the X-ray source at any time is marked as y_SAnd the distance between the detection unit and the center of the view field in the y-axis direction is recorded as y_D. Omega is the included angle between any projection ray and the x-y plane, gamma is the included angle between any two rays with the same omega angle, alpha is the included angle between any two rays with the same omega angle and the projection central ray on the x-y plane, S_OFor detecting the distance of the object to the trajectory of the X-ray source, S_DIs the distance from the X-ray source track to the flat panel detector track, beta is the included angle between the ray and the projection center ray, t is the distance from the projection center to the ray, and t is the element [ -R, R]And R is the field radius. From the imaging model, one can obtain:

R＝S_O·tanγ (2)

(1) straight line scanning process

The X-ray source and the detector perform a relatively parallel linear scan on a plate-like object to be detected along a linear direction (the y direction is taken as an example in the figure). Determining the equivalent projection angle theta (shown in figure 4) of the linear scanning, theta epsilon [0, pi), and the total scanning stroke L of the X-ray source_SComprises the following steps:

meanwhile, the total movement stroke L of the flat panel detector_DComprises the following steps:

line scanningIt can be divided into isometric scan and equiangular scan, and in order to ensure the uniformity of the projected data, equiangular scan is adopted here, as shown in fig. 4. And setting the number p of required sampling points before scanning, wherein the sampling unit angle of the linear equiangular scanning is as follows: and delta theta is theta/p. Further, the X-ray source step distance Δ y_SiComprises the following steps:

step distance delta y of flat panel detector_DiComprises the following steps:

the system amplification ratio K is:

(2) circular scanning process

The X-ray source and the flat panel detector move to the designated positions, and 360-degree projection data are acquired by detecting the indexing movement of the object. Determining a circular scan declination(see fig. 4), the X-ray source moves by a distance L_S1Comprises the following steps:

similarly, the moving distance L of the flat panel detector_D1Comprises the following steps:

if the number of sampling points in the circular scan is p, the indexing angle Δ a of the rotation of the detection object during sampling is 2 pi/p.

The system is adjusted by S_O、S_DThe FOV is adjusted.

Compared with the linear scanning CL and the circumferential scanning CL, the linear-circumferential combined scanning CL imaging method provided by the invention has the following characteristics: the circular scanning can obtain equal CL image resolving power in all directions, and the linear scanning can carry out linear interpolation on the data set which is missing in the Fourier space by the circular scanning.

2. Image reconstruction algorithm

The CL rotary scanning mechanism differs from the CT in that the center line of the CL ray beam does not coincide with the central axis of rotation of the object, and there is a deflection angle smaller than 90 °, and the CL technique is essentially a finite angle CT technique for non-coaxial scanning. The CL can be regarded as a generalized case of the CT, and the image reconstruction algorithm suitable for the CT is also suitable for CL image reconstruction, and because the incomplete projection data causes a limited angle artifact in the reconstructed image from the source, the image reconstruction algorithm is especially important for improving the quality of the CL reconstructed image. The CL image reconstruction algorithms are mainly classified into four categories: shift and Add (SAA), analytic methods, iterative methods, and maximum likelihood statistical methods. The SAA obtains a required focusing plane through translation and direct superposition projection, and has the defects of fast calculation and poor reconstruction quality; the analytical method (FBP/FDK) has high reconstruction speed and good reconstruction quality, and has the defect that a high-quality image can be reconstructed only by complete data; the iterative method (ART/SART/SIRT) can recover images from a small amount of data with low signal-to-noise ratio, and can utilize regularization and prior information, and has the defect of higher calculation cost; the maximum likelihood statistical method can obtain high-quality images for incomplete data, but the calculation cost is high.

In order to verify the effectiveness of the method provided by the invention, an iterative reconstruction algorithm is considered. The classical iterative reconstruction algorithm is ART, i.e. an algebraic reconstruction technique, which corrects the value of each voxel point by adding a correction term during the iterative computation of image reconstruction, and adopts a ray-by-ray update mode, wherein each ray is computed, and all voxel values related to the ray are updated once; the SART algorithm, namely a joint algebraic reconstruction technology, updates the intermediate solution by a joint correction term method under a specific projection angle. A joint correction term, namely a correction term generated by all rays passing through a certain voxel under a specific projection angle; the SIRT algorithm, i.e., the simultaneous iterative reconstruction technique, updates the intermediate solution at all angles by means of a joint correction term, where the joint correction term is a correction term generated by all rays passing through a voxel at all angles.

The method mainly comprises the following steps:

step 1: determining the scanning parameters theta,K；

Step 2: linear scanning, namely determining the movement strokes of an X-ray source and a detector according to the angle theta, wherein the X-ray source and the detector do relative parallel linear movement, the X-ray source emits X rays in the process, the flat panel detector receives the X rays attenuated by a detection object and acquires a group of projection data information with a certain included angle;

step 3: performing circular scanning, namely moving an X-ray source and a detector to specified positions according to a deflection angle phi, wherein a detected object is driven by a wheel disc to perform circular indexing motion in the whole scanning process, the X-ray source emits X rays, and a flat panel detector receives the X rays attenuated by the detected object and acquires a group of projection data information with a certain deflection angle;

step 4: image reconstruction was performed using an iterative SIRT algorithm using the projection data acquired at Step2 and Step 3.

The system modeling may be modeled as a linear matrix equation AX ═ b, whereFor the projection measurement matrix, M is the total amount of data (or total number of rays);for reconstructing the object, N is the total number of voxel points; a ═ a_mn) For the system measurement matrix, M1.., M, N1.., N.

The SIRT algorithm for image reconstruction specifically comprises the following steps:

step 41: calculating a correction term of an equation corresponding to the first ray for each pixel point, and storing the correction term in an array; calculating a correction term of an equation corresponding to the second ray for each voxel point, and adding the correction term into the array; until the correction term of the equation corresponding to the last ray to each voxel point is calculated and added into the array, one-time iterative updating utilizes all projection equations, and the updating processing of the iterative solution under all projection angles is completed;

step 42: step41 is applied to the projection data acquired at Step2 and Step3 until the reconstructed image meets certain criteria.

The update processing of the iterative solution under one projection angle is specifically as follows:

the iterative formula isWherein λ^kIs a relaxation factor for suppressing overcorrection, k is the number of iterations, i is 1. j 1.. N, N is the total number of voxel points, p_iIs the projection value of the ith ray,is an estimate of the ith ray, a_ijIs a projection coefficient reflecting the contribution of the jth voxel point to the ith ray; the iterative process specifically comprises:

input projection data p_iAnd an initial value is given to the user,whereinRepresenting an initial value of a jth voxel;

calculating estimated projection values of all rays:

thirdly, calculating the error of the image to be processed,

and (4) repeating the steps S2 and S3 until the error calculation of all equations under all projection angles is completed, and performing cumulative summation:

calculating the correction value of the jth unknown quantity:

sixthly, carrying out the correction,performing one correction on all voxel points of the reconstructed image to complete one iteration;

and seventhly, taking the result of the iteration as a temporary solution, wherein k is k +1, and repeating the steps from the step two to the step six until the criterion requirement is met.

The above-mentioned modifications are made only for those voxels traversed by the ray, for which the ray does not traverse, a_ijNot greater than 0, so C_j0 corresponds to no correction.

3. In order to verify the feasibility of the method, the simulated reality test is carried out on an MATLAB platform by using a simulated printed circuit board. In the experiment, single-segment linear scanning, orthogonal scanning and linear-circumferential combined scanning are respectively carried out on a simulation die body of 256 × 64 in a cone beam scanning mode, and an image is reconstructed by using an SIRT algorithm in the experimental process. The reconstruction result is shown in fig. 5, and the simulation experiment parameters are shown in table 1 below:

TABLE 1 simulation experiment parameters

Fig. 5 is a diagram showing the reconstruction result of the 32 th layer of the analog printed circuit board phantom, wherein (a) the diagram is an original image, (b) the diagram is a single-segment linear scanning reconstructed image, (c) the diagram is an orthogonal linear scanning reconstructed image, (d) the diagram is a circular scanning reconstructed image, and (e) the diagram is a linear-circular combined scanning reconstructed image. Comparing the four images, (b) the images have a plurality of aliasing artifacts, and the reconstruction effect is the worst; (c) aliasing artifacts exist at the edge part of the image, but the image can be basically reconstructed; (d) no apparent aliasing artifact appears in the graph, but the edge information in the original graph is lost, as indicated by the arrow of the graph; (e) the visual effect of the map is best, no aliasing artifacts are evident, and the edge information remains intact.

In order to investigate the performance of each scanning method, the vertical center waveform of the reconstructed image of layer 32 in the three-dimensional reconstruction result graph is further compared, as shown in fig. 6, the waveforms of each scanning method are compared to find that: the wave form graph result of the orthogonal linear scanning CL reconstruction result is slightly better than that of the linear scanning CL, the wave form graph of the circumferential scanning CL reconstruction result is the worst, and the coincidence degree of the wave form graph of the linear-circumferential combined CL reconstruction result and the wave form graph of the original image is the highest.

As can be seen from simulation experiments, compared with a single-segment linear scanning CL imaging method, the orthogonal linear scanning CL imaging method obtains data in the transverse direction and the longitudinal direction, the projection data volume is richer, and the image reconstruction result is better. The circular scan CL imaging method effectively suppresses aliasing artifacts, but loses some detail information compared to other CL scan imaging methods. Compared with the orthogonal linear scanning CL imaging method and the circular scanning CL imaging method, the linear-circular combined scanning CL imaging method has richer data volume, almost no aliasing artifact can be seen in a reconstructed image, and detailed information is rich. The linear-circumferential combined scanning CL imaging method solves the problems of data loss, longitudinal spatial resolution loss and the like of the traditional scanning CL imaging method, and the feasibility of the method is verified through theoretical simulation imaging.

4. Designing a combined scanning CL imaging system

The system mainly comprises the following components: the system comprises a plurality of mechanical system assemblies, a computer control system, an X-ray source and a flat panel detector, and can meet the detection requirements of the plate-shaped component, such as the operations of adjusting the size of the FOV, adjusting the magnification ratio and the like. As shown in FIG. 7, the mechanical structure of the system is simple, and the system is generally divided into an upper layer structure, a middle layer structure and a lower layer structure. The flat panel detector is arranged above the integral structure, has two moving freedom degrees in the y direction and the z direction, and receives the attenuated X-rays; the detection object is positioned on the middle layer and is arranged on the object stage, the object stage is fixed with the wheel disc, the object stage is driven by the wheel disc to rotate around the axis z, namely, the object can rotate 360 degrees on the plane of the object stage; the X-ray source is arranged below the whole structure, has two moving degrees of freedom in the y direction and the z direction, and emits cone beam X-rays upwards. M1 is the flat panel detector moving along the z direction, M2 is the flat panel detector moving along the y direction, M3 is the detection object rotating around the z axis, M4 is the X ray source moving along the y direction, and M5 is the X ray source moving along the z direction.

In order to reduce the motion error of projection data acquisition in the motion process, a high-precision servo motor is selected by a system motor. The outermost periphery is provided with a radiation protection device to prevent X-rays from harming testers in the operation process of the equipment.

Finally, the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting, although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions may be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions, and all that should be covered by the claims of the present invention.