# gregfreeman/matlabPyrTools

### Subversion checkout URL

You can clone with
or
.
Fetching contributors…

Cannot retrieve contributors at this time

executable file 44 lines (38 sloc) 1.19 KB
 % imStats(IM1,IM2) % % Report image (matrix) statistics. % When called on a single image IM1, report min, max, mean, stdev, skew, % and kurtosis (4th moment about the mean, divided by squared variance) % % When called on two images (IM1 and IM2), report min, max, mean, % stdev of the difference, and also SNR (relative to IM1). % Eero Simoncelli, 6/96. function [] = imStats(im1,im2) if (~isreal(im1)) error('Args must be real-valued matrices'); end if (exist('im2') == 1) difference = im1 - im2; [mn,mx] = range2(difference); mean = mean2(difference); v = var2(difference,mean); if (v < realmin) snr = Inf; else snr = 10 * log10(var2(im1)/v); end fprintf(1, 'Difference statistics:\n'); fprintf(1, ' Range: [%c, %c]\n',mn,mx); fprintf(1, ' Mean: %f, Stdev (rmse): %f, SNR (dB): %f\n',... mean,sqrt(v),snr); else [mn,mx] = range2(im1); mean = mean2(im1); var = var2(im1,mean); stdev = sqrt(real(var))+sqrt(imag(var)); sk = skew2(im1, mean, stdev^2); kurt = kurt2(im1, mean, stdev^2); fprintf(1, 'Image statistics:\n'); fprintf(1, ' Range: [%f, %f]\n',mn,mx); fprintf(1, ' Mean: %f, Stdev: %f, Skew: %f, Kurt: %f\n',mean,stdev,sk,kurt); end
Something went wrong with that request. Please try again.