Device for analyzing occlusion pressure, program for analyzing occlusion pressure, and method for analyzing occlusion pressure

Abstract

An object is to provide a device for analyzing occlusion pressure with which it is possible to measure occlusion pressure with a high positional accuracy by using data obtained by versatile methods, including an operation device, wherein the operation device calculates a formula of a first dental arch based on an occlusion pressure data, calculates a formula of a second dental arch based on an occlusal contact positional data, calculates a transformation matrix that approximates the formula of the first dental arch to the formula of the second dental arch, carries out an operation to transform positional information included in the occlusion pressure data to a new positional information with the transformation matrix, and outputs an occlusion pressure data including the new positional information.

Claims

1. A device for analyzing occlusion pressure comprising: at least one memory storing program code; and at least one hardware processor configured to execute the program code, the program code including: generating code configured to cause the at least one hardware processor to generate a formula of a first dental arch based on an occlusion pressure data, and generate a formula of a second dental arch based on an occlusal contact positional data, generate a transformation matrix that approximates the formula of the first dental arch to the formula of the second dental arch; transforming code configured to cause the at least one hardware processor to transform positional information included in the occlusion pressure data to a new positional information based on the transformation matrix; and output code configured to cause the at least one hardware processor to output, based on the transformation, an occlusion pressure data including the new positional information, wherein outputting the occlusion pressure data including the new positional information comprises overlaying and displaying of the occlusion pressure data including the new positional information to an image of at least a row of teeth.

2. A non-transitory computer readable medium comprising instructions which, when executed, cause at least one hardware processor to analyze occlusion pressure by: generating a formula of a first dental arch based on an occlusion pressure data; generating a formula of a second dental arch based on an occlusal contact positional data; generating a transformation matrix that approximates the formula of the first dental arch to the formula of the second dental arch; and transforming positional information included in the occlusion pressure data to a new positional information based on the transformation matrix.

3. A method for analyzing occlusion pressure performed by at least one processor and comprising: generating a formula of a first dental arch based on an occlusion pressure data; generating a formula of a second dental arch based on an occlusal contact positional data; generating a transformation matrix that approximates the formula of the first dental arch to the formula of the second dental arch; and transforming positional information included in the occlusion pressure data to a new positional information based on the transformation matrix.

Description

BRIEF DESCRIPTION OF DRAWINGS

(1) FIG. 1 is a view to explain a flow of an occlusion pressure analysis method S10.

(2) FIGS. 2A and 2B include views schematically showing a curve obtained with the formula A of a first dental arch (FIG. 2A) and a curve obtained with the formula B of a second dental arch (FIG. 2B).

(3) FIG. 3 is a view to explain that the positional information (X, Y) of occlusion pressure data is transformed by a transformation matrix to obtain a new positional information (X′, Y′).

(4) FIG. 4 is a view to explain the structure of a device for analyzing occlusion pressure 10.

DETAILED DESCRIPTION OF EMBODIMENTS

(5) FIG. 1 shows a flow of an occlusion pressure analysis method S10 according to one embodiment. As seen from FIG. 1, the occlusion pressure analysis method S10 includes the steps S11 to S18. Each step is explained below.

(6) The step S11 is to obtain occlusion pressure data. That is, data of two-dimensional positional information (X, Y) and data of occlusion pressure information (P.sub.XY), of occlusion pressure distributions from the means of detecting occlusion pressure including occlusion pressure distribution information, are obtained.

(7) Here, the means of detecting occlusion pressure is not particularly limited. Examples thereof may include a pressure sensitive paper (pressure sensitive sheet) and a pressure sensitive sensor that detects the positions of occlusal contact with changes of electric properties. Among them, a pressure sensitive paper, which is more versatile than the others and can display occlusion pressure easily, is preferable.

(8) The method for obtaining X, Y, and P.sub.XY from the means of detecting occlusion pressure is not particularly limited. When the means of detecting occlusion pressure is a pressure sensitive paper, examples of the method may include a scanner or camera to read in colors and gradations present on the pressure sensitive paper. If the means of detecting occlusion pressure is a pressure sensitive sensor, examples of the method may include a device to measure electrical properties (voltage, current, resistance value, or variation of them) corresponding to occlusion pressure.

(9) The step S12 is to calculate a formula A of a dental arch (formula A of a first dental arch) with occlusion pressure data as shown in FIG. 2A, from the information data of occlusion pressure generation points (X, Y) among the occlusion pressure data obtained by the step S11.

(10) The calculation of the formula A can be made by regression analysis with the distribution of the occlusion pressure generation points included in the occlusion pressure data. Examples of the regression analysis include least square method.

(11) A more specific example of the calculation is as follows.

(12) The data in which the occlusion pressure (P.sub.XY) is no less than P.sub.0, which is a predetermined threshold value, is extracted, and a regression analysis for the positional information data of the extracted data (X, Y) is carried out, to obtain the following formula of quadratic curve.

(13) Quadratic curve:
A.sub.1X.sup.2+A.sub.2XY+A.sub.3Y.sup.2+A.sub.4X+A.sub.5Y+A.sub.6=0

(14) Here, X and Y are variable numbers and A.sub.1 to A.sub.6 are coefficients obtained by the regression analysis.

(15) The threshold value P.sub.0 may be appropriately determined in view of the performance of the means of detecting occlusion pressure and securing the numbers of data enough to obtain the formula of the quadratic curve.

(16) In addition, the above-described formula of quadratic curve may be a formula of parabola in which rotation is controlled or a formula of ellipse, in view of easy analysis. They are specifically as follows. B.sub.1 to B.sub.3 and C.sub.1 to C.sub.5 in the formulas are coefficients obtained by the regression analysis.

(17) Parabola:
Y=B.sub.1X.sup.2+B.sub.2X+B.sub.3

(18) Ellipse:
C.sub.1X.sup.2+C.sub.2X+C.sub.3Y.sup.2+C.sub.4Y+C.sub.5=0

(19) This parabola has an axis parallel to Y axis, and the ellipse has an axis parallel to X axis and Y axis. Therefore, it is possible to make the calculation easy by reducing the coefficients to be obtained.

(20) In addition, in view of easy analysis, not all of the information data of occlusion pressure generation points of no less than the threshold value P.sub.0, but only a position taking the largest value of each occlusion pressure area that can be extracted at the threshold value P.sub.0, (the area where occlusion pressure generation points connect to each other to be a mass so that an area is formed) or only the data of barycentric position etc. of this area may be used for the regression analysis.

(21) The step 13 is to obtain occlusal contact positional data. That is, it is a step to obtain two-dimensional positional information (x, y) or three-dimensional positional information (x, y, z) where occlusal contact positional distribution is emerged.

(22) Here, the sources of the information where the occlusal contact positional distribution is emerged are not particularly limited, and the following sources may be exemplified. An impression of oral cavity, or two-dimensional images (pictures, scanner images, etc.) of wax taken when occlusion is made in which occlusal contact positions are shown A picture of oral cavity in which occlusal contact positions are shown A two-dimensional image (picture, scanner image, etc.) of a model of row of teeth (plaster model etc.) in which occlusal contact positions are shown Shape data of three-dimensional scan of an oral cavity in which occlusal contact positions are shown Shape data of three-dimensional scan of a model of row of teeth (plaster model etc.) in which occlusion contact positions are shown

(23) Here, the method for showing occlusal contact positions in each picture and data is not particularly limited. For example, the positions in an impression of oral cavity and bitten wax may be identified by shape measurement or light transmission measurement on the portions where became thin by occlusion. The positions in an oral cavity and a model of row of teeth may be identified by image processing with marking of occlusal contact points with a biting paper etc., or with a pressure sensitive paper to which the function of biting paper is added.

(24) As explained above, in the occlusion pressure analysis method S10, the models, pictures, and data used in another treatment and the like may be used as the method for obtaining the occlusal contact positional data, with which general versatility may be increased.

(25) It is preferred that the occlusal contact positional data be obtained from the direction of the normal line of the occlusal plane. This makes it possible to obtain a quadratic curve comparatively similar to the shape of occlusion pressure data, and to increase the correction accuracy.

(26) When pictures etc. obtained from a position different from the direction of the normal line (such as a direction oblique to the normal direction having an angle), the positional information may be identified three-dimensionally (three-dimensional coordinate) by using pictures taken from two or more directions, or by including a jig having a known coordinate in the image, to correct the occlusal contact positional data.

(27) The obtainment method of the data of two-dimensional occlusal contact positional information (x, y) or three-dimensional occlusal contact positional information (x, y, z) where the occlusal contact positional distribution is emerged is not particularly limited, and a three-dimensional scanner may be exemplified.

(28) The step S14 is to calculate a formula B of a dental arch (formula B of a second dental arch) with occlusion pressure data as shown in FIG. 2B, from the distribution of information data of occlusion pressure generation points (x, y) among the occlusal contact positional data obtained by the step S13.

(29) The calculation of the formula B can be made by calculation of the dental arch by regression analysis with the distribution of the occlusal contact positional data. Examples of the regression analysis include least square methods. It is preferred that the kind of the formula be the same as the formula A of the first dental arch calculated in the above step S12. That is, it is preferred that the formula B of the second dental arch be a formula of parabola when the formula A of the first dental arch is a formula of parabola, and that the formula B of the second dental arch be a formula of ellipse when the formula A of the first dental arch is a formula of ellipse. This makes it possible to increase transformation accuracy which is described later.

(30) It is preferred that, when the occlusal contact positional information is three-dimensional, the formula of occlusal plain face or occlusal curved face as shown below be obtained by regression analysis etc. first, then the formula of dental arch projected to the plain face or curved face be obtained in the same way as in the two-dimensional way. Regarding the occlusal curved face, the formula may be decided based on theories such as Monson Spherical Theory and Spee Curve.

(31) Formula of occlusal plain face:
D.sub.1x+D.sub.2y+D.sub.3z+D.sub.4=0

(32) Formula of occlusal curved face (spherical surface):
E.sub.1x.sup.2+E.sub.2x+E.sub.3y.sup.2+E.sub.4y+E.sub.5z.sup.2+E.sub.6z+E.sub.7=0

(33) Here, D.sub.1 to D.sub.4 and E.sub.1 to E.sub.7 are coefficients.

(34) In this case as well, in view of easy analysis, not all of the occlusal contact positional data but only data of such as barycentric position of each occlusal contact area may be used for the regression analysis.

(35) In the step S15, a transformation matrix to approximate the formula A of the first dental arch to the formula B of the second dental arch is obtained by operation. When the measurement of occlusion pressure is carried out with a pressure sensitive paper for example, the formula A of the first dental arch is considered to be different from the actual dental arch since the shape of the pressure sensitive paper is different when occlusion is made and when data is analyzed, as mentioned above. Therefore, the formula A of the first dental arch is transformed in the way to be approximated to the formula B of the second dental arch which is considered to be more similar to the actual dental arch.

(36) As the transformation matrix R.sub.1 when the positional data is three-dimensional, for example the following matrix may be exemplified. This is in the case where the occlusal contact positional information is three-dimensional (x, y, z), and a.sub.11 to a.sub.43 are coefficients.

(37) $\begin{array}{cc}{R}_1*\left(\begin{array}{c}X\\ Y\\ 1\end{array}\right)=\left(\begin{array}{ccc}{a}_{11}& a_{12}& a_{13}\\ a_{21}& a_{22}& a_{23}\\ a_{31}& a_{32}& a_{33}\\ a_{41}& a_{42}& a_{43}\end{array}\right)*\left(\begin{array}{c}X\\ Y\\ 1\end{array}\right)=\left(\begin{array}{c}x\\ y\\ z\\ 1\end{array}\right)& \left[\mathrm{Formula}1\right]\end{array}$

(38) As the transformation matrix R.sub.2 when the positional data is two-dimensional, for example the following matrix may be exemplified. This is in the case where the occlusal contact positional information is two-dimensional (x, y), and b.sub.11 to b.sub.33 are coefficients.

(39) $\begin{array}{cc}{R}_2*\left(\begin{array}{c}X\\ Y\\ 1\end{array}\right)=\left(\begin{array}{ccc}{b}_{11}& b_{12}& b_{13}\\ b_{21}& b_{22}& b_{23}\\ b_{31}& b_{32}& b_{33}\end{array}\right)*\left(\begin{array}{c}X\\ Y\\ 1\end{array}\right)=\left(\begin{array}{c}x\\ y\\ 1\end{array}\right)& \left[\mathrm{Formula}2\right]\end{array}$

(40) Calculation of the coefficients a.sub.11 to a.sub.43 and b.sub.11 to b.sub.33 is not particularly limited, and for example they may be obtained in the following ways. Here, two examples are shown.

(41) In one example, two or more characteristic points that the formulas A and B correspond to each other are extracted, and the above-mentioned coefficients are obtained from a system of equations with the points.

(42) For example, considering the transformation matrix R.sub.2 which is two-dimensional, when the formulas A and B are parabolas, the coefficients b.sub.11 to b.sub.33 that satisfy transformation equation are obtained with the coordinate values (X, Y, x, y) of the characteristic points such as the apex, focal point, directrix of each parabola by solving a system of equations.

(43) On the other hand, when the formulas A and B are ellipses, the coefficients b11 to b33 that satisfy transformation equation are obtained with the coordinate values (X, Y, x, y) of the characteristic points such as the focal point, directrix, center, edge of major axis, edge of minor axis of each ellipse by solving a system of equations.

(44) In this example, together with the above characteristic points, the point closest to the occlusion point positioned on the most distal side of each formulas A and B may also be used as a characteristic point. With this point, not only correction of the shape but also correction of scale may be carried out, which makes it possible to further increase the accuracy of the position of the final occlusion point.

(45) In another example, affine transformation matrixes are obtained to be the transformation matrixes R.sub.1 and R.sub.2. That is, the affine transformation matrixes are obtained by specifying the rotation angle of the axis, magnification/reduction percentage of the axis, the positional shift of the origin of each of the quadric curved of the formulas A and B, to be the transformation matrixes R.sub.1 and R.sub.2.

(46) It is noted in the step S15, after the calculation of the transformation matrixes of R.sub.1 and R.sub.2, fine adjustment may be manually allowed to the transformation matrixes of R.sub.1 and R.sub.2 with consideration of the condition of the actual row of teeth, based on the results of the calculation.

(47) In the step S16, as schematically shown in FIG. 3, the transformation matrixes R.sub.1 and R.sub.2 obtained by the step S15 are applied to each of the positional data of occlusion pressure obtained by the step S11 (data of coordinates (X, Y) shown by o in FIG. 3), to obtain occlusion pressure positional data after transformation (X′, Y′). This new occlusion pressure positional data after transformation is transformed with the transformation matrixes R.sub.1 and R.sub.2. Thus, the data corresponds to the occlusal contact positional data with a good accuracy and shows the positions where occlusion pressure is made with a good accuracy.

(48) This means that the occlusion pressure distribution is corrected to proper positions with a good accuracy with the data obtained by versatile methods (occlusion pressure data, occlusal contact positional data).

(49) The occlusion pressure does not always occur at the all occlusion contact positions, which makes it difficult to individually coordinate the occlusion pressure data and occlusal contact positional data, and the coordination is very troublesome. Comparatively speaking, making transformation with the dental arch like the present embodiment easily makes the coordination with a good accuracy.

(50) In the step S17, the occlusion pressure data (P.sub.xy) among the above-described occlusion pressure data is added to the occlusion contact positional data. With the transformation of the positional information by the above transformation matrixes R.sub.1 and R.sub.2, the new occlusion pressure positions after transformation (X′, Y′) takes the same or similar value as the coordination of the occlusal contact positional data (x, y). Thus, the new occlusion pressure positional data obtained after transformation may be added to the coordination of the occlusal contact positional data as it is.

(51) With this, the positions of the occlusion pressure and the occlusal contact positions are combined with a good accuracy.

(52) It is noted that, after the occlusion pressure position after transformation (X′, Y′) is obtained, the agreement rate of the occlusion pressure position after transformation (X′, Y′) and the coordination of the occlusal contact positional data (x, y) is calculated, and if the agreement rate is low, the calculation may be carried out again with change of at least one of the formulas A, the formula B, the transformation matrix R.sub.1, and the transformation matrix R.sub.2. In this case, regarding the calculation of the agreement rate, it is preferred that the parameter be the data of the occlusion pressure position (X′, Y′) since the occlusion pressure does not always occur at the all occlusal contact positions.

(53) In the step S18, the occlusion pressure including positional information obtained in the step S17 is overlaid with image of row of teeth to be shown. This makes it possible to visually see the positions where occlusion pressure is made and their sizes. As the image of row of teeth, the image used in the step S13 may be used for example.

(54) According to the occlusion pressure analysis method S10 described above, it is possible to correct the occlusion pressure distribution to proper positions with a good accuracy with the data from information obtained from versatile methods, regardless of acquisition methods of the occlusion pressure and occlusal contact positional information which are sources of data.

(55) By using this analysis result, it is possible to carry out occlusion pressure measurement at the accurate occlusion pressure positions.

(56) It is noted that when the occlusion area is changed because of enlargement and contraction of the axis with the occlusion pressure position (X′, Y′) transformed with the transformation matrixes R.sub.1, R.sub.2, the correction may be made to show the area correctly and the correction of the value of the occlusion pressure may be made so that the value of occlusion power at occlusion points (occlusion area x occlusion pressure) becomes correct.

(57) Next, a device for analyzing occlusion pressure 10 according to one embodiment will be explained. FIG. 4 conceptually shows the structure of the device for analyzing occlusion pressure. The device for analyzing occlusion pressure 10 of the present embodiment includes an operation device 11, an input means 12, and a display device 13. The operation device 11 includes an operation means 11a, a RAM 11b, a memorizing means 11c, a receiving means 11d, and an output means 11e. The input means 12 includes a keyboard 12a, a mouse 12b, an external storage device 12c which functions as one storage media, and a communication means 12d.

(58) The operation means 11a is formed from a so-called CPU (central processing unit), connected to each structural member, and can control the members. The operation means 11a also executes various programs stored in the memorizing means 11c etc. that function as a storage media, and functions as a means to carry out calculation included in the above occlusion pressure analysis method S10 and generation of various data based on the programs.

(59) RAM 11b is a structural member functioning as a work area of the operation means 11a and a temporal data memorizing means. RAM 11b may be formed from SRAM, DRAM, flush memory, etc., and is same as known RAMs.

(60) The memorizing means 11c is a structural member functioning as a storage media to save the above-described operation, various programs which are to be bases of the obtainment of data to be generated, and the data. Thus, as one of the programs, the storage means 11c storages a program including each step to execute the occlusion pressure analysis method S10. The memorizing means 11c may also storage various intermediate and final results obtained by the execution of the programs.

(61) The receiving means 11d is a structural member having a function to properly take information from the outside into the operation device 11, to which the input means 12 is to be connected. A so-called input port, input connecter etc. are included in the receiving means 11d.

(62) The output means 11e is a structural member having a function to properly output information that is to be output to the outside, among the obtained results, to which the display device 13 and various devices are to be connected. A so-called output port, output connecter etc. are included in the output means 11e.

(63) The input device 12 includes the keyboard 12a, the mouse 12b, the outer storage device 12c, the transmitting means 12d etc. for example. A known keyboard and mouse may be used for the keyboard 12a and the mouse 12b and explanations thereof are omitted.

(64) The outer storage means 12c is a known memorizing means that can make external connection and can also function as a storage media. Necessary various programs and data may be stored in the outer storage means 12c without particular limitation. For example, programs similar to the programs recorded in the above-described memorizing means 11c may be stored. As the outer storage device 12c, a known device may be used. The examples thereof include CD-ROM, CD-ROM drive, DVD, DVD drive, hard discs, various memories etc. The transmitting means 12d is a means to transmit data from another device to the receiving means 11d directly. The examples of the transmitting means 12d include the occlusion pressure detecting means shown in the step S13, a means to obtain occlusion pressure from the detecting means, and a device to obtain the occlusal contact positional data.

(65) The display device 13 is a device to show operation results of each situation, selections of various functions etc. The operator operates the input means 12 referring to the display device 13 to make the operation device 11 function. As the display device 13, a known display device may be applied. For example, a monitor and the like may be used.

(66) Examples of the operation device 11 include a computer. Normally the computer itself includes each device corresponding to the operation device 11, to which the input means 12 and the display device 13 may be connected.

(67) According to the operation device 11 described above, based on the programs stored in the operation device 11 having each step to carry out the above steps S11 to S18, the occlusion pressure data and the occlusal contact positional data are input to the operation device 11 from the input means 12 as described in the steps S11 and S13, and the above-described operation of each step may be carried out by the operation device 11. Regardless of the acquisition method of the occlusion pressure and the occlusal contact positional information which are sources of data, it is possible to correct the occlusion pressure distribution to proper positions with a good accuracy by means of the data from information obtained by versatile methods.

(68) The device for analyzing occlusion pressure that corrects the occlusion pressure distribution to proper positions with a good accuracy is explained here. This device may be configured to be incorporated in an occlusion pressure measurement device or an occlusal contact position measurement device to obtain occlusion pressure with the results by the device for analyzing occlusion pressure.

(69) Alternatively, the device for analyzing occlusion pressure 10 may be configured to have at least one of the functions of occlusion power measurement, occlusion pressure measurement, and occlusal contact position measurement.

REFERENCES SIGN LIST

(70) 10 device for analyzing occlusion pressure 11 operation device 12 input means 13 display device