BLS-GSM Image Denoising Matlab Toolbox 1.0.3
============================================ 14 May 2004
Latest Revision: 23 Feb 2005
Thank you for using this software.
BLS-GSM stands for "Bayesian Least Squares - Gaussian Scale Mixture".
This toolbox implements the algorithm described in:
J Portilla, V Strela, M Wainwright, E P Simoncelli, "Image Denoising
using Scale Mixtures of Gaussians in the Wavelet Domain," IEEE
Transactions on Image Processing. vol 12, no. 11, pp. 1338-1351,
November 2003
adding some more possibilities for representations, besides the full
steerable pyramid, as explained in the package description below. This
tool implements image denoising assuming that the noise is Gaussian,
and that we know its power spectral density (it does not need to be white).
I have written it using as a starting seed some block-Wiener denoising
Matlab code written by Vasily Strela back in the year 2000.
The package also contains some previously released public domain
software authored by Eero P. Simoncelli, belonging to his MatlabPyrTools
toolbox (available at http://www.cns.nyu.edu/~eero/software.html), and
also some public domain routines taken from Rice Wavelet Toolbox
(http://www.dsp.rice.edu/software/rwt.shtml).
INSTALLATION ============
1) Use MATLAB@ 6.0 version or posterior (it has not been tested with
older versions)
2) Unpack the zip file into a directory in your hard drive, respecting
the original tree of subdirectories included in the pack
3) Make available this directory and all subdirectories to MATLAB@
by including the parent directory in the MATLAB@ path and checking the
appropriated box for including the subdirectories in the path
After the installation, you should be able to run the demo program
("denoi_demo").
PACKAGE DESCRIPTION ===================
MAIN DIRECTORY
* Demo program ("denoi_demo.m", in the main directory)
It shows an example of aplication of the method to a simulated noisy
image. The idea is that the parameters can be changed by the user,
experimenting their influence. This applies especially to the noise
parameters (variance and power spectral density) and the representation
parameters (explained below). It is also possible to modify the model
parameters (size of the neighborhood, including a parent or not, etc.),
although we recommend not to change them (except for, possibly, the
inclusion or not of a "parent" in the neighborhood).
Please, note that, by default, the code uses representation parameters
that do NOT correspond to the representation used in our work IEEE TIP
Nov 2003. The default set of parameters is faster, but does not always
provide the optimal results. For reproducing the results of the referred
article, please uncomment the appropriated lines in the "denoi_demo.m"
program.
* ReadMe.txt (this file)
DENOISING_SUBPROGRAMS DIRECTORY
* Main function for noise removal ("denoi_BLS_GSM.m")
This function takes a noisy image, and the (noise, representation and
model) parameters and applies the referred algorithm for cleaning it.
Besides dealing with issues such as boundary handling (mirror reflection
to avoid boundary artifacts), its main task is performed by a call to
the routines "decomp_recons_X.m", which are described below.
* Converting from the image domain to the wavelet domain and vice versa
("decomp_reconst_X.m")
(The "X" in "denoi_reconst_X.m" represents several options that are
described below). These functions take the noisy image (possibly
extended for boundary handling purposes) and apply the chosen
representation to transform it into the wavelet domain. Once it has been
decomposed into subbands, the latter are processed sequentially (in
couples corresponding to subbands of the same orientations adjacent in
the scale, in case of the model parameter "parent" has been set to "1",
or individually, otherwise), and then recomposed into a single image. In
the current version, 4 main representation options are available
(corresponding to the representation parameter "repres1"):
- Option "w": Orthogonal wavelet. Very fast, but of relatively poor
denoising performance, because of not being translation invariant.
- Option "uw": Undecimated version of orthogonal wavelet. In fact, this
is a highly redundant version of the orthogonal wavelet that avoids
intra-subband aliasing, but still keeps a multiscale structure. (The
lowest level of the pyramid are extended 2 in each dimension, whereas
the rest of the scales are extended 4 times in each dimension.) It
provides a good trade-off between computational cost and denoising
performance (which is very high).
- Option "s": Steerable Pyramid (Simoncelli & Freeman, 1995). It allows
the user to choose an arbitrary number of orientations. The splitting in
oriented subbands is not applied to the very high frequencies, which are
represented in a single subband (high-pass residual). With a moderate
computational cost for a modest number of orientations (4 or less), its
results depend on the type of image, being comparable (or slightly
worse) on average than those obtained with option "uw".
- Option "fs": Full steerable pyramid (described in the article cited
above). Same as previous one, but now also the very high frequencies are
splitted into orientations. It provides very high denoising performance
(for some images slightly better than "uw"), especially with high Nor
(for instance, 8), but it is very demanding computationally.
The relative effect on the denoising performance of using each of these
representations is evaluated in the referred IEEE TIP article. For the
"w" and "uw" representation options, there is also the option of
choosing a particular wavelet. We recommend the use of Daubechies
filters of order from 1 (Haar) to 3 or 4. For more information, please
read the help information included within the functions (or type "help
<function_name>").
* Subband Noise removal function ("denoi_BLS_GSM_band.m")
This function is called by "decomp_reconst_X.m" and implements the
kernel of the algorithm. It is only one, independently from which
representation is applied.
SIMONCELLI_PYR_TOOLS DIRECTORY
A subset with the required functions from Prof. Eero Simoncelli's
"MatlabPyrTools" toolbox (available at
http://www.cns.nyu.edu/~eero/software.html ). It includes the functions
that implement the steerable pyramid decomposition and reconstruction,
and some other functions used for the manipulation of the subbands, in
both its standard and its full version (full steerable pyramid, for the
"fs" option).
ADDED_PYR_TOOLS DIRECTORY
Functions that implement additional utilities, like the DWT code used
for building aliasing-free multi-scale transforms (for the "uw" option),
which was taken from the Rice Wavelet Toolbox
(http://www.dsp.rice.edu/software/rwt.shtml). Some others ({\it
expand.m, shrink.m, snr.m}) are general purpose functions written by
me (to expand and shrink an image in the Fourier domain and to compute
the signal-to-noise ratio in dB, respectively).
IMAGES
Several images in PNG (lossless compressed) format, to test the method.
They include the ones used in the IEEE TIP referred article.
Please, report to javier@decsai.ugr.es any bug or comment.
Javier Portilla
http://decsai.ugr.es/~javier
没有合适的资源?快使用搜索试试~ 我知道了~
BLS_GSM算法,用于图像去噪等
共89个文件
m:47个
png:8个
dll:7个
3星 · 超过75%的资源 需积分: 16 48 下载量 15 浏览量
2011-03-28
14:20:53
上传
评论
收藏 1.41MB RAR 举报
温馨提示
BLS_GSM 算法,很经典的。这是一个软件包,matlab编写,具体使用说明见txt文件
资源推荐
资源详情
资源评论
收起资源包目录
BLS_GSM_denoise.rar (89个子文件)
BLS_GSM_denoise
change_log.txt 1KB
denoising_subprograms
denoi_BLS_GSM_band.m 7KB
decomp_reconst_full.m 2KB
decomp_reconst_W.m 3KB
denoi_BLS_GSM.m 7KB
decomp_reconst.m 2KB
decomp_reconst_WU.m 3KB
denoi_demo.m 3KB
t2.bmp 63KB
Simoncelli_PyrTools
upConv.mexlx 23KB
vector.m 240B
corrDn.mex4 137KB
corrDn.mexsol 29KB
reconWpyr.m 4KB
innerProd.mexsol 4KB
reconSFpyrLevs.m 2KB
rcosFn.m 1KB
upConv.c 6KB
spyrNumBands.m 500B
pointOp.m 1KB
upConv.mexsol 30KB
subMtx.m 441B
showIm.m 6KB
pointOp.mexsol 7KB
range2.m 477B
range2.mexsol 5KB
buildSFpyrLevs.m 2KB
upConv.mex4 139KB
buildWpyr.m 3KB
namedFilter.m 3KB
innerProd.mexglx 7KB
innerProd.c 1KB
reconFullSFpyr2.m 3KB
maxPyrHt.m 628B
corrDn.mexlx 22KB
innerProd.m 415B
range2.c 2KB
corrDn.m 2KB
spyrBand.m 853B
upConv.dll 50KB
pyrBandIndices.m 613B
pointOp.mexlx 6KB
pointOp.dll 34KB
innerProd.mexlx 4KB
steer2HarmMtx.m 2KB
buildSFpyr.m 3KB
range2.dll 33KB
range2.mexlx 5KB
upConv.m 3KB
wpyrHt.m 285B
range2.mex4 10KB
innerProd.dll 41KB
buildFullSFpyr2.m 3KB
pyrBand.m 406B
var2.m 393B
pointOp.c 4KB
pointOp.mex4 12KB
spyrHt.m 321B
corrDn.dll 49KB
reconSFpyr.m 3KB
wpyrBand.m 912B
modulateFlip.m 480B
pixelAxes.m 2KB
shift.m 453B
corrDn.c 5KB
Added_PyrTools
mirdwt.c 3KB
mirdwt.m 5KB
mrdwt.m 5KB
snr.m 248B
buildWUpyr.m 2KB
bound_extension.m 2KB
daubcqf.m 4KB
shrink.m 977B
mrdwt.c 3KB
mrdwt.dll 9KB
reconWUpyr.m 2KB
mirdwt.dll 9KB
expand.m 1KB
change_log.txt.$$$ 1KB
ReadMe.txt 7KB
images
t2.bmp 63KB
barco.png 53KB
peppers256.png 39KB
flinstones.png 171KB
barbara.png 181KB
boat.png 174KB
fingerprint.png 180KB
lena.png 148KB
house.png 34KB
共 89 条
- 1
资源评论
- willsnowliu2014-07-02这个64位的运行不了,而且好像C都已经是dll了,看不到源码
- songyao20092012-05-17运行不起来啊,看来是需要修改的。
- shagirl2012-12-09算法有问题,缺少c语言与matlab的连接文件。
- lcx70282013-08-28不能运行啊,我要能运行的,图片什么格式
silvanmm
- 粉丝: 0
- 资源: 1
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功