1. A method of enhancing final 3D image information density by fusing images stacks acquired at different observation angles, comprising: performing a deconvolution process for each of a plurality of 3D image samples acquired from different observation angles, defined by rotation angles of a revolvable stage, by a processing unit designed to recover the distortion caused by a point spread function, wherein said performing a deconvolution process comprises performing a blind deconvolution to find said point spread function and the reconstructed undistorted 3D optical slices from said 3D image samples; performing a registration process of said deconvoluted plurality of 3D image samples acquired from different observation angles by said processing unit, wherein said performing a registration process comprises generating SIFT (Scale Invariant Feature Transform) matched features of every image of said plurality of 3D image samples; and performing an interpolation process of said registered 3D image samples acquired from different observation angles by said processing unit to output a final 3D image with increased depth resolution in said final 3D image; wherein said performing an interpolation process comprises performing a searching process, and performing an intensity rectification process to fix artifacts generated from an intensity-based interpolation; wherein said performing an intensity rectification process comprises providing a weighting function to maintain continuous intensity changes between upper slice and lower slice and a rectification function to satisfy said continuous intensity changes.