Automatic or semi-automatic segmentation and tracking of artery trees from computed tomography angiography (CTA) is an important step to improve the diagnosis and treatment of artery diseases, but it still remains a significant challenging problem. In this paper, we present an artery extraction method to address the challenge. The proposed method consists of two steps: (1) a geometric moments based tracking to secure a rough centerline, and (2) a fully automatic generalized cylinder structure-based snake method to refine the centerlines and estimate the radii of the arteries. In this method, a new line direction based on first and second order geometric moments is adopted while both gradient and intensity information are used in the snake model to improve the accuracy. The approach has been evaluated on synthetic images as well as 8 clinical coronary CTA images with 32 coronary arteries. Our method achieves 94.7% overlap tracking ability within an average distance inside the vessel of 0.36 mm.