An apparatus adapted for evaluating optical characteristics of an optical system based on an image forming position of a point image that is formed from a point light source through the optical system, the apparatus comprising: a point image producing unit (1 to 4; 4b) having a plurality of the point light sources positioned in a vicinity of an object surface of the optical system and arranged in a planar form, and adapted to form a plurality of the point images through the optical system from the point light sources; an imaging unit adapted to image a plurality of the point images to produce a point image distribution image; a moving unit (5, 12) adapted to change a relative distance (i) between the point light source and the optical system or (ii) between the optical system and the imaging unit in an optical axis direction of the optical system; an imaging control unit adapted to control the moving unit to change the relative distance (i) between the point light source and the optical system, or (ii) between the optical system and the imaging unit and adapted to control the imaging unit to image the point image distribution image in every change of the relative distance (i) between the point light source and the optical system or (ii) between the optical system and the imaging unit; an image axis shape calculation unit adapted to perform processing of detecting an image forming position in each different relative distance of the point image based on image information of a plurality of the point image distribution images imaged at the different relative distance, characterized in that the processing includes calculating a regression expression concerning the image forming position and the relative distance for each of the point images to obtain a shape parameter of the image axis determined by the regression expression in each of the plurality of point images; and an evaluation value calculation unit adapted to fits a characteristic distribution model function representing distribution of the optical characteristics to the shape parameter in each of the plurality of point images, and adapted to calculate an evaluation value of the optical characteristics based on the fitting characteristics distribution model function. a point image producing unit (1 to 4; 4b) having a plurality of the point light sources positioned in a vicinity of an object surface of the optical system and arranged in a planar form, and adapted to form a plurality of the point images through the optical system from the point light sources; an imaging unit adapted to image a plurality of the point images to produce a point image distribution image; a moving unit (5, 12) adapted to change a relative distance (i) between the point light source and the optical system or (ii) between the optical system and the imaging unit (i) between the point light source and the optical system or (ii) between the optical system and the imaging unit in an optical axis direction of the optical system; an imaging control unit adapted to control the moving unit to change the relative distance (i) between the point light source and the optical system, or (ii) between the optical system and the imaging unit and adapted to control the imaging unit to image the point image distribution image in every change of the relative distance (i) between the point light source and the optical system or (ii) between the optical system and the imaging unit; an image axis shape calculation unit adapted to perform processing of detecting an image forming position in each different relative distance of the point image based on image information of a plurality of the point image distribution images imaged at the different relative distance, characterized in that the processing includes calculating a regression expression concerning the image forming position and the relative distance for each of the point images to obtain a shape parameter of the image axis determined by the regression expression in each of the plurality of point images; and an evaluation value calculation unit adapted to fits a characteristic distribution model function representing distribution of the optical characteristics to the shape parameter in each of the plurality of point images, and adapted to calculate an evaluation value of the optical characteristics based on the fitting characteristics distribution model function.