if blur:
            _prj = blur_edges(prj_copy, rin, rout)
            _sim = blur_edges(sim, rin, rout)
        else:
            _prj = prj_copy