The recent application of Fourier Based Iterative Reconstruction Method (FIRM) has made it possible to achieve high-quality 2D images from a fan beam Computed Tomography (CT) scan with a limited number of projections in a fast manner. The proposed methodology in this article is designed to provide 3D Radon space in linogram fashion to facilitate the use of FIRM with cone beam projections (CBP) for the reconstruction of 3D images in a sparse view angles Cone Beam CT (CBCT). For this reason, in the first phase, the 3D Radon space is generated using CBP data after discretization and optimization of the famous Grangeat's formula. The method used in this process involves fast Pseudo Polar Fourier transform (PPFT) based on 2D and 3D Discrete Radon Transformation (DRT) algorithms with no wraparound effects. In the second phase, we describe reconstruction of the objects with available Radon values, using direct inverse of 3D PPFT. The method presented in this section eliminates noises caused by interpolation from polar to Cartesian space and exhibits no thorn, V-shaped and wrinkle artifacts. This method reduces the complexity to for images of size n × n × n The Cone to Radon conversion (Cone2Radon) Toolbox in the first phase and MATLAB/ Python toolbox in the second phase were tested on three digital phantoms and experiments demonstrate fast and accurate cone beam image reconstruction due to proposed.