We aim to segment a 3D CT volume data into a set of organs and solve the segmentation problem as a labeling problem, in which each label represents an organ. We build a higher-order Markov random field (MRF) an solve the energy minimization problem by graph cuts.