CN115670497B - Digital oral panorama reconstruction method based on automatic focusing - Google Patents

Digital oral panorama reconstruction method based on automatic focusing Download PDF

Info

Publication number
CN115670497B
CN115670497B CN202211136767.8A CN202211136767A CN115670497B CN 115670497 B CN115670497 B CN 115670497B CN 202211136767 A CN202211136767 A CN 202211136767A CN 115670497 B CN115670497 B CN 115670497B
Authority
CN
China
Prior art keywords
focusing
coordinates
layer
new
point
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.)
Active
Application number
CN202211136767.8A
Other languages
Chinese (zh)
Other versions
CN115670497A (en
Inventor
张倩
张康平
孙宇
吴宏新
王亚杰
张文宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BEIJING LANGSHI INSTRUMENT CO LTD
Original Assignee
BEIJING LANGSHI INSTRUMENT CO LTD
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by BEIJING LANGSHI INSTRUMENT CO LTD filed Critical BEIJING LANGSHI INSTRUMENT CO LTD
Priority to CN202211136767.8A priority Critical patent/CN115670497B/en
Publication of CN115670497A publication Critical patent/CN115670497A/en
Application granted granted Critical
Publication of CN115670497B publication Critical patent/CN115670497B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)

Abstract

The invention relates to an automatic focusing-based digital oral panorama reconstruction method, which comprises the following steps: s1, scanning around a target object according to a scanning track to obtain multi-frame projection images; s2, according to the preset coordinate positions of the focusing surfaces, the positions of the focusing surfaces are adjusted by shifting along the beam direction, so that the coordinates of a plurality of groups of focusing surfaces are obtained; s3, reconstructing a multi-layer panoramic image according to the coordinates of a plurality of groups of focusing planes; s4, performing filtering operation on the reconstructed multi-layer panoramic image; step S5, calculating definition values of the panoramic images, and searching a definition layer of each part of the teeth based on the calculated definition values; s6, calculating offset of the coordinate of the preset focusing plane according to the layer sequence number of the clear layer where each point on each part of the teeth is located, and obtaining a new focusing plane coordinate suitable for the trend of the teeth of the patient; s7, correcting the coordinates of the new focusing surface; and S8, reconstructing based on the new focusing plane coordinates to obtain a clear image.

Description

Digital oral panorama reconstruction method based on automatic focusing
Technical Field
The invention relates to the technical field of medical images, in particular to a digital oral panorama reconstruction method based on automatic focusing.
Background
The oral full-view film, also called a curved surface broken layer film, is an X-ray photographic image obtained by rotating a narrow-strip X-ray detector along the outside of the cheek of a patient; the oral cavity full view film can display the images of the maxilla, the mandible, the temporomandibular joint, the maxillary sinus and the full mouth teeth on one image at the same time, and has a very important role in the diagnosis of dental diseases in the dental diagnosis.
The digital panoramic reconstruction is carried out by rotating around a patient, collecting a plurality of projection data under each angle, and then reconstructing an image at any given focusing point according to the calculation of the coordinate relation; the current method for obtaining the full-view film basically comprises the steps of determining a dental arch-shaped focusing surface in advance according to priori knowledge, and then reconstructing the focusing surface according to projection data to obtain a panoramic image. This approach has problems: the fixed focusing surface cannot be suitable for each patient, when the teeth of the user are uneven or the inclination angle of the teeth is large, the teeth are imaged by the existing method, and each tooth or each structure of the teeth cannot be focused at the same time, so that some teeth are unclear in the image, which is not beneficial to diagnosis of doctors.
Disclosure of Invention
In order to solve the technical problems, the invention provides a method for automatically judging the focusing position of each tissue in an imaging range and generating a clear image.
Firstly, determining an adjusting range of a focusing layer, then reconstructing a plurality of panoramic pictures, finding out focusing positions of each tissue structure in a mode of image definition judgment, and adjusting preset focusing plane coordinates; then, a new focusing surface suitable for the dental arch trend of a specific patient is regenerated according to the adjusted offset, and a final image is reconstructed.
The invention relates to an automatic focusing-based digital oral panorama reconstruction method, which comprises the following steps:
s1, scanning around a target object according to a scanning track to obtain multi-frame projection images;
s2, according to the preset coordinate positions of the focusing surfaces, the positions of the focusing surfaces are adjusted by shifting along the beam direction, so that the coordinates of a plurality of groups of focusing surfaces are obtained;
s3, reconstructing a multi-layer panoramic image according to the coordinates of a plurality of groups of focusing planes;
s4, performing filtering operation on the reconstructed multi-layer panoramic image;
step S5, calculating definition values of the panoramic images, and searching a definition layer of each part of the teeth based on the calculated definition values;
step S6, calculating the offset of the coordinate of the preset focusing plane according to the layer sequence number of the clear layer where each point on each part of the teeth is located, and adding the offset to the coordinate point corresponding to the preset focusing plane to obtain a new focusing plane coordinate suitable for the trend of the teeth of the patient;
s7, correcting the coordinates of the new focusing surface, eliminating some abnormal values caused by noise, and enabling the new focusing surface to be continuous and free of jump;
and S8, reconstructing based on the new focusing plane coordinates to obtain a clear image.
Further, in the step S2, according to a preset coordinate position of the focusing surface, the position of each focusing surface is adjusted by shifting along the normal direction of the scanning track, so as to obtain the coordinate of the multi-focusing surface, which specifically includes the following steps:
let P be any point on the focal plane, its three-dimensional coordinates be P (x, y, z), and when the position of each focal plane is adjusted according to the offset sigma, the new coordinates P '(x', y ', z') after adjustment are:
x′=x+σ x
y′=y+σ y
z′=z
wherein sigma x Sum sigma y The x, y direction offsets are shown, respectively, and the values are related to the set offset sigma and the current beam angle theta:
σ x =σ*cosθ
σ y =σ*sinθ
further, the step 5 of calculating the sharpness values of the panoramic images, and searching the sharpness layer where each tooth is located based on the calculated sharpness values specifically includes the following steps:
calculating each point f in multiple panoramic images i Gradient values of (x, y):
wherein i is more than or equal to 1 and less than or equal to N, i represents an ith Zhang Quanjing image, N is the number of pre-reconstructed panoramic images, and G is a corresponding point f i Gradient values of (x, y);
and searching a layer with the maximum gradient value in the plurality of panoramic images, namely a clear layer, recording layer sequence numbers and converting the offset corresponding to the layer.
Further, assume that the preset focal plane is shifted forward and backward along the normal line of the track for M times, and N images are reconstructed, where n=2xm+1, i ranges from [1,2 x m+1], and assuming that the difference of the shift amounts of the coordinates of two adjacent images is 1, σ ranges from [ -M, +m ]; and (3) obtaining the sigma corresponding to the point (x, y) when the gradient value is maximum based on the calculation.
Further, in the step 7, the coordinates of the new focusing surface are corrected, and abnormal values caused by noise are removed, so that the new focusing surface is continuous and has no jump;
let P be any point on the focal plane, its three-dimensional coordinates be P (x, y, z), the new coordinates after adjustment be P '(x', y ', z'), a region of m x n is selected around it, and the mean values of the x, y and z coordinates (x mean ,y mean ,z mean ) Comparing (x ', y ', z ') with the correction:
wherein threshold is a threshold.
The beneficial effects are that:
according to the digital oral panorama reconstruction method based on automatic focusing, imaging adjustment of a plurality of layers of focusing planes can be automatically carried out according to the offset, and the definition is automatically judged and the focusing planes are obtained based on a plurality of panoramic images obtained by the plurality of layers of focusing planes; and regenerating a new focusing surface suitable for the trend of the dental arch of the current patient, and reconstructing to obtain a final image, thereby ensuring that the imaging of each tooth is clear.
Drawings
Fig. 1: a method flow chart of the present invention;
fig. 2: defining a schematic diagram of a target and a scanning orbit coordinate system;
fig. 3: reconstructing an obtained image by using a preset focusing surface;
fig. 4: the obtained reconstructed image is automatically focused by the algorithm of the invention.
Detailed Description
The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are only some embodiments of the present invention, but not all embodiments, and all other embodiments obtained by those skilled in the art without the inventive effort based on the embodiments of the present invention are within the scope of protection of the present invention.
According to an embodiment of the present invention, an automatic focusing-based digital oral panorama reconstruction method is provided, as shown in fig. 1, comprising the following steps:
s1, scanning around a target object according to a scanning track to obtain multi-frame projection images;
s2, according to the preset coordinate positions of the focusing surfaces, the positions of the focusing surfaces are adjusted by shifting along the beam direction, so that the coordinates of a plurality of groups of focusing surfaces are obtained;
the definition of the direction is shown in fig. 2, the direction in which the face faces is in an X direction, the direction in which the head top extends upwards is in a Z direction, the direction in the left side of the face is in a Y direction, the Y direction is perpendicular to a plane formed by X, Z, and a X, Y, Z axis forms a vertical coordinate system;
let P be any point on the focal plane, its coordinates be denoted as P (x, y, z), then the new coordinates P '(x', y ', z') after adjustment are:
x′=x+σ x
y′=y+σ y
z′=z
wherein sigma x Sum sigma y The x, y direction offsets are shown, respectively, and the values are related to the set offset sigma and the current beam angle theta:
σ x =σ*cosθ
σ y =σ*sinθ。
s3, reconstructing a multi-layer panoramic image according to the coordinates of a plurality of groups of focusing planes;
s4, performing filtering operation on the reconstructed panoramic images to eliminate the influence of noise on follow-up images;
step S5, calculating definition values of the panoramic images, and searching a definition layer of each part of the teeth based on the calculated definition values;
according to the embodiment of the invention, the specific implementation mode is as follows:
calculating each point f in the image i The gradient values of (x, y), as in equation 1,
wherein i is more than or equal to 1 and less than or equal to N, i represents an ith Zhang Quanjing image, N is the number of pre-reconstructed panoramic images, and G is the gradient value of a corresponding point;
searching a layer with the maximum gradient value at the point in a plurality of groups of images, namely a clear layer, recording layer sequence numbers and converting the offset corresponding to the layer;
according to one embodiment of the present invention, for example, 15 shifts are respectively performed on the preset focus plane forward and backward, 31 images are reconstructed, where n=31, i ranges from [1,31], and if the difference of the shift amounts of the coordinates of two adjacent images is 1, σ ranges from [ -15,15]; if the gradient value is maximum when the calculated point (x, y) is i=20, the corresponding offset sigma=4 at the moment;
step S6, calculating the offset of the coordinate of the preset focusing plane according to the layer sequence number of the clear layer where each point on each part of the teeth is located, and adding the offset to the coordinate point corresponding to the preset focusing plane to obtain a new focusing plane coordinate suitable for the trend of the teeth of the patient;
and S7, reconstructing based on the new focusing plane coordinates to obtain a clear image.
Fig. 3 is a view of an image reconstructed with a preset focus plane, and fig. 4 is a view of a reconstructed image obtained by auto-focusing with the algorithm of the present invention. By using the digital oral panorama reconstruction method based on automatic focusing, provided by the invention, the imaging adjustment of the multi-layer focusing surface can be automatically carried out according to the offset without manual operation of a user, and the definition is automatically judged and the focusing surface is obtained based on a plurality of panoramic images obtained by the multi-layer focusing surface; and regenerating a new focusing surface suitable for the trend of the dental arch of the current patient, and reconstructing to obtain a final image, thereby ensuring that the imaging of each tooth is clear.
According to an embodiment of the present invention, step 7, correcting the coordinates of the new focal plane, and eliminating some abnormal values caused by noise, so that the new focal plane is continuous and has no jump;
let P be any point on the focal plane, its three-dimensional coordinates be P (x, y, z), the new coordinates after adjustment be P '(x', y ', z'), a region of m x n is selected around it, and the mean values of the x, y and z coordinates (x mean ,y mean ,z mean ) Comparing (x ', y ', z ') with the correction:
wherein threshold is a threshold.
While the foregoing has been described in relation to illustrative embodiments thereof, so as to facilitate the understanding of the present invention by those skilled in the art, it should be understood that the present invention is not limited to the scope of the embodiments, but is to be construed as limited to the spirit and scope of the invention as defined and defined by the appended claims, as long as various changes are apparent to those skilled in the art, all within the scope of which the invention is defined by the appended claims.

Claims (1)

1. An automatic focusing-based digital oral panorama reconstruction method is characterized by comprising the following steps:
s1, scanning around a target object according to a scanning track to obtain multi-frame projection images;
s2, according to the preset coordinate positions of the focusing surfaces, the positions of the focusing surfaces are adjusted by shifting along the beam direction, so that the coordinates of a plurality of groups of focusing surfaces are obtained; let P be any point on the focal plane, its three-dimensional coordinates be P (x, y, z), and when the position of each focal plane is adjusted according to the offset sigma, the new coordinates P '(x', y ', z') after adjustment are:
x′=x+σ x
y′=y+σ y
z′=z
wherein sigma x Sum sigma y The x, y direction offsets are shown, respectively, and the values are related to the set offset sigma and the current beam angle theta:
σ x =σ*cosθ
σ y =σ*sinθ;
s3, reconstructing a multi-layer panoramic image according to the coordinates of a plurality of groups of focusing planes;
s4, performing filtering operation on the reconstructed multi-layer panoramic image;
step S5, calculating definition values of the panoramic images, and searching a definition layer of each part of the teeth based on the calculated definition values;
step S6, calculating the offset of the coordinate of the preset focusing plane according to the layer sequence number of the clear layer where each point on each part of the teeth is located, and adding the offset to the coordinate point corresponding to the preset focusing plane to obtain a new focusing plane coordinate suitable for the trend of the teeth of the patient;
s7, correcting the coordinates of the new focusing surface, eliminating some abnormal values caused by noise, and enabling the new focusing surface to be continuous and free of jump; assuming that the preset focus plane is shifted forward and backward for M times along the normal line of the track respectively, reconstructing N images, wherein N=2×M+1, and assuming that the offset difference of coordinates of two adjacent images is 1, the range of sigma is [ -M, +M ]; based on the calculation, obtaining sigma corresponding to the point (x, y) when the gradient value is maximum;
step S5, performing sharpness calculation on a plurality of panoramic images, and searching a sharpness layer where each tooth is located based on the calculated sharpness values, specifically including the following steps:
calculating each point f in multiple panoramic images i Gradient values of (x, y):
wherein i is more than or equal to 1 and less than or equal to N, i represents an ith Zhang Quanjing image, N is the number of pre-reconstructed panoramic images, and G is a corresponding point f i Gradient values of (x, y);
searching a layer with the maximum gradient value in the point in the plurality of panoramic images, namely a clear layer, recording layer sequence numbers and converting the offset corresponding to the layer;
step S7, correcting the coordinates of the new focusing surface, eliminating abnormal values caused by noise, and enabling the new focusing surface to be continuous without jump; the method specifically comprises the following steps:
let P be any point on the focal plane, its three-dimensional coordinates be P (x, y, z), the new coordinates after adjustment be P '(x', y ', z'), a region of m x n is selected around it, and the mean values of the x, y and z coordinates (x mean ,y mean ,z mean ) Comparing (x ', y ', z ') with the correction:
wherein threshold is a threshold;
and S8, reconstructing based on the new focusing plane coordinates to obtain a clear image.
CN202211136767.8A 2022-09-19 2022-09-19 Digital oral panorama reconstruction method based on automatic focusing Active CN115670497B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211136767.8A CN115670497B (en) 2022-09-19 2022-09-19 Digital oral panorama reconstruction method based on automatic focusing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211136767.8A CN115670497B (en) 2022-09-19 2022-09-19 Digital oral panorama reconstruction method based on automatic focusing

Publications (2)

Publication Number Publication Date
CN115670497A CN115670497A (en) 2023-02-03
CN115670497B true CN115670497B (en) 2023-08-15

Family

ID=85062515

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211136767.8A Active CN115670497B (en) 2022-09-19 2022-09-19 Digital oral panorama reconstruction method based on automatic focusing

Country Status (1)

Country Link
CN (1) CN115670497B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117078662A (en) * 2023-10-11 2023-11-17 杭州睿影科技有限公司 Detection method and device for laminated battery, image processing equipment and medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105741239A (en) * 2014-12-11 2016-07-06 合肥美亚光电技术股份有限公司 Method and apparatus for generating tooth panorama image, and panorama machine for shooting teeth
CN208910300U (en) * 2017-11-28 2019-05-31 合肥美亚光电技术股份有限公司 Arch form curve debugging tool for oral cavity panoramic apparatus
CN114159090A (en) * 2021-12-07 2022-03-11 深圳市菲森科技有限公司 Dental panoramic shooting device and automatic correction method

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7336763B2 (en) * 2005-05-02 2008-02-26 Oy Ajat Ltd Dental extra-oral x-ray imaging system and method
US7676022B2 (en) * 2005-05-02 2010-03-09 Oy Ajat Ltd. Extra-oral digital panoramic dental x-ray imaging system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105741239A (en) * 2014-12-11 2016-07-06 合肥美亚光电技术股份有限公司 Method and apparatus for generating tooth panorama image, and panorama machine for shooting teeth
CN208910300U (en) * 2017-11-28 2019-05-31 合肥美亚光电技术股份有限公司 Arch form curve debugging tool for oral cavity panoramic apparatus
CN114159090A (en) * 2021-12-07 2022-03-11 深圳市菲森科技有限公司 Dental panoramic shooting device and automatic correction method

Also Published As

Publication number Publication date
CN115670497A (en) 2023-02-03

Similar Documents

Publication Publication Date Title
CN111223156B (en) Metal artifact eliminating method for dental cone beam CT system
US7715525B2 (en) Single sensor multi-functional dental extra-oral x-ray imaging system and method
JP5554248B2 (en) Creating tomographic images
EP3239923B1 (en) Method and apparatus for generating tooth panoramic image, and panoramic machine for photographing teeth
Du et al. A convolutional neural network based auto-positioning method for dental arch in rotational panoramic radiography
CN115670497B (en) Digital oral panorama reconstruction method based on automatic focusing
US20240065816A1 (en) Device and method for matching three-dimensional intra-oral scan data using computed tomography image
Gröndahl et al. Influence of variations in projection geometry on the detectability of periodontal bone lesions.
KR101502494B1 (en) X-ray photographic device, image processing method thereof, and program
JP2022518106A (en) Methods and equipment for generating 2D panoramic images
CN111833244B (en) Dental panoramic image generation method and device and computer readable storage medium
KR102302245B1 (en) Apparatus and Method for Eliminating Panoramic Dual and Ghost Image using Deep Learning
CN110619646B (en) Single tooth extraction method based on panorama
JP4527200B1 (en) Ear insertion support mechanism and cephalometric radiography apparatus
CN114972239A (en) Method and system for generating oral panoramic image based on oral CBCT (cone beam computed tomography) image
CN115908601A (en) Digital oral cavity panoramic reconstruction method based on manual interaction
CN117541672A (en) Method for automatically synthesizing panoramic image from oral cavity CBCT
JP2006239049A (en) X-ray ct apparatus
JP2005040506A (en) Method and apparatus for processing x-ray image
CN115736970A (en) Method for generating multilayer dental film image
CN114159090A (en) Dental panoramic shooting device and automatic correction method
KR102248191B1 (en) Method and apparatus for correcting cone-beam artifacts in a cone-beam computed tomography image, cone-beam computed tomography apparatus including the same
JP7187217B2 (en) Medical image processing device, X-ray diagnostic device and medical image processing program
Sohmura et al. A novel method of removing artifacts because of metallic dental restorations in 3‐D CT images of jaw bone
US20240338865A1 (en) Optimization of extraoral panoramic images through model-based prior knowledge of the patient's jaw arch form

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant