Method for reconstructing an object by X-ray tomography including a calculation of regularization terms.

Abstract

The invention relates to a method for reconstructing an object by X-ray tomography using X-ray sources and detectors. The method comprises the following steps: obtaining (E0) the positions of the X-ray sources and detectors; pre-calculating (E1) a set of directions; calculating (E2) a tomographic incompleteness vector map; calculating (E3) a tomographic incompleteness map of the metrics; calculating (E4) regularization terms from the tomographic incompleteness vector map and the tomographic incompleteness map of the metrics; and iteratively solving (E5) a regularized tomographic reconstruction problem using the calculated regularization terms. Figure for the abstract: Fig. 2

Description

Title of the invention: Method for reconstructing an object by X-ray tomography comprising a calculation of regularization terms.

[0001] The invention falls within the field of X-ray imaging, and more particularly in the field of X-ray tomographic reconstruction. It relates to a method for reconstructing an object by X-ray tomography comprising a calculation of regularization terms.

[0002] The invention finds its application in the field of security, in particular for the search of baggage taking place in sensitive sites (airports, stations, etc.), during events or on any occasion of baggage inspection. It also finds its application in the medical field, for example when using X-ray medical imaging equipment. More generally, the invention finds its application when it is desired to reconstruct an object by tomography.

[0003] X-ray tomography provides a two-dimensional or three-dimensional image from a series of tomographic projections via a reconstruction algorithm. The reconstruction algorithms used in tomography have weaknesses when the measured data are highly noisy and / or missing. The reconstruction algorithm then faces the lack of data and can then lead to very inconsistent solutions.

[0004] In order to be able to reconstruct an object, the reconstruction algorithm then requires introducing, in addition to input data, a priori on the object to be reconstructed. The a priori on the object are introduced via a regularization thus giving rise to a regularized reconstruction algorithm.

[0005] Regularized reconstruction algorithms thus solve an inverse problem from projection data. This family of algorithms requires the adjustment of regularization terms. These regularization terms aim to provide a compromise between the information in the data and the a priori on the object to be reconstructed. Regularization makes it possible to introduce an a priori on the statistical distribution of the imaged object and / or take into account the acquisition geometry. This regularization depends on parameters called hyper-parameters here. The hyper-parameters are calculated before the object reconstruction process, as detailed later.

[0006] Most of the time, regularization terms are constructed to limit the number of hyper-parameters to be adjusted. Indeed, the greater the number of hyper-parameters, the more difficult it is to adjust them manually.

[0007] To automatically adjust hyper-parameters, some approaches rely on on cross-validation methods [Kohavi, 1995], or more recently on machine learning methods [Shen, 2018], [Ding, 2021]. In 2021, work [Zhang, 2021] showed that a Directional Total Variation type regularization is well suited for limited angle reconstruction problems. This regularization nevertheless requires tuning more hyper-parameters than a conventional Total Variation [Rudin, 1992].

[0008] The methods currently used to adjust the regularization terms require having reference images, and / or having a large number of diversified data. While it is possible to obtain the reference images by simulations, these methods remain difficult to obtain under experimental conditions.

[0009] Furthermore, ideally, the choice of regularization terms should integrate the geometry of the scanner, i.e. the position of the source(s) and the detector(s).

[0010] 3D computed tomography theory has some tools to provide local information about the object to be reconstructed according to the scanner geometry. Assuming untruncated projections, [Tuy, 1983] gives a condition to check whether a continuous X-ray source trajectory is sufficient to reconstruct the object. The condition can be written as follows: "every plane that intersects the imaged object must intersect the source trajectory at least once." In practice, all scanners are limited by a finite number of source positions and therefore Tuy's condition is never met. We speak of tomographic incompleteness when Tuy's condition is not met.

[0011] [Fig.l] shows examples of source-detector trajectories T of a scanner used for the reconstruction of an object. Example a) shows a circular source-detector trajectory. Such a trajectory does not allow an exact reconstruction of the object outside the plane of the trajectory. Indeed, a plane intersecting the object and which is parallel for example to the circle of the trajectory T does not intersect said source-detector trajectory. A helical geometry of the trajectory T, as shown in example b), allows for an exact reconstruction to be obtained. Indeed, whatever the plane taken, it necessarily intersects the helix of the trajectory T.

[0012] Several measures have been studied to quantify the impact of tomographic completeness. [Metzler, Bowsher, & Jaszczak, 2003], [Lin & Meikle, 2010] and [Liu, et al., 2012] numerically evaluate the percentage of planes that are intersected by the source trajectory by sampling the unit sphere. However, this type of measurement is only possible for a continuous source trajectory. However, in practice, acquisitions are made at discrete points along the source trajectory. Thus, planes necessarily pass between the sources. Furthermore, the authors of these publications indicate that it is difficult to predict reconstruction quality from this metric.

[0013] [Stopp, Winne, Jank, & Keeve, 2010] proposed to measure locally the average of the minimal angles, that is, the angles defined by the planes passing through the voxel under study and the associated optimal source. [Clackdoyle & Noo, 2020] also quantify tomographic incompleteness to measure how far a voxel is from respecting the Tuy condition in a given direction. These methods can be applied to discrete source trajectories.

[0014] The prior art methods do not exploit the theoretical local reconstruction conditions of tomography. The adjustment obtained is identical for any object to be reconstructed, without taking into account local sampling differences. Indeed, in the case of a regularization that would be adapted locally, the number of parameters or terms to be adjusted would be too large. For example, a non-local regularization of the Total Directional Variation type has three hyper-parameters to adjust, whereas a local total directional variation applied to a volume of one million voxels would require adjusting three million hyper-parameters. This makes the selection of hyper-parameters impossible with current methods.

[0015] It should be noted that in this document, voxel means either a pixel or a voxel depending on whether a two-dimensional or three-dimensional reconstruction is considered.

[0016] The invention aims to propose for the reconstruction of an object by X-ray tomography, a voxel-by-voxel regularization adapting to each acquisition system or scanner to which it is applied while being inexpensive in data.

[0017] The invention proposes for this purpose a method for reconstructing an object by X-ray tomography using X-ray sources and detectors, said object belonging to an imaged space defined by a set of voxels. The method comprises the following steps: - obtaining the positions of X-ray sources and detectors; - pre-calculation of a set of directions; - calculation of a tomographic incompleteness vector map comprising for each voxel a vector in a direction v* among said set of directions; - calculation of a tomographic incompleteness map of the metrics defining for each voxel a loss of information according to the direction v*; - calculation of regularization terms from the tomographic incompleteness vector map and the tomographic incompleteness map of metrics; and - iterative resolution of a tomo- reconstruction problem graph regularized using the calculated regularization terms.

[0018] The solution provided by the invention thus makes it possible to exploit the information from the tomographic incompleteness to adjust the terms of a regularization, for example by adjusting the hyper-parameters of a local Total Variation type regularization. The objective is to introduce a voxel-by-voxel specific regularization which will exploit the directional information as well as the regularization strength (i.e. the weight given to the regularization) estimated by the tomographic incompleteness.

[0019] Since the incompleteness information is local and specific to each acquisition system or scanner, this leads to a regularization that adapts to each acquisition system. The advantages obtained by the present invention are: - taking into account the geometry of the scanner; - introduction of local regularization allowing more precise reconstruction; - possibility of reconstructing a better quality image, particularly for all acquisition geometries which lead to image reconstruction problems with missing data (limited angle problems, ill-posed inverse problems, ill-conditioned inverse problems); - possibility of lifting the lock on hyper-parameter adjustment for local Total Directional Variation regularization by using an analytical criterion exploiting the geometry of the acquisition system without the slightest need to generate experimental data, simulations or reconstructions, nor human intervention or artificial intelligence; - invention applicable to known geometries as well as new geometries (multi-source stationary, etc.); - invention applicable to all X-ray tomography systems, in the fields of medical imaging, non-destructive testing or even baggage screening in airports; - invention opening new perspectives for new reconstruction methods in X-ray tomography (dynamic reconstruction with 3D Total Directional Variation + local time).

[0020] Particularly convenient preferred features of the reconstruction method according to the invention are presented below.

[0021] The calculation of the tomographic incompleteness vector map includes for each voxel of the imaged space the sub-steps of: - selecting from said sources a source for each direction from said set of directions; and - choice for each selected source of the direction v*, the loss of information defined by the tomographic incompleteness map of the metrics corresponding to the loss of information according to the direction chosen v* by said detector associated with said selected source.

[0022]

[0023]

[0024]

[0025]

[0026]

[0027]

[0028]

[0029]

[0030]

[0031]

[0032] The step of calculating the tomographic incompleteness vector map includes a preliminary sub-step of selecting, from among said sources, sources projecting the voxel of the object on their associated detectors, the source being selected from the sources projecting the object voxel onto their associated detectors. The selection of said source is carried out by calculating a directional local tomographic incompleteness v) G R+ according to the expression Z(xv)=min _ 1, where pe R3, is the projection of a source Z 1. «i Zî f on the plane defined by the point x of the voxel and the normal direction v. The chosen direction is chosen based on the local directional tomographic incompleteness / ( x, v ) • The chosen direction v* is chosen from the calculation of a norm of several directions for which the local directional tomographic incompleteness [(x, v) is calculated. The norm is the Lp norm for pe [1;+oo]. Preferably, p is equal to +oo. The chosen direction v* is chosen from the calculation of a weighted average of several directions for which the directional local tomographic incompleteness / (x, v) is calculated. Regularization terms are calculated for a regularization function It f II reg SUiVallte: He f II regu “ Jwxiie ->j - - / where 1 corresponds to the finite differences in different directions, j, k, 1 cor- corresponding to the indexing for each voxel, and d to the direction of the finite difference. The regularization function || f || . is of the direct total variation type ' ^local tional. Other features and advantages of the invention will become apparent in the description below with reference to the appended drawings, given as non-limiting examples: - [Fig.l] represents examples of source-detector trajectories of a scanner used for the reconstruction of an object; - [Fig.2] represents a diagram of a reconstruction process of a object by X-ray tomography according to one embodiment of the invention - [Fig.3] represents an example of calculation of a steering assembly; - [Fig.4] represents an example of calculation of tomographic incompleteness directional local; and - [Fig.5] shows examples of 3D maps of tomographic incompleteness.

[0033] [Fig. 2] represents a method for reconstructing an object imaged by X-ray tomography according to one embodiment. The reconstruction algorithm used in the invention for reconstructing the object is a regularized reconstruction algorithm. This method firstly involves a calculation of regularization terms, the steps of which are represented, according to an exemplary embodiment, in [Fig. 2],

[0034] X-ray tomography is performed using a scanner. The scanner includes X-ray sources and detectors. We denote {5^ , 5 e cj d7, , d ;V] g r3x2V the set of positions of source-detector pairs of the scanner geometry. There are N source-detector pairs. Each source has an associated detector.

[0035] The calculation of regularization terms and the reconstruction of the object are implemented by computer.

[0036] The reconstruction method comprises a step of obtaining E0 the positions of the X-ray sources and detectors and of defining an imaged space Q containing the object to be reconstructed. The object to be reconstructed is not taken into account in the method. In particular, whatever the shape of the object, this has no influence on the calculation of the incompleteness and therefore on the regularization. Only the dimensions of the imaged space Q containing the object to be reconstructed are taken into account in the method.

[0037] The imaged space Q is meshed into several voxels. The imaged space Q can have any type of shape and dimensions. The imaged space Q is for example a cube with dimensions 256mm x 256mm x 256mm.

[0038] The reconstruction method also comprises a step of pre-calculating a set of directions. The pre-calculation step makes it possible to discretize the directions in space. A finite number of directions in space is thus predefined. The greater the number of pre-calculated or predefined directions, the more precise the calculated incompleteness direction. However, the greater the number of pre-calculated directions, the greater the number of calculations. The number of directions is chosen so as to describe the entire space well while limiting the calculations to a reasonable number.

[0039] The pre-calculation step consists for example of calculating a set of homogeneous unit vectors on the sphere S2, as illustrated in [Fig.3]. For example, the Fibonacci network method [Gonzalez, 2009] which places the points on a spiral uniformly spherical is used. Iteratively, each new point is placed between the largest gap of the previous points.

[0040] The reconstruction method then comprises steps E2, E3 of calculating tomographic incompleteness maps. In particular, the method comprises a step E2 of calculating a vector map of tomographic incompleteness and a step E3 of calculating a tomographic incompleteness map of the metrics.

[0041] According to an exemplary embodiment, the calculation step E2 of the tomographic incompleteness vector map comprises, for each voxel of the imaged space or point x G £2, a first sub-step E20 of selection of the sources projecting the voxel of the imaged space onto their associated detectors. In other words, the sources such that the direction (5., x) intersects the detector dt are selected. The set of sources projecting the voxel of the imaged space £2 onto their associated detectors is noted [MS ..., S_n}.

[0042] The number of sources is thus reduced. For each voxel of the imaged space £2, only the sources that manage to see said voxel are selected.

[0043] The calculation step E2 of the tomographic incompleteness vector map comprises for each voxel of the imaged space or point x G £2, a second sub-step E21 of selection of the source, called the best source, best seeing the voxel.

[0044] The selection of the best source is carried out by calculating a directional local tomographic incompleteness I ( xv ) e R' for a point x G £2 and a pre-calculated direction ye S2- The directional local tomographic incompleteness is calculated according to 1 expression T / -, . (Ihr / dl. i I(x, v)=minU—7-: / =1,2. ...,n [ H ' rn • p. G R3, is the projection of a source on the plane defined by the pointx of the voxel and the normal direction v: p = gj - ( ( Sj - x ) ■ y ) V- This expression is defined by [Clackdoyle & Noo, 2020].

[0045] The above calculation of the directional local tomographic incompleteness I(x, v) evaluates the minimum tangent of the angles defined by the plane (x, v) and the X-ray lines, i.e. the lines passing through x and the source points along the trajectory. The directional local tomographic incompleteness I(x, y ) thus evaluates to what extent a voxel or pointx is reconstructable in the direction v.

[0046] The larger the value Z(x, y), the further point A is from satisfying the local Tuy condition and lacks tomographic information in the v direction. Indeed, the smaller the ratio Ik-pJI of I(x, y), the more this means that the source is “approaching” the normal plane. I (x, v) is 0 when the source is in the normal plane. Conversely, the more the v direction is aligned with the source lines seeing point x, the higher the ratio IIvpJI of / (x, y) is large, the more difficult it is to satisfy Tuy's condition local and less well we can see the pointx in the plane defined by the codirection v (or normal plane).

[0047] By calculating the minimum of the ratio II v / 'Ji, the direction which is most orthogonal to the Ikpjl source-point vector is searched. The best source is thus obtained.

[0048] The best source is here selected from the sources {.Sj, ..., projecting the voxel of the imaged space Q on their associated detectors. The first sub-step of selection E20 of the sources projecting the voxel of the imaged space Q on their associated detectors is optional. According to an exemplary embodiment, the set of sources {Sp S?, ..., is considered. The selection of the best source is carried out among the set of sources , ,£_«}.

[0049] Figure 4 shows an illustration of the calculation of the directional local tomographic incompleteness / (x, v) for the selection of the best source for the reconstruction of a point x. For the sake of simplification, only two sources si, s2 are represented. We note the angle between the line passing through the source si and the point x, and the plane orthogonal to v. The directional local tomographic incompleteness / (x, v) is equal to tanQp). The calculation of the tangent tan(^^ for the source s2 gives 0.65. The directional local tomographic incompleteness / (x. v) is therefore 0.28 here.

[0050] According to another exemplary embodiment, the source selected in step E21 may be a source different from the best source. For example, the second best source may be selected.

[0051] The calculation step E2 of the tomographic incompleteness vector map comprises for each voxel of the imaged space or point X e Q, a third sub-step of choice E22 for the best selected sources of a direction v*.

[0052] The chosen direction v* can be chosen or calculated in different ways. The chosen direction v* is chosen or calculated according to the directional local tomographic incompleteness I ( x, v ) •

[0053] The chosen direction v* can be chosen from the calculation of a norm. The chosen direction v* is for example chosen from the calculation of the norm of several directions (vl, ...., vn), said directions (called considered directions) corresponding to the set of pre-calculated directions or to a subset of said set of pre-calculated directions. The norm is for example the Lp norm of the considered directions, with pe [1; +oo].

[0054] A preferred example is to choose the worst direction v*. The worst direction v* corresponds to the direction in which each voxel is least well seen by the detector associated with the best selected source.

[0055] The worst direction y* e S2 at point XGQ is calculated by the expression: I \ . In other words, the worst direction v* corresponds to the v 1 x = argmax {I(x. v)} \ / veS2 Lp norm of the directions considered for which p tends towards +oo.

[0056] The worst direction v* is thus the one which maximizes the local directional tomographic incompleteness I( x, v) • H this is the direction in which, even by choosing the best source, we obtain the worst possible result, i.e. we see the point x the least well.

[0057] In a simplified example for the explanation in which the set of directions comprises three directions {vl, v2, v3}, the directional local tomographic incompleteness I (x, v) is calculated for each of the three directions. We thus obtain three values ​​I (x, v 1), Z (x, v2), I (x, v3) • The worst direction v* corresponds to the direction among the three directions {v 1, v2, v3] having the greatest directional local tomographic incompleteness value.

[0058] According to another exemplary embodiment, the chosen direction v* can be chosen by calculating a weighted average of the directions considered. The directional local tomographic incompleteness is calculated for each of the directions considered. A weight is associated with each of the directions considered, the weight being a function of the local incompleteness of the direction considered. The weight is all the more important as the value of the local incompleteness of the direction considered is large.

[0059] Thus, in the same simplified example according to which a set of directions comprising three directions { vl, v2, v3} is considered, the directional local tomographic incompleteness Z(x, v) is calculated for each of the three directions. We thus obtain three values ​​I ( x, v 1 ), Z ( x, v2 ), I ( X, v3 ) • Depending on the values ​​obtained Z ( x, v 1 ), I ( x, v2 ), I ( X, v3 ), weights p 1, p2, p3 are respectively associated with each of the three directions { vl, v2, v3} (weight pl is associated with Z(x, vl ), weight p2 with Z ( x, v2 ), and weight p3 with Z ( x, v3 )) • The chosen direction is calculated by calculating the weighted average of the three directions { vl, v2, v3}.

[0060] Whatever the way of calculating or choosing the direction v*, at each point " of the imaged space û, we obtain a vector giving U the chosen direction. We thus obtain the 3D vector map of the tomographic incompleteness of the directions v*: R3 -» S2, ie. the 3D spatial distribution of the tomographic incompleteness in the imaged space Q.

[0061] The calculation step E3 of the tomographic incompleteness map of the metrics is carried out after this choice of the direction v*. The tomographic incompleteness of the metrics is a function of the local directional tomographic incompleteness I ( x, v ) • In other words, the 3D map of the tomographic incompleteness of the metrics IM: R3 R+ is defined by IM ( x ) = h(I ( x, v ) ), where h is a function.

[0062] The 3D map of the tomographic incompleteness of the JM metrics: R3 -» R+ is for example defined by: IM(x) = I(x, v*(.x))- In other words, the tomographic incompleteness of the IM metrics for a given voxel corresponds to the incompleteness function applied to the chosen direction v*.

[0063] The tomographic incompleteness map of the metrics defines for each voxel a loss of information according to the direction chosen v* by the detector associated with the best source selected. The tomographic incompleteness of the metrics provides information to what extent the pointx of the imaged space Q is not clearly seen in the chosen direction U. The greater the value of IM(x), the more the pointx is incomplete in the sense of the tomographic reconstruction conditions.

[0064] The 3D map of the tomographic incompleteness of the IM metrics can be defined by another function. For example, the 3D map of the tomographic incompleteness of the IM metrics: R3 -» R+ can be defined by: IM(x) = I 2 (x, v*(X)), or by IM ( x ) = tan / (xv* ( x ) )•

[0065] We thus obtain two 3D maps of tomographic incompleteness directly linked to the intrinsic geometry of the scanner: - the tomographic incompleteness vector map indicates at each xed point £2, a direction v*(x) in which there is a lack of information (for example the worst direction v*(*) in which we have the least information on this point in terms of projection); and - the tomographic incompleteness map of the metrics quantifies this lack of information with IM( x ) •

[0066] Figure 5 shows examples of 3D maps of tomographic incompleteness where the imaged space £2 is a cube. For each point x G £2, the maps display a cone with: the chosen direction v*(*) and a value corresponding to the lack of information IM(x ) • For the sake of readability of the representation, the mesh of the imaged space £2 is here 3A3 (i.e. 27 cones). A finer mesh, for example 9A9 is preferred.

[0067] The variation of the lack of information IM(x^) is represented by means of point clouds, said points being arranged on the cones. The higher the value of IM ( x ), the greater the density of points.

[0068] In a) of [Fig.5], a case is shown where the source-detector trajectory T is circular. The figure shows that only the points on the plane of the circle of the trajectory respect the Tuy condition. The value of the tomographic incompleteness increases as one moves away from the plane of the circle. In other words, the voxel is less and less reconstructable tomographically as one moves away from the plane of the circle.

[0069] In b) of [Fig.5], a case is shown where the source-detector trajectory T is helical. The figure shows that the tomographic incompleteness is more homogeneous and very low at all points. All directions (even the worst direction) remain satisfactory for tomographic reconstruction.

[0070] The reconstruction method comprises a step E4 of calculating the regularization terms from the tomographic incompleteness vector map and the tomographic incompleteness map of the metrics. The information v*(x) and IM( x) are used to adjust the regularization terms of the regularization function.

[0071] Regularization makes it possible to give a statistical a priori on the imaged object and / or take into account the acquisition geometry of the scanner.

[0072] The regularization is carried out here by an iterative approach. In this approach, a large-dimensional minimization problem is solved, iteratively, for the tomographic reconstruction. The problem here is to find an approximation y of f such that:

[0073] f^argmin || |P + f / || f || / 2 A,a <lr

[0074] where A is the tomographic projection matrix that allows the transition from the object domain to the projection domain. In other words, matrix A is the transformation that allows the transition from image space to projection data space, b are the acquired data (also called sinogram). || Af -b || is called the data attachment term. The data attachment term ensures that the solution remains close to the measurements made. || f || is called the regulation term & twcde rization of which a is a parameter allowing to adjust the importance of the regularization with respect to the data attachment term. The regularization used in the example presented is of the local total variation type. The regularization term can be noted in this case || f || . local

[0075] The above equation imposes a regularity condition to control the divergences of the reconstruction algorithm due to inconsistencies in the input data. The regularity condition constitutes a priori information.

[0076] The data attachment term can be for example || Af -b || • This is a data attachment assuming that the acquisition noise is Gaussian. The attachment to the data can of course be different.

[0077] The regularization function || f || is notably the following:

[0078] il / • hv ,, , v ? * \ F ÿTwhere (u)corresponds to the dif- || / || finite references in different directions (e.g. the three directions (x, y, z)) and j, k, 1 correspond to the indexing for each voxel. g is a function which can for example be / * \ । * i qÇy*d ' / ) = ; 1i h ^CS finite differences are for example calculated in the manner following (subtraction of two neighboring voxels): [°079] = W (u) w = = L ~ ' f jji,IA

[0080] Of course, finite differences exploiting a different number of neighboring voxels are conceivable.

[0081] The terms IMjaj and allow weights to be introduced into the regularization. These weights are adapted in the directions of space (index d) as well as for each voxel (indices j, k, 1). gives in which direction we want to regularize more or less. 7Mjk i gives the desired regularization strength.

[0082] Let us take the example of a voxel for which there is a difference between the terms v*, Vy, v? according to each of the three directions x, y, z. If the directional local tomographic incompleteness J( x, v ) of the voxel is low, although there is a directional effect (the terms yj, y*, v* being different), the overall regularization strength is set to be low since the incompleteness / (x, v ) of the voxel is low. Conversely, for a voxel that is very poorly reconstructable (i.e. for which the directional local tomographic incompleteness I( x, v) is high), the overall regularization strength IM^j is set to be high and the term v^jkJ allows weighting in the different directions.

[0083] The regularization carried out here is therefore not only directional but also local, that is to say specific to each voxel.

[0084] The reconstruction method according to the invention further comprises a step E5 of solving the minimization problem using calculated regularization terms.

[0085] The invention proposes a criterion for automatic or semi-automatic adjustment of local regularization terms. This criterion is based on tomo incompleteness graphs calculated from the geometry of the scanner. The regularization according to the invention is particularly suitable in the case of a directional total variation type regularization. However, any reconstruction algorithm with a regularization term (or terms) that can adapt locally is possible.

[0086] The main benefit is better image reconstruction quality both visually and quantitatively.

[0087] The invention further makes it possible to both take into account the geometry of the scanner and to construct a regularization analytically. Moreover, this regularization works for all scanner geometries. Finally, the a priori on the object used in the tomographic reconstruction algorithm is adjusted locally, which makes it more accurate.

[0088] The invention makes it possible to introduce a voxel-by-voxel specific regularization which will exploit the directional information as well as the regularization strength. The regularization strength makes it possible to define the weight of the regularization with respect to the attachment to the data which is guided by the data. For a strong incompleteness, a strong weight is given to the regularization and a weaker weight to the data. Conversely, for a weak incompleteness (meaning that the object is highly reconstructable only with the acquired data), a strong weight is given to the attachment to the data while a weaker weight is given to the regularization.

[0089] Bibliography

[0090] [Clackdoyle & Noo, 2020]: Clackdoyle, R., & Noo, F. (2020, January). Quantification of Tomographic Incompleteness in Cone-Beam Reconstruction. 4, 63-80. doi:10.1109 / trpms.2019.2918222.

[0091] [Ding, 2021] : Ding, Q. N. (2021). Deep learning with adaptive hyper-parameters for low-dose CT image reconstruction. IEEE Transactions on Computational Imaging, 7, 648-660.

[0092] [Gonzalez, 2009] : Gonzalez, A. (2009, November). Measurement of Areas on a Sphere Using Fibonacci and Latitude-Longitude Lattices. Mathematical Geosciences, 42, 49-64. doi:10.1007 / sll004-009-9257-x.

[0093] [Kohavi, 1995] : Kohavi, R. (1995). A study of cross-validation and bootstrap for accuracy estimation andmodel sélection. 14, 1137-1145.

[0094] [Lin & Meikle, 2010] : Lin, J., & Meikle, S. R. (2010, October). Truncated pinhole SPECT: Sufficient sampling criteria and applications. IEEE Nuclear Science Symposuim: Medical Imaging Conférence. IEEE, doi: 10.1109 / nssmic.2010.5874140.

[0095] [Liu, et al., 2012] : Liu, B., Bennett, J., Wang, G., Man, B. D., Zeng, K., Yin, Z.,. . . Yu, H. (2012, April). Completeness map évaluation demonstrated with candidate next-generation cardiac CT architectures. Medical Physics, 39, 2405-2416. doi: 10.1118 / 1.3700172.

[0096] [Metzler, Bowsher, & Jaszczak, 2003] : Metzler, S. D., Bowsher, J. E., & Jaszczak, R. J. (2003, October). Geometrical similarities of the Orlov and Tuy sampling criteria and a numerical algorithm for assessing sampling completeness. IEEE Transactions on Nuclear Science, 50, 1550-1555. doi:10.1109 / tns.2003.817385.

[0097] [Rudin, 1992] : Rudin, L. I. (1992). Nonlinear total variation based noise removal al-gorithms. Physica D: nonlinearphenomena, 60.1(4), 259-268.

[0098] [Shen, 2018] : Shen, C. G. (2018). Intelligent parameter tuning in optimization-based itérative CT reconstruction via deep reinforcement leaming. IEEE Transactions on Medical Imaging, 37.6, 1430-1439.

[0099] [Stopp, Winne, Jank, & Keeve, 2010] : Stopp, F., Winne, C., Jank, E., & Keeve, E. (2010, February). Quality évaluation of image recording strategies for limited angle to-mography. Tsinghua Science and Technology, 15, 25-29. doi: 10.1016 / s 1007-0214(10)70004-3.

[0100] [Tuy, 1983] : Tuy, H. K. (1983, June). An Inversion Formula for Cone-Beam Reconstruction. SI AM Journal on Applied Mathematics, 43, 546-552. doi: 10.1137 / 0143035.

[0101] [Zhang, 2021] : Zhang, Z. C. (2021). Directional-TV algorithm for image reconstruction from limited-angular-range data. Medical image analysis, 70.

Claims

1.

2.

3. Claims Method of reconstructing an object by X-ray tomography using sources (sp,y?, 5^) and detectors (JdN) of X-rays, said object belonging to an imaged space (Q) defined by a set of voxels, said method comprising the following steps: - obtaining (EO) the positions of the X-ray sources and detectors; - pre-calculation (El) of a set of directions; - calculation (E2) of a tomographic incompleteness vector map comprising for each voxel a vector in a direction v* among said set of directions; - calculation (E3) of a tomographic incompleteness map of the IM(x) metrics defining for each voxel a loss of information according to the direction v*; - calculation (E4) of regularization terms from the tomographic incompleteness vector map and the tomographic incompleteness map of the metrics; and - resolution (E5), iteratively, of a tomographic reconstruction problem regularized using the calculated regularization terms. Reconstruction method according to claim 1, in which the calculation (E2) of the tomographic incompleteness vector map comprises for each voxel of the imaged space (Q) the sub-steps of: - selection (E21) from among said sources, %) of a source for each direction among said set of directions; and - choice (E22) for each selected source of direction v* the loss of information defined by the tomographic incompleteness map of the metrics IM(x) corresponding to the loss of information according to the direction chosen v* by said detector associated with said selected source. The reconstruction method of claim 2, wherein the step of calculation (E2) of the tomographic incompleteness vector map comprises a preliminary sub-step of selection (E20), among said sources (5^ ... 5 5^), of the sources projecting the voxel of the object onto their associated detectors, the source being selected among the sources (Sj, S^, ..., projecting the voxel of the object onto their associated detectors d2, ...,.

4. Reconstruction method according to claim 2 or claim 3, wherein the selection (E21) of said source is carried out by calculating a directional local tomographic incompleteness / (xv) Œ R+ according to the expression r , , . flIvAil • , ~ 1, where ng R3 Z(xv)=minU—^: / = 1,2, , is the projection of a source on the plane defined by the pointx of the voxel and the normal direction v.

5. A reconstruction method according to claim 4, wherein said chosen direction is chosen as a function of the directional local tomographic incompleteness I ( x, v ) •

6. A reconstruction method according to claim 4 or claim 5, wherein said chosen direction v* is chosen from the calculation of a norm of several directions (vl, ..vn) for which the directional local tomographic incompleteness I (x, v) is calculated.

7. A reconstruction method according to claim 6, wherein said norm is the Lp norm for pe [1;+oo], p preferably being equal to +00.

8. A reconstruction method according to claim 4 or claim 5, wherein said chosen direction v* is chosen from the calculation of a weighted average of several directions (vl, ...., vn) for which the directional local tomographic incompleteness I(x, v) is calculated.

9. Reconstruction method according to one of claims 1 to 8, in which the regularization terms are calculated for a regularization function || / 'U , following: focal Il f II = where (uj corresponds to the finite differences in different directions, j, k, 1 corresponding to the indexing for each voxel, and d to the direction of the finite difference.

10. Reconstruction method according to claim 9, in which the regularization function || / || . { is of the total directional variation type.

Images