%
% GetDistortion.m
%
% function D = GetDistortion(Itrue, Iest)
%

function D = GetDistortion(Itrue, Iest)

D = mean(mean((Itrue-Iest).^2));