Abstract: In a three-dimensional image shading method using a volume rendering algorithm, the amount of light incident to a view plane from voxels of a predetermined value or above is determined by modulating the amount of light incident to the voxels by a function of sin .THETA., where .THETA. is the angle of the surfaces of CT slices to a perpendicular line drawn from an arbitrarily selected origin of an orthogonal three-dimensional coordinate system on the view plane.