24小时热门版块排行榜    

查看: 1389  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

whq_nj

金虫 (小有名气)

[求助] 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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖


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.

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
3楼2011-11-17 19:00:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

libralibra

至尊木虫 (著名写手)

骠骑将军


jjdg(金币+1): 感谢参与 2011-11-18 01:04:46
blkproc(S,[512 512],[3 3],fun);
这样调用可以返回需要的图片大小的,你把源图像S的size传入即可,这里指定是512,好像lena的宽高都是512
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
4楼2011-11-17 22:31:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-25 01:47 by BZKMTicpDhFj
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 6/300 2026-08-25 00:11 by BZKMTicpDhFj
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见