|Publication number||US7783093 B2|
|Application number||US 11/437,524|
|Publication date||Aug 24, 2010|
|Filing date||May 20, 2006|
|Priority date||May 20, 2005|
|Also published as||EP1899892A2, EP1899892A4, EP1899921A2, EP1899921A4, US7460638, US7545909, US7756247, US7817830, US7873192, US20060274881, US20060274926, US20060275740, US20070014461, US20090168954, WO2006127416A2, WO2006127416A3, WO2007021346A2, WO2007021346A3, WO2007030163A2, WO2007030163A3|
|Publication number||11437524, 437524, US 7783093 B2, US 7783093B2, US-B2-7783093, US7783093 B2, US7783093B2|
|Inventors||Arun Singh, Edward Marandola, Omid Kia, Uwe Mundry|
|Original Assignee||Imaging Sciences International Llc|
|Export Citation||BiBTeX, EndNote, RefMan|
|Patent Citations (26), Non-Patent Citations (8), Referenced by (6), Classifications (9), Legal Events (6)|
|External Links: USPTO, USPTO Assignment, Espacenet|
This application claims the benefit of U.S. Provisional Patent Application No. 60/682,971, filed May 20, 2005, which is imported herein by reference in its entirety.
This application is related to the following United States Patent Applications which are filed on even date herewith and which are incorporated herein by reference: application Ser. No. 11/437,523, filed on May 20, 2006, allowed on Jun. 24, 2010, entitled LOCATION OF ELONGATED OBJECT; application Ser. No. 11/437,525, filed on May 20, 2006, now U.S. Pat. No. 7,460,638, entitled LOCATION OF FOCAL PLANE; and application Ser. No. 11/437,526, filed on May 20, 2006, now U.S. Pat. No. 7,545,909, entitled GENERATING PANORAMIC VIEWS OF THE JAW USING A FIXED REFERENCE POINT AND CROSSING POINTS OF MULTIPLE RAYS.
The invention relates to locating a curved object in a three-dimensional array of data, and especially, but not exclusively, to locating an unevenly curved structure in a tomographic imaging dataset. The invention has particular application to locating the maxilla, the mandible, or both in a dataset of part of the head of a human or other mammal.
In certain forms of dental medicine and surgery, a “panoramic” image of the jaw is used to examine the jaw, for example, for monitoring of dental health and condition, diagnosis, and planning of prosthetic and other surgical procedures. The panoramic image of the jaw, like a panoramic photograph, depicts the jaw as if it were imaged onto an imaginary approximately cylindrical sheet with the axis of the sheet upright, and the sheet were then unrolled into a flat form. However, the human jaw is not a perfect circular arc, so the “cylindrical” shape of the imaginary sheet is not exactly circular.
A set of three-dimensional data relating to a property of an object that varies over space within the object may be obtained in various ways. For example, an x-ray image of a target may be obtained by placing the target between a source of x-rays and a detector of the x-rays. In a computed tomography (CT) system, a series of x-ray images of a target are taken with the direction from the source to the detector differently oriented relative to the target. From these images, a three-dimensional representation of the density of x-ray absorbing material in the target may be reconstructed. Other methods of generating a three-dimensional dataset are known, including magnetic resonance imaging, or may be developed hereafter.
From the three-dimensional data, a desired section or “slice” may be generated, including a curved slice. For example, a slice curving along the jaw, corresponding to a panoramic view of the jaw may be generated, provided that the position of the jaw within the three-dimensional dataset is known. It has previously been proposed to display on a visual display unit (VDU) a horizontal section through the mandible or maxilla, and for a human user to identify, for example by controlling a cursor on the VDU or by using a stylus on a touch-sensitive VDU, enough points on the mandible or maxilla for the curve of the jaw to be interpolated. However, many dentists and oral surgeons are not skilled computer operators.
There is therefore a hitherto unfulfilled need for a system by which the arc of the jaw can be accurately identified in a tomographic dataset without relying on the skill of a human operator, so that panoramic and similar images can be reliably automatically generated.
According to one embodiment of the invention, there is provided a method and system for locating in a three dimensional data array an arcuate object having an axial extent. Slices of data generally transverse to the axial extent of the object are selected. Rays generally radially of the arcuate object are selected within the slices. Crossing points where the rays cross the boundaries of the arcuate object are located. The position of the arcuate object is determined from the positions of the located points.
According to a preferred embodiment of the invention, the data array is a tomographic dataset of a human head or part thereof, and the curved object is the upper or lower jaw.
By locating the mandible or maxilla in the three-dimensional tomographic dataset, and exploiting the fact that the jaw is continuous along its length and has relatively sharp boundaries, the curve of the jaw can be identified with greater reliability than is attainable by most human operators without special skills or great effort.
The invention also provides computer software arranged to generate an image in accordance with the method of the invention, and computer-readable media containing such software. The software may be written to run on an otherwise conventional computer processing tomographic data.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory and are intended to provide further explanation of the invention as claimed.
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention.
In the drawings:
Reference will now be made in detail to various embodiments of the present invention, examples of which are illustrated in the accompanying drawings.
Referring to the drawings, and initially to
The imaging of the patient's head and calculation of the spatial distribution may be carried out by methods and apparatus already known in the art and, in the interests of conciseness, are not further described here. Suitable apparatus is available commercially, for example, the i-CAT Cone Beam 3-D Dental Imaging System from Imaging Sciences International of Hatfield, Pa.
In step 106, an initial slice 202 through the tomographic dataset is selected. As shown in
In step 108, a reference point 204 is selected. The reference point 204 is preferably on, or close to, the centerplane of the head, and approximately at the center of the dental arch. Because the dental arch is not an arc of a circle, and may not be symmetrical from side to side, very precise centering of the reference point 204 is not required, and may not be meaningfully possible. The same reference point 204 will be used in subsequent processing of other slices. Therefore, the reference point 204 may be treated as the intersection of a vertical reference line with the slice 202, and either the reference line/point 204 or the slice 202 may be selected first, or they may be selected in parallel. For subsequent slices 202, the reference point 204 is automatically selected as a point on the reference line through the reference point in the first slice.
In step 110, a probability count is set to zero.
In step 112, a ray 206 in the slice 202 extending radially away from the center point 204 through the jaw 200 is selected.
In step 114, the ray 206 is inspected to identify the points 208, 210 where the ray crosses the inner face 212 and the outer face 214 of the jaw 200. These faces 212, 214 are typically recognized as ramps in the density profile along the ray, where the density rises rapidly towards the jaw 200. The bone of the jaws typically has a markedly higher density than the surrounding soft tissue. The crossing points 208, 210 may be taken as points where the ramp crosses a selected density threshold. If the ramps are not clearly defined, for example, if the steepness or height of either or both ramp is less than a preset threshold, the ray 206 may be discarded.
If crossing points 208, 210 are identified, in step 116 the crossing points may be subjected to a quality test. For example, the ray may be discarded if the identified crossing points 208, 210 are not approximately where the jaw 200 is expected to be, are too close together or too far apart, or are out of line with the crossing points 208, 210 in neighboring rays.
If the ray passes the quality test, in step 117 the probability count is increased by one, and the position where the ray 206 crosses the jaw is recorded. In one embodiment, only the center of the jaw, which may be taken as a point 216 midway between the crossing points 208, 210, is recorded. In another embodiment, the crossing points 208, 210 are recorded instead, or in addition. If the ray 206 is discarded either at step 114 because adequate ramps are not found or at step 116 because the quality test is failed, then step 117 is bypassed.
The process then returns to step 112 to select another ray 206 in the same slice. A selected number of rays 206 evenly spaced around the arc of the jaw are selected. When it is determined in step 118 that all of the rays 206 for the slice 202 have been processed, the probability count for the slice 202 is recorded, and the process returns to step 106 to select another slice 202.
The process then repeats for each of a stack of slices 202 over a part of the tomographic dataset where, based on the normal positioning of the headrest 36 relative to the x-ray source 28 and detector 30, the jaw 200 is likely to be present. It is not usually necessary to process every voxel slice of the original dataset. Preferably, a reasonable number of slices, which may be evenly spaced or may be more closely spaced in areas of expected high probability, are processed. The initial slice 202 may be the top or bottom slice of the stack of slices, with the process working systematically down or up the stack. However, that is not necessary, and the initial slice 202 may be an arbitrary slice.
Various methods may be used to define the upper and lower bounds of the stack of slices 202 used in steps 106 to 120. A bound may be set as the upper and/or lower edge of the field of the original scan of part of the patient's head. A bound may be set manually using a vertical display of the dataset. The upper or lower edge of the mandible, or the lower edge of the maxilla, may be easily found by detecting the point where high density tissue (usually teeth or bone) abruptly disappears. Other suitable methods may also be used.
When it is determined in step 120 that all of the slices 202 have been processed, in step 122 the recorded probability counts for the different slices 202 are inspected. As shown by the curve 220 in
In step 124, the slice 202 with the highest probability count is selected, and the center points 216 of the selected slice are selected, if the center points were calculated and recorded in step 117, or are calculated, if only the crossing points 108, 110 were recorded in step 117. If the peak of the probability curve 220 has a “plateau,” where more than one neighboring slice has the identical highest probability count, any one of those slices can be chosen arbitrarily. All of those slices are likely to lie in a “good” region of the jaw. A contour line 222 representing the dental arch is then constructed through the points 216. Where a bimodal distribution is found, separate contour lines 222 may be constructed for the maxilla and the mandible.
In step 126, a synthesized panoramic slice formed from columns of voxels along the contour line 222 is then generated and presented to the user. The columns of voxels may be perpendicular to the plane of the slices 202, or may be at an angle. Where separate mandibular and maxillary contour lines 222 have been constructed, the columns of voxels may be slanted or curved to pass through both contour lines 222, giving a hybrid panoramic view that shows both the maxilla and the mandible. The synthesized panoramic slice may be one or more voxels thick. For example, the thickness of the panoramic slice may correspond to a typical value for the actual thickness of the mandible or maxilla perpendicular to the panoramic slice plane.
The thickness of the human mandible and maxilla vary, both from person to person and from point to point within the jaw. Consequently, the synthesized panoramic slice displayed in step 126 may be less than optimal because it may include too much or too little of the thickness of the jaw. If the panoramic slice is too thick, bony structure may overlay or “shade” detail in the interior of the mandible or maxilla, such as the alveolar nerve, or blur details of the root of a tooth. If the panoramic slice is too thin, bony outgrowths, protruding calcification, or displaced teeth may be missed.
Referring now to
Alternatively, if it is desired to view the internal structure of the mandible or maxilla, part of the cortical bone may be virtually “machined” away to generate a slice in which less bone is present to obscure the view. In step 306, the process may select, or get from a user, a figure for the amount of machining, for example, a “percentage of erosion.” The panoramic slice is then generated in step 308, omitting a corresponding part of the cortical bone. For example, where the panoramic slice in step 304 extends from −50% to +50% of the thickness of the jaw, centered on the contour 222 at 0%, a panoramic slice with 30% erosion may extend from −35% to +35%. Other patterns of erosion or machining may be used, for example, removal of a constant thickness of bone.
In step 310, the tomographic data generated as described above with reference to steps 102, 104 in
Where the segmentation process is carried out using a surface as the contour, the final result may be equivalent to the envelope generated in step 302, and may be passed directly to steps 304 and 306. Where the segmentation process is carried out on individual slices 202, the slice contours may be stacked and an envelope interpolated similarly to step 202. Where both an envelope surface from step 302 and a segmentation surface from step 310 are available, the two surfaces may be compared in step 312 to corroborate the identification of the actual jaw surface. Alternatively, portions of the jaw may be located by the process of steps 106 to 120, and portions by segmentation, and combined in step 312 to give a full surface. For example, a three dimensional segmentation process may be better suited to parts of the head where the bone surface is not close to vertical relative to the axis 204.
Also, the segmentation process can identify and distinguish separate bony structures in the same dataset. In
Various modifications and variations can be made in the present invention without departing from the spirit or scope of the invention. For example, the order of the steps may be changed. In
In a further alternative, in step 124 two or more slices 202 with high probability counts may be combined. An average of the positions of the contour lines 222, which may be weighted according to the probability counts of the respective slices 202, may be used, or a consensus curve may be determined by combining only contour lines 222 that agree within a narrow tolerance.
As described with reference to
Where steps 114 and 116 find only one clear crossing point 208 or 210, step 302 may use that one point, although in step 124, a single crossing point 208 or 210 may be unhelpful for constructing the center point 216.
For example, in the processes described above, it is assumed that the number of rays is equal in all slices, and the rays are uniformly spaced. The number of rays may vary from slice to slice, provided that an appropriate allowance is made at steps 122 and 124. The rays may be more closely spaced at parts of the jaw where a small radius of curvature and/or a rapid change in curvature is expected. In the process of
For example, it is not necessary for the reference points 204 of the different slices 202 all to lie on a straight line. However, using consistent reference points 204 typically gives more consistent, and thus more comparable results, and may give a more accurate decision in step 124. In step 302, using consistent reference points 204 gives a more consistent relationship between the arrays of crossing points 208 and 210 and may make the generation of the envelopes computationally simpler. In addition, using a vertical line of points 204 is usually the simplest approach.
Although distinct embodiments have been described, features of different embodiments may be combined. For example, the contour 222 shown in
Thus, it is intended that the present invention cover modifications and variations of this invention provided they come within the scope of the appended claims and their equivalents.
|Cited Patent||Filing date||Publication date||Applicant||Title|
|US4661967||Nov 2, 1984||Apr 28, 1987||Kabushiki Kaisha Morita Seisakusho||Dental radiographic apparatus for photographing entire jaws|
|US4856038||Sep 28, 1987||Aug 8, 1989||Siemens Aktiengesellschaft||Dental x-ray diagnostics installation for producing panoramic exposures of slices in the skull|
|US5195114||Oct 1, 1991||Mar 16, 1993||Kabushiki Kaisha Morita Seisakusho||Digital panoramic radiographic apparatus|
|US5214686||Dec 13, 1991||May 25, 1993||Wake Forest University||Three-dimensional panoramic dental radiography method and apparatus which avoids the subject's spine|
|US5339367||Dec 2, 1992||Aug 16, 1994||National Research Council Of Canada||Identifying curves within a scanned image|
|US5511106||Jun 30, 1994||Apr 23, 1996||Siemens Aktiengesellschaft||X-ray diagnostics installation for producing x-ray exposures of body parts of a patient|
|US5677940||Mar 20, 1996||Oct 14, 1997||J. Morita Manufacturing Corporation||Digital X-ray imaging apparatus|
|US5784429||Jul 2, 1997||Jul 21, 1998||J. Morita Manufacturing Corporation||Dental panoramic X-ray imaging apparatus|
|US6021222||Jan 8, 1998||Feb 1, 2000||Ricoh Co., Ltd.||System and method for the detection of a circle image for pattern recognition|
|US6049584||Aug 3, 1998||Apr 11, 2000||Sirona Dental Systems Gmbh & Co. Kg||X-ray diagnostic apparatus for producing panorama slice exposure of body parts of a patient|
|US6289074||Sep 1, 1999||Sep 11, 2001||J. Morita Manufacturing Corporation||X-ray computed tomography method and system|
|US6493415||Mar 27, 2000||Dec 10, 2002||Nihon University||X-ray computed tomography method and apparatus|
|US6570953||Mar 22, 2000||May 27, 2003||Airona Dental Systems Gmbh||Method for making and reproducing a tomogram of an object, said tomogram pertaining to a section roentgenogram|
|US6780155||Dec 18, 2001||Aug 24, 2004||Koninklijke Philips Electronics||Method and system for ultrasound blood flow imaging and volume flow calculations|
|US6798860||May 17, 2000||Sep 28, 2004||Ge Medical Systems Global Technology Company, Llc||Methods and apparatus for deconvolving imaging data|
|US6859555||Sep 19, 2000||Feb 22, 2005||Siemens Corporate Research, Inc.||Fast dominant circle detection through horizontal and vertical scanning|
|US7039156||Oct 4, 2001||May 2, 2006||Nihon University||Display method and apparatus of x-ray projection image for medical use, x-ray CT apparatus for medical use and recording medium for recording program to achieve the display method|
|US7369691||Mar 19, 2004||May 6, 2008||Kabushiki Kaisha Toshiba||Processor for analyzing tubular structure such as blood vessels|
|US7421059||Apr 10, 2003||Sep 2, 2008||J. Morita Manufacturing Corporation||X-ray computer tomography apparatus|
|US7460638||May 20, 2006||Dec 2, 2008||Imaging Sciences International Llc||Location of focal plane|
|US7499576||Jun 25, 2004||Mar 3, 2009||Palodex Group Oy||Method and system for detecting a fiducial in digital projection images|
|US20040066877 *||Oct 4, 2001||Apr 8, 2004||Yoshinori Arai||Medical x-ray ct image display method, display device, medical x-ray ct device and reocrding medium recording program implementing this display method|
|US20040236210 *||May 22, 2003||Nov 25, 2004||Ming-Dar Tsai||Image analysis method for vertebral compression curvature|
|JP3621146B2||Title not available|
|JP2001061834A||Title not available|
|JPH04144548A||Title not available|
|1||Office Action from U.S. Appl. No. 11/437,523 dated Dec. 29, 2009 (10 pages).|
|2||Office Action from U.S. Appl. No. 11/437,525 dated Sep. 27, 2007 (9 pages).|
|3||Office Action from U.S. Appl. No. 11/437,526 dated Feb. 26, 2008 (8 pages).|
|4||Office Action from U.S. Appl. No. 11/437,526 dated Sep. 22, 2008 (8 pages).|
|5||PCT/US06/19319 International Search Report and Written Opinion dated Nov. 8, 2007 (9 pages).|
|6||PCT/US2006/019320 International Preliminary Report on Patentability dated Mar. 6, 2008 (2 pages).|
|7||PCT/US6/19320 International Search Report and Written Opinion dated Feb. 11, 2008 (10 pages).|
|8||U.S. Appl. No. 11/437,525 Third-Party Submission of Prior Art in Published Application mailed Oct. 3, 2007 (38 pages).|
|Citing Patent||Filing date||Publication date||Applicant||Title|
|US8059873 *||Nov 16, 2005||Nov 15, 2011||Koninklijke Philips Electronics N V Eindhoven||Volume of interest selection|
|US8325874||Jul 30, 2010||Dec 4, 2012||Imaging Sciences International Llc||Panoramic dental imaging using segmentation and a master arch|
|US8548120||Dec 4, 2012||Oct 1, 2013||Dental Imaging Technologies Corporation||Dental imaging using segmentation and an arch|
|US9020192||Mar 13, 2013||Apr 28, 2015||Access Business Group International Llc||Human submental profile measurement|
|US20080123914 *||Nov 16, 2005||May 29, 2008||Koninklijke Philips Electronics, N.V.||Volume of Interest Selection|
|US20110026671 *||Jul 30, 2010||Feb 3, 2011||Hongjian Shi||Panoramic dental imaging using segmentation and a master arch|
|Cooperative Classification||G06T11/003, G06T7/73, A61B6/5258, A61B6/14, G06T2207/30008|
|European Classification||G06T7/00P1, A61B6/14|
|Aug 11, 2006||AS||Assignment|
Owner name: IMAGING SCIENCES INTERNATIONAL, INC., PENNSYLVANIA
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:SINGH, ARUN;MARANDOLA, EDWARD;KIA, OMID;AND OTHERS;REEL/FRAME:018145/0149
Effective date: 20060731
|Jan 19, 2007||AS||Assignment|
Owner name: IMAGING SCIENCES LLC, PENNSYLVANIA
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:IMAGING SCIENCES INTERNATIONAL, INC.;REEL/FRAME:018780/0801
Effective date: 20070102
|Dec 6, 2007||AS||Assignment|
Owner name: IMAGING SCIENCES INTERNATIONAL LLC, PENNSYLVANIA
Free format text: CHANGE OF NAME;ASSIGNOR:IMAGING SCIENCES LLC;REEL/FRAME:020216/0429
Effective date: 20070418
|Apr 4, 2014||REMI||Maintenance fee reminder mailed|
|Aug 24, 2014||LAPS||Lapse for failure to pay maintenance fees|
|Oct 14, 2014||FP||Expired due to failure to pay maintenance fee|
Effective date: 20140824