The study was approved by the Ethics Committee of Shanghai Jiading Central Hospital and the informed consents were waived for the retrospective study and only using the clinical routine images. All methods were performed in accordance with the relevant guidelines and regulations.
Study population
Forty consecutive patients were retrospectively collected from Jiading Central Hospital from April 2021 to May 2021. The patients met the following criteria were included: The patients with suspected or diagnosed coronary artery disease and were performed CCTA examinations using uCT 960+. The exclusion criteria were: (1) Allergy to iodine- contrast media historically; (2) Renal insufficiency (serum creatinine > 15 mg/L); (3) Hyperthyroidism; (4) Pregnancy; (5) Severe respiratory failure or heart failure; (6) The clinical state is unstable and can't follow the breath-holding guidance; (7) No atrial premature beats, ventricular premature beats, or paroxysmal arrhythmia.
Image acquisition and reconstruction
All CCTA examinations were performed using a 320-row detector CT scanner (uCT 960+, United-imaging Healthcare). The scanner’s detector is 160 mm, and the gantry rotation time is 0.25 s. The scan used a prospectively electrocardiography triggering in a single heartbeat with z-axis covered from the inferior carina to the bottom of heart. Iodinated contrast (Iopamidol, 370 mgI/ml) was administered by using a high-pressure syringe at a rate of 5 ml/s via median cubital vein and followed by 25 ml of saline solution. The dose was 0.7 ml/kg of body. The scanning was triggered using a BolusTracking technology, the region of interest (ROI) was in the descending aorta of the four chambers of the heart and the trigger threshold was set at 120 HU, when exceeded it, the CCTA scanning program is automatically started with a delay of 5.2 s.
In reconstruction, the scanner automatically selected the best phase with a best phase selection method (ePhase) according to ECG6. And all images were reconstructed using hybrid iterative reconstruction algorithms in a commercial workstation (United-imaging Healthcare). The filtered function was C-SOFT-BA, and other reconstruction parameters were 200 mm FOV, 512×512 matrix size, 0.5 mm slice thickness and 0.5 mm slice spacing. Besides, a second reconstruction was carried out by adding a motion correction algorithm (CardioCapture, an commercial AI-assisted motion correction algorithm, United-imaging Healthcare) to acquire another CCTA image series 3. The details about CardioCapture could be found in Yan et al.’ paper3. Finally, two CCTA series were obtained, one series was with MC, and the other was not.
Image quality
All images were evaluated by two independent and experienced radiologists who were blinded to clinical information and reconstruction methods. The quality of image was scored using five-point Likert scale rating system. According to the previous study 2, Score 1 indicated very poor image quality and severe motion artifacts and hard to distinguish the CCTA, Score 2 indicated poor image quality and severe motion artifacts and artery could be recognized. Score 3 indicated adequate, moderate artifacts and could be used in clinical diagnosis. Score 4 indicated good, minor artifacts, Score 5 indicated excellent, no motion artifacts. The image with score 3–5 were diagnostic, while score 4–5 were defined excellent.
RCA and PCAT segmentation
The coronary artery was segmented using a pretrained VNet model backend with niftynet7,8, which were trained using 3000 CCTA images and achieved DICE 0.9. RCA was the targeted artery and saved as the ROI, then evaluate the segmentation results visually, if the segmentation was inaccurate, correct it manually. The pericoronary region was segmented using MITK (www.mitk.org, version: v2018.04.2), the specific steps included: dilating the ROI of RCA with the radius of the diameter of the RCA. The expanded part of the RCA was then re-segmented using − 190 - -30 HU threshold method, and the remained region was defined as PCAT.
Radiomics feature extraction
The radiomics features of the PCAT were extracted using pyradiomics software ( pyradiomics.readthedocs.io/en/latest/, version: 3.0.0)9, which complied with IBSI10,11. During feature extraction, the CCTA image was firstly resampled to 1*1*1 mm3, then the gray values were discretized using binWidth with 25. Meanwhile, except for original image, the transformed images using Laplacian of Gaussian (LOG) transformation with sigma 2, 3, 4, 5 and wavelet transformation were also included. Finally, 1218 radiomics features were obtained, which could be classified into 5 classes, shape feature, first-order feature, texture feature, LOG feature and wavelet feature. In these features, ‘original_firstorder_Mean’ was weighted average attenuation of all voxels within PCAT and used an unadjusted FAI, and the ‘original_shape_VoxelVolume’ was the volume of the PCAT. Besides, the volume of RCA was also calculated using the number of voxel within ROI of RCA multiplying the volume of single voxel.
Statistical analysis
Statistical analyses were performed with R software (R Foundation, version 4.0.3, https://cran.r-project.org/), with statistical significance defined by p < 0.05. The relationship between the image quality and heart rate was analyzed using Pearson’s correlation analysis. The image quality score of the CCTA images reconstructed with MC and without MC were compared using χ2 test. The difference of FAI, volumes of PCAT and coronary RCA between 2 image series were compared using paired sample t-test. The pearson’s correlation analyses were carried out to evaluate correlation between the change of FAI, volume of PCAT and the volume of RCA with the image quality change. Each radiomics features was also analyzed using the method mentioned above.