An automatic three-dimensional non-rigid registration scheme is proposed in this paper and applied to thoracic computed tomography (CT) data of patients with stage III non-small cell lung cancer (NSCLC). According to the registration scheme, initially anatomical set of points such as the vertebral spine, the ribs, and shoulder blades are automatically segmented slice by slice from the two CT scans of the same patient in order to serve as interpolant points. Based on these extracted features, a rigid-body transformation is then applied to provide a pre-registration of the data. To establish correspondence between the feature points, the novel application of the self-organizing maps (SOMs) is adopted. In particular, the automatic correspondence of the interpolant points is based on the initialization of the Kohonen neural network model capable to identify 500 corresponding pairs of points approximately in the two CT sets. Then, radial basis functions (RBFs) using the shifted log function is subsequently employed for elastic warping of the image volume, using the correspondence between the interpolant points, as obtained in the previous phase. Quantitative and qualitative results are also presented to validate the performance of the proposed elastic registration scheme resulting in an alignment error of 6 mm, on average, over 15 CT paired datasets. Finally, changes of the tumor volume in respect to each reference dataset are estimated for all patients, which indicate inspiration and/or movement of the patient during acquisition of the data. Thus, the practical implementation of this scheme could provide estimations of lung tumor volumes during radiotherapy treatment planning.