function [noisy_image, rpsnr] = psnr_noise( image, tpsnr )
    rpsnr = 0;
    sigma = 0;
    noisy_image = image;
    while rpsnr < tpsnr(1) || rpsnr > tpsnr(2)
        sigma = sigma + 1;
        noisy_image = image + randn(size(image)) * sigma;
        rpsnr = psnr(uint8(image), uint8(noisy_image), 255);
    end

    noisy_image = double(noisy_image);

end