|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
! N; {& S, @7 ?( r9 V8 N( H3 T* H一、简介9 Y b7 c& r" ^) f+ _! d4 y5 l. L
1 分形插值是根据分形几何自相似性原理和送代函数系统ifs理论,将已知数据插值成具有自相似结构的曲线或曲面,其中每个局部都与整体自相似或统计自相似,因此,分形插值可以有效地避免传统插值方法对相邻插值点间局部变化特征的掩盖。
7 U. z" J k# X5 L5 v8 ]5 t9 o ~1 R" ?绝大多数情况下,图像处理是为了改善图像的视觉质量。因此,如何评价图像的质量是一个十分重要的问题。
: R- C6 [2 O2 m4 R2 y在图像中,有些像素与相邻像素间灰度值存在突变,即存在灰度不连续性。这些具有灰度值突变的像素就是图像中描述对象的轮廓或纹理图像的边缘像素。在图像放大中,对这些具有不连续灰度特性的像素,如果采用常规的插值算法生成新增加的像素,势会使放大图像的轮廓和纹理模糊,降低图像质量。
; }; J! ^' N+ d, q. l4 a f迭代系统函数对于处理规则分形图形有明显的优势,但是对于自然界现象用随机分形处理会更加形象逼真。0 w. ]- d* U: G
2 案例# W, [/ w9 V3 B) J' J% |3 C$ o
已知人脸和狗脸图片,使用分形插值算法调换图中的人脸和狗脸。
" V5 a5 l4 C8 U0 g0 k
) \% H; g; }7 t
8 j4 g9 a: v" E! V* S* c/ \0 @, P( I# r
" T! b2 z0 v% S% S
) p7 N# Q7 U6 U. L) S4 H' s二、源代码
; v$ O! U8 G4 R1 A# E8 O
% L ~: Q5 b9 L! w5 g- clear all
- clc
- function rgbPic = bw2rgb(bwPic)
- bwPicSize = size(bwPic);
- rgbPic = zeros(bwPicSize(1),bwPicSize(2),3);
- I = imread('people.png');
- I2 = imread('dog.png');
- doge = imresize(I2,[size(I,1) size(I,2)]);
- gray = rgb2gray(I);
- YCbCr = rgb2ycbcr(I);
- heigth = size(gray,1);
- width = size(gray,2);
- for i = 1:heigth
- for j = 1:width
- Y = YCbCr(i,j,1);
- Cb = YCbCr(i,j,2);
- Cr = YCbCr(i,j,3);
- if(Y < 80)
- gray(i,j) = 0;
- else
- if(skin(Y,Cb,Cr) == 1)
- gray(i,j) = 255;
- else
- gray(i,j) = 0;
- end
- end
- end
- end
- SE=strel('arbitrary',eye(5));
- gray = imopen(gray,SE);
- gray = imclose(gray,SE);
- a1=(255-bw2rgb(gray))/255;
- next1=immultiply(I,a1);
- a2=bw2rgb(gray)/255;
- next2=immultiply(doge,a2);
- add=imadd(next1,next2);
- [L,num] = bwlabel(gray,8);
- STATS = regionprops(L,'BoundingBox');
- n = 1;
- result = zeros(n,4);
- part = zeros(h,w);
- for i = y:(y+h)
- for j = x:(x+w)
- if bImage(i,j) == 0
- part(i-y+1,j-x+1) = 255;
- else
- part(i-y+1,j-x+1) = 0;
- end
- end
- end
/ b: m3 U; T6 O3 |. P8 s
+ U: J% [. n- }3 m2 T; n: ?
4 C2 v& @7 ~5 |1 c2 V5 R9 M
! t8 c$ x, B9 h; }三、运行结果$ K2 p0 i4 y ~+ v! u6 l/ B; k, N& g
. ~3 J2 y- T' ^- m% p
`$ z0 }$ I8 s
* {/ B# N( {* o& W0 F q* p
|
|