Matlab PSNR和int16数据类型图像

此代码是否适合计算int16数据类型图像的峰值信噪比

function PSNR = PeakSignaltoNoiseRatio(origImg, distImg)

[M N] = size(origImg);
error = origImg - distImg;
MSE = sum(sum(error.^2)) / (M * N);
PSNR=10*log10(255*255/MSE);

---

不要使用

error

作为变量的名称。这是一个Matlab例程

在计算“错误”之前，将int16矩阵转换为例如双精度矩阵：

double（变量）