[Problem] To provide a thorough image reconstruction method for interior CT which uses a priori knowledge which is more practical. [Solution] In this image reconstruction for interior CT, projection data is obtained from particle beams passing through an ROI within a subject being imaged, approximate reconstruction (first step) is performed by a CT image reconstruction method using the obtained projection data, a region where an image value representing a physical quantity is at least piecewise constant or represented at least piecewise by a polynomial equation is specified in the ROI on the basis of the reconstructed CT image, and reconstruction (second step) which is more precise than the first reconstruction step is performed using the position of the specified region where the image value representing a physical quantity is at least piecewise constant or represented at least piecewise by a polynomial equation, and the feature of the physical quantity therein which is at least piecewise constant or represented at least piecewise by a polynomial equation.La présente invention vise à fournir un procédé de reconstruction d'image complète pour tomodensitométrie (CT) intérieure qui utilise une connaissance a priori qui est plus pratique. Dans cette reconstruction d'image pour CT intérieure, des données de projection sont obtenues à partir de faisceaux de particules passant à travers une région d'intérêt (ROI) à l'intérieur d'un sujet étant imagé, une reconstruction approximative (première étape) est réalisée par un procédé de reconstruction d'image CT à l'aide des données de projection obtenues, une région où une valeur d'image représentant une quantité physique est au moins constante par morceaux ou représentée au moins par morceaux par une équation polynomiale est spécifiée dans la ROI sur la base de l'image CT reconstruite, et une reconstruction (seconde étape) qui est plus précise que la première étape de reconstruction est réalisée à l'aide de la position de la r
