paper 125:NSCT——Nonsubsampled contourlet 变换程序(尺度不变性问题研究)
08年,被老板逼得走投无路,xx所得项目看来是实在躲不过去,只好硬着头皮上。开发一款图像处理软件,里面内容很多,不过核心就是NSCT变换。NSCT变换06年的新东西,现在跟风研究的热点。可惜的是作为一种新的变换滤波器长期以来只有matlab程序,没有VC程序。使很多算法开发不能应用到工程中,08暑假,大概花了我一个月吧,恩。。期间电脑还中了次病毒(可恶的机器狗第四代!!),终于把这个软件搞得差不多了,现在回顾一下NSCT的程序包,下次有时间在写点别的。
Nonsubsampled Contourlet变换算法介绍:
对信号的稀疏表示是许多信号处理及应用的基础,2004年Minh N Do、Martin Vetterli提出了一种能够较好表示二维信号的数学工具--Contourlet变换。Contourlet是用金字塔方向滤波器组(PDFB)来将图像分解成不同尺度下的方向子带的。根据PDFB的结构,PDFB是一个拉普拉斯金字塔滤波器Laplacian Pyramid (LP)和一个方向滤波器组的叠加。实验证明,Contourlet变换在图像降噪,纹理,形状的特征提取方面的性能比2-D离散小波变换有了明显的提高。
为了获得平移不变性,本章所用的Nonsubsampled Contourlet变换(NSCT)是基于Nonsubsampled金字塔(NSP)和Nonsubsampled方向滤波器(NSDFB)的一种变换。首先由NSP对输入图像进行塔形分解,分解为高通和低通两个部分,然后由NSDFB将高频子带分解为多个方向子带,低频部分继续进行如上分解。NSCT是一种新型平移不变,多尺度,多方向性的快速变换。
1. Nonsubsampled Pyramid(NSP):
Nonsubsampled Pyramid(NSP)和Contourlet的Laplacian Pyramid(LP)多尺度分析特性不同。图像通过Nonsubsampled Pyramid(NSP)进行多尺度分解,NSP去除了上采样和下采样,减少了采样在滤波器中的失真,获得了平移不变性。NSP为具有平移不变性滤波结构的NSCT多尺度分析,可以得到与LP分解一样的多尺度分析特性。图2.4(a)处分为3个尺度。
2. Nonsubsampled方向滤波器(NSDFB)
Nonsubsampled方向滤波器(NSDFB)是一个双通道的滤波器,将分布在同方向的奇异点合成NSCT的系数。方向滤波器(DFB)是Bamberger and Smith提出的。其通过一个l层的树状结构的分解,有效的将信号分成了 个子带,其频带分割成为锲形。Nonsubsampled DFB(NSDFB)为非采样,减少了采样在滤波器中的失真,获得了平移不变性。并且每个尺度下的方向子图的的大小都和原图同样大小,Contourlet变换为所有子带之和等于原图。NSCT有更多的细节得以保留,变换系数是冗余的。图2.4(b)为三个尺度下对图像频域的分割图,其中每个尺度的方向子带数目以2倍递增,以在1,2,3尺度下的方向子带数目分别为2,4,8个。
NSCT内部变量:
下面列出程序中主要的变量和结构体定义说明:
图像信息结构体:
typedef struct
{ int high; //-----------------------------------------图像高度
int width; //---------------------------------------图像宽度
double max; //--------------------------------------最大像素值
double min; //---------------------------------------最小像素值
double *data; //--------------------------------------像素数据区域
}ImgData;
NSCT变换空间结构体:(用于存放NSCT变换滤波器,尺度、自带空间等基本信息)
typedef struct
{ unsigned int *nlevels;//子带方向数 必须为偶数 sublevels[0]为最高级方向数
sublevels[4]为最低级方向数 为负时为无效
unsigned int clevels; //------------------尺度层数
ImgData ***high;//------------------------高通子带区域
ImgData low;//-----------------------------低通子带区域
ImgData H1;//------------------------------高通滤波器组
ImgData H2; //------------------------------高通滤波器组
ImgData G1; //------------------------------低通滤波器组
ImgData G2; //------------------------------低通滤波器组
ImgData ***filters_dec; //-----------------分解方向滤波器组
ImgData ***filters_rec; //-----------------重构方向滤波器组
BOOL sign;
}nsct_t;
尺度与分解层数滤波器:
const int dlevels[3]={3,3,3};//--------------尺度为三,个尺度方向数为
const int Qunx[4]={0,0,0,0};//--------------------------Q梅花采样矩阵
BMGImageStruct SourceBmp;//---------------------------位图操作结构体
int pfilter_type,dfilter_type; //滤波器类型选择, 两者都默认为,'maxflat'&'dmaxflat7'滤波器
int pralevels; //------------图像金字塔层数(不能大于4,尺度在4以内足以满足一般图像处理需求)
nsct_t nsct; //nsct-------数据计算区域
int shift[2]; //---------------------延迟分量
2.3.6 NSCT内部函数:
下面列出程序中主要的函数定义说明:
void ImgDataM_Var_3m(ImgData *x,float *mean,float *var,float *mom);
//----------------------------------计算nsct子带系数的均值和三阶中心矩
void ImgDataIFFT_2D(complex<double> * pCFData, complex<double> * pCTData, int nWidth, int nHeight);
//----------------------------------2D-FFT反变换
void ImgDataFFT_2D(complex<double> * pCTData, int nWidth, int nHeight, complex<double> * pCFData);
//----------------------------------2D-FFT反变换
void ImgDataFConv2(ImgData *out,ImgData *img,ImgData *mask);
//----------------------------------频域卷积函数
void Nsct_Extend2(ImgData *in,int ru,int rd,int cl,int cr,ImgData *out,const int flag=1);
//----------------------------------nsct镜像扩展函数
void Nsct_Zconv2(ImgData *in,ImgData *out,ImgData *filter,const int *mup);
//----------------------------------时域滤波器卷积函数(快速算法)
void Nsct_Efilter2(ImgData *in,ImgData *out,ImgData *filter);
//----------------------------------边缘扩展的2D滤波
void Nsct_Nssfbdec(ImgData * input,ImgData * out1,ImgData * out2,ImgData * filter1,ImgData * filter2,const int * Qx=Qunx);
//----------------------------------2通道非采样子带分解
void Nsct_Mirror2(ImgData *tomirror);
//----------------------------------镜像变换
void ImgDataConv2(ImgData *out,ImgData *img,ImgData *mask,ImgData *imgunext);
//----------------------------------相关卷积
void Nsct_Upsample2df(ImgData *h0,int lev,ImgData *ho);
//----------------------------------频域2D上采样
void Nsct_Symext(ImgData *imglow,ImgData *upedHi,ImgData *pp );
//----------------------------------2D平衡矩阵生成
void Nsct_Atrousc(ImgData * out,ImgData * signal,ImgData *filter,int (*matrix)[2]);
//----------------------------------包含抽取的滤波器卷积
void Nsct_Nsdfbdec(ImgData *imghigh,int currentlevel);
//----------------------------------非抽取方向子带分解
void Nsct_iTtransform();
//----------------------------------反变换主函数
void Nsct_New(HBITMAP hBitmap,const int *Nlevels=dlevels,const unsigned int Level=3,const int dfilter=1,const int pfilter=1);
//-----------------------------用来初始化原始目标图像,申请nsct计算空间和存储空间
void Nsct_Transform();
//-----------------------------变换主函数
void Nsct_NsfbDec(ImgData *imglow,ImgData *imghigh,int i);
//-----------------------------nsct图像金字塔分解
void NsctDelete();
//-----------------------------nsct析构
void ImgDataFreeMemery(ImgData *image);
//-----------------------------内存空间释放
BOOL ImgDataAllocMemery(ImgData * image);
//-----------------------------内存空间申请
一个示例程序,和显示的样图:
paper 125:NSCT——Nonsubsampled contourlet 变换程序(尺度不变性问题研究)的更多相关文章
- 【尺度不变性】An Analysis of Scale Invariance in Object Detection – SNIP 论文解读
前言 本来想按照惯例来一个overview的,结果看到1篇十分不错而且详细的介绍,因此copy过来,自己在前面大体总结一下论文,细节不做赘述,引用文章讲得很详细,另外这篇paper引用十分详细,如果做 ...
- paper 78:sniff抓包程序片段
#define INTERFACE "eth0"#define MAX_SIZE 65535 int init_raw_socket();int open_promisc(char ...
- iOS 程序进入后台,包含用户上拉快捷菜单导致程序失去活跃的研究
今日在使用某App时候,突然发现上拉菜单.程序视频扔在播放,咦!引起了我的兴趣. 首先,列出两个方法, 第一个方法是AppDelegate的代理.当程序进入后台时候调用 - (void)applica ...
- 程序跳过UAC研究及实现思路(两种方法,现在可能都不行了)
网上很对跳过UAC资料都是说如果让UAC弹出窗体,并没有真正跳过弹窗,这里结合动态提权+计划任务实现真正意义上的跳过UAC弹窗,运行程序的时候可以不出现UAC窗体,并且程序还是以高权限运行. vist ...
- shell复习笔记----入门知识
Unix 简史 UNIX 最初是由贝尔实验室(Bell Telephone Laborataries)的计算机科学研究中心开发的,第一版诞生于1970年--也就是在贝尔实验室退出Multics项目不久 ...
- Fast-RCNN论文总结整理
此篇博客写作思路是一边翻译英文原文一边总结博主在阅读过程中遇到的问题及一些思考,因为博主本人阅读英文论文水平不高,所以还请大家在看此篇博客的过程中带着批判的眼神阅读!小墨镜带好,有什么不对的地方请在留 ...
- paper 65 :尺度不变特征变换匹配算法[转载]
尺度不变特征变换匹配算法 对于初学者,从David G.Lowe的论文到实现,有许多鸿沟,本文帮你跨越.1.SIFT综述 尺度不变特征转换(Scale-invariant feature transf ...
- 多尺度几何分析(Ridgelet、Curvelet、Contourlet、Bandelet、Wedgelet、Beamlet)
稀疏基的讨论已经持续了近一个月了,这次讨论多尺度几何分析.但由于下面讨论的这些变换主要面向图像,而本人现在主要关注于一维信号处理,所以就不对这些变换深入讨论了,这里仅从众参考文献中摘抄整理一些相关内容 ...
- paper 119:[转]图像处理中不适定问题-图像建模与反问题处理
图像处理中不适定问题 作者:肖亮博士 发布时间:09-10-25 图像处理中不适定问题(ill posed problem)或称为反问题(inverse Problem)的研究从20世纪末成为国际上的 ...
随机推荐
- CSipSimple 简单分析
简介 CSipSimple是一款可以在android手机上使用的支持sip的网络电话软件,可以在上面设置使用callda网络电话.连接使用方式最好是使用wifi,或者3g这样上网速度快,打起电话来效果 ...
- spark 2.0 Vector toBreeze
def toBreeze( _v : Vector ) : BZV[Double] = { _v match { case x : org.apache.spark.mllib.linalg.Dens ...
- IOS网络第四天 -网络文件上传(0923略)
01-NSURLSession02-断点续传 02-文件上传01-基本的上传 03-文件上传03-代码封装 04-文件上传04-获得MIMEType.mp4 05-文件的压缩和解压缩.mp4 06-压 ...
- 通过IP地址屏蔽各种“推广”
事情的起因是这样的:最近老是发现iPhone应用的底部出现各种横条广告,一开始以为是Google的广告推广,所以没管它,但是最近这些广告越来越猖狂,里面的内容越来越垃圾.今天仔细一看,原来不是Goog ...
- Redmi Note3 hennessy 刷机过程记录
本文只是凭记忆,记录大致的步骤,提供线索. 准备 刷机包和supersu刷机包,到xiaomi.eu上下载, 如果是稳定版可能有锁bootloader,需要到下载解锁软件.开发版无锁 刷入recove ...
- myBatis的一对多查询,主要利用resultMap实现一次查询多个结果集
日常开发中有这中场景,一个用户有多个角色,一个角色又有多个菜单,想查出一个用户的所有菜单.除了常见的关联查询之外,更使用的应该是利用myBatis的resultMap来实现一次查询出多个结果集,缺点: ...
- 关于type erasure
哇,好久没有写blog了,再不写的话,blog的秘密都要忘记了,嘿嘿. 最近在试着参与一个开源项目,名字叫avim(A Vibrate IM),别想多了哟.地址是:https://github.com ...
- buyexpressv6
<script type="text/javascript"> var is_enabled = -1; var checkSubmitFlg = false; var ...
- centos下JDK的卸载与安装
linux是自带JDK的,但是它自带的JDK是openJDK,我们如果需要安装ant之类的软件,使用这个JDK是不行的.所以我们需要卸载linux下自带的JDK,并安装我们准备的JDK. JDK的卸载 ...
- .Net深复制、浅复制
在.Net,大家都知道引用类型的深复制.浅复制吧. ,一般int等值类型是值类型(复制时是直接传值),一般的类(List<T>,Class)是引用类型(复制时传地址),默认是浅复制.若ob ...