24小时热门版块排行榜    

查看: 2938  |  回复: 8

ggldggld2010

银虫 (小有名气)

[交流] 【求助】Weickert的各向异性扩散方程的滤波方法的matlab程序代码 已有4人参与

对于Weickert的各向异性扩散方程对图像进行滤波,我的困难在于边界条件的处理,所以难以实现,请大家帮忙如何实现,最好给个算例的matlab代码。非常感谢!
回复此楼

» 猜你喜欢

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

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

qshjdax

银虫 (小有名气)


ggldggld2010(金币+2):谢谢参与
别忘那个基于Biweight Estimator误差模型的扩散系数
2楼2010-04-01 12:46:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ggldggld2010

银虫 (小有名气)

谢谢!我想问一下边界条件如何离散化。
3楼2010-04-01 22:37:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cds0564

荣誉版主 (知名作家)

*﹀╮'梦幻cd¨∞


ggldggld2010(金币+2):谢谢参与
ggldggld2010(金币+3):谢谢!能不能说具体些。 2010-04-13 16:07
图像边界的话,通常镜面反射的比较多。
海纳百川,有容乃大
4楼2010-04-13 01:37:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cds0564

荣誉版主 (知名作家)

*﹀╮'梦幻cd¨∞

★ ★
无奈的云(金币+2):谢谢参与解答!很详细,用心了~ 2010-04-14 08:46
ggldggld2010(金币+9):非常感谢!你的解答很详细! 2010-04-14 16:29
边界的延拓方法通常有零延拓、周期延拓、对称延拓、光滑常数延拓、平滑延拓。为了便于书写,以一维的为例,设信号为S1,S2,...,Sn.
1. 零延拓    在原信号两端补零。也即   ..., 0, S1, S2, ..., ..., Sn, 0, 0, ...
2. 周期延拓  将信号看成周期信号,即S(k+n)=Sk。  ..., S(n-1), Sn, S1, S2, ..., Sn, S1, S2, ....
3. 对称延拓   也即镜面反射,等效于对原信号在边界上做对称折叠。 ...., S3, S2, S1, S2,..., Sn, S(n-1), S(n-2), ...
4. 光滑常数延拓    原信号两端加常数。..., S1, S1, S1, S2, S3, ... , Sn, Sn, Sn, ...
5. 平滑延拓   在原信号两端用线性插值法补充采样值。
本人觉得LZ用对称延拓就可以,如果是显示差分格式的话,好像延拓一个像素的量就可以了吧。
以上说的供LZ参考,希望对你有用。
海纳百川,有容乃大
5楼2010-04-13 23:56:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ggldggld2010

银虫 (小有名气)

请教各位:我的问题中给出了方程、初始条件和边界条件,方程的离散化方案(差分格式)有了,只是不知如何把边界条件=0离散化,其中D是结构张量。
6楼2010-04-15 17:57:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cctt126

金虫 (正式写手)


ggldggld2010(金币+2):谢谢参与
研究图像处理的都用什么软件?
当然是matlab首选.
为什么?
因为matlab是一款氛围最好的算法研究软件.
理由是matlab的官方网站的file exchange.
上面有Weickert的各向异性滤波的matlab代码.
http://www.mathworks.com/matlabc ... ce-enhancing-filter
不仅如此,其他更多更好的图像处理代码也有.

[ Last edited by cctt126 on 2010-5-8 at 13:52 ]
数学,英语,编程,没一样让我省心。。。
7楼2010-05-08 13:48:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

string121

金虫 (小有名气)


ggldggld2010(金币+2):谢谢参与
引用回帖:
Originally posted by ggldggld2010 at 2010-04-15 17:57:55:
请教各位:我的问题中给出了方程、初始条件和边界条件,方程的离散化方案(差分格式)有了,只是不知如何把边界条件=0离散化,其中D是结构张量。

是不是解u_t=div(D gradient u)?
条件较弱的边界条件应是:=0,n是图像区域边界外法向(由变分原理可以导出)。
用 gradient u=0 能保证上述条件成立(其实你的=0也只是同样能保证=0成立)。图像一般是矩形区域,gradient u=0 与Nemann边值条件就等价了。如果用中心差分,将u对称延拓一个像素即可。这时图像边缘像素满足:u_x=0,u_ y=0,即:gradient u =0.

=0数值计算不好处理。

[ Last edited by string121 on 2010-5-8 at 18:05 ]
8楼2010-05-08 18:01:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ggldggld2010

银虫 (小有名气)

谢谢你们!
9楼2010-06-10 16:14:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ggldggld2010 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见