This paper describes a method based on an energy minimizing deformable model applied to the 3D biomechanical modeling of a set of organs considered as regions of interest (ROI) for radiotherapy. The initial model consists of a quadratic surface that is deformed to the exact contour of the ROI by means of the physical properties of a mass-spring system. The exact contour of each ROI is first obtained using a geodesic active contour model. The ROI is then parameterized by the vibration modes resulting from the deformation process. Once each structure has been defined

the method provides a 3D global model including the whole set of ROIs. This model allows one to describe statistically the most significant variations among its structures. Statistical ROI variations among a set of patients or through time can be analyzed. Experimental results are presented using the pelvic zone to simulate anatomical variations among structures and its application in radiotherapy treatment planning.