Global Navigation Satellite System (GNSS) is now widely used for continuous ionospheric observations. Three-dimensional computerized ionospheric tomography (3DCIT) is an important tool for the reconstruction of electron density distributions in the ionosphere through effective use of the GNSS data. More specifically, the 3DCIT technique is able to resolve the three-dimensional electron density distributions over the reconstructed area based on the GNSS slant total electron content (STEC) observations. We present an Improved Constrained Simultaneous Iterative Reconstruction Technique (ICSIRT) algorithm that differs from the traditional ionospheric tomography methods in 3 ways. First, the ICSIRT computes the electron density corrections based on the product of the intercept and electron density within voxels so that the assignment of corrections at different heights becomes more reasonable. Second, an Inverse Distance Weighted (IDW) interpolation is used to restrict the electron density values in the voxels not traversed by GNSS rays, thereby ensuring the smoothness of the reconstructed region. Also, to improve the reconstruction accuracy around the HmF2 (the peak height of the F2 layer) altitude, a multiresolution grid is adopted in the vertical direction, with a 10-km resolution from 200 to 420 km and a 50-km resolution at other altitudes. The new algorithm has been applied to the GNSS data over the European and North American regions in different case studies that involve different seasonal conditions as well as a major storm. In the European region experiment, reconstruction results show that the new ICSIRT algorithm can effectively improve the reconstruction of the GNSS data. The electron density profiles retrieved from ICSIRT are much closer to the ionosonde observations than those from its predecessor, namely, the Constrained Simultaneous Iteration Reconstruction Technique (CSIRT). The reconstruction accuracy is significantly improved. In the North American region experiment, the electron density profiles in ICSIRT results show better agreement with incoherent scatter radar observations than CSIRT, even for the topside profiles.