A computer-implemented method for autonomous segmentation of contrast-filled coronary artery vessels, the method comprising the following steps: receiving (101) a CT scan volume representing a 3D volume of a region of anatomy that includes a pericardium; preprocessing (102) the CT scan volume to output a preprocessed scan volume; converting (103) the CT scan volume to three sets of two-dimensional slices, wherein the first set is arranged along the axial plane, the second set is arranged along the sagittal plane and the third set is arranged along the coronal plane; extracting (104) a region of interest (ROI) by by autonomous segmentation of the heart region as outlined by the pericardium, by means of three individually trained ROI extraction convolutional neural networks (CNN), each trained to process a particular one of the three sets of two-dimensional slices to output a mask denoting a heart region as delineated by the pericardium; combining (105) the preprocessed scan volume with the mask to obtain a masked volume; converting (106) the masked volume to three groups of sets of two-dimensional masked slices, wherein the first group is arranged along the axial plane, the second group is arranged along the sagittal plane and the third group is arranged along the coronal plane and each group includes at least three sets, wherein the first set corresponds to the principal plane of the set and at least two other sets are tilted with respect to the principal plane; and performing autonomous coronary vessel segmentation (107) by autonomous segmentation of the sets of the two-dimensional masked slices by means of three individually trained segmentation convolutional neural networks (CNN), each trained to process a particular one of the sets of the two-dimensional masked slices to output a mask denoting the coronary vessels.