The present invention relates generally to a method of designing an x-ray emitter panel.
Medical imaging with X-rays has been common for the past century. Planar radiology is the most common medical imaging modality. At the heart of the approach lies an x-ray source, comprising vacuum-tubes capable of generating a single cone or fan beam of x-rays over a wide range of energies and currents. However, the imaging geometries possible with these substantially point-like sources is limited; that is, the source must be placed a significant distance away from the object to be imaged in order that the cone or fan covers sufficient area. A minimum value of this distance (usually called the Source to Object Distance, SOD, or stand-off distance) can be determined trivially by trigonometric means using values of the opening angle of the x-ray of the cone and a desired coverage width of the object to be imaged.
In practice, the minimum distance to a point-like source is usually set by the so-called “skin safe distance”. Essentially, the peak dose from a fan or cone source is at the entrance to the body (typically the skin) and hence the minimum distance is restricted to that which would cause an excessive dose at the entry point.
When a single, substantially point-like source of x-rays is used, magnification is governed by the ratio of the Source to Object Distance (SOD) to the Source to Image Distance (SID).
In order to take multiple images of an object from a variety of directions, a mechanical gantry is typically used to move a single source along a sequence of locations. An alternative has been proposed in which multiple sources are placed at fixed locations around an object, each source being selectively activated. Such a system, would allow for a shorter overall period of image capture; however, this system is prohibitive due to the cost of the sources, and their relative bulk would afford only a limited number of viewing angles.
A common alternative approach to producing multiple x-ray sources is to produce multiple cathodes (or emitters) in a ‘source array’ or ‘emitter array’. In particular, Field Enhanced Emitter (FEE) arrays, sometimes referred to as Field Emitter Arrays (FEAs), such as Spindt arrays, may be used in x-ray tubes and serve as an advanced cathode. For instance, an array of moderate field enhancement tips may operate at high voltages as emitters for x-ray production. Cathodes produced from carbon nanotubes (CNTs) may offer extreme field enhancement and hence may allow for control of electron emission at low voltages. In all cases, such FEAs allow for multiple sources of x-rays from an extended or distributed source. In the case of flat-panel sources, the size of the arrays can be large and allow for significant displacement from a source on one corner of the array to a source on the opposite corner.
These approaches allow for tomosynthesis, but impose more complex geometric constraints. For instance, each emitter in the array produces its own x-ray cone (referred to herein as a ‘conelet’). In order for complete coverage of an object, there must be a certain amount of overlap of the conelets; however, the image formed by the object will include multiple images or shadows (such as double images) due to illumination of features from multiple angles.
Prior methods of addressing these issues have included the use of so-called anti-scatter grids, which also find use in the formation of conventional x-rays images. However, while these devices serve to limit the acceptance angle of x-rays to the detector, they also limit the information available for a given dose. Therefore, it is desirable to preserve all the available information, and separate the overlaps in software (and/or use a variable number of conelets to obtain more information about an object). With sufficient prior knowledge of aspects of an object, it is possible to separate spatiotemporally overlapping conelets; however, detailed knowledge of the aspects of an object are not always available. Hardware methods such as structured light can further assist in the software image reconstruction; however, in general, and in medical imaging where prior knowledge is often limited in particular, it is usually desirable to have temporally separated images, and then use other reconstruction methods to render a 3D model of the object.
A wide range of methods have been considered in the literature for image reconstruction. With both deterministic emitter and detector locations, it is possible to apply powerful methods to determine 3D information about the object with only the most rudimentary assumptions. The various methods are often referred to as compressed sensing, basis pursuit and sparse inverse problems. Regardless of the specific approach, voxelization of the space to be imaged is typically used to define the smallest image volume. Iterative methods may then be employed to determine the set of voxel values which best fit the object (and reduce a value function to a minimum). The voxelization of the image volume as well as the information that can be recovered are largely determined by the emitter geometry and the SOD and SID. According to a first aspect of the present invention, there is provided a method of designing an x-ray emitter panel for use as a distributed x-ray source, the x-ray emitter panel for use with an x-ray detector panel, the method comprising the steps of
The step of choosing may comprise selecting a desired value and/or determining an actual value possessed by components intended to be used in manufacture of the designed x-ray emitter panel, for instance by testing.
The method may further comprise selecting an array pattern for placement of x-ray emitters. The array pattern may comprise a grid of triangles, for instance equilateral triangles. Alternatively or additionally, the array pattern may be a grid of squares, rectangles and/or hexagons. The array pattern may be a substantially pseudorandom array pattern.
The array pattern may be a combination of array patterns. For instance, the array pattern may be a combination of two, three, four or five overlapping array patterns. The array patterns may have respective pitch scales that differ between array patterns.
A may be between approximately 1 and 20, in particular between approximately 1 and 16. The choice of value of A may depend on the geometry of the system, and in particular the type of array. In particular, A may be related to the emitter density on the panel. The emitter density of a panel may be proportional to the inverse of the square of the pitch scale r, and A may be proportional to the same constant of proportionality. That is, if the emitter density of a panel may be calculated as:
then:
A∝B
In particular, for a triangular array B may be approximately equal to 8√{square root over (3)}, and for a square array B may be approximately equal to 1.
More particularly, A may be defined as follows:
For a triangular array, A may be between approximately 12 and 21, in particular between approximately 13 and 18, more particularly between approximately 14 and 17. A may be between approximately 15 and 16, in particular A may be between approximately 15.5 and 16, more particularly A may be approximately 15.7. A may be approximately 15.72.
For a square array, A may be between approximately 0.5 and 3, in particular between approximately 0.8 and 2, more particularly between approximately 1 and 1.5. A may be between approximately 1.1 and 1.3, in particular A may be between approximately 1.1 and 1.2, more particularly A may be approximately 1.14. A may be approximately 1.136.
Other array patterns are envisaged with corresponding values of A and/or B as would be understood by the skilled person.
ρdet may be a predetermined average (e.g. mean) density of detectors in the detector panel. ρdet may have a value between approximately 103 m−2 and 1010 m−2, in particular between approximately 107 m−2 and 5×109 m−2, more particularly between approximately 5×107 m−2 and 2×109 m−2, for instance, approximately 5.1×107 m−2, 108 m−2, or 1.1×109 m−2.
F may be between approximately 0.05 m2 and 0.3 m2, in particular between approximately 0.1 m2 and 0.2 m2, more particularly between approximately 0.12 m2 and 0.18 m2, for instance approximately 0.12 m2, 0.16 m2, 0.17 m2 or 0.18 m2. For instance, the emitter panel may have dimensions of 0.3 m, 0.4 m, and/or 0.42 m.
Emin may be between approximately 10 and 1500, in particular between approximately 500 and 1200, more particularly between approximately 800 and 1100, for instance approximately 1000.
Etot may be between approximately 1×1012 and 1×1017 or may be between approximately 1×1013 and 1×1016, in particular between approximately 1×1014 and 1×1015, more particularly between approximately 2×1014 and 8×1014, for instance approximately 6×1014.
ηbre may be between approximately 0.001 and 0.1, in particular between approximately 0.005 and 0.05, more particularly between approximately 0.008 and 0.03, for instance approximately 0.01.
The pitch scale may be referred to as the emitter scale and may be between approximately 0.01 m and 0.1 m, in particular between approximately 0.02 m and 0.09 m, more particularly between approximately 0.03 m and 0.08 m, for instance approximately 0.0100 m, 0.0182 m, 0.0212 m, 0.0363 m or 0.0872 m.
The emitter pitch may be statistically defined as the average distance between adjacent emitters; that is, the distance between successive emitters in the array. In a regular array, the distance between successive emitters in the array is invariant; however, in a pseudorandom arrangement, the average may be taken. In some arrays, the emitter pitch and/or the pitch scale may vary across the array, particularly at the edges. The emitter pitch may be between approximately 1 and 4 times the value of r; in particular between approximately 1 and 3.5 times, more particularly between approximately 1.5 and 3 times; for instance, between approximately 2 and 2.5 times. The emitter pitch may be defined as r (for a square grid),
(for a triangular grid), 1.5r, 2r, 2.5r, 3r or 3.5r, and may be between approximately 0.0001 m or 0.001 m and 0.05 m, in particular between approximately 0.005 and 0.02 m, more particularly approximately 0.01 m.
Solving the inequality may comprise finding an approximate solution.
Solving the inequality may comprise selecting a minimum value of the pitch scale r which satisfies the inequality.
Solving the inequality may comprise solving the equation:
Solving the equation may comprise applying Newton's method.
The method may further comprise the step of selecting a collimation angle, α, that is less than or equal to (or simply less than) twice the arctangent of the ratio of the selected pitch scale r to the maximum emitter-detector panel separation Dmax; for example:
The method may further comprise the step of selecting a collimation angle, α, that is substantially equal to twice the arctangent of the ratio of the selected pitch scale r to the maximum emitter-detector panel separation Dmax; that is:
α may be between approximately 5° and 70°, in particular between approximately 10° and 40°, more particularly between approximately 15° and 30°, for instance approximately 20°, 20.6°, 23° and/or 24°.
The method may further comprise the step of selecting a collimation angle, α, that satisfies the inequality:
in which:
Mdesign may be determined by solving the further equation:
in which:
As can be appreciated by referring to
D=d+δ
The above relations can similarly be applied to maximal, actual and design variables as well, as would be understood by the skilled person. For instance, in actual use (during x-ray imaging), the variables may differ from the design variables. Thus, the method may further comprise a relation between the collimation angle, α, and the other variables in use satisfies the inequality:
in which:
Mactual may be determined by solving the further equation:
in which:
A lower limit for α may be given by
It will be appreciated that the above relations are not all independent.
Dmax may be determined based on an imaging modality. The method may further comprise determining Dmax.
Dmax may be determined by solving the further equation:
in which:
dmax is a predetermined maximum thickness of tissue placed between the emitter panel and the detector panel; and
Mmin is a predetermined minimum value of a factor representative of a multiplicity of overlapping conelets from adjacent x-ray emitters on a given part of the tissue, wherein Mmin has a value between 1 and 4.
Dmax may be between approximately 0.05 m and 2 m, in particular between approximately 0.07 m and 1 m, more particularly between approximately 0.08 m and 0.5 m, for instance 0.1 m, 0.2 m, 0.4 m, 0.48 m and/or 0.5 m. It may be between approximately 0.2 m and 0.4 m.
dmax may be between approximately 0.01 m and 1 m, in particular 0.05 m and 0.4 m, more particularly between approximately 0.1 m and 0.3 m, for instance 0.05 m or 0.3 m.
Mmin may be determined based on a consideration of the specific image reconstruction approach used and the desired speed of imaging, where here the speed may include the complete process of image acquisition and processing. The method may further comprise determining Mmin.
Mmin may be determined based on a minimum value of Mactual.
In general, a larger Mmin corresponds to more information available for image reconstruction and a longer period being required to obtain the image. Mmin may be between approximately 1 and 3.9, in particular between approximately 1 and 3.7, more particularly between approximately 1 and 3.6, for instance approximately 1.5, 2, 3 or 3.6.
The above and other characteristics, features and advantages of the present invention will become apparent from the following detailed description, taken in conjunction with the accompanying drawings, which illustrate, by way of example, the principles of the invention. This description is given for the sake of example only, without limiting the scope of the invention. The reference figures quoted below refer to the attached drawings.
The present invention will be described with respect to certain drawings but the invention is not limited thereto but only by the claims. The drawings described are only schematic and are non-limiting. Each drawing may not include all of the features of the invention and therefore should not necessarily be considered to be an embodiment of the invention. In the drawings, the size of some of the elements may be exaggerated and not drawn to scale for illustrative purposes. The dimensions and the relative dimensions do not correspond to actual reductions to practice of the invention.
Furthermore, the terms first, second, third and the like in the description and in the claims, are used for distinguishing between similar elements and not necessarily for describing a sequence, either temporally, spatially, in ranking or in any other manner. It is to be understood that the terms so used are interchangeable under appropriate circumstances and that operation is capable in other sequences than described or illustrated herein.
Moreover, the terms top, bottom, over, under and the like in the description and the claims are used for descriptive purposes and not necessarily for describing relative positions. It is to be understood that the terms so used are interchangeable under appropriate circumstances and that operation is capable in other orientations than described or illustrated herein.
It is to be noticed that the term “comprising”, used in the claims, should not be interpreted as being restricted to the means listed thereafter; it does not exclude other elements or steps. It is thus to be interpreted as specifying the presence of the stated features, integers, steps or components as referred to, but does not preclude the presence or addition of one or more other features, integers, steps or components, or groups thereof. Thus, the scope of the expression “a device comprising means A and B” should not be limited to devices consisting only of components A and B. It means that with respect to the present invention, the only relevant components of the device are A and B.
Similarly, it is to be noticed that the term “connected”, used in the description, should not be interpreted as being restricted to direct connections only. Thus, the scope of the expression “a device A connected to a device B” should not be limited to devices or systems wherein an output of device A is directly connected to an input of device B. It means that there exists a path between an output of A and an input of B which may be a path including other devices or means. “Connected” may mean that two or more elements are either in direct physical or electrical contact, or that two or more elements are not in direct contact with each other but yet still co-operate or interact with each other.
Reference throughout this specification to “an embodiment” or “an aspect” means that a particular feature, structure or characteristic described in connection with the embodiment or aspect is included in at least one embodiment or aspect of the present invention. Thus, appearances of the phrases “in one embodiment”, “in an embodiment”, or “in an aspect” in various places throughout this specification are not necessarily all referring to the same embodiment or aspect, but may refer to different embodiments or aspects. Furthermore, the particular features, structures or characteristics of any embodiment or aspect of the invention may be combined in any suitable manner, as would be apparent to one of ordinary skill in the art from this disclosure, in one or more embodiments or aspects.
Similarly, it should be appreciated that in the description various features of the invention are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the disclosure and aiding in the understanding of one or more of the various inventive aspects. This method of disclosure, however, is not to be interpreted as reflecting an intention that the claimed invention requires more features than are expressly recited in each claim. Moreover, the description of any individual drawing or aspect should not necessarily be considered to be an embodiment of the invention. Rather, as the following claims reflect, inventive aspects lie in fewer than all features of a single foregoing disclosed embodiment. Thus, the claims following the detailed description are hereby expressly incorporated into this detailed description, with each claim standing on its own as a separate embodiment of this invention.
Furthermore, while some embodiments described herein include some features included in other embodiments, combinations of features of different embodiments are meant to be within the scope of the invention, and form yet further embodiments, as will be understood by those skilled in the art. For example, in the following claims, any of the claimed embodiments can be used in any combination.
In the description provided herein, numerous specific details are set forth. However, it is understood that embodiments of the invention may be practised without these specific details. In other instances, well-known methods, structures and techniques have not been shown in detail in order not to obscure an understanding of this description.
In the discussion of the invention, unless stated to the contrary, the disclosure of alternative values for the upper or lower limit of the permitted range of a parameter, coupled with an indication that one of said values is more highly preferred than the other, is to be construed as an implied statement that each intermediate value of said parameter, lying between the more preferred and the less preferred of said alternatives, is itself preferred to said less preferred value and also to each value lying between said less preferred value and said intermediate value.
The use of the term “at least one” may mean only one in certain circumstances.
The principles of the invention will now be described by a detailed description of at least one drawing relating to exemplary features of the invention. It is clear that other arrangements can be configured according to the knowledge of persons skilled in the art without departing from the underlying concept or technical teaching of the invention, the invention being limited only by the terms of the appended claims.
and denotes the set of integers such that the defined points fit on a given panel. This pattern is shifted such that the panel is covered homogenously by 48 exposures enumerated by the formula
f+4(g−1)+16(h−1), (f=1, . . . , 4; g=1, . . . , 4; h=1, . . . , 3)
The centres of circles that are fired simultaneously in exposure (f, g, h) are given by:
The first 16 exposures (corresponding to h=1) are obtained by shifting the parent pattern to the nodes obtained by bisecting the grid of equilateral triangles twice.
The second and third group of 16 exposures are centred where the first set of exposures left holes (regions not covered by radiation). There are twice as many holes as disks in any given exposure, which leads to the three sets of 16.
Note that all emitters except those near the boundary of the panel are equidistant to their six nearest neighbours, the distance being
We call this distance emitter pitch, while we refer to r as the emitter scale. The emitter scale also has the interpretation as radius of the non-overlapping disks of radiation that reach the detector panel simultaneously in any given exposure. r may be chosen such that these disks are just touching.
Number | Date | Country | Kind |
---|---|---|---|
1418391.7 | Oct 2014 | GB | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/IB2015/057792 | 10/12/2015 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2016/059535 | 4/21/2016 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
6181764 | Solomon | Jan 2001 | B1 |
6233309 | Baptist | May 2001 | B1 |
6259765 | Baptist | Jul 2001 | B1 |
6333968 | Whitlock | Dec 2001 | B1 |
6456691 | Takahashi | Sep 2002 | B2 |
6553096 | Zhou | Apr 2003 | B1 |
6556656 | Hess | Apr 2003 | B2 |
6674837 | Taskar | Jan 2004 | B1 |
6760407 | Price | Jul 2004 | B2 |
6975703 | Wilson | Dec 2005 | B2 |
7023950 | Annis | Apr 2006 | B1 |
7082182 | Zhou | Jul 2006 | B2 |
7103138 | Pelc | Sep 2006 | B2 |
7145981 | Pelc | Dec 2006 | B2 |
7192031 | Dunham | Mar 2007 | B2 |
7295651 | Delgado | Nov 2007 | B2 |
7399987 | Viscor | Jul 2008 | B1 |
7873146 | Okunuki | Jan 2011 | B2 |
7940888 | Tsujii | May 2011 | B2 |
7945015 | Tsujii | May 2011 | B2 |
7976218 | Vermilyea | Jul 2011 | B2 |
7991114 | Okunuki | Aug 2011 | B2 |
7991120 | Okunuki | Aug 2011 | B2 |
8229074 | Mahapatra | Jul 2012 | B2 |
8306184 | Chang | Nov 2012 | B2 |
8447011 | Ohta | May 2013 | B2 |
8488742 | Tsujii | Jul 2013 | B2 |
8699657 | Baeumer | Apr 2014 | B2 |
8774351 | Funk | Jul 2014 | B2 |
9053890 | Sun | Jun 2015 | B2 |
9217719 | Lowell | Dec 2015 | B2 |
9281157 | Utsumi | Mar 2016 | B2 |
9398677 | Tang | Jul 2016 | B2 |
9478385 | Kim | Oct 2016 | B2 |
9490099 | Mackie | Nov 2016 | B2 |
9506877 | Koh | Nov 2016 | B2 |
9517043 | Tamura | Dec 2016 | B2 |
9603575 | Koh | Mar 2017 | B2 |
9653247 | Tang | May 2017 | B2 |
9653251 | Tang | May 2017 | B2 |
9719947 | Yun | Aug 2017 | B2 |
9734979 | Tang | Aug 2017 | B2 |
9748069 | Kim | Aug 2017 | B2 |
9761404 | Tang | Sep 2017 | B2 |
9782136 | Zhou | Oct 2017 | B2 |
9833205 | Sugie | Dec 2017 | B2 |
9897557 | Kim | Feb 2018 | B2 |
9922793 | Hori | Mar 2018 | B2 |
10014148 | Tang | Jul 2018 | B2 |
10165993 | Kim | Jan 2019 | B2 |
10217598 | Evans | Feb 2019 | B2 |
10242836 | Iida | Mar 2019 | B2 |
10269527 | Kenmotsu | Apr 2019 | B2 |
20090232270 | Okunuki | Sep 2009 | A1 |
20090232272 | Tsujii | Sep 2009 | A1 |
20100098218 | Vermilyea | Apr 2010 | A1 |
20100266097 | Okunuki | Oct 2010 | A1 |
20100329416 | Tsujii | Dec 2010 | A1 |
20110038465 | Mahapatra | Feb 2011 | A1 |
20110280367 | Baeumer | Nov 2011 | A1 |
20120051510 | Ohta | Mar 2012 | A1 |
20140192955 | Lowell | Jul 2014 | A1 |
Number | Date | Country |
---|---|---|
2016059535 | Apr 2016 | WO |
Number | Date | Country | |
---|---|---|---|
20170245814 A1 | Aug 2017 | US |