基于GPU和分块技术的巨幅影像快速傅里叶变换算法研究(英文)

快速傅里叶变换(FFT)是遥感影像处理的基础方法,随着高光谱、高空间和高时间分辨率遥感影像获取能力的提升,如何利用快速傅里叶变换技术快速有效地处理巨幅遥感影像是当前遥感影像处理技术中的重要环节和研究热点。傅里叶变换算法FFT是基本的图像处理算法之一,该算法可进行遥感影像的条带噪声去除、影像压缩和影像配准处理等多种用途。CUFFT函数库是NVIDIA公司提供的基于GPU的FFT算法库,FFTW是由MIT科学实验室计算机组在PC平台上开发的基于CPU的FFT算法,是目前在基于CPU的运行速度最快的FFT算法函数库,这两种实现共有的问题是当可用内存或显存的容量小于图像容量时,就会出现内存或显存溢出。针对这种问题,提出了一种基于GPU和分块技术的巨幅遥感影像快速傅里叶变换(huge remote fast Fourier transform,HRFFT)算法。通过对CUDA的CUFFT函数库中的FFT算法进行改进,解决了巨幅图像内存或显存溢出的问题,并结合HJ-1A卫星的CCD影像,通过实验与其他算法进行了对比,证明了该方法的合理性。在实际应用中,利用本文提出的HRFFT算法,改善了影像处理的效果,提高了遥感影像的质量,同时加快了影像处理的速度,节省了计算时间,取得了较好的效果。     

基于JASMIN框架的FFT并行解法器及其应用

为解决并行应用程序使用FFTW(Fastest Fourier Transform in the West)并行软件包所面临的计算规模难以扩展、数据结构变动大、实现不同数据结构间通信难度大、接口不确定等问题,在JASMIN框架内设计实现FFT并行解法器.该解法器封装了数据分布存储、数据通信等并行计算细节,通过重新分布存储数据,调用一维FFT变换实现高维FFT的并行计算,并提供规范接口,支撑用户简便地实现FFT的并行计算.数值测试表明,该解法器具有很好的并行性能.该解法器已应用于激光等离子体成丝不稳定性的数值模拟并行程序,它在2048个处理器上的并行效率可达80%以上.     

基于傅里叶变换光谱技术的Zoom-FFT算法研究

在分析经典的Zoom-FFT算法基础上,提出一种基于傅里叶变换光谱技术的Zoom-FFT算法,用matlab仿真常规FFT算法和Zoom-FFT算法,对不同采样步长的干涉条纹进行数据处理,通过反演出的光谱曲线图和原始光谱曲线图可以看出:采样步长小于20 m时,FFT和Zoom-FFT算法都可以反演出光谱;而当采样步长大于20 m且小于33.3 m时,FFT算法未能反演出光谱,而Zoom-FFT算法仍然可以反演出光谱。     

基于二维傅里叶变换的单帧干涉图相位提取方法

针对大口径光学元件干涉测试过程中,测试装置和干涉腔长较大,气流扰动和环境振动对移相测试过程产生影响等问题,采用一种基于二维傅里叶变换的单帧干涉图处理方法,只需要对一幅空间载频干涉条纹图进行处理即可获得待测相位,具有抗振测试的优点。对该方法的基本原理和算法过程进行分析,并对近红外大口径移相平面干涉仪中600 mm口径的光学平晶进行了面形测试。实验结果表明:采用该方法所得波面峰谷值(PV)为0.112,波面均方根值(RMS)为0.014,与移相算法所得波面数据相比,波面峰谷值偏差不到(1/500);波面均方根值(RMS) 偏差几乎为零。     

Physical modeling of multiphase flow via lattice Boltzmann method: Numerical effects, equation of state and boundary conditions

The aims of the present paper are threefold. First, we further study the fast Fourier transform thermal lattice Boltzmann (FFT-TLB) model for van der Waals (VDW) fluids proposed in Phys. Rev. E, 2011, 84(4): 046715. We analyze the merits of the FFT approach over the traditional finite difference scheme and investigate the effects of smoothing factors on accuracy and stability in detail. Second, we incorporate the VDW equation of state with flexible parameters into the FFT-TLB model. As a result, the revised model may be used to handle multiphase flows with various critical densities and temperatures. Third, we design appropriate boundary conditions for systems with solid walls. These improvements, from the views of numerics and physics, significantly extend the application scope of the model in science and engineering.     

A fast numerical method for a natural boundary integral equation for the Helmholtz equation

A Neumann boundary value problem of the Helmholtz equation in the exterior circular domain is reduced into an equivalent natural boundary integral equation. Using our trigonometric wavelets and the Galerkin method, the obtained stiffness matrix is symmetrical and circulant, which lead us to a fast numerical method based on fast Fourier transform. Furthermore, we do not need to compute the entries of the stiffness matrix. Especially, our method is also efficient when the wave number k in the Helmholtz equation is very large.     

Noise source identification of diesel engine based on variational mode decomposition and robust independent component analysis

Through separating and identifying the noise sources of diesel engine, each independent noise obtained can be used as reference for the noise reduction, condition monitoring and fault diagnosis. In the noise source identification of diesel engine, the combustion noise and the piston slap noise are found to be overlapped in time domain and frequency domain. So it is intricate to accurately separate them. Therefore, the noise source identification method which is based on variational mode decomposition (VMD), robust independent component analysis (RobustICA) and continuous wavelet transform (CWT) is proposed. In the test, a 6-cylinder diesel engine was tested in a semi-anechoic chamber. The lead wrapped method was adopted to wrap No. 1–5 cylinders so as to isolate the interference noises, only the No. 6 cylinder part was bared. The single channel noise signal of cylinder head was measured. Then the variational mode decomposition algorithm is utilized to decompose the noise signal into several variational mode components. The RobustICA algorithm is adopted to extract the independent components. Finally, the continuous wavelet transform and the prior knowledge of diesel engine are applied to further identify the separated results. The results show that by using the proposed method to separate and identify the radiation noise of the cylinder head of the diesel engine, the independent components obtained are respectively the combustion noise and the piston slap noise. Comparing with the EEMD-RobustICA-CWT method, each independent noises obtained through the proposed method are more accurate and pure with less other interference noises.     

对称r-循环矩阵的快速算法和并行算法

借助于快速付立叶变换(FFT),给出了一种判断对称r-循环线性系统是否有解的快速算法,并且在有解的情况下求出其解,该算法的计算复杂度为O(nlogn),且具有很好的并行性,若使用n台处理机并行处理该算法则只需要O(logn)步.当r=0时,对称r-循环矩阵变成一个上三角型Hankel矩阵,我们也给出了此类矩阵求逆的一种算法.最后将该算法推广到线性同余系统,其运算量仅为O(nlogn).     

基于FPGA的数字相位计在煤矿提升机激光测距系统中的应用

介绍了用于矿井提升机红外激光位置跟踪系统中的一种新型数字鉴相器的研究方法。该方法是基于现场可编程门阵列(FPGA)器件,并运用快速傅里叶变换(FFT)算法完成的。阐述了在实时高速测量场合运用FPGA器件的优势,并给出了这种新型数字鉴相器的基本原理和具体实现过程。结果表明:数字相位计在应用于矿井提升机的位置跟踪时,不仅可以保证测量精度,而且在很大程度上简化了电路的复杂程度和消耗,具有良好的应用价值。     

Discussions about FFT-based two-step phase-shifting algorithm

An FFT-based two-step phase-shifting (TPS) algorithm is described in detail and implemented by use of experimental interferograms. This algorithm has been proposed to solve the TPS problem with random phase shift except π. By comparison with the visibility-function-based TPS algorithm, it proves that the FFT-based algorithm has obvious advantages in phase extracting. Meanwhile, we present a π-phase-shift supplement to the TPS algorithm, which combines the two interferograms and demodulates the phase map by locating the extrema of the combined fringes after removing the respective backgrounds. So combining this method and FFT-based one, one could really implement the TPS with random phase shift. Whereafter, we systematically compare the TPS with single-interferogram analysis algorithm and conventional three-step phase-shifting one. The results demonstrate that the FFT-based TPS algorithm has a satisfactory accuracy. At last, based on the polarizing interferometry, a schematic setup of two-channel TPS interferometer with random phase shift is suggested to implement the simultaneous collection of interferograms.