nsct变化工具箱

NSCT方法是由传统的contourlet变化改进而来，contourlet变化是用轮廓段的基结构来对图像的直线奇异和曲线奇异进行逼近检测，但其融合后的图像不具有平移不变性，没有很好的消除混频现象以及吉布斯现象。而本文提出的NSCT不但保留了contourlet变化的多尺度，多方向，各向异性等优点，在图像分解时采用来非下采样形式剪切波变换能够很好的避免图像由于分解与重构带来的细节丢失，更重要的是分解后的图像与原图像大小相同，因此能够更好的完整描述图像的方向性和特征。

/******************************************************************
* atrousc.c -  Written by Arthur Cunha. This routine builds up on 
*               zconv2D_OS.c written by Jason Laska
*
* Inputs:   x - A 2D signal
*           h - 2D filter
*           m - separable upsampling matrix
*         
* Outputs:  y - 2D result of convolution with filter 
*           upsampled by a m only the ‘valid‘ part is returned.
*           Similar to conv2（xh‘valid‘） where h is the upsampled
*           filter.
*  
*          
*
* Usage:    y = zconv2D_O（xhm）;
*
* Notes:    This function does not actually upsample the filter 
*           it computes the convolution as if the filter had been 
*           upsampled. This is the ultimate optimized version.
*           Further optimized for separable （diagonal） upsampling matrices.
*
* This is a MEX-FILE for matlab
*
/********************************************************/

#include “mex.h“
#include 

//Constants for matlab interfacing
#define OUT     plhs[0]
#define SIGNAL  prhs[0] //flip and shift
#define FILTER  prhs[1] //stationary
#define MMATRIX prhs[2]

//MACRO for converting positions to linear
#define LINPOS（rowcolcollen） （row*collen）+col


void mexFunction（int nlhs mxArray *plhs[] int nrhs const mxArray *prhs[]）
{
    //Declarations
   double *FArray*SArray*outArray*M;
/* FArray   - Filter coefficients
   SArray   - Signal coefficients
   outArray - Output coefficients
   M        - upsampling matrix 	*/
   int SColLengthSRowLengthFColLengthFRowLengthO_SColLengthO_SRowLength;
   int SFColLengthSFRowLength;
   int n1n2l1l2k1k2f1f2 kk2 kk1;
   double sum;   
   int M0M3sM0sM3;

    //Get the input sizes
    SColLength = mxGetM（SIGNAL）; 
    SRowLength = mxGetN（SIGNAL）;
    FColLength = mxGetM（FILTER）; 
    FRowLength = mxGetN（FILTER）;
    
    SFColLength = FColLength-1;
    SFRowLength = FRowLength-1;
    

	//Get The Data
    FArray = mxGetPr（FILTER）;
    SArray = mxGetPr（SIGNAL）;
    M = mxGetPr（MMATRIX）;
    M0 = （int）M[0];    
    M3 = （int）M[3];   
    sM0 = M0-1;
    sM3 = M3-1;
    

	// Corrected Lengths

	O_SColLength = SColLength - M0*FColLength + 1;
	O_SRowLength = SRowLength - M3*FRowLength + 1;
	
	
    //Make output size and Allocate out vector
    
 
    OUT      = mxCreateDoubleMatrix（O_SColLength O_SRowLength mxREAL）; 
    outArray = mxGetPr（OUT）;	//outArray is new vector
 
	/* Convoluyion loop */

    for （n1=0;n1		for （n2=0;n2			sum=0;		    
		    kk1 = n1 + sM0;;
			for （k1=0;k1  			    kk2 = n2 + sM3;
				for （k2=0;k2					 f1 = SFRowLength - k1; /* flipped index */
					 f2 = SFColLength - k2;  		
					 sum+= FArray[LINPOS（f1f2FColLength）] * SArray[LINPOS（kk1kk2SColLength）];					
					 kk2+=M3;
				}
				kk1+=M0;
			} 
		    outArray[LINPOS（n1n2O_SColLength）] = sum;
		}
	}

    return;
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2969  2004-10-17 16:42  nsct_toolbox\atrousc.c

     文件       7168  2004-10-17 16:45  nsct_toolbox\atrousc.dll

     文件       8664  2005-01-24 21:51  nsct_toolbox\atrousc.mexmac

     文件       1234  2005-11-03 09:10  nsct_toolbox\atrousdec.m

     文件      19406  2004-09-09 15:01  nsct_toolbox\atrousfilters.m

     文件       1012  2004-10-17 16:45  nsct_toolbox\atrousrec.m

     文件       3214  2005-01-21 15:59  nsct_toolbox\decdemo.m

     文件       1495  2004-12-22 22:48  nsct_toolbox\dfbdecdemo.m

     文件      13234  2004-10-18 12:28  nsct_toolbox\dfilters.m

     文件       2607  2004-12-22 21:26  nsct_toolbox\dmaxflat.m

     文件       1103  2003-04-10 11:30  nsct_toolbox\efilter2.m

     文件       1861  2003-11-06 00:25  nsct_toolbox\extend2.m

     文件        856  2003-04-10 11:30  nsct_toolbox\ld2quin.m

     文件        662  2003-11-05 10:17  nsct_toolbox\ldfilter.m

     文件        854  2003-11-05 10:19  nsct_toolbox\mctrans.m

     文件        775  2003-04-10 11:30  nsct_toolbox\modulate2.m

     文件       4662  2005-11-03 10:47  nsct_toolbox\nsctdec.m

     文件       3257  2005-11-03 10:33  nsct_toolbox\nsctrec.m

     文件       5046  2005-01-20 23:13  nsct_toolbox\nsdfbdec.m

     文件       4931  2005-01-20 23:14  nsct_toolbox\nsdfbrec.m

     文件        889  2005-12-19 11:12  nsct_toolbox\nsfbdec.m

     文件        821  2004-10-21 07:11  nsct_toolbox\nsfbrec.m

     文件       2424  2005-01-20 23:51  nsct_toolbox\nssfbdec.m

     文件       2785  2005-01-20 23:50  nsct_toolbox\nssfbrec.m

     文件       1552  2004-08-07 18:43  nsct_toolbox\parafilters.m

     文件       1161  2003-11-07 19:58  nsct_toolbox\qupz.m

     文件       1557  2005-01-20 23:45  nsct_toolbox\README.txt

     文件       1981  2003-04-10 11:30  nsct_toolbox\resampz.m

     文件       1001  2005-01-20 23:27  nsct_toolbox\shownsct.m

     文件        888  2005-12-19 11:07  nsct_toolbox\symext.m

............此处省略13个文件信息