Method and apparatus for judging deformation of vertebral body

Abstract: A method for judging deformation of a vertebral body, which comprises measuring a central length (a), rear brim length (c), and front brim length (d) from a profile X-ray image of a vertebral body to be judged, determining ratios of c/d, a/c and a/d, and classifying deformation of the vertebral body to be judged by using indices, c, d, c/d, a/c and aid is provided. Furthermore, a method is provided for judging deformation of a vertebral body, which comprises obtaining at least two indices by measuring a central length (a), vertebral body width (b), rear brim length (c), and front brim length (d) from a profile X-ray image of a vertebral body, analyzing these indices by a discriminant function, and classifying a deformation of the vertebral body from results of the analysis.

Description: The present invention relates to a method and apparatus for judging deformation of a vertebral body. More particularly, the present invention relates to a method and apparatus for judging deformation of a vertebral body which is one of the bones forming the spinal column. A judgment of the presence of vertebral body deformation accompanied by osteoporosis, as well as a classification of the deformation, is very important for grasping the progress of osteoporosis and for confirmation of the effect of therapy. In the prior art, the presence of deformation of a vertebral body which is one of the bones forming the spinal column has been judged by a physician by visual observation of the profile X-ray images of thoracic vertebrae or lumbar vertebrae. However, compared with a fracture of the so called long bones such as the femur, tibia, radius, ulna, etc., a judgment of vertebral deformation is difficult since this is usually a wedge-shaped deformation, pressure fracture, depressed fracture, etc., and any judgment inevitably involves individual differences. Also, even if an attempt is made to monitor changes over a period of time, for confirmation of the effect of therapy, because the deformation cannot be quantified, progress of the deformation, if any, can not be easily determined. EP-A 0069229 discloses a method for reducing the effect of patient movement in X-ray diagnosis, in which different images are derived from a video signal containing a sequence of images of a patient under examination, first and subsequent images contained in the signal being stored and used to form a secondary mask which is then subtracted from images contained in the video signal. A system for the analysis of spinal deformities is known from Medical & Biomedical Engineering & Computing 1982, 20 , 715-726 in which a deformity in the curvature of the spine as a whole is measured by the digitization of X-ray data. Also, although a method has been reported in which the ratio (a/d) of the central length (a) to the front brim length (d) of the third lumbar vertebrae is determined as the index (Lumbar Spine Score) for the degree of bone atrophy, or a change in the longitudinal and lateral bone beams of the third lumbar vertebrae is observed as the index for the degree of bone atrophy (severity classified by Jikei University), or a judgment of the fracture of a vertebral body is made by measurement of the central length (a), front brim length (d) and rear brim length (c) of a vertebral body [The New England Journal of Medicine, vol. 306, 446 (February 25, 1982)], a method for classifying the type of deformation of a vertebral body has not heretofore been known. Accordingly, the objects of the present invention are to eliminate the above-mentioned disadvantages of the prior art and to provide a method for judging deformation of a vertebral body in a subject, characterised in that the method comprises the steps of: (i) measuring a central length (a), a rear brim length (c) and a front brim length (d) from a plurality of profile X-ray images obtained from vertebral bodies having substantially no deformation at a portion corresponding to a portion of said vertebral body to be judged, storing the measurements in a memory, and determining values of c/d, a/c and a/d, the average values c and d of the rear brim lengths (c) and front brim lengths (d), respectively, and the standard deviations σ c and σ d of the values c and d, respectively: (ii) measuring a, c, and d from a plurality of the profile X-ray images obtained from vertebral bodies having at least one deformation, storing the measurements in said memory, and determining values of c/d, a/c, a/d, and c , d , σ c , and σ d ; (iii) preparing a judgement standard from the data obtained in steps (i) and (ii) for determining the nature of the deformity according to a classification consisting of "no deformation" vertebrae, "wedge-shaped" vertebrae, "inverse wedge-shaped" vertebrae, "fish" vertebrae, and "flat" vertebrae; (iv) measuring a, c, and d from a profile X-ray image of said vertebral body to be judged, storing the measurements in said memory, and determining c/d, a/c, and a/d; and (v) classifying deformation of the vertebral body to be judged from c/d in the case of "wedge-shaped" and "inverse wedge-shaped" vertebrae, from c , d , σ c , σ d , and c/d in the case of "flat" vertebrae, and from c , d , σ c , σ d , a/c, a/d, and c/d in the case of the "no deformation" and "fish" vertebrae, in order to judge deformation and to monitor the effect of therapy on said vertebral body to be judged. In accordance with the present invention, there is further provided apparatus for judging deformation of a vertebral body according to the method of claim 1, the apparatus being defined in claim 7. The present invention will be better understood from the description set forth below with reference to the accompanying drawings, in which: Figure 1 shows a profile X-ray image of a vertebral body; Figure 2 shows the relationships of the respective vertebral body deformations and front brim length (d), rear brim length (c), etc., according to the judgment method of the present invention; Figure 3 shows an example of vertebral body deformation; Figure 4 shows the front brim length, central length, rear brim length, and vertebral body width in a profile X-ray image of a vertebral body; Figure 5 is a block diagram illustrating the basic structure of apparatus for judging deformation of a vertebral body by using an X-ray image of the vertebral body; Figure 6 is a block diagram illustrating apparatus for judging deformation of a vertebral body in the X-ray image of the vertebral body to be judged by using a discriminating standard obtained mainly from X-ray images of vertebral bodies having no substantial deformation, in which a means for inputting information relating to X-ray images of vertebral bodies having deformation caused therein is further required; Figure 7 is a block diagram illustrating apparatus for judging deformation of a vertebral body by using, as a criteria, a discriminating standard obtained mainly by using a discriminating function; Figure 8 shows apparatus for judging deformation of a vertebral body comprising a digitizer for inputting the points 1 to 6 in Figure 1 or the points 1 to 8 in Figure 4 from the X-ray image to a computer; a personal computer for classifying the presence or absence of the deformation and types of the deformation by measuring the lengths a , c , and d in Figure 1 or a , b , c and d in Figure 4, followed by the image treatment; a memory including software for effecting the abovementioned treatment; and a printer for printing out, for example, the measured values or the types of deformation; and Figure 9 shows apparatus for judging deformation of a vertebral body comprising a TV camera for reading the image of each vertebral body from the X-ray image; a personal computer for classifying the presence or absence of the deformation and the types of the deformation by measuring the length a , c , and d in Figure 1 or the lengths a , b , c , and d in Figure 4 from the read X-ray image; a memory including software for registering the above images and effecting the abovementioned treatment. The present inventors have made an intensive study of a method of objective evaluation of vertebral body deformation, and found that the presence of vertebral body deformation can be objectively judged according to the judgment standard in which lengths of the vertebral body and their ratios are combined by accurately measuring the lengths of a front brim, rear brim, central portion from a profile X-ray image (i.e., an image on a film taken from the side by X-ray radiation) of a vertebral body and determining the ratios of the respective portions of the vertebral body, and further found that the type of vertebral body deformation can be judged according to the judgment standard in which lengths of the vertebral body and their ratios are combined, and the progress over a period of time of the vertebral body deformation can be judged from the change in the deformation type as well as a change in the lengths of the vertebral body, and thus accomplished the present invention. In the prior art, the profile X-ray image of vertebral body was visually observed by a physician when judging the vertebral body deformation, but according to the present invention, the lengths of the respective portions of the vertebral body are measured, judgment standards are prepared with reference to a judgment by a physician, and the presence of vertebral body deformation as well as a classification of the deformation type are objectively performed following the judgment standards. The judgment method of the present invention is now described in more detail. (i) First, from a profile X-ray of a vertebral body without deformation, the central length (a), rear brim length (c), and front brim length (d) of each vertebral body are measured to determine an average value ( c ) of c , an average value ( d ) of d , and a standard deviation (σ c , σ d ). The above measurement is used to obtain standard values which become the judgment standard for the central length (a), rear brim length (c), and front brim length (d), etc., determined from a patient whose vertebral body deformation is to be judged. The vertebral body without deformation to be measured is preferably a vertebral body of a person of from 50 to 75 years old, in view of the age of the patient whose vertebral body deformation is to be judged. Also, it is preferable to perform measurements separately for men and women, and it is desirable to perform measurements of about 50 vertebral bodies, or more, for both men and women. When making a profile X-ray image of a vertebral body, the thoracic vertebrae from the first thoracic vertebrae to the twelfth thoracic vertebrae, lumbar vertebrae from the first lumbar vertebrae to the fifth lumbar vertebrae should be included, giving a total of 17 vertebral bodies, and therefore, for example, it is advantageous to take the photograph of the profile X-ray image with the eighth thoracic vertebrae as the center separately from the profile X-ray image with the third lumbar vertebrae as the center. In the profile X-ray image with the eighth thoracic vertebrae as the center, the first thoracic vertebrae cannot be accurately measured in most cases, and further, the second and the third thoracic vertebrae are unclear in many cases. However, since the frequency of vertebral body deformation is not great in these vertebral bodies, it is sufficient if measurement can be made from the fourth thoracic vertebrae or the fifth thoracic vertebrae. The center length (a), rear brim length (c), and front brim length (d) of a vertebral body are specifically as shown in Fig. 1. For determination of these values, for the profile X-ray image of a vertebral body, the central length (a), rear brim length (c), and front brim length (d) may be measured for each vertebral body with a scale or slide calipers, but the six points of 1 to 6 shown in Fig. 1 can be input through a digitizer into a computer to measure the central length (a) (3-4), rear brim length (c) (5-6), and front brim length (d) (1-2), and to calculate the ratios of c/d, a/c, a/d. Also, for example, the degree of blackening may be recorded with a TV camera, etc., and the lengths of the respective portions of the vertebral body may be automatically measured by image processing. (ii) Subsequently, from the profile X-ray image of the vertebral body of the patient whose vertebral body deformation is to be judged, as in (i), the center length (a), rear brim length (c), and front brim length (d) of each vertebral body are measured, and further, the ratios of c/d, a/c and a/d are determined by calculation. These measurements can be conducted exactly as in the case of the above (i). (iii) Next, a judgment of vertebral body deformation can be made from a , c , d , c/d, a/c and a/d, for example, as follows. That is, (A) when c/d, c , d , a/c, a/d satisfy the following conditions, the judgment can be "no deformation" (type N): (a) 0.7 < c/d < 1.4 ; (b) at least one of c ≧ c - 2σ and d ≧ d - 1.5σ is satisfied; (c) at least one of a/c > 0.8 and a/d > 0.8 is satisfied. As the judgment standard for such "no deformation", the ratio c/d of the rear brim length to the front brim length must be between 0.7 and 1.4 as the first condition (a). For, as described below, if c/d becomes 1.4 or more, the front brim portion becomes a deformed wedge-shaped vertebrae, and if c/d becomes 0.7 or less, the rear brim portion becomes a deformed inverse wedge-shaped vertebrae, which in practice is very rare. Next, as the second condition (b), at least one of rear brim length c and front brim length d must be greater than c - 2σ and d - 1.5σ when the average values of the rear brim length c and front brim length d of vertebral body without deformation determined in (i) are made c and d , respectively. In a typical healthy vertebral body, a , c , d all have a high value, namely a ≧ a - 2σ a , c ≧ c - 2σ c , d ≧ d - 1.5σ d , but even when the vertebral body is slightly compressed to become a ≧ a - 2σ a , provided that at least one of the conditions of c ≧ c - 2σ c or d ≧ d - 1.5σ d , can be satisfied, no clear vertebral body deformation is recognized and judgment of "no deformation" may be made. As the third condition, at least one of a/c and a/d must be more than 0.8. For, if both of a/c and a/d are 0.8 or less, namely the center length a is much smaller, than the rear brim length c and front brim length d , a fish vertebrae condition as described below is determined. Thus, for the above reasons, when the above conditions of (a), (b) and (c) are satisfied, a judgment of "no deformation" can be made. The relationship of the conditions of the above (a), (b), and (c), or the respective conditions as described below with the respective types of vertebral body deformation, inclusive of "no deformation", are shown in Fig. 2. (B) When c/d ≧ 1.4 , the judgment can be a "wedge-shaped vertebrae" (type I). In the wedge-shaped vertebrae, deformation has occurred at the front brim portion as shown in Fig. 3, and the front brim length (d) has become much smaller than the rear brim length (c). For example, when the front brim length is 3/4 of the rear brim length ( c/d = 1.33 ), a wedge-shaped vertebral body deformation can be easily recognized, but in some cases, a clear wedge-shaped vertebrae cannot be easily recognized. On the other hand, when the front brim length is 2/3 of the rear brim length ( c/d = 1.5 ), it can be clearly recognized as a wedge-shaped vertebrae, and therefore, c/d ≧ 1.4 is made the judgment standard for wedge-shaped vertebrae. In addition to the so called wedge-shaped vertebrae, the wedge-shaped vertebrae (type I) is inclusive of deformations with shortened front brim length, like a wedge-shaped vertebrae, as a result of, for example, an upper brim pressure fracture, upper brim depressed fracture, lower brim pressure fracture, lower brim depressed fracture, and as described below, even in the case of a vertebral body which is judged as a flat vertebrae because a < a - 2σ a , c < c - 2σ c , d < d - 1.5σ d , if there is a remarkable deformation at the front brim portion satisfying c/d ≧ 1.4 , it can be judged to be a wedge-shaped vertebrae. (C) When c/d ≦ 0.7 , the judgment can be an "inverse wedge-shaped vertebrae" (type IV). The inverse wedge-shaped vertebrae may be defined as a vertebral body wherein deformation has occurred at the rear brim portion, and thus the rear brim length (c) has become smaller than the front brim length (d), contrary to the wedge-shaped vertebrae wherein deformation has occurred at the front brim portion, and the front length (d) has become smaller than the rear brim length (c) (see Fig. 3), but in practice, such a vertebral body does not substantially exist. In the fifth lumbar vertebrae, as shown in Example 1, the rear brim length (c) is smaller than the front brim length (d), with the average value of c/d being lower than 1, e.g., 0.95, but in some cases, an inverse wedge-shaped vertebrae cannot be clearly recognized even when the rear brim length (c) is 3/4 of the front brim length (d) ( c/d = 0.75 ), and therefore, a judgment of an inverse wedge-shaped vertebrae (type IV) can be made when c/d ≦ 0.7 . (D) When c , d , c/d satisfy the conditions shown below, the judgment can be a "flat vertebrae" (type III): (a) c < c - 2σ c ; (b) d < d - 1.5σ d ; (c) 0.7 < c/d < 1.4 . In the first condition (a) and the second condition (b), the rear brim length (c) and the front brim length (d) are both smaller than c - 2σ c and d - 1.5σ d when the average values of the vertebral body without deformation are made c and d ; namely, c < c - 2σ c and d < d - 1.5σ d . The flat vertebrae is a vertebral body having a relatively uniform deformation at the front brim portion, the central portion, and the rear brim portion under pressure (see Fig. 3); namely, all of the front brim length (d), the central length (a) and the rear brim length (c) are smaller. As shown by the average values of the lengths of the respective portions of the vertebral body without deformation of Example 1, since the central length (a) is smaller than the front brim length (d) or the rear brim length (c) even in the vertebral body without deformation, the specific feature of the flat vertebrae resides in a particularly small front brim length (d) and rear brim length (c), and the vertebral body wherein c < c - 2σ c and d < d - 1.5σ d is judged to be a flat vertebrae (type III). The conditions of -2σ c for c and 1.5σ d for d are set because c is greater than d ( c > d ) as shown in Example 1, and thus c and d become approximately equal values by setting such conditions, to satisfy the conditions for a flat vertebrae. The second condition (c) is 0.7 < c/d < 1.4 . For, even where a < a - 2σ a , c < c - 2σ c , d < d - 1.5σ d , deformation at the front brim portion is particularly marked, and where c/d ≧ 1.4 , it can be judged to be a wedge-shaped vertebrae, while in the case of a marked deformation at the rear brim portion when c/d ≦ 0.7 , it is judged to be an inverse wedge-shaped vertebrae. (E) When c/d, c , d , a/c, a/d satisfy the following conditions, the judgment can be a "fish vertebrae" (type II): (a) 0.7 < c/d < 1.4 ; (b) at least one of c ≧ c - 2σ c and d ≧ d - 1.5σ d is satisfied; (c) a/c ≦ 0.8 and a/d ≦ 0.8 . "Fish vertebrae" is a vertebrae where a depressed fracture or a pressure fracture has occurred at the central portion, whereby the central length (a) is particularly smaller than the front brim length (d) and rear brim length (c) (see Fig. 3). Accordingly, among the so called "no deformation" classifications, except for wedge-shaped vertebrae, flat vertebrae, and inverse wedge-shaped vertebrae, the central length (a) is particularly smaller. When the central length (a) is smaller than the front brim length (d) and the rear brim length (c), thus satisfying a/c < 0.9 and a/d < 0.9 , deformation such as fish vertebrae can be clearly recognized, but in some cases a clear fish vertebrae cannot be recognized, and therefore, a/c < 0.8 and a/c < 0.8 can be made the judgment standards for fish vertebrae. In addition to the so called fish vertebrae, the fish vertebrae (type II) is also inclusive of deformations with a shortened central length like fish vertebrae, as a result of an upper brim pressure fracture, upper brim depressed fracture, lower brim pressure fracture, lower brim depressed fracture, etc. The relationships between these judgment standards for "no deformation", "wedge-shaped vertebrae", "inverse wedge-shaped vertebrae", "flat vertebrae" and "fish vertebrae", are shown in Fig. 2. In the judgment as described above, explanation has been made by selecting the upper limit as 1.4 and the lower limit as 0.7 as the judgment standard for c/d as a preferable example, but these values can be selected freely, for example, in the case of the upper limit, from the range of 1.25 to 1.55, preferably from 1.33 to 1.5, more preferably from 1.4 to 1.45. The lower limit can be selected freely from the range of 0.8 to 0.6, preferably from 0.75 to 0.65. As the judgment standard for a/c and a/d, 0.8 was selected in the case of "no deformation" and "fish vertebrae", but this value can be freely selected from the range of 0.65 to 0.9, preferably from 0.7 to 0.85, more preferably from 0.75 to 0.8. As the judgment standard for c and d , c - 2σ c and d - 1.5σ d were selected in the case of "no deformation", "flat vertebrae", and "fish vertebrae" but the judgment standard for c can be freely selected from the range of c - 1.0σ c to c - 2.5σ c , preferably c - 1.25σ c to c - 2.25σ c , more preferably c - 1.5σ c to c - 2.0σ c . The judgment standard for d also can be similarly selected from the range of d - 1.0σ d to d - 2.5σ d , preferably d - 1.25σ d to d - 2.25σ d , more preferably d - 1.5σ d to d - 2.0σ d . As mentioned above, although the coefficients β₁ to β₆ can be previously set in the present invention, β₁ to β₆ may optionally be determined from a discriminating function capable of discriminating two groups with two variants or two groups with one variant by using values obtained by measuring two of a , b , c , d , and e with respect to at least five profile X-ray image in each of two types among four types (i.e., wedge-shaped vertabrae, fish vertabrae flat vertebral, and no deformation). As typical axamples, β₁ is z₅ obtained by the following discriminant function using the ratio c/d derived from c and d, which are obtained by measuring at least 5 profile X-ray images in each of two types selected from "wedge-shaped vertebrae" and "no deformation": (i) When the standard deviation σ₁ of c/d of at least 5 profileX-ray image of "wedge-type vertebrae" is substantially the same asthe standard deviation σ₂ of c/d of at least 5 profile X-ray images of "no deformation". wherein µ₁ and µ₂ are averages of c/d relating to the profile X-ray images of "wedge-shaped vertebrae" and "no deformation", respectively, π₁ and π₂ are numbers of the X-ray images of "wedge-type vertebrae" and "no deformation", respectively, µ is an average of µ ₁ and µ ₂ , and σ is an average of σ₁ and σ₂ ; and (ii) When σ₁ and σ₂ are substantially different, wherein µ₁ and µ₂ are average of c/d relating to the profile X-ray images of "wedge-shaped vertebrae" and "no deformation", respectively, π₁ and π₂ are numbers of the X-ray images of "wedge-type vertebrae" and "no deformation", respectively. Furthermore, the present inventors have made an intensive study of a method of objective evaluation of vertebral body deformation, and found that the presence of a vertebral body deformation as well as the deformation type can be judged objectively by use of a discriminant function by accurately measuring the lengths of, for example, the front brim (d), rear brim (c), central portion (a), and the vertebral body width (b) from a profile X-ray image of a vertebral body, and further found that the progress over a period of time of the vertebral body deformation can be judged from the change in the type of said deformation as well as the change in the lengths of the vertebral body, and thus accomplished the present invention. In the prior art, vertebral body deformation has been judged by a physician by visual observation of the profile X-ray image of a vertebral body, but according to the method of the present invention, the lengths of the respective portions of a vertebral body are measured and for the vertebral bodies judged by a physician to be "no deformation", "wedge-shaped vertebrae", "fish vertebrae", and "flat vertebrae", discrimination is made by discriminant functions to determine discriminant formulae for the respective types of vertebral body deformations, and thereafter, the presence of vertebral body deformation and the type of deformation are objectively judged by this discriminant formulae for a patient whose vertebral body deformation is to be judged. The types of deformed vertebral bodies are classified in the present invention generally into wedge-shaped vertebrae, fish vertebrae, and flat vertebrae. In the wedge-shaped vertebrae (type I), deformation has occurred at the front brim portion as shown in Fig. 4, and the front brim length has become particularly smaller than the rear brim length, including in addition to the so called wedge-shaped vertebrae deformation with a shortened front brim length like wedge-shaped vertebrae as the result of upper brim pressure fracture, an upper brim depressed fracture, lower brim pressure fracture, and lower brim depressed fracture, etc. In the fish vertebrae (type II), the depressed fracture or pressure fracture has occurred at the central portion, whereby the central length has become particularly smaller than the front brim length and rear brim length (see Fig. 4), including in addition to the so called fish vertebrae deformations with a shortened center length like fish vertebrae as the result of upper brim pressure fracture, an upper brim depressed fracture, lower brim pressure fracture, and lower brim depressed fracture, etc. In the flat vertebrae (type III), the front brim portion, central portion, and rear brim portion are deformed relatively uniformly under pressure (see Fig. 4); namely, all of the front rear length, central length, and rear brim length are smaller. In addition, there is the inverse wedge-shaped vertebrae (type IV). This type is defined as a vertebral body wherein deformation has occurred at the rear brim portion, and thus the rear brim length is smaller than the front brim length, contrary to the wedge-shaped vertebrae where deformation has occurred at the front brim portion, and the front brim portion is smaller than the rear brim length (see Fig. 4). In practice, however, such a vertebrae does not substantially exist. Accordingly, it is sufficient if the four types, wedge-shaped vertebrae, fish vertebrae, flat vertebrae and "no deformation", can be classified. The judgment method of the second aspect of the present invention is now described in more detail. (i) First, the vertebral body, central length (a), vertebral body width (b), rear brim length (c), and front brim length (d) of each of the vertebrae judged by a physician as "no deformation", wedge-shaped vertebrae, fish vertebrae and flat vertebrae, are measured from a profile X-ray of the vertebral body. The above measurement is used to determine the discriminant formula for discriminating the presence of vertebral body deformation as well as the deformation type. Accordingly, the vertebral bodies of these respective types should be preferably vertebral bodies of persons generally 50 to 75 years old, in view of the age of the patient whose vertebral body deformation is to be judged. Also, it is preferable to perform measurements separately for men and women, and it is desirable to perform measurements of about 50 or more vertebral bodies for both men and women. However, since wedge-shaped, vertebrae, fish vertebrae and flat vertebrae do not widely occur, the discriminant formula for discrimination of the vertebral body deformation type according to the method of the present invention can be determined if at least 5 each of these deformed vertebral bodies can be measured. Also, from the first to twelfth thoracic vertebrae and the first to the fifth lumbar vertebrae, the size differs for each vertebral body, and therefore, it is necessary to measure each vertebral body and determine each deformation type. When making a profile X-ray image of a vertebral body, the thoracic vertebrae from the first thoracic vertebrae to the twelfth thoracic vertebrae, lumbar vertebrae from the first lumbar vertebrae to the fifth lumbar vertebrae should be included, giving a total of 17 vertebral bodies, and therefore, for example, it is advantageous to take the photograph of the profile X-ray image with the eighth thoracic vertebrae as the center separately from the profile X-ray image with the third lumbar vertebrae as the center. In the profile X-ray image with the eighth thoracic vertebrae as the center, the first thoracic vertebrae cannot be accurately measured in most cases, and further, the second and the third thoracic vertebrae are unclear in many cases. However, since the frequency of vertebral body deformation is not great in these vertebral bodies, it is sufficient if measurement can be made from the fourth thoracic vertebrae or the fifth thoracic vertebrae. The center length (a), width (b), rear brim length (c), and front brim length (d) of a vertebral body are specifically as shown in Fig. 4. For determination of these values, from the profile X-ray image of a vertebral body, the central length (a), width (b), rear brim length (c), and front brim length (d) may be measured for each vertebral body with a scale or slide calipers, but the eight points of 1 to 8 shown in Fig. 4 can be input through a digitizer into a computer to measure the central length (a) (3-4), width (b) (7-8), rear brim length (c) (5-6), and front brim length (d) (1-2). Also, for example, the degree of blackening may be recorded with a TV camera, etc., and the lengths of the respective portions of the vertebral body may be automatically measured by image processing. (ii) Next, discriminant formula for discriminating "no deformation", wedge-shaped vertebrae, fish vertebrae, and flat vertebrae is determined. Such a discriminating analysis is described in, for example, Chuichi Okuno et al "Multivariate Analytical Method" (published by Nikkagiren); Chuichi Okuno et al "Overall Multivariate Analytical Method" (published by Nikkagiren); Koichi Sugiyama "Introduction to Multivariate Data Analysis" (published by Asakura Shoten), and others, and calculation can be made by referring to these literatures. However, by inputting the central length (a), vertebral body width (b), rear brim length (c), and front brim length (d) of the respective deformation types determined in (i) by a commercially available computer program (e.g., BMD multivariate analysis program from IBM Co.), the discriminant formula for discriminating "no deformation", wedge-shaped vertebrae, fish vertebrae, and flat vertebrae can be simply obtained. Such discriment formulae are specifically described as follows. That is, the following formulae and the vertebral body is judged from the maximum value of z₁ to z₄: z₁ = a₁₀ + a₁₁x₁ + a₁₂x₂ + a₁₃x₃ + a₁₄x₄ z₂ = a₂₀ + a₂₁x₁ + a₂₂x₂ + a₂₃x₃ + a₂₄x₄ z₃ = a₃₀ + a₃₁x₁ + a₃₂x₂ + a₃₃x₃ + a₃₄x₄ z₄ = a₄₀ + a₄₁x₁ + a₄₂x₂ + a₄₃x₃ + a₄₄x₄ wherein z₁ , z₂ , z₃ , and z₄ are discriminant values of "wedge-shaped vertebrae", "fish vertebrae", "flat vertebrae", and "no deformation", respectively, x₁ , x₂ , x₃ , and x₄ are variables corresponding to a central length, vertebral body width, rear brim length, and front brim length, respectively, a₁₀ , a₂₀ , a₃₀ , and a₄₀ are constants, and a₁₁ to a₄₄ are coefficients. This discriminant formula must be determined for each vertebral body from the first thoracic vertebrae to the fifth lumbar vertebrae, generally from the fourth thoracic vertebrae to the fifth lumbar vertebrae. (iii) Next, from the profile X-ray image of the vertebral body of the patient whose vertebral body deformation is to be judged, the central length (a), vertebral body width (b), rear brim length (c), and front brim length (d) are measured for each vertebral body as described in (i), and by substituting a , b , c and d in the discriminant formula determined in (ii), the discriminant values for "no deformation", wedge-shape vertebrae, fish vertebrae, and flat vertebrae are determined respectively, and the type which has the highest numerical value is judged to be the vertebral body deformation type. For example, when the vertebral body deformation type of the third lumbar vertebrae of the patient is to be judged, the respective measured values of a , b , c , and d are substituted in the discriminant formula of the third lumbar vertebrae determined in (ii) as a matter of course. Having described above the method for judging, the present invention is not limited to the above classification of the four vertebral deformation types of "no deformation", wedge-shaped vertebrae, fish vertebrae and flat vertebrae from the four measured values of central length (a), vertebral body width (b), rear brim length (c), and front brim length (d), but it is also possible to classify the vertebral body deformation types by use of the discriminant function as described below. (A) Method for classifying the three vertebral body deformation types from three or four measured values. For example, in the first to the eleventh thoracic vertebrae, when fish vertebrae does not substantially exist, the case occurs wherein classification is required for only the three types of "no deformation", wedge-shaped vertebrae, and flat vertebrae, or the case when sufficient measured values of a , b , c , and d of fish vertebrae cannot be obtained for determining the discriminant formula even if a classification of the four types including fish vertebrae is desired. In such a case, it is possible to determine the discriminant formula for classifying the three types of "no deformation", wedge-shaped vertebrae, and flat vertebrae from the four measured values of the central length (a), vertebral body width (b), rear brim length (c), and front brim length (d), or from three of these four measured values, for example, three measured values of a , c , and d . (B) Method for classifying two vertebral body deformation types from two to four measured values. Fish vertebrae and flat vertebrae can be judged relatively easily even when visually observed, and rarely occur, and therefore, in some cases they are difficult to judge with the use of a discriminant formula. Accordingly, when it is desired to classify only "no deformation" and wedge-shaped vertebrae by use of a discriminant formula, three, for example, a , c and d , or two, for example, measured values of c and d , can be used to determine the discriminant formula. Also, in such a case, it is possible to determine the discriminant value for classifying "no deformation" and wedge-shaped vertebrae by using of a single variate discriminant function of c/d. Specifically speaking, the discriminant function is represented by the following formula in which the ratio (c/d) of c and d measured from at least 5 profile X-ray images in each of two types of :wedge-shaped vertebrae" and "no deformation" and the ratio (c/d) of less than z₅ obtained from the vertebral body to be judged is judged as "no deformation" and the ratio (c/d) equal to or larger than z₅ is judged as "wedge-shaped vertebrae". (i) When the standard deviation σ₁ of c/d of at least 5 profile X-ray image of "wedge-type vertebrae" is substantially the same asthe standard deviation σ₂ of c/d of at least 5 profile X-ray images of "no deformation". wherein µ₁ and µ₂ are averages of c/d relating to the profile X-ray images of "wedge-shaped vertebrae" and "no deformation", respectively, π₁ and π₂ are numbers of the X-ray images of "wedge-type vertebrae" and "no deformation", respectively, µ is an average of µ ₁ and µ ₂ , and σ is an average of σ₁ and σ₂ ; and (ii) When σ₁ and σ₂ are substantially different, wherein µ₁ and µ₂ are averages of c/d relating to the profile X-ray images of "wedge-shaped vertebrae" and "no deformation", respectively, π₁ and π₂ are numbers of the X-ray images of "wedge-type vertebrae" and "no deformation", respectively. The apparatus for judging the deformation of vertebral body according to the present invention is used for the practice of the above-mentioned method of the present invention. This apparatus is basically composed of an imputting means for at least two indices of a , b , c , and d , an arithmetic means for effecting necessary operations for the descrimination, a discriminating means by using the calculation results, a means for inputting the discriminant function and/or standard necessary for the calculation and discrimination, and an output means for outputting the judgment result. As the inputting means, a rule means, a digitizer means, an image processing means in which an image inputting means such as a TV camera is used. In the preferred embodiment of the present apparatus, c/d, a/c, and a/d are determined by the arithmetic means, and c , d , σ c , σ d , and β₁ - β₆ are inputted as discriminating standards by the discriminating standard inputting means, and the discriminating means can discriminate the results by these standards. Furthermore, the present apparatus can be further provided with a means for obtained the discriminating standards from a number of the profile X-ray images having no substantial deformation. Furthermore, the present apparatus can be optionally provided with an arithmetic means capable of discriminating two groups with two variants, which are used for determining the indices β₁ to β₆ , or of discriminating two groups with one variant. In this case, the data of at least two types of the four types of the vertebral body, i.e., "wedge-shaped vertebrae", "fish vertebrae", "flat vertebrae", and "no deformation" should be input. I another preferred embodiment of the present apparatus, the arithmetic means is connected to a discriminating standard input means for inputting the discriminant function, whereby the arithmetic means calculates the discriminating values with said discriminant function and the discriminating means determine the type of the vertebral body based on the maximum discriminating value. In this case, the present apparatus is provided with a function capable of discriminating the four types of the vertebral body, "wedge-shaped vertebrae", "fish vertebrae", "flat vertebrae", and "no deformation" from a one-dimensional formula containing four variables a , b , c , and d . Furthermore, there is also the apparatus which is capable of judging the two types of the vertebral body with one variant. According to the first aspect of the judgment method and the apparatus of the present invention, the type of vertebral body deformation can be objectively evaluated, and also the change in the said deformation type and the progress over a period of time of the vertebral body deformation can be judged. Also, the present method and apparatus is very useful for determining the progress of bone disease such as osteoporosis, etc., and for confirmation of the effect of therapy. According to the second aspect of the judgment method and apparatus of the present invention, the type of vertebral body deformation can be objectively evaluated, and further the change in the type of deformation type and the progress over a period of time of the vertebral body deformation can be judged. Also, the present method and apparatus is very useful for determining the progress of bone disease such as osteoporosis, etc., and for confirmation of the effect of therapy. Examples The present invention will now be further illustrated by, but is by no means limited to, the following Examples. Example 1 From the profile X-ray images of thoracic and lumbar vertebrae, with the eighth thoracic vertebrae and the third lumbar vertebrae as the center, of 50 to 75 years old women, the front brim length (d), central length (a), rear brim length (c), and vertebral body width (b) were measured for each vertebral body, and for the vertebral bodies classified by a physician as "no deformation", the average values ( x ) and standard deviations (σ) for the respective vertebral bodies of the measured values of a , b , c , and d, and the calculated values of c/d, are as shown in Table 1. Example 2 From the profile X-ray image of the thoracic and lumber vertebrae of a 70 years old osteoporosis patient (woman), the center length (a), vertebral body width (b), rear brim length (c), and front brim length (d) were measured for each vertebral body, and c/d, a/c and a/d were calculated. The results of the judgment of the presence of deformation of each vertebral body, and the type of change, are shown in Table 2. The judgment was made according to the judgment standard shown in Table 2. The same standards were also used in Examples 3 and 4. Example 3 From the profile X-ray image, with the third lumbar vertebrae as the center, of ten osteoporosis patients (women) 59 to 76 years old, the central length (a), vertebral body width (b), rear brim length (c), and front brim length (d) of the third lumber vertebrae were measured and c/d, a/c, and a/d were calculated. The results of the judgment of the presence of deformation, and the type of deformation, are shown in Table 3. Example 4 From the profile X-ray image, with the eighth thoracic vertebrae and the third lumbar vertebrae as the center, of a 67 years old osteoporosis patient (woman), the center length (a), vertebral body width (b), rear brim length (c), and front brim length (d) from the sixth thoracic vertebrae to the fifth lumbar vertebrae were measured, and c/d, a/c, and a/d were calculated to judge the presence of deformation and the type of deformation. Six months later, the profile X-ray image of the same site was photographed again to judge the presence of deformation and the type of deformation according to the same method, and at the same time, the change in the ratios of a , b , c , and d were calculated. The results are shown in Table 4. The second lumbar vertebrae changed from a "no deformation" classification to a type III (flat vertebrae) within six months. Although not showing a clear deformation, the conditions of the front brim and rear brim of the seventh thoracic vertebrae, and the center of the tenth thoracic vertebrae, tend to have worsened. Also, the conditions of the centers of the tenth thoracic vertebrae and the third lumbar vertebrae had further worsened. On the other hand, substantially no change has occurred in the sixth, eighth, ninth, eleventh, and twelfth thoracic vertebrae, or in the fourth and fifth lumbar vertebrae. Example 5 For 186 vertebral bodies, including 6 vertebral bodies judged as wedge-shaped vertebrae, 25 vertebral bodies judged as fish vertebrae, eight vertebral bodies judged as flat vertebrae and 147 vertebral bodies judged as "no deformation" according to visual observation by a physician of a profile X-ray image, with the third lumbar vertebrae as the center, of osteoporosis patients (women), eight points of 1 to 8 shown in Fig. 4 were input through a digitizer into a computer to measure the central length (a), vertebral body width (b), rear brim length (c), and front brim length (d), and calculation was performed by a computer by inputting the above a , b , c , and d into a multivariate discriminant analytical program in the BMD multivariate analytical program of IBM Co. to obtain the coefficients and constants of the discriminant formula as shown in Table 5. The discriminant value for each vertebral body was calculated by substituting the measured values of a , b , c , and d , into the above four discriminant formulae, and the type having the maximum numerical value was judged to be the vertebral body deformation type. The coincidence of the vertebral body deformation type thus judged according to the discriminant formula with the vertebral body deformation type according to the physician's judgment is shown in Table 6, indicating a good coincidence between both judgments. The coincidence ratio of the respective types is 94.6%. Example 6 From a profile X-ray image, with the third lumbar vertebrae as the center, of 10 osteoporosis patients (women), the central length (a), vertebral body width (b), rear brim length (c), and front brim length (d) of the third lumbar vertebrae were measured and substituted in the discriminant formula determined in Example 5 (coefficients and constants are shown in Table 5) to calculate the respective discriminant values, and the type having the maximum numerical value was judged to be the vertebral body deformation type. The results are shown in Table 7. Example 7 From vertebral body profile X-ray images with the eighth thoracic vertebrae and the third thoracic vertebrae as the centers thereof, respectively, according to a physician's judgment, as shown in Table 8, from the sixth thoracic vertebrae to the third thoracic vertebrae, 257 vertebral bodies of wedge-shaped vertebrae and 2,157 vertebral bodies of "no deformation" were selected, and the central length (a), vertebral width (b), rear brim length (c), and front brim length (d) of each vertebral body were measured as in Example 5, and in addition, the c/d was calculated. For discriminating wedge-shaped vertebrae from "no deformation" with one variate of c/d, the value of c/d for each vertebral body was input to the one variate discriminant analytical program in the BMD multivariate analytical program of IBM Co. to obtain the discriminant values shown in Table 8. The coincidence with the physicians judgment when discrimination was made between wedge-shaped vertebrae-and "no deformation" by use of these discriminant values is shown in Table 8. For all vertebral bodies, the coincidence ratio of "no deformation" was 98.15%, and the coincidence ratio of wedge-shaped vertebrae 87.16%, thus giving good results. Example 8 From the profile X-ray images, with the eighth thoracic vertebrae as the center, of 10 osteoporosis patients (women), as described in Example 5, the center length (a), vertebral body width (b), rear brim length (c), and front brim length (d) of the eighth thoracic vertebrae were measured, and in addition, the c/d was calculated. The results are shown in Table 9. Since these vertebral bodies were neither clearly fish vertebrae nor flat vertebrae, a discrimination between wedge-shaped vertebrae and "no deformation" by the c/d determined in Example 7 was conducted. The results are shown in the right hand column in Table 9. No. Vertebrae body measured values (mm) Deformation type a b c d c/d 1 20.0 22.5 24.7 16.5 1.50 Wedge-shaped vertebrae 2 21.6 30.2 24.4 21.3 1.14 No deformation 3 17.5 32.9 21.7 13.6 1.60 Wedge-shaped vertebrae 4 17.7 38.3 29.6 15.6 1.90 " 5 22.0 28.8 23.8 19.2 1.24 No deformation 6 12.9 38.2 21.8 13.2 1.64 Wedge-shaped vertebrae 7 21.0 28.3 24.5 21.2 1.15 No deformation 8 18.9 32.0 22.4 18.5 1.21 " 9 19.8 30.8 24.3 18.7 1.30 " 10 19.3 31.3 22.2 15.6 1.43 Wedge-shaped vertebrae

Claims: A method for judging deformation of a vertebral body in a subject, comprising the steps of: (i) measuring a central length a, a rear brim length c and a front brim length d from a plurality of profile X-ray images obtained from vertebral bodies having substantially no deformation at a portion corresponding to a portion of said vertebral body to be judged, storing the measurements in a memory, and determining values of c/d, a/c and a/d, the average values c and d of the rear brim lengths c and front brim lengths d, respectively, and the standard deviations σ c and σ d of the values c and d, respectively: (ii) measuring a, c, and d from a plurality of profile X-ray images obtained from vertebral bodies having at least one deformation, storing the measurements in said memory, and determining values of c/d, a/c, a/d, and c , d , σ c , and σ d ; (iii) preparing a judgement standard from the data obtained in steps (i) and (ii) for determining the nature of the deformity according to a classification consisting of "no deformation" vertebrae, "wedge-shaped" vertebrae, "inverse wedge-shaped" vertebrae, "fish" vertebrae, and "flat" vertebrae; (iv) measuring a, c, and d from a profile X-ray image of said vertebral body to be judged, storing the measurements in said memory, and determining c/d, a/c, and a/d; and (v) classifying deformation of the vertebral body to be judged from c/d in the case of "wedge-shaped" and "inverse wedge-shaped" vertebrae, from c , d , σ c , σ d , and c/d in the case of "flat" vertebrae, and from c , d , σ c , σ d , a/c, a/d, and c/d in the case of the "no deformation" and "fish" vertebrae, in order to judge deformation and to monitor the effect of therapy on said vertebral body to be judged. A method according to claim 1, wherein said step (v) is carried out according to the following sub-steps: (i) c/d ≧ β₁ ... wedge-shaped vertebrae (or anterior wedge fracture) (ii) c/d ≦ β₂ ... inverse wedge-shaped vertebrae (or inverse anterior wedge fracture) (iii) c < c - β₃σ c , d < d - β₄σ d , and β₂ < c/d < β₁ ... flat vertebrae (or compression fracture) (iv) β₂ < c/d < β₁ , a/c ≦ β₅ , a/d ≦ β₆ and c ≧ c - β₃σc and/or d ≧ d - β₄σ d ... fish vertebrae (or central biconcave fracture) (v) β₂ < c/d β₁ , c ≧ c - β₃σ c and/or d ≧ d - β₄σ d , and a/c > β₅ and/or a/d > β₆ ... No deformation (or no fracture), wherein β₁ to β₆ are selected as follows: 1.25 ≦ β₁ ≦ 1.55, 0.6 ≦ β₂ ≦ 0.8, 1 ≦ β₃ ≦ 2.5, 1 ≦ β₄ ≦ 2.5, 0.65 ≦ β₅ ≦ 0.9, and 0.65 ≦ β₆ ≦ 0.9. A method as claimed in claim 2, wherein β₁ to β₆ are 1.33 ≦ β₁ ≦ 1.5, 0.65 ≦ β₂ ≦ 0.75, 1.25 ≦ β₃ ≦ 2.25, 1.25 ≦ β₄ ≦ 2.25, 0.7 ≦ β₅ ≦ 0.85, and 0.7 ≦ β₆ ≦ 0.85. A method as claimed in claim 3, wherein β₁ is 1.4, β₂ 0.7, β₃ 2, β₄ 1.5, β₅ 0.8, and β₆ 0.8. A method as claimed in any of claims 2 to 4, wherein at least one of β₁ to β₆ are determined by a discriminant function capable of discriminating two groups with one or two variables by using two of a, b, c, and d or the ratio thereof, which are obtained by measuring at least 5 profile X-ray images in each of two types selected from four types, "wedge-shaped vertebrae", "fish vertebrae", "flat vertebrae", and "no deformation". A method as claimed in any preceding claim, wherein the measurement is carried out by a rule means, a digitizer means, or a means for automatically measuring the lengths from image treatment by memorizing the X-ray image with an X-ray cage automatic input means. Apparatus for judging deformation of a vertebral body according to the method of claim 1, the apparatus comprising (i) an input means for inputting values of at least central length a, rear brim length c, and front brim length d from a profile X-ray image of a vertebral body to be judged, (ii) an arithmetic means for performing the calculation of c/d, a/c, a/d, the average values c , d of c and d and the standard deviations σ c and σ d of the values c and d, respectively, necessary for determination of the nature of the deformity according to a classification consisting of "wedge-shaped vertebrae", "inverse wedge-shaped vertebrae", "fish vertebrae", "flat vertebrae", and "no deformation" by using said input values, (iii) a means for discriminating "wedge-shaped" and "inverse wedge-shaped" vertebrae from c/d, "flat" vertebrae from c , d , σ c , σ d , and c/d, and "no deformation" and "fish" vertebrae from c , d , σ c , σ d , a/c, a/d, and c/d, (iv) a means for inputting the judgement standard necessary for the said calculation and discrimination, and (v) an output means for outputting the judgement result. Use of the apparatus according to claim 7, for carrying out the method of any one of claims 2-6.

Inventors: AOKI CHOJU, UOTANI YASUHIRO, YAMASHITA GENTARO, AOKI, CHOJU, UOTANI, YASUHIRO, YAMASHITA, GENTARO
Applicants: TEIJIN LTD, TEIJIN LIMITED
Application year: 1987
Main classification: A61B5
IPC codes: G06F19