CN112700487B - Method and system for acquiring measuring scale marks in skull lateral position slice - Google Patents
Method and system for acquiring measuring scale marks in skull lateral position slice Download PDFInfo
- Publication number
- CN112700487B CN112700487B CN202011642265.3A CN202011642265A CN112700487B CN 112700487 B CN112700487 B CN 112700487B CN 202011642265 A CN202011642265 A CN 202011642265A CN 112700487 B CN112700487 B CN 112700487B
- Authority
- CN
- China
- Prior art keywords
- scale
- skull
- measurement
- side position
- acquiring
- 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
Links
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The application provides a detection method and a detection system for a measuring scale in a skull lateral position slice, comprising the following steps: determining a picture of a skull measurement scale in the skull side position slice according to the skull mark points; performing digital processing on the picture of the skull measurement scale according to a first preset algorithm, and determining the relative position information of the skull measurement scale in the skull side position slice; performing image processing on the skull measurement staff gauge graph according to a second preset algorithm to obtain a scale clearance of the skull measurement staff gauge and a scale phase of the staff gauge in the skull side position slice; according to the scale gap of the skull measurement scale and the scale phase of the scale, and combining the length average value of the human upper jaw to calculate the scale of the skull measurement scale. The application solves the problems that the staff gauge distance of the skull side position slice needs to be manually participated, the distance data has errors and the efficiency is lower in the prior art.
Description
Technical Field
The application belongs to the technical field of tooth correction, and particularly relates to a method and a system for acquiring a measuring scale in a skull lateral position slice, electronic equipment and a computer storage medium.
Background
The skull side position slice is one of the criteria used by doctors to judge whether the patient is odontopathy or osseous convexity, and uses standard measuring software to measure related data and compare the data with the standard value, so that the diagnosis and treatment design of the odontopathy are more accurate. In the process of calculating the relevant measured value, besides the mark points of the skull side position slice, the scale on the skull side position slice needs to be detected so as to accurately and quantitatively give the measured data.
At present, the scale of most skull side position plates needs to manually point out two points on the scale and give the distance represented by the two points on the scale, so that the corresponding actual distance between pixels in the image can be judged. The judging mode has a very large short plate in the automatic analysis of the skull side position slice, needs manual participation and has lower efficiency.
Based on the above, the present application provides a technical solution to the above technical problems.
Disclosure of Invention
The application aims to overcome the defects in the prior art, provides a measuring scale acquisition system, electronic equipment and a computer storage medium in a skull side position slice, and solves the problems that the scale distance of the skull side position slice in the prior art needs to be manually participated, the data has errors and the efficiency is lower.
The technical scheme provided by the application is as follows:
a method of acquiring a scale of a measurement scale in a cranial flap, comprising:
determining a picture of a skull measurement scale in the skull side position slice according to the skull mark points; performing digital processing on the picture of the skull measurement scale according to a first preset algorithm, and determining the relative position information of the skull measurement scale in the skull side position slice; performing image processing on the skull measurement staff gauge graph according to a second preset algorithm to obtain a scale clearance of the skull measurement staff gauge and a scale phase of the staff gauge in the skull side position slice; according to the scale gap of the skull measurement scale and the scale phase of the scale, and combining the length average value of the human upper jaw to calculate the scale value of the skull measurement scale.
Further preferably, determining the relative position information of the skull measurement scale in the skull side position slice comprises: dividing the picture of the skull measurement scale into a plurality of sub-blocks according to a preset dividing algorithm;
and carrying out Fourier transform processing on each sub-image block, carrying out marking arrangement on the sub-image blocks with the same frequency in the transformation information of each sub-image block after the Fourier transform processing, and extracting the relative position information of the skull measuring scale in the skull side position piece.
Further preferably, dividing into sub-blocks comprises:
L max ＝L K ＝min(N 1 ,N 2 )/N B
L min ＝L 1 ＝min(h,L max /K),
L n ＝(L max /L min ) (n-1)/(K-1) L min )；N 1 ,N 2 the length of the length and width of the picture (i.e. pixels) representing the skull measurement scale, respectivelyLength of (d) a); n (N) B Representing the number of preset divided sub-blocks on the side of minimum length, L n The size of the nth sub-picture block is obtained through logarithmic size; h is the minimum block size of the sub-picture block; k refers to the number of L values to be given, and the length L of the sub-block n-1 With a value of L or less n The method comprises the steps of carrying out a first treatment on the surface of the The smallest part of interest in image processing can be assumed to be a square of h x h.
Further preferably, determining the relative position information of the skull measurement scale in the skull side position slice comprises: setting a picture of a skull measurement scale in a preset traversal area, and acquiring pose information of the preset traversal area in a skull side position slice; the pose information includes: inclination angle and relative position information of the skull measurement scale in the skull patch.
Further preferably, the method comprises: generating a traversing line corresponding to each angle value in a preset angle range,
acquiring pixel values of a skull measurement scale map on the generated traversing line; carrying out Fourier transform on the pixel values of the obtained skull measurement scale image to obtain Fourier transform images of different positions under different angles; and acquiring traversing line information corresponding to the highest amplitude point in the generated Fourier transform graph, and extracting the pose information of the skull measurement scale graph from the traversing line information.
Further preferably, the method further comprises: and acquiring linear position information in the skull side position slice through a Hough transformation algorithm, and setting the acquired linear position information as the relative position information of the skull measuring scale in the skull side position slice.
Further preferably, acquiring the scale gap of the scale and the scale phase of the scale includes:
acquiring the inclination angle of the skull measurement scale in the skull patch by using a power density function;
acquiring a straight line perpendicular to the scale mark of the ruler according to the inclination angle, and extracting a pixel value of the straight line;
acquiring the fundamental frequency of the relative position information of the skull measurement scale by accumulating the average normalized square error function through the pixel values of the pixels on the extracted straight line, and acquiring the scale clearance value of the skull measurement scale through the fundamental frequency;
and performing least square fitting on the skull measurement staff drawings to obtain a sine model of the staff, and obtaining the scale phase of the staff through the sine model.
Further preferably, the inclination angle of the measurement scale includes:
the coordinates of the point where the power spectral density is the greatest are (v, u).
Further preferably, the power density of each sub-tile comprises: PSDF (ω) =f { Y B }F * {Y B Z, F is Y B Is a frequency spectrum of (c).
Further preferably, obtaining the fundamental frequency value includes:
autocorrelation ACF function: is the similarity between the function value y (j) and the function value y (j+T) after the self translation T; t is the y translation number.
Further preferably, the method comprises: a sine model of the graduated scale is obtained by least squares fitting,
(Y B (x)-M(x,α)) 2 ，
m- -sine model F- -Y B The frequency spectrum, 2pi/lambda is atEach pixel period in the ruler direction; lambda is the number of pixels between the scale marks. Alpha 1 ,α 2 Respectively representing the magnitudes of sin and cos functions, α= { α 1 ,α 2 The } is obtained by iteration.
A measuring scale acquiring system in a skull side position slice, which can execute the method for acquiring the measuring scale in the skull measurement, comprising a measuring scale image acquiring module, a measuring scale image acquiring module and a measuring scale image acquiring module, wherein the measuring scale image acquiring module is used for acquiring the image of the skull measuring scale according to the skull mark points in the skull side position slice;
the scale image position acquisition module is used for carrying out digital processing on the picture of the skull measurement scale acquired by the measurement scale image acquisition module according to a first preset algorithm, and determining the relative position information of the skull measurement scale in the skull side position piece;
the scale image processing module is used for carrying out image processing on the position information of the skull measurement scale acquired by the scale image position acquisition module according to a second preset algorithm to acquire a scale gap of the skull measurement scale and a scale phase of the scale in the skull side position piece;
the scale acquisition module is used for calculating the scale of the skull measurement scale according to the scale gap of the skull measurement scale and the scale phase of the scale, which are acquired in the scale image processing module, and combining the length average value of the upper jaw of the human body.
An electronic device comprising a processor and a memory, the processor executing computer instructions stored by the memory, causing the electronic device to perform the method of obtaining a measurement scale in a cranial-lateral plate of any one of the preceding claims.
A computer storage medium comprising computer instructions which, when run on an electronic device, cause the electronic device to perform the method of measuring scale acquisition in a skull side-slice as claimed in any preceding claim.
The application provides a method and a system for acquiring a measuring scale in a skull lateral position slice, electronic equipment and a computer storage medium, which can bring at least one of the following beneficial effects:
according to the application, the position of the staff gauge is roughly positioned in the skull side position plate, the accurate position of the staff gauge is determined, the scale clearance of the staff gauge and the scale phase of the staff gauge are further calculated in an iterative mode, the staff gauge is realized through a computer language, and the staff gauge is displayed in a data form, so that the problem of low manual participation efficiency in the prior art is solved.
According to the application, the staff scale is obtained by taking the average value of the upper jaw length of a human body in different tooth replacement stages in a standard database and taking the average value as a reference to finish the staff scale obtaining in various modes such as Fourier transform, hough transform and the like, so that the automatic obtaining is realized, manual participation is liberated, the data is more refined, an accurate data basis is further provided for invisible orthodontic, and meanwhile, the working efficiency is improved; and particularly, accurate data information is given for the difficult and complicated cases, so that a perfect experience is provided for the patients in need.
Drawings
The above features, technical features, advantages and implementation thereof will be further described in the following detailed description of preferred embodiments with reference to the accompanying drawings in a clearly understandable manner.
FIG. 1 is a flow chart of an embodiment of a method for obtaining a measurement scale in a cranial flap of the present application;
FIG. 2 is a scale drawing of a measurement scale in a cranial flap of the present application;
FIG. 3 is another flowchart of an embodiment of a method for obtaining a measurement scale in a cranial lateral plate according to the present application;
FIG. 4 is another flow chart of an embodiment of a method for obtaining a measurement scale in a cranial lateral plate according to the present application;
FIG. 5 is a view showing the position of the maxillary length of the cranial flap of the present application;
FIG. 6 is a segmented sub-segment of the inventive skull side slice after measurement scale map image processing;
FIG. 7 is a graph of transformed data of a measurement scale map in a cranial lateral plate of the present application;
FIG. 8 is a scale map of the measuring scale map image processed in the cranial lateral plate of the present application;
FIG. 9 is a block diagram of an embodiment of a measurement scale acquisition system in a cranial flap of the present application;
FIG. 10 is a graph traversal of the measurement scale algorithm in the cranial flap of the present application;
fig. 11 is a structural view of the electronic device of the present application.
Detailed Description
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the following description will explain the specific embodiments of the present application with reference to the accompanying drawings. It is evident that the drawings in the following description are only examples of the application, from which other drawings and other embodiments can be obtained by a person skilled in the art without inventive effort.
The method comprises the steps of obtaining distance information between mark points of a skull lateral plate, comparing the distance information with a standard value, judging the oral state of a correction patient through definite numerical information, facilitating a doctor to formulate a correction scheme of the patient, specifically obtaining corresponding scale distances according to pixel values and scale distances by obtaining scale scales in the lateral plate, wherein the scale distances are used for analyzing the pixel/scale values to determine the distance between two points in the medial plate; in the prior art, the measurement is manually operated or finished by virtue of the experience of doctors, so that the digital age of the existing large environment is not satisfied, and the automatic diagnosis and treatment cannot be satisfied. Based on this, the present application provides a technical solution for solving the technical problems, referring to the following examples.
Referring to fig. 1 to 11, the present application provides an embodiment of a method for detecting a measuring scale in a skull side-position plate, comprising the steps of:
step S100, determining a picture of a skull measurement scale in a skull lateral position slice according to a skull mark point;
according to the rule of the head side position plate shooting, the position of the scale in the figure is referred to the position of the upper right corner shown in fig. 2, the position of the scale is generally in front of the nasion point (N) when the head side position plate shooting is carried out, and the position of the scale of the side position plate facing to the right is referred to the left side of the nasion point in combination with fig. 5, namely the position of the scale of the side position plate facing to the left is referred to the left side of the nasion point, the picture with the scale is cut according to the coordinates of the nasion point, and the picture with the scale can be obtained based on the world coordinate system, and the picture with the corresponding scale is obtained based on the head side position plate.
In the application, the nasion point is used as one of the cranium marking points, which can be obtained by automatic searching by a method which is not limited to random forest, and a plurality of marking points are arranged in the cranium. The selection of the mark points is not limited in the application.
Step S200, carrying out digital processing on the picture of the skull measurement scale according to a first preset algorithm, and determining the relative position information of the skull measurement scale in the skull side position slice;
referring to fig. 6, the acquired skull measurement scale is digitized, that is, converted into digital information, and the positional relationship thereof is located, and relevant scale information is obtained according to the positional relationship.
The following embodiments are provided for the first preset algorithm to obtain the relative position information of the skull measurement scale in the skull lateral position slice:
referring to fig. 3 and 6, a first embodiment is as follows:
step S211, dividing the picture of the skull measurement scale into a plurality of sub-blocks according to a preset dividing algorithm;
step S212, carrying out Fourier transform processing on each sub-image block;
step S213, marking and sorting sub-blocks with the same frequency in the sub-blocks after the Fourier transform processing;
step S214 extracts the relative position information of the skull measurement scale in the skull side position slice.
Specifically, in order to accurately position the scale, the obtained picture of the skull measurement scale can be divided into different blocks, the size of the blocks is L, fourier transform is performed on each small block, the small blocks with the same frequency characteristic are marked, and the small block area is the scale position needing accurate positioning. As shown in fig. 6 below, a block with the same fourier characteristics is circled by a rectangular block, which is a fourier calculation result graph after a certain rough localization map is divided into a block, and the relative position of the skull measurement scale is obtained after the same-frequency block is sorted.
The specific implementation mode of dividing the sub-blocks into a plurality of sub-blocks according to a preset dividing algorithm is as follows:
different block sizes L may estimate different scale spacing distances d (L). For the input image, { L 1 ,L 2 ,...L K Continuing to automatically select blocks containing part of the scale image, K refers to the number of L values that need to be given. Maximum square edge L (L max ) And the smallest edge (L min ) Are set to be the fraction times of the smallest edge of the input image;
L max ＝L K ＝min(N 1 ,N 2 )/N B
L min ＝L 1 ＝min(h,L max /K),
L n ＝(L max /L min ) (n -1)/(K-1)L min )；
N 1 ,N 2 the length of the length and width of the picture (i.e. the length of the pixel) representing the skull measurement scale, respectively; n (N) B Representing the number of preset divided sub-blocks on the side of minimum length, L n The size of the nth sub-picture block is obtained through logarithmic size; h is the minimum block size of the sub-picture block; k refers to the number of L values to be given, and the length L of the sub-block n-1 With a value of L or less n The method comprises the steps of carrying out a first treatment on the surface of the The smallest part of interest in image processing can be assumed to be a square of h x h.
From the smallest L 1 Initially, let i' th L, i.e. L i Then the image is divided into (N) 1 /L i )*(N 2 /L i ) A block for processing the signal of each small block until at least one B is found min Successive dice have the same frequency characteristics, B min For the minimum number of consecutive blocks of the same frequency, the selection is stopped, i.e. the selected block length is L i . If the ith does not satisfy the conditional feature, then the (i+1) th L is selected. Through numerical experiments, it is found that B min Is the minimumIs sufficient to eliminate the number of most binary relations (candidate points simply selected due to noise signals). Naturally, the larger will reduce the false positive data (background environment). At the same time, however, if B min Too large, the number of false negatives will likewise be reduced.
The second embodiment is as follows:
step S221, setting the picture of the skull measurement scale in a preset traversal area;
step S222, acquiring pose information of a preset traversal region in a skull side position slice; the pose information includes: the angle of inclination and position information (i.e., relative to position information) of the skull measurement scale in the skull flap.
Referring to fig. 10, the traversal process is: generating a traversing line corresponding to each angle value in a preset angle range, and acquiring pixel values of a skull measurement scale map on the generated traversing line; carrying out Fourier transform on the pixel values of the obtained skull measurement scale map to obtain Fourier transform scale maps with different amplitudes; and acquiring traversing line information corresponding to the highest point of the amplitude in the generated Fourier transform scale graphs with different amplitudes, and extracting pose information of the skull measurement scale graph from the traversing line information.
The position and the angle of the scale are searched by using a traversing method by utilizing the periodicity rule of the scale marks. In fig. 10, θ is the included angle between the traversing straight line and the vertical direction, the value range of θ is assumed to be [ -10 °,10 ° ] according to the rule of thumb, each traversing line takes the corresponding pixel value, and performs a fourier transform, so that a fourier transform graph is generated under each θ, as shown in fig. 7 below, and the traversing line corresponding to the point with the highest amplitude in fig. 7 can be the line where the scale mark is located. And extracting the relative position information of the graduated scale from the function of the line where the graduated scale is located.
The third embodiment mode is as follows:
referring to fig. 8, the position information of the straight line is obtained from the skull side position slice by the Hough transformation algorithm, and the obtained position information of the straight line is set as the relative position information of the skull measuring scale in the skull side position slice.
Step S300, performing image processing on the skull measurement staff gauge graph according to a second preset algorithm, and acquiring a scale gap of the skull measurement staff gauge and a scale phase of the staff gauge in the skull side position slice;
specifically, in the accurate positioning of the relative position information of the scale, the characteristic that the autocorrelation function is affected by noise relatively low can be utilized to calculate the cumulative average normalized square difference function (CSDF Cumulative Mean Normalised Difference Functions) to calculate the fundamental frequency, and after the fundamental frequency is calculated, the sine model of the scale is obtained through least square fitting. The scale gap of the scale can be obtained through the fundamental frequency.
The sine model of the scale is obtained by least square fitting, which comprises the following steps:
(Y B (x)-M(x；α)) 2 ，--(1)，Y B (x) A rule grid function for cutting;
m- -sine model F- -is the spectrum of YB, 2pi/λ is each pixel period in the ruler direction; lambda is the image pixel. Is obtained by iteration, alpha 1 α 2 Respectively representing the magnitudes of sin and cos functions, α= { α 1 ,α 2 }。
Step S310, acquiring the inclination angle of the skull measurement scale in the skull patch by using a power density function; in the first position acquisition mode of the present application, the scale may be divided into small squares, the power density (spectrum of the autocorrelation function) is obtained for the small squares, and the coordinates of the point with the maximum power spectral density are obtained in the power spectral density map, and the direction of the scale:
measuring the inclination angle of the scale includes:the coordinates of the point with the maximum power spectral density are (v, u);
power density of each sub-tile: PSDF (ω) =f { Y B }F * {Y B }--(4)
Step S320, obtaining a straight line perpendicular to the scale mark of the ruler according to the inclination angle direction, and extracting a pixel value of the straight line; the straight line equation includes: (x-a)/(b-y) =tan θ— (5) assuming that the relative position of the scale is determined to be (a, b), (in the case of a plot of the truncated scale, the coordinates are the center point position of the screenshot), the tilt angle of the scale is θ.
Step S330, the pixel value of the pixel on the extracted straight line is used for obtaining the fundamental frequency of the relative position information of the skull measurement scale through the cumulative average normalized square error function, and the scale gap value of the skull measurement scale is obtained through the fundamental frequency; the obtaining the fundamental frequency value comprises:
autocorrelation ACF function: is the similarity between the function value y (j) and the function value y (j+T) after the self translation T; t is a y translation number; the autocorrelation function ACF in the application: the degree of correlation between values of the same sequence at different times is reflected. y acquires only a series of values, which are autocorrelation-processed, and the T value corresponds to a value obtained by delaying the series by T. At each j, an ACF value is obtained, and a CSDF value is obtained. W is a window and is typically a windowing process for the signal processing, which may be referred to herein as the length of the truncated data. j is the jth data point.
Step S340, performing least square fitting on the skull measurement scale to obtain a sinusoidal model of the scale, and obtaining the scale phase of the scale through the sinusoidal model.
(Y B (x)-M(x,α)) 2 ，
M-sine model F-is Y B Is used for the spectrum of the (c),for each pixel period in the ruler direction; lambda is the number of pixels between the scale marks; alpha 1 ,α 2 Respectively representing the magnitudes of sin and cos functions, α= { α 1 ,α 2 The } is obtained by iteration.
Step S400, according to the scale gap of the skull measurement scale and the scale phase of the scale, and combining the length average value of the human upper jaw, obtaining the scale phase of the skull measurement scale, and calculating the scale of the skull measurement scale.
S＝[(L ANP /g)*λ]--(7)
S- -scale value of skull measurement staff L ANP -human upper jaw length average value, g number of pixels corresponding to human upper jaw length average value length, lambda-number of pixels between scale gaps of skull measurement scale, []Is rounded.
Referring to fig. 5, the maxillary length (ANS-Ptm): the distance between the pterygoid lobe and the anterior nasal spine and the foot drop on the FH plane represents the length of the upper jaw.
Referring to fig. 9, the present application further provides an embodiment of a system for acquiring a measurement scale in a skull lateral slice, wherein the method for acquiring a measurement scale in skull measurement according to any of the above embodiments may include: a measurement scale image acquisition module 100 for determining a picture of a skull measurement scale in the skull side position slice according to the skull mark points; the scale image position acquisition module 200 is used for carrying out digital processing on the picture of the skull measurement scale acquired by the measurement scale image acquisition module according to a first preset algorithm, and determining the relative position information of the skull measurement scale in the skull side position slice; the scale image processing module 300 performs image processing on the position information of the skull measurement scale acquired by the scale image position acquisition module according to a second preset algorithm to acquire a scale gap of the skull measurement scale and a scale phase of the scale in the skull side position slice; the scale obtaining module 400 is used for calculating the scale of the skull measurement scale according to the scale gap of the skull measurement scale and the scale phase of the scale, which are obtained in the scale image processing module, and combining the length average value of the upper jaw of the human body.
Specifically, the embodiment of the system for acquiring the measuring scale in the skull side position slice of the present application may execute the embodiment of the method for acquiring the measuring scale in the skull side position slice, which is not described herein.
According to the application, the position of the staff gauge is roughly positioned in the skull side position plate, the accurate position of the staff gauge is determined, the scale clearance of the staff gauge and the scale phase of the staff gauge are further calculated in an iterative mode, the staff gauge is realized through a computer language, and the staff gauge is displayed in a data form, so that the problem of low manual participation efficiency in the prior art is solved.
According to the application, the staff scale is obtained by taking the average value of the upper jaw length of a human body in different tooth replacement stages in a standard database and taking the average value as a reference to finish the staff scale obtaining in various modes such as Fourier transform, hough transform and the like, so that the automatic obtaining is realized, manual participation is liberated, the data is more refined, an accurate data basis is further provided for invisible orthodontic, and meanwhile, the working efficiency is improved; and particularly, accurate data information is given for the difficult and complicated cases, so that a perfect experience is provided for the patients in need.
The block diagram of the electronic device 1000 may be a tablet computer, a notebook computer or a desktop computer as shown in fig. 11. Electronic device 1000 may also be referred to by other names of portable terminals, laptop terminals, desktop terminals, etc.
The electronic device 1000 has a processor 1001 and a memory 1002, where the memory 1002 stores a computer program, and when the processor 1001 runs the computer program in the memory 1002, the method for detecting a measurement scale in a skull side bit slice is implemented according to an embodiment.
The processor 1001 may include one or more processing cores, such as a 4-core processor, an 8-core processor, and so on. The processor 1001 may be implemented in at least one hardware form of DSP (Digital Signal Processing ), FPGA (Field-Programmable Gate Array, field programmable gate array), PLA (Programmable Logic Array ). The processor 1001 may also include a main processor, which is a processor for processing data in an awake state, also referred to as a CPU (Central Processing Unit ), and a coprocessor; a coprocessor is a low-power processor for processing data in a standby state.
In some embodiments, the processor 1001 may integrate a GPU (Graphics Processing Unit, image processor) for rendering and drawing of content required to be displayed by the display screen. In some embodiments, the processor 1001 may also include an AI (Artificial Intelligence ) processor for processing computing operations related to machine learning.
Memory 1002 may include one or more computer-readable storage media, which may be non-transitory. Memory 1002 may also include high-speed random access memory, as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. In some embodiments, a non-transitory computer readable storage medium in memory 1002 is configured to store at least one instruction, at least one program, code set, or instruction set for execution by processor 1001 to implement a method of measuring scale detection in a skull side bit slice according to an embodiment of the present application.
In some embodiments, the electronic device 1000 further includes: peripheral interface 1003 and peripheral devices. The processor 1001, the memory 1002, and the peripheral interface 1003 may be connected by a bus or signal line. Peripheral devices may be connected to peripheral device interface 1003 via buses, signal lines, or a circuit board.
In particular to the present embodiment, the peripheral devices may include an intraoral scanner 1004 and a 3D printing device 1005. The processor 1001 obtains the digital dental model in the patient's mouth through the intraoral scanner 1004, the processor 1001 obtains the digital dental model collected by the intraoral scanner 1004 through a program command in the process of executing a computer program, obtains the gingival deformation parameters through executing the measuring scale detection method in the skull side slice according to the embodiment, designs the shell-shaped dental appliance according to the obtained gingival deformation parameters, transmits the data information corresponding to the designed digital shell-shaped dental appliance model to the 3D printing equipment 1005, and directly prints and prepares the shell-shaped dental appliance through the 3D printing equipment 1005.
The present embodiment also provides a computer-readable storage medium, which may be a nonvolatile computer-readable storage medium, and may also be a volatile computer-readable storage medium. The computer readable storage medium has instructions stored therein that, when executed on a computer, cause the computer to perform a method of measuring scale detection in a skull side-slice according to an embodiment.
The modules of the embodiments may be stored in a computer readable storage medium if implemented as software functional modules and sold or used as a stand-alone product. Based on such understanding, the technical solution of the present application may be embodied essentially or in part or all of the technical solution or in part in the form of a software product stored in a storage medium, including instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-only memory (ROM), a random access memory (Random access memory, RAM), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
It should be noted that the above embodiments can be freely combined as needed. The foregoing is merely a preferred embodiment of the present application and it should be noted that modifications and adaptations to those skilled in the art may be made without departing from the principles of the present application, which are intended to be comprehended within the scope of the present application.
Claims (13)
1. A method of acquiring a scale of a measurement scale in a cranial flap, comprising:
determining a picture of a skull measurement scale in the skull side position slice according to the skull mark points;
performing digital processing on the picture of the skull measurement scale according to a first preset algorithm, and determining the relative position information of the skull measurement scale in the skull side position slice;
performing image processing on the skull measurement staff gauge graph according to a second preset algorithm to obtain a scale clearance of the skull measurement staff gauge and a scale phase of the staff gauge in the skull side position slice;
according to the scale gap of the skull measurement scale and the scale phase of the scale, and combining the length average value of the human upper jaw to calculate the scale value of the skull measurement scale,
the acquiring of the scale clearance of the skull measurement scale and the scale phase of the scale in the skull side position slice comprises,
acquiring an inclination angle in a skull patch of the skull measurement scale by using a power density function, calculating coordinates of a point with the maximum power spectrum density by using the power density function, and acquiring the inclination angle based on the coordinates of the point with the maximum power spectrum density;
acquiring a straight line perpendicular to the scale mark of the ruler according to the inclination angle, and extracting pixel values of all pixels on the straight line;
acquiring the fundamental frequency of the relative position information of the skull measurement scale by accumulating the average normalized square error function through the pixel values of the pixels on the extracted straight line, and acquiring the scale clearance value of the skull measurement scale through the fundamental frequency;
and performing least square fitting on the skull measurement staff drawings to obtain a sine model of the staff, and obtaining the scale phase of the staff through the number of pixels between the scales of the staff obtained by the sine model.
2. The method for acquiring the scale of the measuring scale in the skull side position slice according to claim 1, wherein determining the relative position information of the skull measuring scale in the skull side position slice comprises:
dividing the picture of the skull measurement scale into a plurality of sub-blocks according to a preset dividing algorithm;
and carrying out Fourier transform processing on each sub-image block, carrying out marking arrangement on the sub-image blocks with the same frequency in the transformation information of each sub-image block after the Fourier transform processing, and extracting the relative position information of the skull measuring scale in the skull side position piece.
3. The method of obtaining a measurement scale mark in a skull side position slice according to claim 2, the dividing into a plurality of sub-blocks comprising:
L max ＝L K ＝min(N 1 ,N 2 )/N B
L min ＝L 1 ＝min(h,L max /K)
L n ＝(L max /L min ) (n-1)/(K-1) L min
N 1 ,N 2 length of length and width of the picture representing the skull measurement scale, respectively; n (N) B Representing the number of preset divided sub-blocks on the side of minimum length, L n The size of the sub-picture block is the nth preferred; h is the minimum block size of the sub-picture block; k refers to the number of L values to be given, and the length L of the sub-block n-1 With a value of L or less n The method comprises the steps of carrying out a first treatment on the surface of the The smallest part of interest in image processing is the h x h block.
4. The method for acquiring the scale of the measuring scale in the skull side position slice according to claim 1, wherein determining the relative position information of the skull measuring scale in the skull side position slice comprises:
setting a picture of a skull measurement scale in a preset traversal area, and acquiring pose information of the preset traversal area in a skull side position slice; the pose information includes: inclination angle and relative position information of the skull measurement scale in the skull patch.
5. The method for obtaining a measurement scale in a skull side position slice according to claim 4, comprising: generating a traversing line corresponding to each angle value in a preset angle range,
acquiring pixel values of a skull measurement scale map on the generated traversing line;
carrying out Fourier transform on the pixel values of the obtained skull measurement scale image to obtain Fourier transform images at different positions under different angles;
and acquiring traversing line information corresponding to the highest amplitude point in the generated Fourier transform graph, and extracting the pose information of the skull measurement scale graph from the traversing line information.
6. The method for obtaining a measurement scale mark in a skull side position plate according to claim 1, further comprising: and acquiring the linear position information in the skull side position slice through a Hough transformation algorithm, and setting the relative position information of the skull measuring scale in the skull side position slice according to the acquired linear position information.
7. The method for acquiring the scale of the measuring scale in the skull side position plate according to claim 1, wherein the inclination angle of the measuring scale comprises:
the coordinates of the point where the power spectral density is the greatest are (v, u).
8. The method of obtaining a measurement scale in a cranial flap according to claim 7, wherein the power density of each sub-tile comprises: PSDF (w) =f { Y } B }F * {Y B Z, F is Y B Is a frequency spectrum of (c).
9. The method for obtaining a measurement scale mark in a skull side position slice according to claim 8, wherein obtaining the fundamental frequency value comprises:
wherein, the relevant ACF function: is the similarity between the function value y (j) and the function value y (j+T) after the self translation T; t is the y translation number.
10. The method for obtaining a measurement scale mark in a skull side position plate according to claim 9, comprising: by least squares fitting (Y B (x)-M(x,α)) 2 A sinusoidal model of the graduated scale is obtained,
M(x,y；θ R ,λ)＝α 1 cos[F(x,y；θ R ,λ)]+α 2 sin[F(x,y；θ R ,λ)]
wherein M represents a sinusoidal model F represented as Y B 2 pi/lambda is the frequency spectrum of each pixel period in the ruler direction; lambda is the number of pixels between the scale marks; alpha 1 ，α 2 Respectively representing the magnitudes of sin and cos functions, α= { α 1 ，α 2 The } is obtained by iteration.
11. A measuring scale acquisition system in a skull lateral slice, which is capable of performing the method for acquiring measuring scale scales in skull measurement according to any of claims 1-10, comprising,
the measuring scale image acquisition module is used for determining a picture of the skull measuring scale in the skull lateral position slice according to the skull mark points;
the scale image position acquisition module is used for carrying out digital processing on the picture of the skull measurement scale acquired by the measurement scale image acquisition module according to a first preset algorithm, and determining the relative position information of the skull measurement scale in the skull side position piece;
the scale image processing module is used for carrying out image processing on the position information of the skull measurement scale acquired by the scale image position acquisition module according to a second preset algorithm to acquire a scale gap of the skull measurement scale and a scale phase of the scale in the skull side position piece;
the scale acquisition module is used for calculating the scale of the skull measurement scale according to the scale gap of the skull measurement scale and the scale phase of the scale, which are acquired in the scale image processing module, and combining the length average value of the upper jaw of the human body.
12. An electronic device comprising a processor and a memory, the processor executing computer instructions stored by the memory, causing the electronic device to perform the method of measuring scale acquisition in a skull side-piece of any of claims 1-10.
13. A computer storage medium comprising computer instructions which, when run on an electronic device, cause the electronic device to perform the method of measuring scale acquisition in a skull side-piece of any of claims 1 to 10.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011642265.3A CN112700487B (en) | 2020-12-31 | 2020-12-31 | Method and system for acquiring measuring scale marks in skull lateral position slice |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011642265.3A CN112700487B (en) | 2020-12-31 | 2020-12-31 | Method and system for acquiring measuring scale marks in skull lateral position slice |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112700487A CN112700487A (en) | 2021-04-23 |
| CN112700487B true CN112700487B (en) | 2023-09-15 |
Family
ID=75514162
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011642265.3A Active CN112700487B (en) | 2020-12-31 | 2020-12-31 | Method and system for acquiring measuring scale marks in skull lateral position slice |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112700487B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104105442A (en) * | 2012-03-01 | 2014-10-15 | 株式会社日立医疗器械 | Medical image display device and medical image display method |
| EP3261056A1 (en) * | 2016-06-24 | 2017-12-27 | Justus-Liebig-Universität Gießen | Method for measurement and ranking of cranial shape |
| CN110059670A (en) * | 2019-04-29 | 2019-07-26 | 杭州雅智医疗技术有限公司 | Human body Head And Face, limb activity angle and body appearance non-contact measurement method and equipment |
| CN110169782A (en) * | 2019-06-27 | 2019-08-27 | 北京大学第三医院（北京大学第三临床医学院） | A kind of cephalometry method of craniofacial skeleton structure |
| CN111643088A (en) * | 2019-03-04 | 2020-09-11 | 刘帅 | Measuring scale for rapidly measuring soft and hard tissue characteristics of lateral skull tablet |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10007988B2 (en) * | 2014-02-28 | 2018-06-26 | The Board Of Regents Of The Nevada System Of Higher Education On Behalf Of The University Of Nevada, Las Vegas | Systems and methods for approximating the soft tissue profile of the skull of an unknown subject |
-
2020
- 2020-12-31 CN CN202011642265.3A patent/CN112700487B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104105442A (en) * | 2012-03-01 | 2014-10-15 | 株式会社日立医疗器械 | Medical image display device and medical image display method |
| EP3261056A1 (en) * | 2016-06-24 | 2017-12-27 | Justus-Liebig-Universität Gießen | Method for measurement and ranking of cranial shape |
| CN111643088A (en) * | 2019-03-04 | 2020-09-11 | 刘帅 | Measuring scale for rapidly measuring soft and hard tissue characteristics of lateral skull tablet |
| CN110059670A (en) * | 2019-04-29 | 2019-07-26 | 杭州雅智医疗技术有限公司 | Human body Head And Face, limb activity angle and body appearance non-contact measurement method and equipment |
| CN110169782A (en) * | 2019-06-27 | 2019-08-27 | 北京大学第三医院（北京大学第三临床医学院） | A kind of cephalometry method of craniofacial skeleton structure |
Non-Patent Citations (1)
| Title |
|---|
| X线头影测量计算机自动识别系统的建立;何艳,邵金陵,张辉,吴刚,吴莎;中国美容医学(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112700487A (en) | 2021-04-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Damas et al. | Forensic identification by computer-aided craniofacial superimposition: a survey | |
| Wang et al. | Automatic analysis of lateral cephalograms based on multiresolution decision tree regression voting | |
| Grau et al. | Automatic localization of cephalometric landmarks | |
| Guyomarc'h et al. | Anthropological facial approximation in three dimensions (AFA 3D): Computer‐assisted estimation of the facial morphology using geometric morphometrics | |
| CN109615642B (en) | Automatic organ-at-risk delineation method and device in radiotherapy plan | |
| CN107833248B (en) | Medical image scanning method and medical imaging equipment | |
| KR101090375B1 (en) | Ct image auto analysis method, recordable medium and apparatus for automatically calculating quantitative assessment index of chest-wall deformity based on automized initialization | |
| JP4504377B2 (en) | Automatic optimal surface determination for heart related acquisition | |
| US20150235104A1 (en) | Method and system for registration/alignment of 3d digital models | |
| CN112037200A (en) | Method for automatically identifying anatomical features and reconstructing model in medical image | |
| CN108888284A (en) | Image adjusting method, device and equipment, storage medium | |
| CN115797730B (en) | Model training method and device and head shadow measurement key point positioning method and device | |
| CN115222937A (en) | Method and device for detecting scoliosis | |
| CN112790879A (en) | Tooth axis coordinate system construction method and system of tooth model | |
| Lee et al. | Variation within physical and digital craniometrics | |
| WO2022160731A1 (en) | Image processing method and apparatus, electronic device, storage medium, and program | |
| CN113241155B (en) | Method and system for acquiring mark points in skull side position slice | |
| JP5640280B2 (en) | Osteoporosis diagnosis support device and osteoporosis diagnosis support program | |
| CN112700487B (en) | Method and system for acquiring measuring scale marks in skull lateral position slice | |
| CN111223117A (en) | Tongue surface image segmentation device, method and computer storage medium | |
| CN116158774B (en) | Methods, apparatus, and computer storage media for measuring tibial plateau caster angle | |
| LEITÃO | de Menezes | |
| da Silva et al. | A systematic review of photogrammetry as a reliable methodology in gender identification of human skull | |
| US20220051428A1 (en) | Systems and methods for assessment of nasal deviation and asymmetry | |
| US20210201512A1 (en) | Method and apparatus for registering live medical image with anatomical model |
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 | ||
| CB02 | Change of applicant information |
Address after: 201210 Pudong New Area, Shanghai, 2305, A, 122, north gate, two floor. Applicant after: Zhengya Dental Technology (Shanghai) Co.,Ltd. Address before: 201210 Pudong New Area, Shanghai, 2305, A, 122, north gate, two floor. Applicant before: SHANGHAI SMARTEE DENTI-TECHNOLOGY Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |