|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
源代码
6 {0 b$ K, l: Z a! b" w
$ b" w7 B6 d8 r* `0 o``
8 ]# s3 R/ I5 j9 R ?clear all* I9 j0 p% B3 P( o' a9 T
clc- R/ I' `9 C& U! ^
$ x) m1 R$ i9 E7 ~6 S6 u
im = imread(‘1.bmp’);8 b2 s( i# q& }- B3 x1 ?
/ o4 |& |* S5 v1 c7 @7 c Iquality = CCF(im)
?' j' e' v: y& L; qfunction [quality] = CCF(im)
% r6 q( g+ H' V9 @6 e5 Y4 a6 P0 z& K, Y* _3 S- n1 v- Y( V
%-------------CCF_coloRFulness -----------------
9 M. `9 T) S( Y$ BimColor = double(im);
7 f9 E8 V) q& ?" U% e9 e
- o: Q% H& p3 V4 F: v. H5 I' ]2 g( p% [: t0 {% y6 z
4 }8 m3 |3 v* h* r
- 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);1 T( l1 a$ J+ |( K+ ?
9 ~# H2 R* I U( c {
5 e+ A7 a9 j7 ?$ H& ^/ g%-------------CCF_contrast----------------------, c- @7 B! }2 y& @
" `# @0 H2 J/ P# f( t, l6 M! X) u- im1=rgb2gray(im);
- CCF_contrast = CCFcontrast(im1);3 u. r3 X! p9 V' d9 l
9 ]$ L" v& _; `
, O3 Z: e1 p. X- O( x( ]
%-------------CCF_FADE--------------------------9 E. S W3 W j5 Z- z. n c$ l
# ~3 T. i, k9 w. W7 t, K4 q3 o- CCF_FADE = FADE(im);: T, N1 h3 ]2 q0 J
! K: @2 w4 a% E- ]; ~4 @, p: h. A, L; d, ` d5 e4 e0 G3 S/ c
%------------ normalization ------------------% K% Z# P& P- ~- n
7 W+ j9 r/ |: l) l- CCF_colorfulness=mapminmax(CCF_colorfulness,1,10);
- CCF_contrast=mapminmax(CCF_contrast,1,10);
- CCF_FADE=10-mapminmax(CCF_FADE,1,10);
; [; F- F" [9 H% D
4 L# @# Y1 q; h
6 i8 V _, d2 T/ ?! I4 S% ------------calculate image quality with coefficients---------------------
6 `1 O; B% n7 p1 _4 k' I8 n
4 F" O [1 {5 E( u; L- c=[0.17593 0.61759 0.33988 ];
- quality = c(1)*(CCF_colorfulness) + c(2)*(CCF_contrast) + c(3)*(CCF_FADE)
9 `. F# W& h* j3 L0 E" b
0 q! x9 U8 X( y- U- P. |/ x5 q
5 o/ Y# C: h$ Y) b' Cend- u( N/ L7 E5 J- ^; R1 u
% r- K* ^; o0 v/ P2 L0 r
1 X, z3 W% O% j9 W* U" r |
|
/1 
发表于 2021-4-22 17:28
收藏
淘帖
支持！
反对！