|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
源代码4 _, p3 B4 i7 Z/ t5 D& y9 k
9 g5 ]( v! Q- W! }``4 ?9 X2 k; s: v* N9 Y' b! H
clear all+ N5 p7 {9 r; _2 B Q4 ~
clc
7 [* X- G# V9 J
* { w! A9 T" E0 q" I6 p3 Aim = imread(‘1.bmp’);
: }! W. H/ o' y# P1 f4 i0 @* ]% s/ U! b+ I6 S @# S7 e# u+ X2 w9 ]" |
quality = CCF(im)
: E0 E: L& }& l1 j$ Y' lfunction [quality] = CCF(im)
; ^8 L- k5 J C. U) r/ R( E. Z
$ m. p$ Z% y6 d/ h3 Z( b k9 g' S: N+ K%-------------CCF_coloRFulness -----------------# l6 [6 t. ?+ I6 I; j: y" m
imColor = double(im); [( F9 g$ M1 s) ?/ _! B
* u7 l/ J5 A. K Q/ H
* N+ o( k q/ d! a2 \8 O6 i6 W! Q
1 i, W' `2 O2 c& f- R = imColor(:,:,1);
- G = imColor(:,:,2);
- B = imColor(:,:,3);
- RR = log(R+0.00001) - mean2(log(R+0.00001));
- GG = log(G+0.00001) - mean2(log(G+0.00001));
- BB = log(B+0.00001) - mean2(log(B+0.00001));
- alpha = RR-GG;
- beta = 0.5*(RR+GG)-BB;
- mu_alpha=mean(mean(alpha));
- mu_beta=mean(mean(beta));
- var_alpha=var(var(alpha));
- var_beta=var(var(beta));
- CCF_colorfulness=1000*((sqrt(var_alpha+var_beta)+0.3*sqrt(mu_alpha*mu_alpha+mu_beta*mu_beta))/85.59);
" Z4 o* Q! f: b , P: s0 U# G# k V' s+ g; h
+ B) a; w4 X/ |+ H# d/ N4 g1 [%-------------CCF_contrast----------------------
b7 f3 z" {' \! C4 U
7 X2 q o V8 O" h. ]- im1=rgb2gray(im);
- CCF_contrast = CCFcontrast(im1);- v1 ~$ j7 t$ f9 _; O* Q, k
- k& w; A, N0 g4 e1 u/ B6 d; S
- I! D0 D, |9 w+ ?; q%-------------CCF_FADE--------------------------
; r- u1 I& B) O6 X$ i% C+ h& z" ]# X8 ^: S
- CCF_FADE = FADE(im);: Y. k* L7 c6 E2 d Y
0 }+ o$ n0 ]: x+ A8 d0 _# {. l4 Z; f9 d5 i1 m$ U
%------------ normalization ------------------
7 T% o, l) ?' j6 ]: W6 h4 D, P' P; ^ u" H- V2 w1 d" O
- CCF_colorfulness=mapminmax(CCF_colorfulness,1,10);
- CCF_contrast=mapminmax(CCF_contrast,1,10);
- CCF_FADE=10-mapminmax(CCF_FADE,1,10);
- D" V3 O& {; s: k 1 C( X& s! F* z5 [( N
* r* C t5 i4 S3 j: P1 X' w, ]- J
% ------------calculate image quality with coefficients---------------------
+ k! c8 v- u F- R0 O6 V
( g' P g. c% v9 {1 A* Y$ Q- c=[0.17593 0.61759 0.33988 ];
- quality = c(1)*(CCF_colorfulness) + c(2)*(CCF_contrast) + c(3)*(CCF_FADE)5 S( C8 d( ]; {7 J
, {. X8 x9 a6 ]; z$ ^* B2 a* u# k8 I8 x* _$ C
end
! J2 q! }- r2 E. W. x. l2 F; ~3 R4 D8 ^2 B
- f9 e. h& Q( a, M; C/ ~$ c: x/ A
|
|