ivan@117: function [RMSErn, RMSEcd, rn_im, cd_im] = SMALL_vmrse_type2(orig, corr, recon)
ivan@128: %%  Implementation of VectorRMSE type2
ivan@117: % 
ivan@117: %   Input:
ivan@117: %       - Original image
ivan@117: %       - Corrupted image
ivan@117: %       - Reconstructed Image
ivan@117: %
ivan@117: %   Output:
ivan@117: %       - RMSErn - RMSE from residual noise (noise not completely removed)
ivan@117: %       - RMSEcd - RMSE from collateral distortion - excessive filtering
ivan@117: %       - rn_im  - image of residual noise
ivan@117: %       - cd_im  - image of collateral distortion
ivan@117: %
ivan@117: %   F. Russo, "New Method for Performance Evaluation of Grayscale Image
ivan@117: %   Denoising filters", IEEE Signal Processing Letters, vol. 17, no. 5,
ivan@117: %   pp.417-420, May 2010
ivan@117: 
ivan@117: %   Centre for Digital Music, Queen Mary, University of London.
ivan@117: %   This file copyright 2011 Ivan Damnjanovic.
ivan@117: %
ivan@117: %   This program is free software; you can redistribute it and/or
ivan@117: %   modify it under the terms of the GNU General Public License as
ivan@117: %   published by the Free Software Foundation; either version 2 of the
ivan@117: %   License, or (at your option) any later version.  See the file
ivan@117: %   COPYING included with this distribution for more information.
ivan@117: %%
ivan@117: 
ivan@117:     recon_int = round(recon);
ivan@117:     
ivan@117:     RN1 = ((orig<recon_int)&(recon_int<=corr));
ivan@117:     RN2 = ((orig>recon_int)&(recon_int>=corr));
ivan@117:     CD1 = ((orig<recon_int)&(recon_int>corr));
ivan@117:     CD2 = ((orig>recon_int)&(recon_int<corr));
ivan@117:     
ivan@117:     RMSErn = sqrt(sum(sum((RN1+RN2).*(orig-recon).^2)))/512;
ivan@117:     RMSEcd = sqrt(sum(sum((CD1+CD2).*(orig-recon).^2)))/512;
ivan@117:     rn_im=RN1+RN2;
ivan@117:     cd_im=CD1+CD2;
ivan@117:     
ivan@117: end
ivan@117: