The 2D+1 method for an efficient modeling of collimator blur in 3D SPECT imaging