# Developing of a Mathematical Model to Perform Measurements of Axial Vertebral Rotation on Computer-Aided and Automated Diagnosis Systems, Using Raimondi’s Method

Authors: José Hurtado-Aviles; Joaquín Roca-González; Konstantsin Sergeevich Kurochka; Jose Manuel Sanz-Mengibar; Fernando Santonja-Medina
Journal: Radiology Research and Practice
 (2021)

---

## Abstract

Introduction. Axial vertebral rotation (AVR) is a basic parameter in the study of idiopathic scoliosis and on physical two-dimensional images. Raimondi’s tables are the most used method in the quantification of AVR. The development of computing technologies has enabled the creation of computer-aided or automated diagnosis systems (CADx) with which measurement on medical images can be carried out more quickly, simply, and with less intra and interobserver variabilities than manual methods. Although there are several publications dealing with the measurement of AVR in CADx systems, none of them provides information on the equation or algorithm used for the measurement applying Raimondi’s method. Goal. The aim of this work is to perform a mathematical modelling of the data contained in Raimondi’s tables that enable the Raimondi method to be used in digital medical images more precisely and in a more exact manner. Methods. Data from Raimondi’s tables were tabulated on a first step. After this, each column of Raimondi’s tables containing values corresponding to vertebral body width (D) were adjusted to a curve determined by AVR = f (d). Third, representative values of each rotation divided by D were obtained through the equation of each column D. In a fourth step, a regression line was fitted to the data in each row, and from its equation, the mean value of the D/d distribution is calculated (value corresponding to the central column, D = 45). Finally, a curve was adjusted to the obtained data using the least squares method. Summary and Conclusion. Our mathematical equation allows the Raimondi method to be used in digital images of any format in a more accurate and simplified approach. This equation can be easily and freely implemented in any CADx system to quantify AVR, providing a more precise use of Raimondi’s method, as well as being used in traditional manual measurement as it is performed with Raimondi tables.


---

## Body

## 1. Introduction

Axial vertebral rotation (AVR) is a basic parameter in the study of idiopathic scoliosis in adolescents [1, 2]. Its measurement is necessary to evaluate the severity of scoliosis and quantify the risk of progression [1, 3–7] for the selection of treatment [3, 4, 8, 9] and for the analysis of orthopaedic surgery [1, 6–8, 10, 11].Two-dimensional medical images, especially anterior-posterior and lateral X-rays of the full spine when standing [12–14], continue to be the method of choice for diagnosing and monitoring scoliosis [15].On physical two-dimensional images, Raimondi’s tables [16–18] are the most used method in the quantification of AVR [19–21]. This method possesses a more than acceptable degree of reproducibility between observers, and it is easy to apply [19, 21, 22].Figure1 shows a schematic description of the anatomical references used to measure AVR in two-dimensional images.Figure 1
A schematic description of the anatomical references used to measure AVR in two-dimensional images.Based on the position of the two closest lateral faces vertebral body points and the two opposite sides of the shadow of the pedicle turned towards the centre of the vertebra in the anteroposterior projection, the software calculates the width of the vertebral body (D) and the distance from the centre of the pedicle to the side of the vertebral body (d). Adapted by permission from Nature/Springer, European Spine Journal, Vrtovec T, Pernus F, and Likar B, a review of methods was conducted for quantitative evaluation of axial vertebral rotation [7].The development of computing technologies has enabled the creation of computer-aided or automated diagnosis systems (CADx) with which manual measurement methods on medical images (for instance, Cobb angle, distances, and wedgings) can be applied more quickly, simply, and with less intra and interobserver variabilities [23–28].The digital X-rays require a software for the clinicians to state a precise diagnosis and treatment, as well as quantify the follow-up changes. These kinds of software are usually not designed to quantify specific AIS parameters, or they are not easily accessible due to their complex use or higher cost.Although there are several publications dealing with the measurement of AVR in CADx systems, none of them provides information on the equation or algorithm used for the measurement applying Raimondi’s method [3, 28–30].The aim of this work is to perform a mathematical modelling based on the data contained in Raimondi’s tables that enables them to be applied on CADx systems, allowing the Raimondi method to be used in digital medical images more precisely and in a more exact manner.Data from Raimondi’s tables were tabulated on a first step. After this, each column of Raimondi’s tables containing values corresponding to vertebral body width (D) were adjusted to a curve determined by AVR = f (d). In this mathematical equation, “d” represents each value of the column D. Third, representative values of each rotation divided by D were obtained through the equation of each column D. A representative value for each AVR was obtained this way. The fourth step was to adjust the data included in each row to a regression line, and this equation was used to calculate the average value of the distribution D/d (middle column value = 45). Finally, a curve was adjusted to the obtained data using the least squares method.The novelties of our work are (a) to provide an adimentional equation to describe in a theoretically way the empirical data tabulated in Raimondi’s tables with more accuracy and less precision error; (b) our mathematical equation allows the Raimondi method to be used in CADx to quantify AVR in medical images without the need for a scale; and (c) this equation can be easily and freely implemented in any CADx system, providing a more precise use of Raimondi’s method than on the traditional printed X-ray.


## 2. Materials and Methods

In Raimondi’s tables, the data pertaining to one same column correspond to a particular vertebral diameter,D, between 20 mm and 70 mm (Table 1). This value is given by the width of the vertebral body D projected on the X-ray image.Table 1
Simplified Raimondi’s tables with data to obtain AVR. Modification of Raimondi. The quantification of the vertebral rotation in the scoliosis. Comparing experiences [32].
2.32.42.52.62.7…7.8AVR = 2°2.62.72.82.93…9AVR = 4°2.933.13.23.4…9.8AVR = 6°3.23.33.43.63.7…10.5AVR = 8°3.53.63.83.94.1…11.7AVR = 10°3.73.844.24.3…12.8AVR = 12°3.944.24.44.6…13.9AVR = 14°4.14.24.54.74.9…15AVR = 16°4.34.54.74.95.2…16AVR = 18°4.54.755.25.5…17.3AVR = 20°4.95.15.45.75.9…18.8AVR = 22°5.35.55.86.16.4…20.1AVR = 24°5.766.36.66.9…21.7AVR = 26°66.46.777.4…23.1AVR = 28°……………………17.518.319.22020.9…61.5AVR = 60°D = 20D = 21D = 22D = 23D = 24…D = 70Within each column, the data correspond to the distance,d, between the centre of the shadow of the pedicle furthest from the edge of the vertebra and said edge, in millimetres.Each pair of values (D, d) corresponds to a value of AVR between 2° and 60°.From the distribution of empirical data in Raimondi’s tables, a mathematical model that describes the information contained in them in an ideal manner is sought.First, the data in each column are analysed separately. Using the least squares method, the equation of a curve,(1)fd|D=x=AVR,is pursued, such that every d,AVR pair is as close as possible to this curve.Using the “Statistix 10” software by Analytical Software, a rational function in the form(2)y=a+bx1+cx+dx2,where a, b, and c are the constants, was chosen as the best option.By adjusting the function model to each set of values in each column of Raimondi’s tables, the adjustment values “pseudoR2” reflected in Table 2 are obtained.Table 2
Adjustment of data in Raimondi’s tables to the function model.
DR2200.9981210.9987220.9988230.9989240.9991250.9992260.9992270.9994280.9994290.9994300.9996310.9996320.9997330.9996340.9997350.9997360.9997370.9997380.9997390.9997400.9997410.9997420.9997430.9997440.9993450.9997460.9996470.9997480.9997490.9996500.9996510.9996520.9996530.9996540.9996550.9996560.9996570.9996580.9996590.9995600.9995610.9994620.9994630.9994640.9994650.9994660.9994670.9994680.9993690.9993700.9993The quality of the adjustment of each curve improves as the width of the vertebra,D, increases, reaching a maximum of R2 = 0.9997. The explanation for this may lie in the difficulty of performing a good measurement when this refers to very small distances (Figure 2).Figure 2
Adjustment of the data in columnD = 20 of Raimondi’s tables to the curve.From the valueD = 48, R2 decreases slightly until at D = 68, we have R2 = 0.9993.It can also be discerned that forD = 44, a value of R2 outside of the collection of values described exists due to an atypical datum (Figure 3).Figure 3
Adjustment of the data in columnD = 44 of Raimondi’s tables to the curve.Residuals of the data of various columns from Raimondi’s tables regarding their respective regression curves can be observed in Figures4(a) and 4(b).Figure 4
Residuals of the data of various columns from Raimondi’s tables regarding their respective regression curves.
(a)(b)When using the Raimondi’s tables on CADx systems, it is necessary for the vertebrae shown in the digital image to be of sufficient size (20≤D≤70) in order to be able to apply the equation [31]. In addition, the dimensions of variables D and d are given in millimetres.All pairs of (D, d) values corresponding to a certain AVR are proportional, meaning they should produce ideally equal D/d values.To avoid having to know the scale of the image to measure AVR, the independent variable will have the formD/d in the calculated equation, such that the value is adimensional.To calculate the equation(3)fD,d=AVR,the data in Raimondi’s tables is replaced by other data obtained in the following manner (Table 3):(i)
For the set of data in each column, the nonlinear regression curve (equation (2)) is obtained(ii)
From the equation of each regression curve, thed values corresponding to the rotations 2°, 4°, 6°, …, 60° are obtained to five decimal places(iii)
The valueD/d is entered to five decimal places in each cell in the tablesTable 3
Raimondi’s tables completed with the results of dividing theD value of each column by the values obtained on the regression curves.
8.154088.270128.362288.429318.46827…8.79462AVR = 2°7.573687.680057.737917.786017.80927…7.89888AVR = 4°7.040257.137397.169037.202687.21409…7.16873AVR = 6°6.548556.6376.648946.671536.67412…6.537AVR = 8°6.094276.174466.171986.186146.18232…5.98232AVR = 10°5.673635.746015.733295.741245.73293…5.49782AVR = 12°5.28345.348475.328885.332345.32102…5.06923AVR = 14°4.920844.979034.955254.955584.94246…4.69047AVR = 16°4.583534.63534.609464.607764.59374…4.34998AVR = 18°4.269394.315174.288954.286084.27189…4.04111AVR = 20°3.976614.016843.991523.988133.97433…3.76408AVR = 22°3.703623.738723.715253.711873.69886…3.51061AVR = 24°3.449043.479453.458493.45553.44357…3.2801AVR = 26°3.211673.237813.21983.217473.20683…3.06784AVR = 28°……………………1.12651.133271.131671.133361.13382…1.11129AVR = 60°20/d21/d22/d23/d24/d…70/dIn Table3, the values of each row are very similar. To find the representative value of each rotation, a regression line is adjusted to the data in each row, and from its equation, the average D/d value of distribution is calculated (value corresponding to the central column, D = 45).Figures5–7 show that as the width of the vertebra increases, the linearity of the point cloud increases. It can also be observed how the general aspect of the data is quite regular, with the exception of those corresponding to D = 60 and D = 70, where the pattern followed loses its regularity. Consequently, each regression line is calculated by omitting some of the values as shown in the figures.Figure 5
Variation of theD/d value corresponding to the AVR of 8° as D varies.Figure 6
Variation of theD/d value corresponding to the AVR of 12° as D varies.Figure 7
Variation of theD/d value corresponding to the AVR of 50° as D varies.Table4 shows the representative D/d values of each AVR calculated.Table 4
RepresentativeD/d values of each AVR.
AVR (°)D/d28.884848.072967.363986.7457106.1983125.7112145.2751164.8817184.5261204.2023223.9068243.6362263.3881283.1603302.9495322.7553342.5753362.4087382.2544402.1112421.9784441.8552461.7409481.6348501.5365521.4452541.3607561.2825581.2102601.1434


## 3. Results

The least squares method is used to adjust a curve to the data of Table4. Using the “Statistix 10” software, a rational function in the form(4)y=a+bx+cx21+dx+ex2,where a, b, c, d, and e are the constants, is chosen. The adjustment obtained with this model is R2 = 1.The equation obtained is(5)AVR=20.22483−330.5077D/d+33.46082D/d21−3.93825D/d−1.322272D/d2.Convergence criterion met after 16 iterations.Residual SS (SSE): 9.722E − 04Residual MS (MSE): 3.889E − 05Standard deviation: 6.236E − 03Degrees of freedom: 25AICc: −294.46PseudoR2: 10000


## 4. Discussion

AVR is a necessary parameter in being able to correctly evaluate and provide a prognosis for monitor and treat scoliosis [1–7].Raimondi’s tables are widely used on physical X-rays due to its notable reliability and ease of use in quantifying AVR [21, 32].We are not aware that until now, a method to calculate ARV by CADx systems has been developed using the Raimondi method.The data contained in the Raimondi’s tables were obtained experimentally from thousands of measurements on educational and anatomical spines [32]. In consequence, Raimondi recognises that the measurements obtained with his method involve systematic errors [19, 32]. These errors are small in comparison with the real value of the measures made, meaning their presence will not vary the decision on the evaluation or treatment of the patient. However, these errors may propagate along with random errors due to a multitude of variables present in the measuring process such as (A) the loss of information when showing a three-dimensional structure on a flat image; (B) the characteristics of the medical image (e.g., existence of noise that hinders the legibility of the image); (C) limitations of the observer (e.g., exactitude in the identification of measurement references points, need for training, and experience); (D) the particular morphological characteristics of the vertebrae (e.g., deformations, the variability of the interpedicular distance between vertebrae); and (E) the position of the patient in the image.On the other hand, medical images today reach the medical professional on a CD, and Raimondi’s tables cannot be applied to these digital images, when they do not contain information about the scale (that is, the image does not have the scale imprinted on it or the DICOM archive entry does not have this information).The obtained equation allows the AVR measurement to be implemented in CADx systems with the Raimondi method, allowing more accurate and precise results to be obtained, due to the elimination of the systematic error contained in the Raimondi tables and the use of auxiliary tools that usually incorporate CADx systems, such as enlargement of regions of interest, application of filters, and automatic detection algorithms.We have not found any previous work that offers an equation or an analogous algorithm to use Raimondi’s method in CADx systems that show it openly. The obtained equation is available to be used freely.The form of the equation obtained also has the advantage that allows the Raimondi method to be used in digital medical images that do not contain information about its scale (for example, images in TIFF format, in DICOM files that do not contain this information or when it is not printed on the same image).The limitations of our study are that the programming code to implement our equation in any software is not provided. Our equation was tested using the C++ programming language and OpenCV library, so we can recommend these tools to implement the obtained equation. Also, our equation does not innovate in Raimondi’s procedure but achieves a linear theoretical model based on its empirical data, improving its validity and reliability. Therefore, the highlight of our study is that this equation can be easily and freely implemented in any CADx system, allowing the use of the Raimondi method in digital medical images, whether or not they contain the scale, in a more exact and precise way.


## 5. Summary and Conclusion

A mathematical equation has been calculated that can be implemented in CADx systems from the mathematical modelling of the empirical data tabulated in Raimondi’s tables.This allows the Raimondi method to be used in digital images of any format more accurately, simply, and quickly, by eliminating the small systematic error existing in the discrete values of the tables and being able to use in conjunction with auxiliary tools in the CADx system such as zooming over regions of interest, algorithms for enhancing medical image structures, adjusting brightness, and contrast.The equation can be easily and freely implemented in any CADx system.



---

*Source: 5523775-2021-02-06.xml*