You can compute the scale space in a different (and more efficient) way: You smooth the image with a Gaussian filter with increasing scale parameter, either by adding `2*ds` or multiplying the scale with a constant factor. Then you take the difference of two subsequent images.