US20180286087A1 - Volume image reconstruction using projection decomposition - Google Patents
Volume image reconstruction using projection decomposition Download PDFInfo
- Publication number
- US20180286087A1 US20180286087A1 US15/472,613 US201715472613A US2018286087A1 US 20180286087 A1 US20180286087 A1 US 20180286087A1 US 201715472613 A US201715472613 A US 201715472613A US 2018286087 A1 US2018286087 A1 US 2018286087A1
- Authority
- US
- United States
- Prior art keywords
- volume
- image
- projection images
- sub
- acquired
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000354 decomposition reaction Methods 0.000 title description 4
- 238000000034 method Methods 0.000 claims abstract description 59
- 238000000638 solvent extraction Methods 0.000 claims description 38
- 238000009877 rendering Methods 0.000 claims description 11
- 230000003750 conditioning effect Effects 0.000 claims description 9
- 238000003384 imaging method Methods 0.000 description 24
- 238000012545 processing Methods 0.000 description 22
- 230000008569 process Effects 0.000 description 15
- 238000007408 cone-beam computed tomography Methods 0.000 description 13
- 238000002591 computed tomography Methods 0.000 description 12
- 238000004422 calculation algorithm Methods 0.000 description 11
- 238000004590 computer program Methods 0.000 description 11
- 238000005192 partition Methods 0.000 description 11
- 238000013459 approach Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 7
- 238000009499 grossing Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 230000005855 radiation Effects 0.000 description 5
- 238000002601 radiography Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000033001 locomotion Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000002059 diagnostic imaging Methods 0.000 description 3
- 239000002184 metal Substances 0.000 description 3
- 238000003325 tomography Methods 0.000 description 3
- 210000003484 anatomy Anatomy 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000005055 memory storage Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 239000012536 storage buffer Substances 0.000 description 2
- 241000220317 Rosa Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 239000000872 buffer Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013506 data mapping Methods 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003707 image sharpening Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G06T5/003—
-
- G06T5/009—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/436—Limited angle
Definitions
- the invention relates generally to the field of diagnostic imaging and in particular to radiographic volume imaging and image reconstruction techniques that are directed to reducing aliasing effects and maintaining image quality while reducing the number of projection images obtained.
- Digital radiographic volume imaging provides three-dimensional (3-D) images that have been reconstructed from a series of 2-D images taken over a succession of angles of the X-ray source relative to the detector.
- Acquisition of the 2-D projection images used for cone beam CT employs a large-area digital detector, such as a digital radiography (DR) detector that is typically used for conventional single projection radiography.
- DR digital radiography
- Computed tomography (CT) systems such as cone beam computed tomography (CBCT) or cone beam CT systems offer considerable promise as one type of diagnostic tool for providing 3-D volume images.
- Cone beam CT systems capture volume data sets using a high frame rate flat panel digital radiography (DR) detector and an X-ray source.
- the X-ray source and detector are typically affixed to a gantry that revolves about the object to be imaged, with the X-ray source directing, from various points along its orbit around the subject, a divergent cone beam of X-rays toward the subject.
- the CBCT system captures projection images throughout the source-detector orbit, for example, with one 2-D projection image at every angular increment of rotation.
- the projections are then reconstructed into a 3D volume image using various techniques.
- FBP filtered back projection
- a factor that affects the quality of volume reconstruction relates to the number of 2-D projection images acquired.
- Projection images are generally obtained at evenly spaced angular increments; having images at a sufficient number of angles helps to provide enough data to minimize or eliminate aliasing effects such as view aliasing artifacts, typically appearing in the form of regularly spaced streaks, and other image processing problems.
- Each projection image requires exposure of the patient.
- 2-D projection image data is advantageous for 3-D image reconstruction, it would be preferable to reduce the number of 2-D projection images that are obtained, thereby reducing exposure risks to the patient.
- volume imaging methods that can obtain sufficient projection image content for accurate volume reconstruction while reducing dosage requirements.
- Embodiments disclosed herein offer methods that can help to reduce patient exposure levels for 2-D projection image acquisition without compromising 3-D volume reconstruction results.
- the present disclosure seeks to improve methods for volume image reconstruction from projection images by a dynamic partitioning of the image space.
- a method for forming an image comprising: reconstructing a full volume image of a subject according to a plurality of X-ray projection images acquired from a digital detector at a plurality of acquisition angles; partitioning the full volume image to form at least a first sub-volume and a second non-overlapping sub-volume; within each partitioned sub-volume, calculating a corresponding set of forward projection images for the sub-volume, wherein members of the corresponding set of forward projection images are computed at angles corresponding to the acquisition angles, and calculating an intermediate set of intermediate forward projection images at angles between the acquisition angles; calculating a weight factor relating to the contribution of each pixel in the plurality of acquired X-ray projection images to each sub-volume at each acquisition angle; forming a plurality of synthesized sub-volume projection images according to the calculated weight factors and the plurality of acquired projection images in each sub-volume; combining the plurality of formed synthesized sub-volume projection images corresponding to each intermediate forward projection angle to form each
- FIG. 1 shows, in schematic form, the scanning activity of a conventional CBCT imaging apparatus.
- FIG. 2 is a logic flow diagram showing a sequence for enhanced 3D reconstruction.
- FIGS. 3A and 3B are schematic diagrams illustrating a detailed sequence for enhanced 3-D volume reconstruction using synthetic images.
- FIGS. 4A and 4B show projections using a ray-trace approach.
- FIG. 5 shows sinogram curves for volume image content.
- FIGS. 6A and 6B show geometric partitioning.
- FIG. 7 shows a logic flow diagram using volume decomposition.
- FIG. 8 shows an exemplary window that defines a search area.
- FIGS. 9A, 9B, and 9C show exemplary geometric partitioning schemes.
- FIG. 10 shows the generation of a forward projection signal data using partitioning according to an embodiment of the present disclosure.
- FIG. 11 illustrates the meaning of boundary matching terms.
- FIG. 12A shows an acquired projection image
- FIG. 12B shows a synthesized projection image generated according to an embodiment of the present disclosure.
- volume image is synonymous with the terms “3-dimensional image” or “3-D image”.
- Embodiments of the present invention are particularly well suited for suppressing the types of metal artifacts that occur in 3-D volume images, including cone-beam computed tomography (CBCT) as well as fan-beam CT images.
- CBCT cone-beam computed tomography
- Partition is used herein to indicate a division of a whole into non-overlapping parts. Partition of an imaged volume segments the volume into a number of non-overlapping portions.
- the mathematical definition of a cylinder includes not only the familiar right circular cylinder, but any number of other shapes whose outer surface can be traced out by moving a straight line parallel to a fixed straight line, wherein the moving straight line intersects a fixed planar closed curve or base.
- rendering is the active process of generating and forming an image for display and generating the pattern of signals needed for displaying it to a user.
- Image data content that is used for rendering can be transformed from a 2-D or 3-D model (or models), typically stored as scene content in some type of scene file, into suitable patterns of light energy that are emitted from a display screen.
- a scene file contains objects in a strictly defined language or data structure, describing aspects of the image content such as geometry, viewpoint, texture, lighting, and shading information as a description of a scene.
- the data contained in the scene content or scene file is passed to a rendering program to be processed and output or streamed to a display driver or graphics processing unit (GPU) for direct presentation on a display or to a digital image or raster graphics image file.
- the digital image data file can alternately be available for presentation on a display.
- rendering provides a transformation that can be considered as analogous to an “artist's rendering” of a scene; different artists working in different media can generate different renderings of the same scene content.
- the same image content can be rendered, for example, on a monochrome display or in color on a full color display.
- modalities is a term of art that refers to types of imaging. Modalities for an imaging system may be conventional X-ray radiography, fluoroscopy or pulsed radiography, tomosynthesis, tomography, ultrasound, MRI, or other types of imaging.
- subject refers to the patient who is being imaged and, in optical terms, can be considered equivalent to the “object” of the corresponding imaging system.
- set refers to a non-empty set, as the concept of a collection of elements or members of a set is widely understood in elementary mathematics.
- subset or “partial subset”, unless otherwise explicitly stated, are used herein to refer to a non-empty proper subset, that is, to a subset of the larger set, having one or more members.
- a subset may comprise the complete set S.
- a “proper subset” of set S is strictly contained in set S and excludes at least one member of set S.
- a “partition of a set” is a grouping of the set's elements into non-empty subsets so that every element is included in one and only one of the subsets. Two sets are “disjoint” when they have no element in common.
- pixel refers to a picture element unit cell containing a photo-conversion circuit and related circuitry for converting incident electromagnetic radiation to an electrical signal.
- pixels for picture image data elements, conventionally used with respect 2-D imaging and image display, and “voxels” for volume image data elements, often used with respect to 3-D imaging, can be used interchangeably.
- the 3-D volume image is itself generated from image data obtained as pixels on a 2-D sensor array and displays as a 2-D image from some angle of view.
- 2-D image processing and image analysis techniques can be applied to the 3-D volume image data.
- techniques described as operating upon pixels may alternately be described as operating upon the 3-D voxel data that is stored and represented in the form of 2-D pixel data for display.
- techniques that operate upon voxel data can also be described as operating upon pixels.
- volume imaging refers to volume radiographic imaging modalities such as, but not limited to, computed tomography (CT) and/or tomosynthesis imaging.
- Volume imaging methods form a volume 3-D image of a subject that can be viewed as a planar slice or plane section taken at a specified depth and angle.
- volume imaging obtains 3-D depth information by changing the relative angle between the X-ray source and the subject for each 2-D projection image that is acquired during scanning.
- the phrase “in signal communication” indicates that two or more devices and/or components are capable of communicating with each other via signals that travel over some type of signal path.
- Signal communication may be wired or wireless.
- the signals may be communication, power, data such as image data, or energy signals.
- the signal paths may include physical, electrical, magnetic, electromagnetic, optical, wired, and/or wireless connections between the first device and/or component and second device and/or component.
- the signal paths may also include additional devices and/or components between the first device and/or component and second device and/or component.
- CBCT imaging apparatus and the imaging algorithms used to obtain 3-D volume images using such systems are well known in the diagnostic imaging art and are, therefore, not described in detail in the present application.
- Some exemplary algorithms and approaches for forming 3-D volume images from the source 2-D images, projection images that are obtained in operation of the CBCT imaging apparatus can be found, for example, in the disclosure of U.S. Pat. No. 5,999,587 entitled “Method of and System for Cone-Beam Tomography Reconstruction” to Ning et al. and U.S. Pat. No.
- Embodiments of the present invention can be readily adapted to the particular geometry of the CBCT or other volume imaging apparatus.
- an extremity imaging apparatus can generate volume images suitable for application of methods described herein.
- Embodiments of the present disclosure provide approaches for reducing patient exposure and reducing artifacts by acquiring only a partial percentage of the N+MX-ray projection images that would otherwise be needed to for artifact-free reconstruction. Then, to obtain an ordered set of projection images that is sufficient to be used for tomographic reconstruction, a number of synthesized projection images M is generated and added to the acquired set of N projection images. This forms a larger set of N+M projection images, effectively providing a reduced angular spacing ⁇ ′ between adjacent projection images in the sequence used for reconstruction. The result is an enhanced 3-D reconstruction that can be substantially artifact-free.
- a computer or other type of dedicated logic processor for obtaining, processing, and storing image data is part of the CT or CBCT system, along with one or more displays for viewing image results.
- a computer-accessible memory is also provided, which may be a memory storage device used for longer term storage, such as a device using magnetic, optical, or other data storage media.
- the computer-accessible memory can comprise an electronic memory such as a random access memory (RAM) that is used for shorter term storage, such as employed to store a computer program having instructions for controlling one or more computers to practice the method according to the present invention.
- RAM random access memory
- FIG. 1 there is shown, in schematic form and using enlarged distances for clarity of description, the activity of a conventional CBCT imaging apparatus for obtaining the individual 2-D images that are used to form a 3-D volume image.
- a cone-beam radiation source 22 directs a cone of radiation toward a subject 20 , such as a patient or other subject.
- a sequence of images is obtained in rapid succession at varying angles about the subject, such as one image at each 1-degree angle increment in a 200-degree orbit.
- a DR detector 24 is moved to different imaging positions about subject 20 in concert with corresponding movement of radiation source 22 .
- FIG. 1 there is shown, in schematic form and using enlarged distances for clarity of description, the activity of a conventional CBCT imaging apparatus for obtaining the individual 2-D images that are used to form a 3-D volume image.
- a cone-beam radiation source 22 directs a cone of radiation toward a subject 20 , such as a patient or other subject.
- a sequence of images is obtained in rapid succession at varying angles about the
- FIG. 1 shows a representative sampling of DR detector 24 positions to illustrate how these images are obtained relative to the position of subject 20 .
- a suitable imaging algorithm such as filtered back projection (FBP) or other conventional technique, is used for generating the 3-D volume image.
- Image acquisition and program execution are performed by a computer 30 or by a networked group of computers 30 that are in image data communication with DR detectors 24 .
- Image processing and storage is performed using a computer-accessible memory 32 .
- the 3-D volume image can be presented on a display 34 .
- FBP is a discrete implementation of an analytic model that assumes that CT transmission measurements are linear functions of the attenuation line integrals along the corresponding primary photon trajectories through the subject and are noiseless.
- CT transmission measurements are linear functions of the attenuation line integrals along the corresponding primary photon trajectories through the subject and are noiseless.
- rendering is the active process of generating and forming an image for display and generating the pattern of signals needed for displaying it to a user.
- Image data content that is used for rendering can be transformed from a 2-D or 3-D image model (or models), typically stored as scene content in some type of scene file, into suitable patterns of light energy that are emitted from a display screen.
- Embodiments of the present disclosure provide approaches for reducing patient exposure and reducing artifacts by acquiring a number N of actual projection images that represents only a partial percentage of the full number (N+M) of X-ray projection images that would otherwise be needed to for artifact-free reconstruction. Then, to obtain an ordered set of projection images that is sufficient to be used for tomographic reconstruction, a number of synthesized projection images M are generated and added to the acquired set of N projection images. This forms a larger set N+M of projection images, effectively providing a reduced angular spacing ⁇ ′ between adjacent projection images in the sequence used for reconstruction. The result is an enhanced 3-D reconstruction that can be substantially artifact-free.
- the number M of synthesized projection images that is generated can be smaller than, equal to, or even larger than the number N of actual X-ray projection images. Examples given following show generation of a single synthesized image M corresponding to an angle midway between the corresponding angles of first and second adjacent X-ray projection images n and (n+1); the image obtained at this midway angle is represented herein using the notation n.5. Thus, in the examples shown here, M ⁇ N. However, multiple synthesized images M can be generated for angles between any two adjacent X-ray projection images n and (n+1). For example, using the same processing approach, two synthesized images M can be generated between any two adjacent X-ray projection images n and (n+1).
- An embodiment of the present disclosure provides a method for improving reconstruction results without additional exposure to the patient by generating synthesized images that are formed using partitioning, forward projection, and an image data mapping sequence.
- the synthetic images can then be used to supplement the image content provided from actual acquired projection images without additional exposure to the patient.
- FIG. 2 is a logic flow diagram showing an overview of the sequence for forming and using synthetic images for enhanced volume reconstruction.
- an ordered set having a number N of actual X-ray projection images ⁇ A0, A1, A2, A3, . . . ⁇ , termed the “acquired” images in the context of the present disclosure is obtained.
- Each X-ray projection image in the set is acquired at one of a corresponding sequence of N capture angles ⁇ as indicated in FIG. 1 .
- the ordered set of N acquired X-ray projection images can be considered ordered according to acquisition angle.
- X-ray projection images can be acquired in any order with respect to angle; the ordered set structure is convenient arrangement for showing the acquisition sequence and for illustrating the processing that follows, used to generate synthesized images.
- every acquired X-ray projection image An has a corresponding capture angle ⁇ n and is adjacent to a previous acquired X-ray projection image (n ⁇ 1) with corresponding capture angle ⁇ n ⁇ 1 and to a subsequent or following acquired X-ray projection image (n+1) with corresponding capture angle ⁇ n+1 .
- Angular spacing between any two adjacent acquired X-ray projection images An in the sequence can be expressed as ⁇ and can be at equal intervals.
- an initial reconstruction step S 210 forms a 3-D volume image using an analytic reconstruction method such as filtered back projection (FBP), or using an iterative reconstruction process, or using some other suitable reconstruction algorithm.
- FBP filtered back projection
- the procedural sequence of FIG. 2 performs a forward projection step S 220 that generates an intermediate set of N+M forward projection images, shown in FIG. 2 and described herein as ⁇ F0, F0.5, F1, F1.5, F2, F2.5, F3 . . . ⁇ at corresponding angles ⁇ 0 ⁇ 0.5 ⁇ , 1 ⁇ , 1.5 ⁇ , 2 ⁇ , 2.5 ⁇ , 3 ⁇ , . . . ⁇ .
- the notation using “0.5” indicates angles between the acquisition angles ⁇ , such as half-way between.
- Forward projection thus generates projection images at each of the angles ⁇ corresponding to the N actual X-ray projections and adds M additional projection images at intermediate angles, such as at angles n.5 ⁇ , half-way between the original angles ⁇ . It should be emphasized that the intermediate angles can be at other than half-way between acquisition angle, such as at 1 ⁇ 3 intervals or at some arbitrary angular offset between the nearest acquired images.
- the process then proceeds to re-build a set of M synthesized projection images ⁇ S0.5, S1.5, S2.5, S3.5, . . . ⁇ using the M added projection images as templates in a synthesized image generation step S 230 .
- a sequence for generating the synthesized images using decomposition or partitioning of the volume is described in more detail subsequently.
- a final reconstruction step S 250 executes, generating an enhanced 3-D volume reconstruction using the original acquired X-ray projection images ⁇ A0, A1, A2, A3, . . . ⁇ supplemented by the set of synthesized images ⁇ S0.5, S1.5, S2.5, S3.5, . . . ⁇ that have been formed using the process described subsequently.
- the final reconstruction step S 240 thus uses the combined set ⁇ A0, S0.5, A1, S1.5, A2, S2.5, A3, S3.5, . . . ⁇ for forming the enhanced 3-D reconstruction.
- the data can then be displayed, rendered for viewing in a display step S 250 .
- rendering relates to the selection, processing, and display of a particular 2-D view of the reconstructed 3-D image scene content.
- FIG. 3A is a schematic diagram that shows supplementing the ordered set of N acquired X-ray projection images An with a set of M synthesized images Sn.
- FIG. 3B is a schematic diagram that show's a combined set of acquired X-ray and synthesized images used for volume image reconstruction.
- volume reconstruction While the process described with reference to FIGS. 2-3B provides suitable image content for volume reconstruction with reduced exposure to the patient, there is room for improvement.
- One inherent difficulty that impacts volume reconstruction relates to relative movement of volume content in different directions.
- FIG. 4A schematically shows, as a top view, a simplified model used to describe the Radon transform and the resulting sinogram that is obtained by combining successive images, each representative of the amount of energy that passes through object O onto detector 24 .
- FIG. 4A shows the transform portion obtained for a projection image acquired at an angle of ⁇ 1 .
- FIG. 4B shows the transform portion that is obtained at angle ⁇ 2 .
- rotation of the radiation source 22 and detector 24 from angle ⁇ 1 to ⁇ 2 changes the relative positions of points A and B in the sequence of projection images that are acquired.
- the source 22 /detector 24 pair rotates from angle ⁇ 1 to ⁇ 2 , the radiation energy for point A appears to move toward the left.
- FIG. 5 shows the changing spatial relationships between two imaged points, such as those described with reference to FIGS. 4A and 4B , in a sinogram view.
- the sinogram can be considered as a stacked collection of projections of a single slice of the subject, viewed from the side.
- the rotation axis is designated as R.
- the sinogram SN 1 for a first object point shows that the first object point is at the front of each projection, and moving in one direction, for angles over the range indicated J 1 and at the back of each projection, moving in the opposite direction, for angles in the range J 2 .
- Sinogram SN 2 for a second object point shows that the second object point is at the rear of each projection, and moving in one direction for angles over the range indicated K 1 and at the back of each projection over the range of angles K 2 . Over the range of angles K 3 , the second object point is once again at the front of each projection.
- At least one embodiment of the present disclosure addresses the difficulties that are inherent to the scan geometry, as described with reference to FIGS. 4A-5 , by using a geometric partitioning of the volume space.
- geometric partitioning the positions of the X-ray source and detector relative to the subject volume are considered and the volume is partitioned to help compensate for the relative movement conflicts noted with reference to FIGS. 4A, 4B, and 5 .
- one partitioning approach divides the volume in half for each relative position of the source and detector. That portion of the imaged volume V that is forward of the central axis of rotation with reference to the source, and is defined by a bisecting plane P that is orthogonal to a central ray of illumination, forms one partial sub-volume. The portion of the imaged volume that lies behind the bisecting plane P is a second partial sub-volume for the projection images at that angle. In the geometric partitioning logic, the partial sub-volume varies at each projection angle.
- Partitioning can also address difficulties with various types of image content.
- Another type of partitioning that can be applied to the volume and used herein relates to voxel content or density data value for the imaged object.
- the partitioning is static, so that the same voxels occupy each sub-volume for every source-to-detector acquisition angle.
- metal content can form a sub-volume for separate forward projection processing.
- Other known components, such as bone or soft tissue can similarly be used to define sub-volumes.
- a static partitioning defines a fixed portion of the imaged volume, which can include any imageable portion of the object space, the cylindrical volume that lies between the source and detector as they orbit the subject.
- a weight can be computed for each calculated signal needed to form the forward projection images Fn.
- a front and rear weight are generated for each pixel:
- the front and rear partitions can then be individually identified and processed using the computed weights.
- partitioning can generate two or more sub-volumes for the imaged content and that similar processing applies for each of two or more sub-volumes.
- the calculated weight values are then used to determine how each pixel value is assigned relative to the partitioning.
- An embodiment of the present disclosure provides a method for improving reconstruction results without additional exposure to the patient using the synthetic image generation process described with reference to FIGS. 2-3B along with volume decomposition or partitioning as described above and with reference to FIGS. 4A-6B .
- the logic flow diagram of FIG. 7 shows a sequence for volume image reconstruction using partitioning.
- an ordered set having a number N actual X-ray projection images ⁇ A0, A1, A2, A3, . . . ⁇ termed the “acquired” images in the context of the present disclosure.
- Each X-ray projection image in the set is acquired at one of a corresponding sequence of N capture angles ⁇ as indicated in FIG. 1 .
- the ordered set of N acquired X-ray projection images can be considered ordered according to acquisition angle.
- X-ray projection images can be acquired in any order with respect to angle; the ordered set structure is convenient arrangement for showing the acquisition sequence and for illustrating the processing that follows, used to generate synthesized images.
- every acquired X-ray projection image An has a corresponding capture angle ⁇ n and is adjacent to a previous acquired X-ray projection image A(n ⁇ 1) with corresponding capture angle ⁇ n ⁇ 1 and to a subsequent or following acquired X-ray projection image A(n+1) with corresponding capture angle ⁇ n+1 .
- Angular spacing between any two adjacent acquired X-ray projection images An in the sequence can be expressed as ⁇ .
- an initial reconstruction step S 710 forms a full 3-D volume image using filtered back projection (FBP), using an iterative reconstruction process, or using some other suitable reconstruction algorithm.
- FBP filtered back projection
- a partitioning step S 720 then applies partitioning logic to the full 3-D volume define a set of two or more sub-volumes ⁇ SV1, SV2, . . . ⁇ .
- the sequence of FIG. 7 shows processing for the complete partition of the volume, with a number p of sub-volumes.
- sub-volumes SVn can be defined using geometry of the scanning operation or using structural or content data from the reconstructed object itself.
- a forward projection step S 730 within each sub-volume SVp, a forward projection step S 730 generates forward projection images at angles ⁇ n corresponding to the acquisition angles used for the actual X-ray projection images.
- step S 730 also generates a set of intermediate forward projection images, at angles between each angle ⁇ n and its respective, adjacent following angle ⁇ n+1 .
- Intermediate forward projection images are represented at angles (n+0.5) midway between preceding and following forward projections; as noted earlier, the intermediate forward projections can be calculated at other angular positions between angle ⁇ n and its respective, adjacent following angle ⁇ n+1 .
- a weighting step S 740 then calculates and applies a weight or weighting factor to each pixel value for the corresponding acquired image data that contributes to the sub-volume.
- the weighting factor is an estimate of the proportionate contribution of each acquired pixel of the sub-volume to the computed forward projection value and is based on the proportion of pixel values in the sub-volume to the total data value obtained for that pixel at the corresponding acquisition angle.
- a synthesized sub-volume projection calculation step S 750 then uses the weighted values from step S 740 to generate a synthetic projection image for each sub-volume accordingly.
- a combination step S 760 combines the synthetic projection image results for the corresponding angle from each sub-volume to form a synthesized projection image for subsequent volume reconstruction, forming a set of M synthesized projection images, using the variable notation described previously.
- a reconstruction step S 770 performs volume reconstruction using the initial set of N acquired projection images (corresponding to the N actual acquisition angles ⁇ 0 ⁇ ; 1 ⁇ , 2 ⁇ , 3 ⁇ , . . . ⁇ plus the M synthesized projection images (corresponding to acquisition angles between the set of N actual acquisition angles).
- a display step S 780 then allows display, storage, and transmission of the volume image data, as well as display, storage, and transmission of one or more of the synthetic projection images formed using the FIG. 7 process.
- FIGS. 9A-9C show partitioning using front and rear segments of the volume image space; however, other partitioning models can be applied.
- FIGS. 9A-9C Some of the alternate partition models are shown in FIGS. 9A-9C .
- FIG. 9A shows additional segments defined within the reconstructed volume.
- FIG. 9B shows a top view of a partition based on relative rate of change for volume content according to distance from rotation axis R.
- FIG. 9C shows a wedged partition of the volume.
- the same processing described for partitioning into front and rear portions, as described with reference to FIGS. 6A and 6B can be used for the geometric partitioning shown in FIGS. 9A-9C or for other geometric partitioning arrangements that define a partition having non-overlapping sub-volumes.
- image conditioning processes can be applied to the volume data that is used for forward projection as well as for any of the generated forward projection or synthesized image data.
- Image conditioning improves the appearance of a projection image or of a synthesized image and can be applied as part of display rendering in order to change image sharpening or blurring, such as to improve image patch identification and processing, for example.
- Image conditioning such as metal artifacts reduction can be applied to the original acquired X-ray projection images or to the volume data that is reconstructed from the original X-ray projection images.
- FIG. 10 shows mapping and combining patch content for a candidate patch from two adjacent acquired X-ray projection images to form patch p on an intermediate image designated Fn.5 using forward projection, in order to form synthesized image Sn.5 at completion of this processing.
- a synthesized image constructed using the method described herein can be formed as a mosaic of patches. Often, boundaries or seams between patches are readily discernable, causing perceptible degradation in the final reconstruction volume.
- One method to correct for seams is to perform a one dimensional smoothing filter across the seams. If the filter is applied uniformly, however, some sharp edges are overly smoothed and clash with the actual images.
- seam smoothing is executed with respect to the forward projected image.
- the synthesized image seams are smoothed by applying the following one directional algorithm vertically and horizontally:
- the approach outlined above can help to dampen the seam boundary adaptively and limit image degradation that can otherwise occur with a global smoothing approach.
- FIG. 11 shows the meaning of boundary matching terms D A , D R , D B , D L , used herein.
- the basic energy relationship E can be expressed as follows:
- Block matching energy can then be formulated as follows:
- E blockmatching ⁇ ( x , y ) ⁇ S ⁇ ( B ) ⁇ ⁇ f left ⁇ ( x - u , y - v ) - f right ⁇ ( x + u , y + v ) ⁇ Equation ⁇ ⁇ ( 2 )
- Embodiments of the present disclosure provide an improvement over earlier interpolation methods, particularly with respect to reduction of artifacts in general, particularly view aliasing artifacts.
- FIG. 12A shows an acquired projection image An of a hand.
- FIG. 12B shows a synthesized image generated according to an embodiment of the present disclosure.
- the synthesized image can be displayed, stored, transmitted, or can be used with acquired projection images for volume reconstruction.
- the present disclosure utilizes a computer program with stored instructions that perform on image data accessed from an electronic memory.
- a computer program of an embodiment of the present invention can be utilized by a suitable, general-purpose computer system, such as a personal computer or workstation.
- a suitable, general-purpose computer system such as a personal computer or workstation.
- many other types of computer systems can be used to execute the computer program of the present invention, including networked processors.
- the computer program for performing the method of the present invention may be stored in a computer readable storage medium.
- This medium may comprise, for example; magnetic storage media such as a magnetic disk such as a hard drive or removable device or magnetic tape; optical storage media such as an optical disc, optical tape, or machine readable bar code; solid state electronic storage devices such as random access memory (RAM), or read only memory (ROM); or any other physical device or medium employed to store a computer program.
- the computer program for performing the method of the present invention may also be stored on computer readable storage medium that is connected to the image processor by way of the internet or other communication medium. Those skilled in the art will readily recognize that the equivalent of such a computer program product may also be constructed in hardware.
- memory can refer to any type of temporary or more enduring data storage workspace used for storing and operating upon image data and accessible to a computer system, including a database.
- the memory could be non-volatile, using, for example, a long-term storage medium such as magnetic or optical storage. Alternately, the memory could be of a more volatile nature, using an electronic circuit, such as random-access memory (RAM) that is used as a temporary buffer or workspace by a microprocessor or other control logic processor device. Displaying an image requires memory storage. Display data, for example, is typically stored in a temporary storage buffer that is directly associated with a display device and is periodically refreshed as needed in order to provide displayed data.
- This temporary storage buffer can also be considered to be a memory, as the term is used in the present disclosure.
- Memory is also used as the data workspace for executing and storing intermediate and final results of calculations and other processing.
- Computer-accessible memory can be volatile, non-volatile, or a hybrid combination of volatile and non-volatile types.
- the computer program product of the present invention may make use of various image manipulation algorithms and processes that are well known. It will be further understood that the computer program product embodiment of the present invention may embody algorithms and processes not specifically shown or described herein that are useful for implementation. Such algorithms and processes may include conventional utilities that are within the ordinary skill of the image processing arts. Additional aspects of such algorithms and systems, and hardware and/or software for producing and otherwise processing the images or co-operating with the computer program product of the present invention, are not specifically shown or described herein and may be selected from such algorithms, systems, hardware, components and elements known in the art.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
- The invention relates generally to the field of diagnostic imaging and in particular to radiographic volume imaging and image reconstruction techniques that are directed to reducing aliasing effects and maintaining image quality while reducing the number of projection images obtained.
- Digital radiographic volume imaging provides three-dimensional (3-D) images that have been reconstructed from a series of 2-D images taken over a succession of angles of the X-ray source relative to the detector. Acquisition of the 2-D projection images used for cone beam CT employs a large-area digital detector, such as a digital radiography (DR) detector that is typically used for conventional single projection radiography.
- Computed tomography (CT) systems, such as cone beam computed tomography (CBCT) or cone beam CT systems offer considerable promise as one type of diagnostic tool for providing 3-D volume images. Cone beam CT systems capture volume data sets using a high frame rate flat panel digital radiography (DR) detector and an X-ray source. The X-ray source and detector are typically affixed to a gantry that revolves about the object to be imaged, with the X-ray source directing, from various points along its orbit around the subject, a divergent cone beam of X-rays toward the subject. The CBCT system captures projection images throughout the source-detector orbit, for example, with one 2-D projection image at every angular increment of rotation. The projections are then reconstructed into a 3D volume image using various techniques. Among the most common methods for reconstructing the 3-D volume image are filtered back projection (FBP) approaches.
- A factor that affects the quality of volume reconstruction relates to the number of 2-D projection images acquired. Projection images are generally obtained at evenly spaced angular increments; having images at a sufficient number of angles helps to provide enough data to minimize or eliminate aliasing effects such as view aliasing artifacts, typically appearing in the form of regularly spaced streaks, and other image processing problems.
- Each projection image, however, requires exposure of the patient. Thus, although having more 2-D projection image data is advantageous for 3-D image reconstruction, it would be preferable to reduce the number of 2-D projection images that are obtained, thereby reducing exposure risks to the patient. In addition, it would be beneficial to reduce the required scan time in order to help reduce image reconstruction problems due to patient motion.
- Thus, there would be advantages to volume imaging methods that can obtain sufficient projection image content for accurate volume reconstruction while reducing dosage requirements.
- It is an object of the present disclosure to advance the art of diagnostic 3-D volume imaging. Embodiments disclosed herein offer methods that can help to reduce patient exposure levels for 2-D projection image acquisition without compromising 3-D volume reconstruction results. In particular, the present disclosure seeks to improve methods for volume image reconstruction from projection images by a dynamic partitioning of the image space.
- These objects are given only by way of illustrative example, and such objects may be exemplary of one or more embodiments of the invention. Other desirable objectives and advantages inherently achieved by the disclosed invention may occur or become apparent to those skilled in the art. The invention is defined by the appended claims.
- According to an embodiment of the present disclosure, there is provided a method for forming an image comprising: reconstructing a full volume image of a subject according to a plurality of X-ray projection images acquired from a digital detector at a plurality of acquisition angles; partitioning the full volume image to form at least a first sub-volume and a second non-overlapping sub-volume; within each partitioned sub-volume, calculating a corresponding set of forward projection images for the sub-volume, wherein members of the corresponding set of forward projection images are computed at angles corresponding to the acquisition angles, and calculating an intermediate set of intermediate forward projection images at angles between the acquisition angles; calculating a weight factor relating to the contribution of each pixel in the plurality of acquired X-ray projection images to each sub-volume at each acquisition angle; forming a plurality of synthesized sub-volume projection images according to the calculated weight factors and the plurality of acquired projection images in each sub-volume; combining the plurality of formed synthesized sub-volume projection images corresponding to each intermediate forward projection angle to form each of a plurality of synthesized projection images for the full volume image; reconstructing a second volume image according to the plurality of acquired X-ray projection images and the plurality of formed synthesized projection images; and displaying, storing, or transmitting a rendering of the reconstructed second volume image.
- The foregoing and other objects, features, and advantages of the invention will be apparent from the following more particular description of the embodiments of the invention, as illustrated in the accompanying drawings. The elements of the drawings are not necessarily to scale relative to each other.
-
FIG. 1 shows, in schematic form, the scanning activity of a conventional CBCT imaging apparatus. -
FIG. 2 is a logic flow diagram showing a sequence for enhanced 3D reconstruction. -
FIGS. 3A and 3B are schematic diagrams illustrating a detailed sequence for enhanced 3-D volume reconstruction using synthetic images. -
FIGS. 4A and 4B show projections using a ray-trace approach. -
FIG. 5 shows sinogram curves for volume image content. -
FIGS. 6A and 6B show geometric partitioning. -
FIG. 7 shows a logic flow diagram using volume decomposition. -
FIG. 8 shows an exemplary window that defines a search area. -
FIGS. 9A, 9B, and 9C show exemplary geometric partitioning schemes. -
FIG. 10 shows the generation of a forward projection signal data using partitioning according to an embodiment of the present disclosure. -
FIG. 11 illustrates the meaning of boundary matching terms. -
FIG. 12A shows an acquired projection image. -
FIG. 12B shows a synthesized projection image generated according to an embodiment of the present disclosure. - The following is a detailed description of preferred embodiments of the invention, reference being made to the drawings in which the same reference numerals identify the same elements of structure in each of the several figures.
- In the drawings and text that follow, like components are designated with like reference numerals, and similar descriptions concerning components and arrangement or interaction of components already described are omitted. Where they are used, the terms “first”, “second”, and so on, do not necessarily denote any ordinal or priority relation, but are simply used to more clearly distinguish one element from another.
- In the context of the present disclosure, the term “volume image” is synonymous with the terms “3-dimensional image” or “3-D image”. Embodiments of the present invention are particularly well suited for suppressing the types of metal artifacts that occur in 3-D volume images, including cone-beam computed tomography (CBCT) as well as fan-beam CT images.
- The term “partition” is used herein to indicate a division of a whole into non-overlapping parts. Partition of an imaged volume segments the volume into a number of non-overlapping portions.
- The mathematical definition of a cylinder includes not only the familiar right circular cylinder, but any number of other shapes whose outer surface can be traced out by moving a straight line parallel to a fixed straight line, wherein the moving straight line intersects a fixed planar closed curve or base.
- In the image processing context of the present disclosure, “rendering” is the active process of generating and forming an image for display and generating the pattern of signals needed for displaying it to a user. Image data content that is used for rendering can be transformed from a 2-D or 3-D model (or models), typically stored as scene content in some type of scene file, into suitable patterns of light energy that are emitted from a display screen. A scene file contains objects in a strictly defined language or data structure, describing aspects of the image content such as geometry, viewpoint, texture, lighting, and shading information as a description of a scene. The data contained in the scene content or scene file is passed to a rendering program to be processed and output or streamed to a display driver or graphics processing unit (GPU) for direct presentation on a display or to a digital image or raster graphics image file. The digital image data file can alternately be available for presentation on a display. In general, the term “rendering” provides a transformation that can be considered as analogous to an “artist's rendering” of a scene; different artists working in different media can generate different renderings of the same scene content. The same image content can be rendered, for example, on a monochrome display or in color on a full color display.
- The term “modality” is a term of art that refers to types of imaging. Modalities for an imaging system may be conventional X-ray radiography, fluoroscopy or pulsed radiography, tomosynthesis, tomography, ultrasound, MRI, or other types of imaging. The term “subject” refers to the patient who is being imaged and, in optical terms, can be considered equivalent to the “object” of the corresponding imaging system.
- The term “set”, as used herein, refers to a non-empty set, as the concept of a collection of elements or members of a set is widely understood in elementary mathematics. The terms “subset” or “partial subset”, unless otherwise explicitly stated, are used herein to refer to a non-empty proper subset, that is, to a subset of the larger set, having one or more members. For a set S, a subset may comprise the complete set S. A “proper subset” of set S, however, is strictly contained in set S and excludes at least one member of set S. A “partition of a set” is a grouping of the set's elements into non-empty subsets so that every element is included in one and only one of the subsets. Two sets are “disjoint” when they have no element in common.
- With respect to an image detector, the term “pixel” refers to a picture element unit cell containing a photo-conversion circuit and related circuitry for converting incident electromagnetic radiation to an electrical signal. For the image processing steps described herein, the terms “pixels” for picture image data elements, conventionally used with respect 2-D imaging and image display, and “voxels” for volume image data elements, often used with respect to 3-D imaging, can be used interchangeably.
- It should be noted that the 3-D volume image is itself generated from image data obtained as pixels on a 2-D sensor array and displays as a 2-D image from some angle of view. Thus, 2-D image processing and image analysis techniques can be applied to the 3-D volume image data. In the description that follows, techniques described as operating upon pixels may alternately be described as operating upon the 3-D voxel data that is stored and represented in the form of 2-D pixel data for display. In the same way, techniques that operate upon voxel data can also be described as operating upon pixels.
- In the context of the present disclosure, “volume imaging” refers to volume radiographic imaging modalities such as, but not limited to, computed tomography (CT) and/or tomosynthesis imaging. Volume imaging methods form a volume 3-D image of a subject that can be viewed as a planar slice or plane section taken at a specified depth and angle. As noted previously, volume imaging obtains 3-D depth information by changing the relative angle between the X-ray source and the subject for each 2-D projection image that is acquired during scanning.
- In the context of the present disclosure, the phrase “in signal communication” indicates that two or more devices and/or components are capable of communicating with each other via signals that travel over some type of signal path. Signal communication may be wired or wireless. The signals may be communication, power, data such as image data, or energy signals. The signal paths may include physical, electrical, magnetic, electromagnetic, optical, wired, and/or wireless connections between the first device and/or component and second device and/or component. The signal paths may also include additional devices and/or components between the first device and/or component and second device and/or component.
- Reference is hereby made to an article by Bertram, Wiegert, Schafer, Rose, and Aach entitled “Directional View Interpolation for Compensation of Sparse Angular Sampling in Cone Beam CT” in IEEE Transactions on Medical Imaging Vol. 28, No. 7, July 2009, pp. 1011-1022, incorporated herein in its entirety.
- Reference is hereby made to U.S. Pat. No. 7,602,879 (Chen et al.) and U.S. Pat. No. 9,224,216 (Zamyatin et al), both of which are incorporated herein in their entirety.
- CBCT imaging apparatus and the imaging algorithms used to obtain 3-D volume images using such systems are well known in the diagnostic imaging art and are, therefore, not described in detail in the present application. Some exemplary algorithms and approaches for forming 3-D volume images from the source 2-D images, projection images that are obtained in operation of the CBCT imaging apparatus can be found, for example, in the disclosure of U.S. Pat. No. 5,999,587 entitled “Method of and System for Cone-Beam Tomography Reconstruction” to Ning et al. and U.S. Pat. No. 5,270,926 entitled “Method and Apparatus for Reconstructing a Three-Dimensional Computerized Tomography (CT) Image of an Object from Incomplete Cone Beam Data” to Tam, both of which are incorporated herein in their entirety.
- Embodiments of the present invention can be readily adapted to the particular geometry of the CBCT or other volume imaging apparatus. In particular, an extremity imaging apparatus can generate volume images suitable for application of methods described herein.
- Embodiments of the present disclosure provide approaches for reducing patient exposure and reducing artifacts by acquiring only a partial percentage of the N+MX-ray projection images that would otherwise be needed to for artifact-free reconstruction. Then, to obtain an ordered set of projection images that is sufficient to be used for tomographic reconstruction, a number of synthesized projection images M is generated and added to the acquired set of N projection images. This forms a larger set of N+M projection images, effectively providing a reduced angular spacing Δα′ between adjacent projection images in the sequence used for reconstruction. The result is an enhanced 3-D reconstruction that can be substantially artifact-free.
- In typical applications, a computer or other type of dedicated logic processor for obtaining, processing, and storing image data is part of the CT or CBCT system, along with one or more displays for viewing image results. A computer-accessible memory is also provided, which may be a memory storage device used for longer term storage, such as a device using magnetic, optical, or other data storage media. In addition, the computer-accessible memory can comprise an electronic memory such as a random access memory (RAM) that is used for shorter term storage, such as employed to store a computer program having instructions for controlling one or more computers to practice the method according to the present invention.
- In order to more fully understand the methods of the present invention and the problems addressed, it is instructive to review principles and terminology used for CBCT image capture and reconstruction. Referring to the perspective view of
FIG. 1 , there is shown, in schematic form and using enlarged distances for clarity of description, the activity of a conventional CBCT imaging apparatus for obtaining the individual 2-D images that are used to form a 3-D volume image. A cone-beam radiation source 22 directs a cone of radiation toward a subject 20, such as a patient or other subject. A sequence of images is obtained in rapid succession at varying angles about the subject, such as one image at each 1-degree angle increment in a 200-degree orbit. ADR detector 24 is moved to different imaging positions about subject 20 in concert with corresponding movement ofradiation source 22.FIG. 1 shows a representative sampling ofDR detector 24 positions to illustrate how these images are obtained relative to the position ofsubject 20. Once the needed 2-D projection images are captured in this sequence, a suitable imaging algorithm, such as filtered back projection (FBP) or other conventional technique, is used for generating the 3-D volume image. Image acquisition and program execution are performed by acomputer 30 or by a networked group ofcomputers 30 that are in image data communication withDR detectors 24. Image processing and storage is performed using a computer-accessible memory 32. The 3-D volume image can be presented on adisplay 34. - FBP is a discrete implementation of an analytic model that assumes that CT transmission measurements are linear functions of the attenuation line integrals along the corresponding primary photon trajectories through the subject and are noiseless. When scanning subjects comprised of anatomically native materials under normal conditions, relatively simple corrections to the raw projection data are sufficient to assure that these assumptions (i.e. linear relationship) are at least approximately true. This treatment allows acquisition and accurate volume reconstruction without visually observable artifacts.
- In the image processing context of the present disclosure, “rendering” is the active process of generating and forming an image for display and generating the pattern of signals needed for displaying it to a user. Image data content that is used for rendering can be transformed from a 2-D or 3-D image model (or models), typically stored as scene content in some type of scene file, into suitable patterns of light energy that are emitted from a display screen.
- Embodiments of the present disclosure provide approaches for reducing patient exposure and reducing artifacts by acquiring a number N of actual projection images that represents only a partial percentage of the full number (N+M) of X-ray projection images that would otherwise be needed to for artifact-free reconstruction. Then, to obtain an ordered set of projection images that is sufficient to be used for tomographic reconstruction, a number of synthesized projection images M are generated and added to the acquired set of N projection images. This forms a larger set N+M of projection images, effectively providing a reduced angular spacing Δα′ between adjacent projection images in the sequence used for reconstruction. The result is an enhanced 3-D reconstruction that can be substantially artifact-free.
- It should be noted that the number M of synthesized projection images that is generated can be smaller than, equal to, or even larger than the number N of actual X-ray projection images. Examples given following show generation of a single synthesized image M corresponding to an angle midway between the corresponding angles of first and second adjacent X-ray projection images n and (n+1); the image obtained at this midway angle is represented herein using the notation n.5. Thus, in the examples shown here, M<N. However, multiple synthesized images M can be generated for angles between any two adjacent X-ray projection images n and (n+1). For example, using the same processing approach, two synthesized images M can be generated between any two adjacent X-ray projection images n and (n+1).
- An embodiment of the present disclosure provides a method for improving reconstruction results without additional exposure to the patient by generating synthesized images that are formed using partitioning, forward projection, and an image data mapping sequence. The synthetic images can then be used to supplement the image content provided from actual acquired projection images without additional exposure to the patient.
-
FIG. 2 is a logic flow diagram showing an overview of the sequence for forming and using synthetic images for enhanced volume reconstruction. In an acquisition step S200, an ordered set having a number N of actual X-ray projection images {A0, A1, A2, A3, . . . }, termed the “acquired” images in the context of the present disclosure, is obtained. Each X-ray projection image in the set is acquired at one of a corresponding sequence of N capture angles α as indicated inFIG. 1 . The ordered set of N acquired X-ray projection images can be considered ordered according to acquisition angle. In practice, X-ray projection images can be acquired in any order with respect to angle; the ordered set structure is convenient arrangement for showing the acquisition sequence and for illustrating the processing that follows, used to generate synthesized images. - In the example process shown in
FIG. 2 , every acquired X-ray projection image An has a corresponding capture angle αn and is adjacent to a previous acquired X-ray projection image (n−1) with corresponding capture angle αn−1 and to a subsequent or following acquired X-ray projection image (n+1) with corresponding capture angle αn+1. Angular spacing between any two adjacent acquired X-ray projection images An in the sequence can be expressed as Δα and can be at equal intervals. From the ordered set of acquired projection images, an initial reconstruction step S210 forms a 3-D volume image using an analytic reconstruction method such as filtered back projection (FBP), or using an iterative reconstruction process, or using some other suitable reconstruction algorithm. - In order to generate an expanded set of projection images for reconstruction, the procedural sequence of
FIG. 2 performs a forward projection step S220 that generates an intermediate set of N+M forward projection images, shown inFIG. 2 and described herein as {F0, F0.5, F1, F1.5, F2, F2.5, F3 . . . } at corresponding angles {0α 0.5α, 1α, 1.5α, 2α, 2.5α, 3α, . . . }. As described previously, the notation using “0.5” indicates angles between the acquisition angles α, such as half-way between. Forward projection thus generates projection images at each of the angles α corresponding to the N actual X-ray projections and adds M additional projection images at intermediate angles, such as at angles n.5α, half-way between the original angles α. It should be emphasized that the intermediate angles can be at other than half-way between acquisition angle, such as at ⅓ intervals or at some arbitrary angular offset between the nearest acquired images. - Continuing with the
FIG. 2 sequence, the process then proceeds to re-build a set of M synthesized projection images {S0.5, S1.5, S2.5, S3.5, . . . } using the M added projection images as templates in a synthesized image generation step S230. A sequence for generating the synthesized images using decomposition or partitioning of the volume is described in more detail subsequently. - To continue the sequence of
FIG. 2 , a final reconstruction step S250 executes, generating an enhanced 3-D volume reconstruction using the original acquired X-ray projection images {A0, A1, A2, A3, . . . } supplemented by the set of synthesized images {S0.5, S1.5, S2.5, S3.5, . . . } that have been formed using the process described subsequently. The final reconstruction step S240 thus uses the combined set {A0, S0.5, A1, S1.5, A2, S2.5, A3, S3.5, . . . } for forming the enhanced 3-D reconstruction. The data can then be displayed, rendered for viewing in a display step S250. For step S250, rendering relates to the selection, processing, and display of a particular 2-D view of the reconstructed 3-D image scene content. -
FIG. 3A is a schematic diagram that shows supplementing the ordered set of N acquired X-ray projection images An with a set of M synthesized images Sn. -
FIG. 3B is a schematic diagram that show's a combined set of acquired X-ray and synthesized images used for volume image reconstruction. - While the process described with reference to
FIGS. 2-3B provides suitable image content for volume reconstruction with reduced exposure to the patient, there is room for improvement. One inherent difficulty that impacts volume reconstruction relates to relative movement of volume content in different directions. -
FIG. 4A schematically shows, as a top view, a simplified model used to describe the Radon transform and the resulting sinogram that is obtained by combining successive images, each representative of the amount of energy that passes through object O ontodetector 24.FIG. 4A shows the transform portion obtained for a projection image acquired at an angle of θ1.FIG. 4B shows the transform portion that is obtained at angle θ2. As shown inFIGS. 4A and 4B , rotation of theradiation source 22 anddetector 24 from angle θ1 to θ2 changes the relative positions of points A and B in the sequence of projection images that are acquired. As thesource 22/detector 24 pair rotates from angle θ1 to θ2, the radiation energy for point A appears to move toward the left. Conversely, as thesource 22 moves over the same angular range, point B appears to move toward the right. At angles near θ2, the image data from points A and B are moving in different directions relative to the projection images captured atdetector 24. This change in spatial relationship for the obtained data, varying with the location of anatomy relative to the rotational geometry of the CT or CBCT scan, can complicate the task of volume image reconstruction, since different portions of the anatomy shift in position in different directions and at different rates from one projection image to the next. -
FIG. 5 shows the changing spatial relationships between two imaged points, such as those described with reference toFIGS. 4A and 4B , in a sinogram view. For this view, the sinogram can be considered as a stacked collection of projections of a single slice of the subject, viewed from the side. The rotation axis is designated as R. The sinogram SN1 for a first object point shows that the first object point is at the front of each projection, and moving in one direction, for angles over the range indicated J1 and at the back of each projection, moving in the opposite direction, for angles in the range J2. Sinogram SN2 for a second object point shows that the second object point is at the rear of each projection, and moving in one direction for angles over the range indicated K1 and at the back of each projection over the range of angles K2. Over the range of angles K3, the second object point is once again at the front of each projection. - At least one embodiment of the present disclosure addresses the difficulties that are inherent to the scan geometry, as described with reference to
FIGS. 4A-5 , by using a geometric partitioning of the volume space. For geometric partitioning, the positions of the X-ray source and detector relative to the subject volume are considered and the volume is partitioned to help compensate for the relative movement conflicts noted with reference toFIGS. 4A, 4B, and 5 . - Referring to the perspective view of
FIG. 6A and top view ofFIG. 6B , one partitioning approach divides the volume in half for each relative position of the source and detector. That portion of the imaged volume V that is forward of the central axis of rotation with reference to the source, and is defined by a bisecting plane P that is orthogonal to a central ray of illumination, forms one partial sub-volume. The portion of the imaged volume that lies behind the bisecting plane P is a second partial sub-volume for the projection images at that angle. In the geometric partitioning logic, the partial sub-volume varies at each projection angle. - Partitioning can also address difficulties with various types of image content. Another type of partitioning that can be applied to the volume and used herein relates to voxel content or density data value for the imaged object. With this object-based partitioning approach, the partitioning is static, so that the same voxels occupy each sub-volume for every source-to-detector acquisition angle. Thus, for example, metal content can form a sub-volume for separate forward projection processing. Other known components, such as bone or soft tissue, can similarly be used to define sub-volumes. A static partitioning defines a fixed portion of the imaged volume, which can include any imageable portion of the object space, the cylindrical volume that lies between the source and detector as they orbit the subject.
- In the partitioning sequence, a weight can be computed for each calculated signal needed to form the forward projection images Fn. According to an embodiment that uses geometric partitioning with two partitions as shown in
FIGS. 6A and 6B , a front and rear weight are generated for each pixel: -
- The front and rear partitions can then be individually identified and processed using the computed weights.
- In order to simplify discussion, the description given herein is directed to embodiments in which only two partitions are generated. However, it should be noted that partitioning can generate two or more sub-volumes for the imaged content and that similar processing applies for each of two or more sub-volumes.
- The calculated weight values are then used to determine how each pixel value is assigned relative to the partitioning.
- Using Partitioning when Generating Synthetic Images
- An embodiment of the present disclosure provides a method for improving reconstruction results without additional exposure to the patient using the synthetic image generation process described with reference to
FIGS. 2-3B along with volume decomposition or partitioning as described above and with reference toFIGS. 4A-6B . - The logic flow diagram of
FIG. 7 shows a sequence for volume image reconstruction using partitioning. In an acquisition step S700, an ordered set having a number N actual X-ray projection images {A0, A1, A2, A3, . . . }, termed the “acquired” images in the context of the present disclosure, is obtained. Each X-ray projection image in the set is acquired at one of a corresponding sequence of N capture angles α as indicated inFIG. 1 . The ordered set of N acquired X-ray projection images can be considered ordered according to acquisition angle. In practice, X-ray projection images can be acquired in any order with respect to angle; the ordered set structure is convenient arrangement for showing the acquisition sequence and for illustrating the processing that follows, used to generate synthesized images. - In the example process shown in
FIG. 7 , every acquired X-ray projection image An has a corresponding capture angle αn and is adjacent to a previous acquired X-ray projection image A(n−1) with corresponding capture angle αn−1 and to a subsequent or following acquired X-ray projection image A(n+1) with corresponding capture angle αn+1. Angular spacing between any two adjacent acquired X-ray projection images An in the sequence can be expressed as Δα. From the ordered set of acquired projection images, an initial reconstruction step S710 forms a full 3-D volume image using filtered back projection (FBP), using an iterative reconstruction process, or using some other suitable reconstruction algorithm. - A partitioning step S720 then applies partitioning logic to the full 3-D volume define a set of two or more sub-volumes {SV1, SV2, . . . }. The sequence of
FIG. 7 shows processing for the complete partition of the volume, with a number p of sub-volumes. As noted previously, sub-volumes SVn can be defined using geometry of the scanning operation or using structural or content data from the reconstructed object itself. - Continuing with the
FIG. 7 sequence, within each sub-volume SVp, a forward projection step S730 generates forward projection images at angles αn corresponding to the acquisition angles used for the actual X-ray projection images. In addition, step S730 also generates a set of intermediate forward projection images, at angles between each angle αn and its respective, adjacent following angle αn+1. Intermediate forward projection images are represented at angles (n+0.5) midway between preceding and following forward projections; as noted earlier, the intermediate forward projections can be calculated at other angular positions between angle αn and its respective, adjacent following angle αn+1. - A weighting step S740 then calculates and applies a weight or weighting factor to each pixel value for the corresponding acquired image data that contributes to the sub-volume. The weighting factor is an estimate of the proportionate contribution of each acquired pixel of the sub-volume to the computed forward projection value and is based on the proportion of pixel values in the sub-volume to the total data value obtained for that pixel at the corresponding acquisition angle. A synthesized sub-volume projection calculation step S750 then uses the weighted values from step S740 to generate a synthetic projection image for each sub-volume accordingly. A combination step S760 combines the synthetic projection image results for the corresponding angle from each sub-volume to form a synthesized projection image for subsequent volume reconstruction, forming a set of M synthesized projection images, using the variable notation described previously. A reconstruction step S770 performs volume reconstruction using the initial set of N acquired projection images (corresponding to the N actual acquisition angles {0α; 1α, 2α, 3α, . . . } plus the M synthesized projection images (corresponding to acquisition angles between the set of N actual acquisition angles). A display step S780 then allows display, storage, and transmission of the volume image data, as well as display, storage, and transmission of one or more of the synthetic projection images formed using the
FIG. 7 process. - By way of illustration, the examples given herein show partitioning using front and rear segments of the volume image space; however, other partitioning models can be applied. Some of the alternate partition models are shown in
FIGS. 9A-9C .FIG. 9A shows additional segments defined within the reconstructed volume.FIG. 9B shows a top view of a partition based on relative rate of change for volume content according to distance from rotation axis R.FIG. 9C shows a wedged partition of the volume. The same processing described for partitioning into front and rear portions, as described with reference toFIGS. 6A and 6B , can be used for the geometric partitioning shown inFIGS. 9A-9C or for other geometric partitioning arrangements that define a partition having non-overlapping sub-volumes. - It should be noted that various image conditioning processes can be applied to the volume data that is used for forward projection as well as for any of the generated forward projection or synthesized image data. Image conditioning improves the appearance of a projection image or of a synthesized image and can be applied as part of display rendering in order to change image sharpening or blurring, such as to improve image patch identification and processing, for example. Image conditioning such as metal artifacts reduction can be applied to the original acquired X-ray projection images or to the volume data that is reconstructed from the original X-ray projection images.
-
FIG. 10 shows mapping and combining patch content for a candidate patch from two adjacent acquired X-ray projection images to form patch p on an intermediate image designated Fn.5 using forward projection, in order to form synthesized image Sn.5 at completion of this processing. - A synthesized image constructed using the method described herein can be formed as a mosaic of patches. Often, boundaries or seams between patches are readily discernable, causing perceptible degradation in the final reconstruction volume. One method to correct for seams is to perform a one dimensional smoothing filter across the seams. If the filter is applied uniformly, however, some sharp edges are overly smoothed and clash with the actual images.
- According to a method of the present disclosure, seam smoothing is executed with respect to the forward projected image. The synthesized image seams are smoothed by applying the following one directional algorithm vertically and horizontally:
- (i) Identify adjacent pixels p1 and p2, each belonging to a different block.
- (ii) Calculate the difference ΔSyn of the two pixels, p1 and p2.
- (iii) Calculate the difference ΔAct, of pixels in the same location in the forward projected image.
- (iii) Compare ΔSyn, and ΔAct. If |ΔSyn|<|ΔAct|, then make no change to pixel values; otherwise add one-third of the difference to p1 and remove one third of the difference from p2.
- The approach outlined above can help to dampen the seam boundary adaptively and limit image degradation that can otherwise occur with a global smoothing approach.
- To help mitigate blocking artifacts, the Applicants can alternately use an improved spatial similarity metric, formed by adding terms that encourage spatial correlation with neighboring patches. This includes two types of terms: a boundary matching energy and a flow field smoothing energy.
FIG. 11 shows the meaning of boundary matching terms DA, DR, DB, DL, used herein. - The basic energy relationship E can be expressed as follows:
-
E=block matching energy+boundary matching energy+flow field smoothing energy Equation (1): - Block matching energy can then be formulated as follows:
-
- A flow field term expresses the consistency of flow direction:
-
- where, as shown in
FIG. 11 : -
- Embodiments of the present disclosure provide an improvement over earlier interpolation methods, particularly with respect to reduction of artifacts in general, particularly view aliasing artifacts.
- By way of example,
FIG. 12A shows an acquired projection image An of a hand.FIG. 12B shows a synthesized image generated according to an embodiment of the present disclosure. The synthesized image can be displayed, stored, transmitted, or can be used with acquired projection images for volume reconstruction. - Consistent with at least one embodiment, the present disclosure utilizes a computer program with stored instructions that perform on image data accessed from an electronic memory. As can be appreciated by those skilled in the image processing arts, a computer program of an embodiment of the present invention can be utilized by a suitable, general-purpose computer system, such as a personal computer or workstation. However, many other types of computer systems can be used to execute the computer program of the present invention, including networked processors. The computer program for performing the method of the present invention may be stored in a computer readable storage medium. This medium may comprise, for example; magnetic storage media such as a magnetic disk such as a hard drive or removable device or magnetic tape; optical storage media such as an optical disc, optical tape, or machine readable bar code; solid state electronic storage devices such as random access memory (RAM), or read only memory (ROM); or any other physical device or medium employed to store a computer program. The computer program for performing the method of the present invention may also be stored on computer readable storage medium that is connected to the image processor by way of the internet or other communication medium. Those skilled in the art will readily recognize that the equivalent of such a computer program product may also be constructed in hardware.
- It should be noted that the term “memory”, equivalent to “computer-accessible memory” in the context of the present disclosure, can refer to any type of temporary or more enduring data storage workspace used for storing and operating upon image data and accessible to a computer system, including a database. The memory could be non-volatile, using, for example, a long-term storage medium such as magnetic or optical storage. Alternately, the memory could be of a more volatile nature, using an electronic circuit, such as random-access memory (RAM) that is used as a temporary buffer or workspace by a microprocessor or other control logic processor device. Displaying an image requires memory storage. Display data, for example, is typically stored in a temporary storage buffer that is directly associated with a display device and is periodically refreshed as needed in order to provide displayed data. This temporary storage buffer can also be considered to be a memory, as the term is used in the present disclosure. Memory is also used as the data workspace for executing and storing intermediate and final results of calculations and other processing. Computer-accessible memory can be volatile, non-volatile, or a hybrid combination of volatile and non-volatile types.
- It will be understood that the computer program product of the present invention may make use of various image manipulation algorithms and processes that are well known. It will be further understood that the computer program product embodiment of the present invention may embody algorithms and processes not specifically shown or described herein that are useful for implementation. Such algorithms and processes may include conventional utilities that are within the ordinary skill of the image processing arts. Additional aspects of such algorithms and systems, and hardware and/or software for producing and otherwise processing the images or co-operating with the computer program product of the present invention, are not specifically shown or described herein and may be selected from such algorithms, systems, hardware, components and elements known in the art.
- The invention has been described in detail with particular reference to a presently preferred embodiment, but it will be understood that variations and modifications can be effected within the spirit and scope of the invention. The presently disclosed embodiments are therefore considered in all respects to be illustrative and not restrictive. The scope of the invention is indicated by the appended claims, and all changes that come within the meaning and range of equivalents thereof are intended to be embraced therein.
Claims (18)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/472,613 US10089758B1 (en) | 2017-03-29 | 2017-03-29 | Volume image reconstruction using projection decomposition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/472,613 US10089758B1 (en) | 2017-03-29 | 2017-03-29 | Volume image reconstruction using projection decomposition |
Publications (2)
Publication Number | Publication Date |
---|---|
US10089758B1 US10089758B1 (en) | 2018-10-02 |
US20180286087A1 true US20180286087A1 (en) | 2018-10-04 |
Family
ID=63638734
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/472,613 Expired - Fee Related US10089758B1 (en) | 2017-03-29 | 2017-03-29 | Volume image reconstruction using projection decomposition |
Country Status (1)
Country | Link |
---|---|
US (1) | US10089758B1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11375964B2 (en) * | 2018-05-28 | 2022-07-05 | Riken | Acquisition method, acquisition device, and control program for tomographic image data by means of angular offset |
US11557038B2 (en) | 2019-10-17 | 2023-01-17 | Electronics And Telecommunications Research Institute | Apparatus and method for x-ray data generation |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023146545A1 (en) * | 2022-01-31 | 2023-08-03 | Hewlett-Packard Development Company, L.P. | Transforming object data models |
CN114712728A (en) * | 2022-03-31 | 2022-07-08 | 西安大医集团股份有限公司 | Projection parameter correction method and computer equipment |
Family Cites Families (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5270926A (en) | 1990-12-21 | 1993-12-14 | General Electric Company | Method and apparatus for reconstructing a three-dimensional computerized tomography (CT) image of an object from incomplete cone beam projection data |
US5999587A (en) | 1997-07-03 | 1999-12-07 | University Of Rochester | Method of and system for cone-beam tomography reconstruction |
US7120283B2 (en) * | 2004-01-12 | 2006-10-10 | Mercury Computer Systems, Inc. | Methods and apparatus for back-projection and forward-projection |
US7602879B2 (en) | 2005-11-02 | 2009-10-13 | Wisconsin Alumni Research Foundation | Method for increasing the resolution of a CT image during image reconstruction |
EP2083390B1 (en) * | 2008-01-24 | 2016-06-01 | Brainlab AG | Method for segmenting a 3D image data set, accompanying computer program product and accompanying system |
JP5334423B2 (en) * | 2008-02-13 | 2013-11-06 | 株式会社日立メディコ | Image processing apparatus, image processing method, and program |
DE102008009128B4 (en) * | 2008-02-14 | 2014-11-06 | Siemens Aktiengesellschaft | Tomosynthetic image reconstruction method and diagnostic device using this method |
US9402590B2 (en) * | 2008-10-15 | 2016-08-02 | Toshiba Medical Systems Corporation | Three-dimensional image processing apparatus and X-ray diagnostic apparatus |
HU0900478D0 (en) * | 2009-07-31 | 2009-09-28 | Holografika Hologrameloeallito | Method and apparatus for displaying 3d images |
DE102009053664A1 (en) * | 2009-11-17 | 2011-05-19 | Ziehm Imaging Gmbh | Method for the empirical determination of a correction function for the correction of radiation hardening and stray radiation effects in projection radiography and in computed tomography |
DE102010013498B4 (en) * | 2010-03-31 | 2016-06-02 | Siemens Aktiengesellschaft | Method for determining three-dimensional volume data, imaging device and data carrier |
CN102793552B (en) * | 2011-05-23 | 2014-05-21 | 北京东方惠尔图像技术有限公司 | Computed tomography (CT) image acquisition device and CT scanning imaging system |
US9659390B2 (en) * | 2011-10-28 | 2017-05-23 | Carestream Health, Inc. | Tomosynthesis reconstruction with rib suppression |
WO2014000279A1 (en) * | 2012-06-29 | 2014-01-03 | 北京东方惠尔图像技术有限公司 | 3d image generation method and device for g-arm x-ray machine and g-arm x-ray machine |
US9224216B2 (en) | 2013-07-31 | 2015-12-29 | Kabushiki Kaisha Toshiba | High density forward projector for spatial resolution improvement for medical imaging systems including computed tomography |
US11045151B2 (en) * | 2013-11-18 | 2021-06-29 | Varian Medical Systems, Inc. | Cone-beam computed tomography imaging devices, systems, and methods |
US10083543B2 (en) * | 2013-12-20 | 2018-09-25 | Carestream Dental Technology Topco Limited | Metal artifacts reduction for cone beam CT using image stacking |
FR3020871B1 (en) * | 2014-05-07 | 2017-08-11 | Areva | DEVICE AND METHOD FOR THREE-DIMENSIONAL RECONSTRUCTION OF A SCENE BY IMAGE ANALYSIS |
US20150348259A1 (en) * | 2014-06-03 | 2015-12-03 | Carestream Health, Inc. | Quantitative method for 3-d bone mineral density visualization and monitoring |
DE102014213408B4 (en) * | 2014-07-10 | 2018-10-25 | Siemens Healthcare Gmbh | Method for determining a three-dimensional model data set of a blood vessel system comprising at least one vessel segment |
JP6009502B2 (en) * | 2014-07-29 | 2016-10-19 | 株式会社ソニー・インタラクティブエンタテインメント | Information processing apparatus and information processing method |
JP6293619B2 (en) * | 2014-08-28 | 2018-03-14 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Image processing method, apparatus, and program |
US9626587B2 (en) * | 2014-08-28 | 2017-04-18 | Toshiba Medical Systems Corporation | Iterative reconstruction scheme for phase contrast tomography |
US9934597B2 (en) * | 2014-09-11 | 2018-04-03 | Carestream Health, Inc. | Metal artifacts reduction in cone beam reconstruction |
US9836861B2 (en) * | 2014-12-12 | 2017-12-05 | Samsung Electronics Co., Ltd. | Tomography apparatus and method of reconstructing tomography image |
US10796463B2 (en) * | 2015-08-05 | 2020-10-06 | General Electric Company | Tomographic imaging for time-sensitive applications |
US10395395B2 (en) * | 2016-03-15 | 2019-08-27 | Carestream Health, Inc. | Virtual projection images for tomosynthesis artifact reduction |
MA44842A (en) * | 2016-05-06 | 2021-05-26 | L3 Security & Detection Systems Inc | SYSTEMS AND METHODS FOR GENERATING PROJECTION IMAGES |
US9965875B2 (en) * | 2016-06-21 | 2018-05-08 | Carestream Health, Inc. | Virtual projection image method |
-
2017
- 2017-03-29 US US15/472,613 patent/US10089758B1/en not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11375964B2 (en) * | 2018-05-28 | 2022-07-05 | Riken | Acquisition method, acquisition device, and control program for tomographic image data by means of angular offset |
US11557038B2 (en) | 2019-10-17 | 2023-01-17 | Electronics And Telecommunications Research Institute | Apparatus and method for x-ray data generation |
Also Published As
Publication number | Publication date |
---|---|
US10089758B1 (en) | 2018-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9245320B2 (en) | Method and system for correcting artifacts in image reconstruction | |
CN107427274B (en) | Tomographic apparatus and method for reconstructing tomographic image thereof | |
US8121374B2 (en) | Bands artifact reduction for cardiac CT imaging | |
EP2030170B1 (en) | Dynamic computed tomography imaging | |
US10789738B2 (en) | Method and apparatus to reduce artifacts in a computed-tomography (CT) image by iterative reconstruction (IR) using a cost function with a de-emphasis operator | |
US9600924B2 (en) | Iterative reconstruction of image data in CT | |
US10083543B2 (en) | Metal artifacts reduction for cone beam CT using image stacking | |
US10089758B1 (en) | Volume image reconstruction using projection decomposition | |
EP2633498B1 (en) | Apparatus and method for hybrid reconstruction of an object from projection data. | |
CN101433464A (en) | X-ray CT apparatus and image processing apparatus | |
US20180089864A1 (en) | Cbct image processing method | |
JP6505513B2 (en) | X-ray computed tomography imaging apparatus and medical image processing apparatus | |
US9965875B2 (en) | Virtual projection image method | |
US7272205B2 (en) | Methods, apparatus, and software to facilitate computing the elements of a forward projection matrix | |
US10966670B2 (en) | Imaging system and method for dual-energy and computed tomography | |
US8379948B2 (en) | Methods and systems for fast iterative reconstruction using separable system models | |
US20080086052A1 (en) | Methods and apparatus for motion compensation | |
CN110073412B (en) | Image noise estimation using alternating negatives | |
US10307114B1 (en) | Iterative volume image reconstruction using synthetic projection images | |
US10395395B2 (en) | Virtual projection images for tomosynthesis artifact reduction | |
US9558569B2 (en) | Method and system for substantially reducing cone beam artifacts based upon image domain differentiation in circular computer tomography (CT) | |
CN117425918A (en) | Method for CT reconstruction | |
Sun | Rigid motion correction for head CT imaging | |
Heverhagen | Physics of computed tomography scanning | |
US20110317894A1 (en) | Projection truncation processing for cbct |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: CARESTREAM HEALTH, INC., NEW YORK Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:RAY, LAWRENCE A.;SIMON, RICHARD A.;VOGELSANG, LEVON O.;SIGNING DATES FROM 20170403 TO 20170404;REEL/FRAME:041927/0228 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
AS | Assignment |
Owner name: CREDIT SUISSE AG, CAYMAN ISLANDS BRANCH, NEW YORK Free format text: SECURITY INTEREST;ASSIGNORS:CARESTREAM HEALTH, INC.;CARESTREAM HEALTH HOLDINGS, INC.;CARESTREAM HEALTH CANADA HOLDINGS, INC.;AND OTHERS;REEL/FRAME:048077/0587 Effective date: 20190114 Owner name: CREDIT SUISSE AG, CAYMAN ISLANDS BRANCH, NEW YORK Free format text: SECURITY INTEREST;ASSIGNORS:CARESTREAM HEALTH, INC.;CARESTREAM HEALTH HOLDINGS, INC.;CARESTREAM HEALTH CANADA HOLDINGS, INC.;AND OTHERS;REEL/FRAME:048077/0529 Effective date: 20190114 |
|
CC | Certificate of correction | ||
FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
AS | Assignment |
Owner name: CARESTREAM HEALTH WORLD HOLDINGS LLC, NEW YORK Free format text: RELEASE OF SECURITY INTEREST IN INTELLECTUAL PROPERTY (FIRST LIEN);ASSIGNOR:CREDIT SUISSE AG, CAYMAN ISLANDS BRANCH;REEL/FRAME:061683/0529 Effective date: 20220930 Owner name: CARESTREAM HEALTH ACQUISITION, LLC, NEW YORK Free format text: RELEASE OF SECURITY INTEREST IN INTELLECTUAL PROPERTY (FIRST LIEN);ASSIGNOR:CREDIT SUISSE AG, CAYMAN ISLANDS BRANCH;REEL/FRAME:061683/0529 Effective date: 20220930 Owner name: CARESTREAM HEALTH CANADA HOLDINGS, INC., NEW YORK Free format text: RELEASE OF SECURITY INTEREST IN INTELLECTUAL PROPERTY (FIRST LIEN);ASSIGNOR:CREDIT SUISSE AG, CAYMAN ISLANDS BRANCH;REEL/FRAME:061683/0529 Effective date: 20220930 Owner name: CARESTREAM HEALTH HOLDINGS, INC., NEW YORK Free format text: RELEASE OF SECURITY INTEREST IN INTELLECTUAL PROPERTY (FIRST LIEN);ASSIGNOR:CREDIT SUISSE AG, CAYMAN ISLANDS BRANCH;REEL/FRAME:061683/0529 Effective date: 20220930 Owner name: CARESTREAM HEALTH, INC., NEW YORK Free format text: RELEASE OF SECURITY INTEREST IN INTELLECTUAL PROPERTY (FIRST LIEN);ASSIGNOR:CREDIT SUISSE AG, CAYMAN ISLANDS BRANCH;REEL/FRAME:061683/0529 Effective date: 20220930 Owner name: CARESTREAM HEALTH WORLD HOLDINGS LLC, NEW YORK Free format text: RELEASE OF SECURITY INTEREST IN INTELLECTUAL PROPERTY (SECOND LIEN);ASSIGNOR:CREDIT SUISSE AG, CAYMAN ISLANDS BRANCH;REEL/FRAME:061683/0681 Effective date: 20220930 Owner name: CARESTREAM HEALTH ACQUISITION, LLC, NEW YORK Free format text: RELEASE OF SECURITY INTEREST IN INTELLECTUAL PROPERTY (SECOND LIEN);ASSIGNOR:CREDIT SUISSE AG, CAYMAN ISLANDS BRANCH;REEL/FRAME:061683/0681 Effective date: 20220930 Owner name: CARESTREAM HEALTH CANADA HOLDINGS, INC., NEW YORK Free format text: RELEASE OF SECURITY INTEREST IN INTELLECTUAL PROPERTY (SECOND LIEN);ASSIGNOR:CREDIT SUISSE AG, CAYMAN ISLANDS BRANCH;REEL/FRAME:061683/0681 Effective date: 20220930 Owner name: CARESTREAM HEALTH HOLDINGS, INC., NEW YORK Free format text: RELEASE OF SECURITY INTEREST IN INTELLECTUAL PROPERTY (SECOND LIEN);ASSIGNOR:CREDIT SUISSE AG, CAYMAN ISLANDS BRANCH;REEL/FRAME:061683/0681 Effective date: 20220930 Owner name: CARESTREAM HEALTH, INC., NEW YORK Free format text: RELEASE OF SECURITY INTEREST IN INTELLECTUAL PROPERTY (SECOND LIEN);ASSIGNOR:CREDIT SUISSE AG, CAYMAN ISLANDS BRANCH;REEL/FRAME:061683/0681 Effective date: 20220930 |
|
LAPS | Lapse for failure to pay maintenance fees |
Free format text: PATENT EXPIRED FOR FAILURE TO PAY MAINTENANCE FEES (ORIGINAL EVENT CODE: EXP.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20221002 |