Abstract-Dynamic single photon emission computed tomography (SPECT) has demonstrated the potential to quantitatively estimate physiological parameters in the brain and the heart. The generalized linear least square (GLLS) method is a well-established method for solving linear compartment models with fast computational speed. However, the high level of noise intrinsic in the SPECT data leads to reliability and instability problems of GLLS for generating parametric images. An integrated method is proposed to restrict the noise in both the temporal and spatial domains to estimate multiple parametric images for dynamic SPECT. This method comprises three steps which are optimum image sampling schedule in the projection space, cluster analysis applied postreconstruction and parametric image generation with GLLS.