24小时热门版块排行榜

返回列表

whq_nj

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 588
散金: 15
帖子: 216
在线: 66小时
虫号: 874419
注册: 2009-10-16
性别: GG
专业: 模式识别

[求助] matlaB图像边界问题，简化程序

tic
clear all
S=imread('lena.png');
S=double(S);
S=Normalized(S);
[p,q]=size(S);
average=zeros(p,q);
rme=zeros(p,q);
beta=zeros(p,q);
heta=zeros(p,q);
u=zeros(p,q);

for i0=2:1:p-1 %此处开始图像循环
for i1=2:1:q-1
      V=[S(i0-1,i1-1) S(i0-1,i1) S(i0-1,i1+1);
            S(i0,i1-1) S(i0,i1) S(i0,i1+1);
            S(i0+1,i1-1) S(i0+1,i1) S(i0+1,i1+1)];
      average(i0,i1)=mean(V(: ));%求窗口均值

%%%下面小循环为了求当前窗口的标准差
su=0;
for i=1:1:3
         for j=1:1:3
         su=(V(i,j)-average(i0,i1))^2+su;
         end
end
rme(i0,i1)=sqrt(su/9);%标准差
%%%%下面是最终想要求三个矩阵（应用均值和标准差）
      beta(i0,i1)=(rme(i0,i1)/average(i0,i1));
      u(i0,i1)=log(1/average(i0,i1));
      heta(i0,i1)=average(i0,i1);
end
end

function [normalized_matrix]=Normalized(matrix)
input_matrix=abs(matrix);
Max_input=max(input_matrix(: ));
Min_input=min(input_matrix(: ));
min_matrix=ones(size(input_matrix)).*Min_input;
normalized_matrix=(input_matrix-min_matrix)./(Max_input-Min_input+eps);

Elapsed time is 5.915782 seconds.

（V是一个3x3的窗口，用此窗口遍历整幅图像，求窗口内图像的均值和标准差；应用标准差和均值求三个矩阵beta，u和heta）

需要解决的问题：
1最后求得的图像均值，标准差和最后三个矩阵边界像素都为0。这个边界问题怎么解决？

2运算时间太长，我想这样窗口遍历图像的算法应该有更加高效的程序，精简程序或其他高效程序？

请高手指点，或可以站内和我联系，如果可以解决可以追加BB，呵呵[ Last edited by jjdg on 2011-11-16 at 01:49 ]

回复此楼

» 猜你喜欢

垃圾破二本职称评审标准已经有18人回复
职称评审没过，求安慰已经有53人回复
毕业后当辅导员了，天天各种学生超烦已经有5人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复
EST投稿状态问题已经有7人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2011-11-15 20:20:20

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

whq_nj

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 588
散金: 15
帖子: 216
在线: 66小时
虫号: 874419
注册: 2009-10-16
性别: GG
专业: 模式识别

自己顶一下

回复此楼

2楼2011-11-16 12:46:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

程序强帖: 40
应助: 817 (博后)
金币: 12914.1
红花: 64
帖子: 2238
在线: 287.3小时
虫号: 696514
注册: 2009-02-05
专业: 计算机软件

【答案】应助回帖

★
whq_nj(金币+20): 非常感谢及时回复，程序效率提高了不少，可beta，u和rme三个矩阵大小和源图像不一致，希望您再帮我看看 2011-11-17 22:22:27
余泽成(金币+1): 谢谢参与应助！ 2011-11-19 14:28:20

运行了5次,速度提高一倍以上

CODE:

 % 清空变量

clear all;

% blkproc后续变成blockproc了,但输入不再是矩阵,而是struct,调用略有差别

% 这里不显示这个warning

warning off;

% 计时

tic;

% 读图

S = imread('lena.png');

% 归一化

S = double(S);

abss = abs(S);

S = (S-min(abss(:)))/(max(abss(:))-min(abss(:)));

% 均值

funmean = @mean2;

heta = blkproc(S,[3,3],funmean);

% 标准差

funstd = @std2;

rme = blkproc(S,[3,3],funstd);

% 标准差/均值

beta = rme./heta;

% log结果

u = log(1./heta);

% 打印处理时间

toc;

运行5次,最小2.5s不到,跟分块处理有关

CODE:

 
Elapsed time is 2.481932 seconds.

Elapsed time is 2.462658 seconds.

Elapsed time is 2.448344 seconds.

Elapsed time is 2.451187 seconds.

Elapsed time is 2.454321 seconds.