24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jjgy_2002

木虫 (著名写手)

应助: 16 (小学生)
金币: 1285.3
散金: 1556
红花: 4
帖子: 1947
在线: 984.2小时
虫号: 646930
注册: 2008-11-05
专业: 信号理论与信号处理

[交流] 【求助】因不能追加赏币重新发帖：急求Tchebichef矩不变量或Krawtchouk矩不变量代码已有3人参与

因不能追加赏币重新发帖：
150金币急求Tchebichef矩不变量或Krawtchouk矩不变量代码，matlab,c,c++的都行。麻烦anyuezhiji和其他虫友辛苦辛苦给出正确结果代码，本人不胜感激，谢谢啦～～～

原帖见下面链接：
http://muchong.com/bbs/viewthread.php?tid=1957402&fpage=1
相关资料：
已上传T不变矩的资料，‘Image Analysis by Tchebichef Moments’这篇为经典的T矩文献，但里面没有涉及不变矩；其他几篇涉及T不变矩，最后一篇‘Image Analysis by Krawtchouk Moments’关于K不变矩的经典文献，与T矩类似，如能把K不变矩的代码编写出来并利用K不变矩对图像进行重构也行。能把K不变矩代码编好，T不变矩的就更没问题。下面是文献链接：
1http://d.namipan.com/d/(very%20g ... 5b0da3b37027c510300
2http://d.namipan.com/d/(muse%20s ... 361e9c2c8dbcecf0300
3http://d.namipan.com/d/Image%20A ... 385a3c14271526d0c00

[ Last edited by jjgy_2002 on 2010-5-2 at 10:35 ]

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有55人回复
最近几年招的学生写论文不引自己组发的文章已经有5人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复

1楼 2010-05-02 10:22:07

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

lijinfeng042

木虫 (小有名气)

Matlab

仿真EPI: 2
应助: 1 (幼儿园)
金币: 2156.1
散金: 115
帖子: 291
在线: 31.5小时
虫号: 1019062
注册: 2010-05-15
性别: GG
专业: 通信理论与系统

★ ★ ★ ★
zzuwangshilei(金币+2):多谢专家解释 2010-06-06 15:59:10
jjgy_2002(金币+50):非常感谢lijinfeng042的辛勤付出！我需求的是k不变矩，并利用其进行图像重构。 2010-06-16 12:08:46
robert2020(金币+2):辛苦了！ 2010-06-16 12:36:24

引用回帖:

Originally posted by lijinfeng042 at 2010-06-05 22:27:42:

我有疑问啊想了两天了呵呵卡住就是 -N 那不是pochhammer函数没意义? 求解答毕竟我不是专业的哦不知理解是否有误

呵呵到这里终结了剩下的问题我就不懂了我没学过图形处理抱歉了
经画图结果和你提供的文献Image Analysis by Krawtchouk Moments.pdf对比应该是正确的
先上程序

CODE:

%关于不动点Krawtchouk

%2010-06-05 by lijinfeng042

%hypergeometric functions (2F1)

%hypergeom([a,b],c,z) is the Gauss hypergeometric function 2F1(a,b;c;z).

%K1=Kn(x;p,N)=hypergeom([-n,-x],-N,1/p)

%p1=ρn(n;p,N)=((p-1)/p)^n*n!/(gamma(-N+n)/gamma(-N))

%w1=ωn(x;p,N)=binomial(N,x)*p^x*(1-q)^(N-x))

%K=Kn(x;p,N)*sqrt(ωn(x;p,N)/ρn(n;p,N))

%  warming;0

function K=Kraw(n,x,p,N)

if x> N 

     error('MATLAB:nchoosek:KOutOfRange','Kraw(n,x,p,N) x must be an integer between 0 and N.'); 

 end 

K1=hypergeom([-n,-x],-N,1/p);

p1=gamma(n+1)*((p-1)/p)^n/proohc(-N,n);

w1=nchoosek(N,x)*p^x*(1-p)^(N-x);

K=K1*sqrt(w1/p1);

end

CODE:

function prooh=proohc(a,n)

%给出 Pochhammer 符号函数 (a)_n。

%只针对a,n是整数的情形 

%2010-06-05 by lijinfeng042

    if a+n>a

      prooh=prod(a:a+n-1);

    else

         prooh=prod(a+n,a);

    end

end

CODE:

%作图

c=[];x=1:100;

for i=x;

    y=Kraw(1,i,0.3,100);

    c=[c y];

end

line(x,c)

CODE:

%Krawtchouk moments

%Kmn=symsum(symsum(K(n,x,p1,N-1)*K(m,y,p2,M-1)*f(x,y),0,M-1),0,N-1)

%Krawtchouk Moment Invariants

%Mmn=symsum(symsum(x^n*y^m*f(x,y),0,M-1),0,N-1)

function K=Krawtchouk_Moment(n,m,x,y,p1,p2,N,M,f)

% 对于一幅大小为N*M的数字图像f，用f(x,y)

% 表示坐标（x，y）点像素的灰阶，则该数字图像的n+m阶Krawtchouk矩为：Kmn

%---------------云是风的梦  2010-06-06

syms MM NN

ff=inline(f,'x','y');

Kmn=symsum(symsum(Kraw(n,x,p1,N-1)*Kraw(m,y,p2,M-1)*feval(ff,x,y),MM,0,M-1),NN,0,N-1);%双重求和

K=vpa(Kmn,10);%显示小数结果

end