This video demonstrates the use of a novel graphical tool for measuring the spatially weighted calcium score (SWCS), an alternative to the Agatston score, for quantifying coronary artery calcification. The graphical tool computes SWCS based on image data from gated cardiac computed tomography and user-defined paths of the coronary arteries.
The current standard for measuring coronary artery calcification to determine the extent of atherosclerosis is by calculating the Agatston score from computed tomography (CT). However, the Agatston score disregards pixel values less than 130 Hounsfield Units (HU) and calcium regions less than 1 mm2. Due to this thresholding, the score is not sensitive to small, weakly attenuating regions of calcium deposition and may not detect nascent micro-calcification. A recently proposed metric called the spatially weighted calcium score (SWCS) also utilizes CT but does not include a threshold for HU and does not require elevated signals in contiguous pixels. Thus, the SWCS is sensitive to weakly attenuating, smaller calcium deposits and may improve the measurement of coronary heart disease risk. Currently, the SWCS is underutilized owing to the added computational complexity. To promote translation of the SWCS into clinical research and reliable, repeatable computation of the score, the aim of this study was to develop a semi-automatic graphical tool that calculates both the SWCS and the Agatston score. The program requires gated cardiac CT scans with a calcium hydroxyapatite phantom in the field of view. The phantom allows for deriving a weighting function, from which each pixel's weight is adjusted, allowing for the mitigation of signal variations and variability between scans. With all three anatomical views visible simultaneously, the user traces the course of the four main coronary arteries by placing points or regions of interest. Features such as scroll-to-zoom, double-click to delete, and brightness/contrast adjustment, along with written guidance at every step, make the program user-friendly and easy to use. Once tracing the arteries is complete, the program generates reports, which include the scores and snapshots of any visible calcium. The SWCS may reveal the presence of subclinical disease, which may be used for early intervention and lifestyle changes.
Measuring the amount of calcium within arteries using computed tomography (CT) is an established way to assess the severity of coronary atherosclerosis. Knowing and quantifying the extent of atherosclerosis is key to determining the risk of future coronary heart disease1,2,3,4. The most common way of measuring calcium in the coronary arteries is using the Agatston score5. However, part of the Agatston score calculation relies on the intensity of the chosen pixels, measured in Hounsfield Units (HU). Any pixels less than 130 HU are not accounted for in the calculation. Similarly, calcifications with an area less than 1 mm2 are not considered. Due to these thresholds, the Agatston score is not sensitive to small, weakly attenuating foci of calcification, which may still be important in revealing the presence of subclinical disease6.
A previously described metric called the spatially weighted calcium score (SWCS) was proposed to assess the risk of atherosclerotic plaque in patients with low levels of calcification7. Unlike the Agatston score, the SWCS does not use signal thresholding to reduce the impact of image noise. Instead, it makes use of a phantom-an object with known concentrations of calcium hydroxyapatite (CHA) placed on the participant such that it is in the scan's field of view. Here, a phantom with 0 mg/mL, 50 mg/mL, 100 mg/mL, and 200 mg/mL CHA was used during development; however, in the current implementation of the graphical tool, only the 0 mg/mL and 100 mg/mL sections are required. The phantom is used to create a scan-specific weighting function, which is then used to weigh each of the user-selected pixels as well as its neighbors. Pixels with neighboring pixels that have a high attenuation level are given more weight than ones surrounded by pixels with lower attenuation levels. This process makes the SWCS tolerant to noise and comparable from scan to scan8. The SWCS is continuous and produces a score even when there are low levels of calcification, allowing for quantification of the extent of atherosclerosis when the Agatston score is zero. By allowing the evaluation of micro-calcification even when the Agatston score is zero, the SWCS may be important in revealing the presence of subclinical disease. This may allow a better understanding of the genetic, environmental, and other risk factors in atherosclerosis9,10. A previous study, which examined individuals with an Agatston score of zero at baseline and non-zero at a follow-up approximately 15 years later, observed that those with a higher SWCS at baseline had a higher coronary heart disease (CHD) event rate. The predictive power of the SWCS is especially important in younger populations, where the detection and monitoring of residual risk over a long term may be helpful6.
Presented here is a semi-automatic tool for calculating the SWCS along with the Agatston score. The tool utilizes a graphical user interface running on a compatible programming language. The user is able to interact with the images to generate a final series of reports, which include the two calcium scores. To start, the user selects a case, or a series of Digital Imaging and Communications in Medicine (DICOM) files, to input into the program. These images must be breath-held, electrocardiogram-gated CT scans, acquired only during diastole to avoid respiratory and cardiac motion. While the program is operational with any cardiac CT images, to produce meaningful results, the source images should meet the minimum clinical calcium scoring guidelines11,12. For reference, a slice thickness of 3 mm, peak tube voltage of 100 kVp, average CT dose index-vol of 1.19 mGy, and image resolution of 512 x 512 pixels are used in the study here. Any images that are not 512 x 512 pixels are resampled in the program automatically to ensure adequate and consistent resolution of small areas of calcification. Once the images are loaded, the user is able to see them in the axial, sagittal, and coronal views. One may then adjust the brightness and contrast of the images for better visualization before selecting the 0 mg/mL and 100 mg/mL sections of the phantom. Next, the user can trace each of the four coronary arteries-left anterior descending (LAD), left coronary artery (LCA), left circumflex (LCX), and right coronary artery (RCA)-by placing either a point, a region of interest (ROI), or a combination of both to allow for a thorough selection of an artery's pixels regardless of how the artery appears in the axial plane. The user may delete and replace or redraw points and ROIs as needed. Clicking the SWCS button generates the final reports. Cases are auto-saved so that images, along with the points and ROIs, can be reloaded at a later time. Written instructions are also available at every point while using the program, making the program easy to use.
This study was conducted with approval of the Mount Sinai Institutional Review Board (HS-20-01011), and all subjects gave written informed consent.
1. Preparation before starting the protocol
Figure 1: Format of main project folder. This figure shows how the project's main folder should be structured and formatted for proper use of the program. Please click here to view a larger version of this figure.
2. Launching the Program
Figure 2: Initial program window. The program, when initially launched, has the buttons laid out along with an art image. Please click here to view a larger version of this figure.
Figure 3: Graphical user interface (GUI). Once images are loaded in, the program's GUI shows three anatomical views of the images along with crosshairs on each view, representing the cursor. Please click here to view a larger version of this figure.
3. Analyzing coronary artery calcification
Figure 4: Draw ROI feature. When the Draw ROI option is chosen, a pop-up of the current axial slice appears. The yellow shows an ROI that was previously drawn on this slice. Please click here to view a larger version of this figure.
4. Accessing the results
The representative results shown in this section display what successful use of the program entails. Here, a patient with an Agatston score greater than zero is used as an example. As discussed earlier, the results within a patient's metadata folder will have spreadsheets in the form of CSV files, images in the form of PNG files, and reports in the form of PDF files, as shown in Figure 5. The number of PNG files differs from case to case, since only snapshots of selected pixels with noticeable calcium (HU > 130) are included. There are also images of the weighting function, the phantom, and the trajectory of points/ROIs for each artery in 3D space. These images appear in the reports. There is one report for each artery analyzed. The reports give the SWCS and the Agatston score, the weighting function, point/ROI trajectory, and any snapshots of the phantom and artery. A report with the total scores is also included and only has the scores and the weighting function. Figure 6 shows what the first page of the LCX report looks like for this case, while Figure 7 shows representative images of the points'/ROIs' trajectory graph, the phantom snapshot, and a noticeable calcium snapshot.
To validate the program's calculation of the Agatston score, a validation study was carried out comparing the program's output to that of commercially available software. A total of 10 cases known to have calcium in the coronary arteries were analyzed by two image analysts separately on both the program and the commercial software. Cases with calcium present in the arteries were used to avoid Agatston scores of zero, which would not be useful for comparison purposes. The total Agatston score (sum of the Agatston scores from each artery) from both tools was collected for the 10 cases and analyzed on a Bland-Altman plot (Figure 8). The 95% confidence interval was ± 17 percentage points from the mean.
Figure 5: Contents of results folder. The metadata folder for a given patient has the shown CSV, PNG, and PDF files if the program is used correctly. Please click here to view a larger version of this figure.
Figure 6: LCX report. This example shows what the first page of a report should look like. The SWCS and Agatston score are displayed in red, along with the extent of calcification-the number of slices included in the Agatston score calculation. The phantom-derived weighting function is also displayed, which shows a given pixel's weight according to its attenuation level. Please click here to view a larger version of this figure.
Figure 7: Various PNGs. The report for each artery includes A) a graph showing the trajectory of the labeled points/ROIs, B) a snapshot of the phantom, and C) one or more snapshots of noticeable calcium, if any. Please click here to view a larger version of this figure.
Figure 8: Validation of program Agatston score. This Bland-Altman plot shows the percent difference between the Agatston score obtained from the program versus that obtained from the commercial software for 10 cases known to have calcium in one or more coronary arteries. Please click here to view a larger version of this figure.
Supplementary File 1: SWCS code. This file contains the program's code that is to be run for SWCS measurement. Please click here to download this File.
Supplementary File 2: Lamprocapnos spectabilis. This is the art image displayed on the main program window when initially launched. Please click here to download this File.
While the protocol for this program is relatively easy to follow, there are a few critical steps that are necessary for successful use and reliable results. Before starting, it is important to make sure the patient data that will be used in this program is anonymized to ensure patient confidentiality. The initial formatting and naming of the project’s main folder must be correct for the program to recognize where to pull and place data. Incorrect naming and/or placement of folders, especially the Meta_Data folder, leads to errors in the program. Not including the cover art image in the project’s folder also leads to the inability of the program to run, since it specifically looks for the image. It is also critical to make sure the true center of each phantom is chosen by checking all three views. This ensures that accurate points are pulled from not only the slice where the point is placed, but also the slices above and below. Placing a point too high up or too far down the phantom can result in points of empty space, or air, being used in the weighting function calculation. Finally, it is important to close the program window after each case is done running. To analyze another patient’s images, the program is relaunched by clicking the Run button. This ensures proper image viewing in the program window.
Since the majority of the method presented is software-based, troubleshooting mainly involves checking the inputs for the program. As stated before, the folder formatting and naming are critical and should be the first thing to check when running into errors. Another simple check is making sure the cover art is in the correct location. Also, one should be sure that only DICOM files are inputted into the program; any other file type within the patient’s original data folder will lead to errors. Another, less common reason for errors in the program is not having the correct toolboxes downloaded for the chosen programming language necessary for DICOM processing and some mathematical calculations. For errors that are not explained here, it is helpful to use the software’s help center to explain errors found in the command window.
While this method is effective and successful in obtaining an accurate SWCS, there are some limitations to the program. The reliance on folder structure and setup limits how a user can store project data up until the final reports are generated. This may take some adjusting if the user is not used to the folder structure that is required. Another limitation falls within the program itself. The ability to only place single points or free-form ROIs and having to label all slices for each artery limits how quickly arteries can be traced. The act of having to close each pop-up window after drawing an ROI also adds to the amount of time spent analyzing each case. However, despite these limitations, this method for generating the SWCS is effective and easy to learn.
The presented method is significant due to its novel nature. While the method for calculating the SWCS has been laid out thoroughly by others7, a program that calculates both the SWCS and Agatston score semi-automatically does not exist currently. The fact that this program calculates both scores saves time for the user by cutting out the extra step of using another program to get the Agatston score. As the importance of quantifying low levels of calcification continues to grow6, the need for a program that can generate the SWCS will also grow. This program will mainly be helpful to the field of cardiology, as the SWCS helps better understand the risk factors associated with atherosclerosis.
In conclusion, a novel tool has been implemented to calculate the SWCS and Agatston score, with the Agatston score validated against an independent tool. The tool will allow for robust calculation of the SWCS in future studies by multiple users to further the understanding and detection of subclinical coronary heart disease.
The authors have nothing to disclose.
This work was supported by NIH grant R01ES029967.
Calcium Hydroxyapatite | Sigma-Aldrich | 289396-100G | Suspended in EpoxAcast 690 resin for phantom creation |
Clinical Cardiac CT Scanner | Siemens | SOMATOM Force Dual Source CT | Used for the source images; Any cardiac CT will be sufficient |
EpoxAcast 690 | Smooth-On | 03641 | Used for phantom creation |
MATLAB | Mathworks | R2019a | Requires Image Processing Toolbox and Statistics and Machine Learning Toolbox; Any version compatible with and able to run version R2019a scripts is sufficient |
Standard Computer | N/A | N/A | macOS or Windows operating system |
syngo.via | Siemens | VB60A_HF04 | Commercial software used for computing Agatston score for validation study |