function SNR = mComputeSNR( original_image, reconstructed_image );

signal = mean( mean( original_image.^2 ) );
noise  = mean( mean( ( original_image - reconstructed_image ).^2 ) );
  
SNR = 10*log10( signal/noise );