Classifications

• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T7/00—Image analysis
  • G06T7/10—Segmentation; Edge detection
  • G06T7/11—Region-based segmentation
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T7/00—Image analysis
  • G06T7/10—Segmentation; Edge detection
  • G06T7/12—Edge-based segmentation
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
  • A61B6/03—Computed tomography [CT]
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
  • A61B6/4035—Arrangements for generating radiation specially adapted for radiation diagnosis the source being combined with a filter or grating
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/48—Diagnostic techniques
  • A61B6/482—Diagnostic techniques involving multiple energy imaging
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
  • A61B6/51—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for dentistry
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
  • A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
  • A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
  • A61B6/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T12/00—Tomographic reconstruction from projections
  • G06T12/10—Image preprocessing, e.g. calibration, positioning of sources or scatter correction
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
  • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
  • A61B6/5252—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data removing objects from field of view, e.g. removing patient table from a CT image
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; 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 OR CALCULATING; 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/10116—X-ray image
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T2207/00—Indexing scheme for image analysis or image enhancement
  • G06T2207/30—Subject of image; Context of image processing
  • G06T2207/30004—Biomedical image processing
  • G06T2207/30036—Dental; Teeth

Definitions

• the present invention relates to the field of X-ray imaging and, more particularly, to using multi-energy X-ray scans to segment teeth in a reconstructed image.
• Teeth segmentation identifies the voxels that belong or correspond to teeth in a three dimensional (3D) reconstructed image of an X-ray scan. More specifically, teeth segmentation may identify a part of the image that comprises teeth, identify individual teeth in the image, and identify parts of a tooth in an image. Different dental applications require different levels of segmentation and it is highly desirable that teeth segmentation be as automatic as possible, requiring little or no human interaction.
• teeth segmentation in reconstructions of X-ray scans is very difficult.
• Currently it is not possible to segment teeth fully automatically in all cases. The reason for this is generally two-fold.
• First it is difficult to distinguish between tooth roots and surrounding alveolar bone because they have similar material composition.
• Second, reconstructions have artifacts due to beam hardening, the presence of metal, and scatter which cause material of uniform material composition to appear non-uniform in a reconstructed image.
• the present invention comprises methods for producing a three-dimensional (3D), segmented representation of one or more teeth using multi-energy X-ray spectra and/or a multi-energy X-ray scanner.
• tooth segmentation may be performed automatically or semi-automatically for images which are reconstructed from a multi-energy X-ray scan.
• the results of the tooth segmentation may be represented in a number of ways.
• the tooth segmentation results may be represented as an image map which identifies voxels which are within a tooth.
• the tooth segmentation results may be represented in the form of a three- dimensional (3D) grid or any other representation of a three-dimensional (3D) spatial region.
• the example embodiments herein describe the present invention in connection with a dual X-ray spectrum scanner and dual X-ray spectra, but other example embodiments include the use of multiple X-ray spectrum scanners and more than two X-ray spectra. Therefore, the scope of the present invention is not limited to a dual X-ray spectrum scanner or dual X-ray spectra.
• a dual energy scan can be performed by changing the source voltage and/or filtration of the X-ray source during the scan (fast switching) or performing two separate scans with different source voltage and/or filtration. Alternatively, a dual energy scan may be performed simultaneously with two different sources and detectors.
• One example embodiment of the present invention includes and uses an energy discriminating photon counting detector with at least two energy bins.
• the measured data of an X-ray scan is the exposure value of each pixel of the detector.
• the exposure value is related to the X-ray attenuation as the photons travel along a line from the source to the detector.
• the measured data is generally corrected for X-ray source non uniformity and detector response (flat field correction) and detector defects before it is used for image reconstruction.
• the measured data at all detector pixels is often referred to as an X-ray projection because it is a radiographic projection of an object onto the detector.
• a scan consists of a series of projections at different source and detector locations. Often the source and detector move about an axis-of-rotation (AOR). The patient is positioned so that the AOR is located at the center of a region-of-interest (ROI). For dental applications, the ROI is usually within the dental arch.
• two sets of projections are collected. The two projection sets may be for the same or different X-ray paths (source/detector locations).
• the methods of the present invention provide the ability to process the scan data before or during the tooth segmentation process so that segmentation can be performed with little or no human intervention. This is accomplished at least in part and alone or in combination by reconstructing the dual energy data in a way that increases the contrast between a tooth and other material in the scanned object such as soft tissue and bone, by reducing artifacts that are caused by the change in X-ray spectrum as it propagates through the scanned material (beam hardening), by reducing artifacts that are due to photon starvation and beam hardening caused by the presence of metal and other dense material, and by reducing artifacts that are caused by X-ray scatter.
• the methods of the present invention include combining the measured scan data at the two X-ray spectra or combining the reconstructions of the scan data at the two X-ray spectra in order to provide an image which is better suited for tooth segmentation than separate reconstructions at each of the two X-ray spectra. Furthermore, the methods provide an ability to incorporate the dual energy scan data into the teeth segmentation process so as to optimize the utility of the dual-spectral data.
• Figs. 1A, 1B and 1C display pictorial images of a slice of a reconstruction of a tooth along with a contour that corresponds to the outline of a segmented region.
• Figs. 2A and 2B display pictorial images of a slice of a reconstruction of a tooth having a metal filling.
• Fig. 3 displays a schematic diagram of an x-ray scanner positioned relative to a patient.
• Fig. 4 displays a flowchart representation of a method for tooth segmentation in accordance with a first example embodiment of the present invention.
• Fig. 5 displays a flowchart representation of a method for tooth segmentation in accordance with a second example embodiment of the present invention.
• Fig. 6 displays a flowchart representation of a method for tooth segmentation in accordance with a third example embodiment of the present invention.
• Fig. 7 displays a flowchart representation of a method for evaluating the quality of tooth segmentation according to an example embodiment of the present invention.
• Fig. 8 displays a flowchart representation of a method for determining if under segmentation is present in the tooth segmentation results.
• Figs. 1A, 1B, and 1C display images 100, 110, 120 illustrating two problems that are solved by the methods of the present invention.
• Image 100 is a slice of a reconstruction of a cone beam scan of a dental arch. Tooth root 102 and surrounding bone 104 appear identical in this image.
• Image 110 shows the result of segmenting tooth root 102. Because of the inability to distinguish between root 102 and bone 104, the segmentation fails and the segmented region includes not only root 102, but also bone 104 and roots of an adjacent tooth 106.
• Image 120 in Fig. 1C shows a slice of a reconstruction with tooth 122 along with contour 124 which is the outline of the segmented region.
• Region 126 of the tooth is missing from the segmented region because of imaging artifacts which cause tooth 122 to appear non- uniform in the reconstruction.
• the image artifact may be caused by beam hardening.
• image 200 is a slice of a reconstruction having a tooth 202 which contains a metal filling 204.
• the dark areas in the tooth 206 are an artifact which is caused by the metal filling.
• Image 220 shows a three-dimensional (3D) representation of the results of segmenting tooth 202 and adjacent teeth.
• the segmentation of the tooth 222 is missing at least one root because of the metal artifacts present in the reconstruction.
• Fig. 3 displays an X-ray scanner. X-rays from source 300 pass through collimator
• the filter 310 modifies the X-ray energy spectrum and can be used, along with modification of the source's voltage, to choose the X-ray spectrum.
• the X-rays pass through dental arch region-of-interest (ROI) 308 in the patient's head 304 and are incident on detector 306.
• ROI region-of-interest
• the detector 306 is an energy discriminating photon counting detector
• the dual energy scans are actually a single scan. Otherwise, the voltage of source 300 and filter 310 is changed within a single scan or by performing two scans.
• the essential outcome of a dual energy scan are two sets of projections for different X-ray spectra which can be used to reconstruct a three dimensional (3D) image of a ROI.
• Fig. 4 One example embodiment of this invention is shown in Fig. 4.
• the dual energy scan is described at a low energy scan 400 and a high energy scan 402. This means that the average X-ray photon energy of scan 400 is lower than scan 402.
• scan 400 is the photon count in the low energy bin
• scan 402 is the photon count in the high energy bin.
• the low energy scan data 404 and high energy scan data 406 are combined in step 408.
• the purpose of the step 408 is to combine the low and high energy scan data so that when the data is reconstructed in step 410, the reconstruction has reduced artifacts and increased material contrast.
• the low energy a h and high energy an scan data may be combined using a polynomial function
• step 408 the low and high data may be combined in several different ways.
• the data is combined to enhance the contrast between tooth roots and surrounding alveolar bone.
• the data may be combined in another way to enhance the contrast between tooth and soft tissue such as the surrounding gum.
• pi and pi correspond to line integrals of material density for two basis materials.
• Preferred basis materials for image decomposition are soft tissue and hydroxyapatite, although other materials can be used.
• contrast between different materials is not only related to the difference in code values of the materials in the reconstruction, but also to the variation and noise in the code values of each material.
• One measure of the contrast between two materials is the Mahalanobis distance between the distribution of code values of the materials.
• the combined scan data 408 is used in step 410 to create reconstructions that are artifact reduced and preferably artifact free.
• this reconstruction is a virtual monochromatic reconstruction meaning that it appears as if it is reconstructed from a scan using a monochromatic X-ray source.
• Such a reconstruction is free of beam hardening artifacts.
• the monochromatic energy can be set to maximize the ability to differentiate between materials such as tooth, bone, and soft tissue to enable the subsequent segmentation step 412.
• step 412 one or more teeth in the reconstruction are segmented. This means that each tooth is distinguished from surrounding bone and tissue and from other teeth. This may also include segmenting individual parts of a tooth including crown, enamel, dentin, neck, pulp, and root.
• This step may use any image segmentation method including neural nets, clustering, active contours, snakes, thresholding, and level sets.
• the result of this step is a three- dimensional (3D) representation of teeth 414 which may take the form of a three-dimensional (3D) image mask, a surface map, a mesh, or any other means of representing a region in space.
• Fig. 5 shows another example embodiment of the present invention. This example embodiment of the invention is most appropriate when the low and high energy scans correspond to different X-ray paths through the object.
• the low energy scan 500 produces low energy scan data 504 and high energy scan 502 produces high energy scan data 506.
• the low energy scan data is reconstructed in step 508 and the high energy scan data in step 510.
• the low and high energy reconstructions are combined in order to facilitate the tooth segmentation step 514 which results in a three-dimensional (3D) representation of teeth 516.
• the methods of the present invention use multi-energy scans to improve tooth segmentation, even in the case of truncated projections, by including a way to evaluate the quality of tooth segmentation and to feedback the results into the step in which scan data or reconstructions at two or more energies is combined so that the processing of the scan data and/or reconstruction can be modified in order to facilitate tooth segmentation.
• a low energy scan 600 and high energy scan 602 are performed to generate low energy 604 and high energy 606 scan data.
• the combined scan data is processed in step 608 and reconstructed in step 610.
• the quality of tooth segmentation in step 612 is evaluated in step 613 and the results are input to step 608 in which the scan data is reprocessed in order to improve the segmentation results.
• Step 613 can take many different forms. Two example embodiments are described in detail below, but the essence of this step is to provide a measure of teeth segmentation quality. Step 613 may include several quality measures.
• Fig. 7 illustrates an image uniformity quality evaluation method that determines if the teeth in a reconstruction are being properly segmented. If not, steps 608 and 610 are modified, for example, to improve beam hardening and metal artifact correction.
• Fig. 7 displays an example of the steps that occur within step 613. These steps are based on the fact that teeth are usually convex in shape. If the segmentation results are concave this indicates that the dual energy scan data was not sufficiently processed to remove artifacts.
• step 700 the degree of convexity of the segmented teeth is calculated. If the convexity is sufficiently low, then the segmentation process is complete and no further processing is necessary. Otherwise, in step 702 concave regions are identified.
• Contour 124 in Fig. 1B corresponds to an example of a concave region which shows over-segmentation. A concave region may also indicate under-segmentation in which multiple teeth are segmented as a single tooth.
• step 704 the code value distribution of one or more reconstructions inside and outside the segmented region are evaluated.
• the reconstruction is a virtual monochromatic reconstruction.
• a difference in code value distributions may indicate that scan data processing in step 608 was not sufficient to produce a virtual monochromatic reconstruction which is completely free of beam hardening artifacts.
• step 706 it is determined if additional scan data processing is necessary. It is possible that convex segmented regions correspond to locations at which a tooth is forming into multiple roots. It is a part of this step to distinguish between concavity due to insufficient artifact removal and variation in tooth shape.
• the reconstruction code values can take several forms.
• the code values may be X-ray attenuation coefficients in units of cm 1 .
• the code values may be in Hounsfield units.
• the code values may not measure a physical property of the scanned object, but are nevertheless useful for tooth segmentation.
• Fig. 8 displays another set of steps that can occur within and form step 613, possibly in parallel with the processing steps in Fig. 7.
• the steps in Fig. 8 are directed at reducing the under- segmentation problem that is illustrated in image 110 in Fig. 1.
• step 800 the tooth segmentation in adjacent axial slices are compared. A large change in segmentation, for example, as measured by the Sorensen-Dice coefficient may indicate that tooth segmentation may be extending into surrounding bone or multiple teeth are segmented as one.
• step 802 regions of under- segmentation are identified.
• step 806 it is determined if additional processing is necessary. If so, steps 608, 610, and 612 are repeated in such a way to boost material differentiation.

Landscapes

• Engineering & Computer Science (AREA)
• Health & Medical Sciences (AREA)
• Life Sciences & Earth Sciences (AREA)
• Medical Informatics (AREA)
• Physics & Mathematics (AREA)
• Biomedical Technology (AREA)
• Molecular Biology (AREA)
• Biophysics (AREA)
• Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
• Optics & Photonics (AREA)
• Pathology (AREA)
• Radiology & Medical Imaging (AREA)
• Veterinary Medicine (AREA)
• Heart & Thoracic Surgery (AREA)
• High Energy & Nuclear Physics (AREA)
• Surgery (AREA)
• Animal Behavior & Ethology (AREA)
• General Health & Medical Sciences (AREA)
• Public Health (AREA)
• Computer Vision & Pattern Recognition (AREA)
• General Physics & Mathematics (AREA)
• Theoretical Computer Science (AREA)
• Dentistry (AREA)
• Oral & Maxillofacial Surgery (AREA)
• Apparatus For Radiation Diagnosis (AREA)

Priority Applications (3)

Applications Claiming Priority (2)

Publications (1)

Family

ID=68425296

Family Applications (1)

Country Status (4)

Cited By (1)

Families Citing this family (1)

Family Cites Families (12)

• 2019
  • 2019-10-11 JP JP2021545258A patent/JP2022508691A/ja active Pending
  • 2019-10-11 US US17/284,493 patent/US20210343020A1/en not_active Abandoned
  • 2019-10-11 EP EP19797453.8A patent/EP3864619A1/en not_active Withdrawn
  • 2019-10-11 WO PCT/US2019/055760 patent/WO2020077166A1/en not_active Ceased

Non-Patent Citations (3)

Also Published As

Similar Documents

Legal Events