波导环境多目标时延-多普勒解卷成像

有源声呐感兴趣的参量是目标距离和径向速度,它们无法直接观测得到,需要通过估计而获得。利用波导多路径环境多目标时延-多普勒模型,可以导出采样互模糊度函数均值是发射信号自模糊度函数与广义目标反射性密度函数的两维卷积,其中广义目标反射性密度函数为信道扩展函数与目标反射性密度函数的两维卷积。依据信息理论最小Csiszar鉴别准则,可导出R-L (Richardson-Lucy)迭代解卷算法,对采样互模糊度函数均值进行两维迭代解卷积,消除发射信号和信道引入的模糊,序贯地实现时延-多普勒两维像的估计,进而获得多目标的时延和多普勒参量估计。仿真结果和海上实验数据分析验证了R-L解卷算法的可行性和有效性,较之常规的匹配滤波和维纳滤波算法,R-L算法有效地提高了时延和多普勒估计的分辨力和精度。      

Minimum entropy deconvolution of pulse-echo signals acquired from attenuative layered media.

In this article deconvolution of ultrasonic pulse-echo data acquired from attenuative layered media is considered. The problem is divided in two subproblems: treating the sparse reflection sequence caused by the layered structure of the media and treating the frequency-dependent attenuation. The first subproblem is solved by means of joint maximum a posteriori estimation of the assumed zero mean, white, nonstationary reflection sequence and its corresponding sequence of unknown standard deviations. This approach leads to an algorithm that seeks minimum entropy solutions for the reflection sequence and therefore the algorithm serves as a novel link between the classical Wiener filter and methods for sparse or minimum entropy deconvolution. The second subproblem is solved by introducing a new signal processing-oriented, linear discrete-time model for frequency-dependent attenuation in isotropic and homogeneous media. The deconvolution algorithm is tested using simulated data and its performance for real normal incidence pulse-echo data from a composite material is also demonstrated. The results show that the algorithm, in combination with the attenuation model, yields estimates that reveal the internal structure of the composite and, thus, simplify the interpretation of the ultrasonic data.     

Performance of laser-ultrasonic F-SAFT imaging

The resolution and signal-to-noise ratio of laser-ultrasonics to detect small and buried defects can be greatly enhanced by using the synthetic aperture focusing technique (SAFT). Originally developed in the time domain, SAFT can also be implemented in the frequency domain (F-SAFT) using the angular spectrum approach for a significant reduction in processing time. In this paper, an F-SAFT based data processing method especially adapted to laser-ultrasonic data is presented. This method allows for further significant improvements towards laser-ultrasonic imaging of small defects. It includes temporal deconvolution of the waveform data, control for an optimal aperture and frequency bandwidth as well as spatial interpolation of the subsurface images. All the above operations are well adapted to the frequency domain calculations and embedded in the F-SAFT data processing. Also, the aperture control and spatial interpolation allow a reduction of sampling requirements to further decrease both inspection and processing times. The above improvements are illustrated using laser-ultrasonic data taken from an aluminum sample with flat-bottom holes.     

稀疏孔径光学系统成像的图像恢复算法研究

提出两种稀疏孔径光学系统成像的图像恢复模型.分析维纳滤波、最小二乘方滤波和极大似然法盲去卷积三种图像恢复算法的适用条件.针对存在噪声干扰的稀疏孔径光学系统,通过实验对比,指出维纳滤波和最小二乘方滤波把相机光学传函当作系统传函,其理论推导能够达到最优.盲去卷积把大气传输函数和相机光学传函作为系统传函进行恢复,其恢复结果优于维纳滤波带入常数K和最小二乘方滤波调整参量结果.     

Experiences with synthetic aperture focusing technique in the field

The detection and evaluation of defects in industrial components relies strongly on ultrasonic inspection techniques. Distance gain size (DG) or reference reflector methods can be improved concerning their localization, signal-to-noise ratio and sizing accuracy by the synthetic aperture focusing technique (SAFT). To obtain a high quality image, parameters like focal probe versus contact technique probe, achieved resolution or features of SAFT images compared with B-scan images are discussed. The implementation of SAFT in a CAD environment allows us to present stacked 2D reconstructions dynamically. On a cladded testblock with half-penny shaped cracks the advantage of combining CAD with SAFT is shown. A 3D SAFT example finalizes the overview of two decades of experience in applying this technique.     

Synthetic aperture imaging using sources with finite aperture: deconvolution of the spatial impulse response

A method for ultrasonic synthetic aperture imaging using finite-sized transducers is introduced that is based on a compact, linear, discrete model of the ultrasonic measurement system developed using matrix formalism. Using this model a time-domain algorithm for deconvolution of the transducer's spatial impulse responses (SIRs) is developed that is based on a minimum mean square error (MMSE) criterion. The algorithm takes the form of a spatiotemporal filter that compensates for the SIRs associated with a finite-sized transducer at every point of the processed image. A major advantage of the proposed method is that it can be used for any transducer, provided that its associated SIRs are known. This is in contrast to the synthetic aperture focusing technique (SAFT), which treats the transducer as a point source. The performance of the method is evaluated with simulations and experiments, performed in water using a linear phased array. The results obtained using the proposed method are compared to those obtained with a classical time-domain SAFT algorithm. For a finite aperture source, it is clearly shown that the resolution obtained using the proposed method is superior to that obtained using the SAFT algorithm.     

Deconvolution of vibroacoustic images using a simulation model based on a three dimensional point spread function

Vibro-acoustography (VA) is a medical imaging method based on the difference-frequency generation produced by the mixture of two focused ultrasound beams. VA has been applied to different problems in medical imaging such as imaging bones, microcalcifications in the breast, mass lesions, and calcified arteries. The obtained images may have a resolution of 0.7–0.8 mm. Current VA systems based on confocal or linear array transducers generate C-scan images at the beam focal plane. Images on the axial plane are also possible, however the system resolution along depth worsens when compared to the lateral one. Typical axial resolution is about 1.0 cm. Furthermore, the elevation resolution of linear array systems is larger than that in lateral direction. This asymmetry degrades C-scan images obtained using linear arrays. The purpose of this article is to study VA image restoration based on a 3D point spread function (PSF) using classical deconvolution algorithms: Wiener, constrained least-squares (CLSs), and geometric mean filters. To assess the filters’ performance on the restored images, we use an image quality index that accounts for correlation loss, luminance and contrast distortion. Results for simulated VA images show that the quality index achieved with the Wiener filter is 0.9 (when the index is 1.0 this indicates perfect restoration). This filter yielded the best result in comparison with the other ones. Moreover, the deconvolution algorithms were applied to an experimental VA image of a phantom composed of three stretched 0.5 mm wires. Experiments were performed using transducer driven at two frequencies, 3075 kHz and 3125 kHz, which resulted in the difference-frequency of 50 kHz. Restorations with the theoretical line spread function (LSF) did not recover sufficient information to identify the wires in the images. However, using an estimated LSF the obtained results displayed enough information to spot the wires in the images. It is demonstrated that the phase of the theoretical and the experimental PSFs are dissimilar. This fact prevents VA image restoration with the current theoretical PSF. This study is a preliminary step towards understanding the restoration of VA images through the application of deconvolution filters.     

三次相位板波前编码系统彩色图像恢复的迭代算法

采用三次相位板进行景深延拓的波前编码系统得到非对称的点扩展函数.为了获得最终清晰的彩色图像.研究了一种基于广义极小残差法(GMREs)的迭代算法,结合Tikhonov规整化方法,并利用多通道处理过程对中间图像进行左卷积恢复.为了消除恢复图像边界的振铃效应,推导了新的光学成像过程数学模型,该模型采用反镜像边界条件并利用直积近似对卷积核进行处理.模拟数据的分析表明,采用多通道处理过程对彩色图片进行恢复时,新的算法在给出精确的反巷积结果的同时能有效地抑制噪声的放大;实验结果显示,较之经典的维纳滤波恢复结果,新算法能够更好的消除边界的振铃和图像边缘的振动波纹.     

Delay-Doppler deconvolution image formation for multiple targets in waveguide environment

The ranges and radial velocities of targets are interesting parameters in the active sonar, which can't be observed directly but can only be estimated. Firstly, by making use of the delay-Doppler model of multi-targets in waveguide multipath environment, one finds that sample cross-ambiguity function is a two-dimensional(2D) convolution of the auto-ambiguity function of the transmitted signal with the generalized target reflectivity density, which is a 2D convolution of the spread function of channel with the reflectivity density as well. Secondly,from the perspective of information theory, an iterative deconvolution algorithm named R-L(Richardson-Lucy) is derived based on minimum Csiszar discrimination criterion. Finally, the blurs caused by both of the transmitted signal and channel are removed by 2D deconvolution of the expectation of sample cross-ambiguity function, 2D image and then parameters of time-delay and Doppler is estimated sequentially. Results of both numerical simulation and sea experimental data processing verify the feasibility and effectiveness of R-L deconvolution algorithm, which effectively improves the resolution and precision of the time delay and Doppler estimation, when compared to the classical match filtering and Wiener filtering.     

基于维纳滤波的迭代盲解卷积

盲解卷积是在两个卷积因子未知的情况下,通过卷积结果来获知卷积因子的。不考虑噪声,针对高斯模糊图像,在初始估计点扩展函数之后,利用维纳滤波的方法进行频域迭代盲解卷积,达到图像恢复的目的。实验表明,恢复的图像纹理比较清晰,边缘有所改善,主观视觉效果显著。该算法提高了分辨率。     

Ultrasonic synthetic aperture focusing using planar pulse-echo transducers

The lateral resolution of digital data from the planar (unfocused) pulse-echo transducers used in conventional ultrasonic inspections can be improved using the synthetic aperture focusing technique (SAFT).

For practical applications it is important to minimize the level of sidelobes (artefacts) introduced by SAFT, without significant loss of resolution. This may be achieved by the inclusion of a suitable aperture-weighting function in the SAFT algorithm, combined with a synthetic aperture size related to the width of the transducer beam-spread.

The properties of the resulting optimized SAFT algorithm are quantified using experimental data from a series of artificial flaws (slots) of different sizes.     

增量维纳滤波法在波前探测解卷积中的应用

一个常规的自适应光学系统通常包含三个重要环节：波前探测、波前校正和波前重构。因此对系统的技术要求非常高,造成系统复杂,成本昂贵。基于哈特曼-夏克波前探测的图像解卷积处理就是“事后”处理的一种,它省去了波前校正环节,使常规的自适应光学系统得到简化,降低了系统成本。其基本原理为对瞬时波面进行短时间曝光探测,同时记录相应的短时间曝光图像,来进行解卷积处理。将增量维纳滤波法应用于基于哈特曼-夏克波前探测的解卷积中,并对室内模拟点源情况下的三组畸变光斑图像数据进行了解卷积恢复处理。结果表明,将增量维纳滤波法应用于基于波前探测的解卷积是完全可行的,在室内模拟点源情况下,恢复的图像可以达到衍射极限分辨力。与维纳滤波相比,它扩展了噪声抑制因子的选取范围,在噪声抑制因子选取不正确的情况下,仍能得到比维纳滤波更好的结果。     

基于微通道板的二维位置灵敏像探测器X射线成像研究

研制了一套基于微通道板的二维位置灵敏像探测器,开展了X射线成像实验研究。对游标卡尺构建的狭缝成像显示,狭缝图像清晰可见,图像的边缘分辨率能达到约500 μm。对组合样品进行了成像测试,在X射线机靶流为1 mA的条件下,像探测器系统在1 ms内即可实现成像,且成像结果清晰。采用中值滤波与锐化、盲反卷积、维纳滤波等方法,对游标卡尺的成像结果进行了初步的图像处理。结果显示,三种方法都能改善图像的质量,其中维纳滤波法效果最优,相比原图边缘分辨率提高了约7%。     

The application of maximum entropy processing to the deconvolution of coupling patterns in NMR

The deconvolution of J-coupling patterns in NMR by iterative maximum entropy processing is demonstrated. Both the in-phase and the antiphase coupling patterns are considered. The deconvolution of the coupling pattern, either for one value of the coupling constant (1D J deconvolution) or for a range of coupling constants (2D J deconvolution) is shown. It is demonstrated that the method can be used for improving the signal-to-noise ratio for known coupling patterns by removing the coupling structure, as well as for extracting coupling constants from an unknown spectrum. Examples are shown both in 1D NMR and in slicewise processing of 2D spectra.     

Single-Channel FMCW-Radar-Based Multi-Passenger Occupancy Detection Inside Vehicle

In this paper, we provide the results of multi-passenger occupancy detection inside a vehicle obtained using a single-channel frequency-modulated continuous-wave radar. The physiological characteristics of the radar signal are analyzed in a time-frequency spectrum, and features are proposed based on these characteristics for multi-passenger occupancy detection. After clutter removal is applied, the spectral power and Wiener entropy are proposed as features to quantify physiological movements arising from breathing and heartbeat. Using the average means of both the power and Wiener entropy at seats 1 and 2, the feature distributions are expressed, and classification is performed. The multi-passenger occupancy detection performance is evaluated using linear discriminant analysis and maximum likelihood estimation. The results indicate that the proposed power and Wiener entropy are effective features for multi-passenger occupancy detection.     

气体自发瑞利-布里渊散射的理论分析及压强反演

气体压力是描述体系状态的重要参数,许多物理、化学性质都与压力有关.传统侵入式的压力测量方法会对气体状态产生干扰,影响测量精度,因此需要一种无扰式的测量方法.本实验测量了压强为2,4和6 atm(1 atm=1.01325×10~5 Pa)下加入气溶胶的N2在90?散射方向的自发瑞利-布里渊散射光谱,利用卷积后的Tenti S6模型对测量光谱进行直接拟合,拟合得到的压强值总体误差小于6.0%,求和归一化的均方根误差总体小于6.5%;利用理想的Tenti S6模型对经维纳滤波器反卷积处理后的测量光谱进行拟合,拟合得到的压强值误差总体小于5.0%,求和归一化的均方根误差总体小于6.0%.通过对两种方法的详细对比,发现压强低于2 atm时,对测量光谱进行反卷积处理在一定程度上可以消除仪器函数的影响,提高测量光谱的准确性,其光谱拟合效果和压强反演精度要优于卷积光谱.而在压强高于2 atm的情况下,卷积光谱的拟合效果和压强反演精度要优于反卷积光谱.     

红外连续激光反射镜热畸变的有限元分析

 基于热传导方程和热弹方程，利用有限元分析方法，就光斑尺寸远大于或接近于基底材料热扩散长度的情况以及反射镜在固定或自由的边界条件下，分别计算了硅和石英两种基底材料多层膜红外连续激光反射镜的最大温升、最大形变及最大热应力，并探讨了它们随光斑尺寸的变化规律。结果表明：在自由边界条件下，反射镜表面的最大轴向位移与光斑半径之间近似为线性关系；而在固定边界条件下，反射镜的最大热应力与光斑半径之间近似为线性关系；反射镜的夹持状态对最大轴向位移及最大热应力的影响随着光斑尺寸的增加而增强；在相同的入射激光光源及相同的边界条件下，硅镜具有较低的温升值或较高的抗热损伤阈值，而石英镜具有较好的抗热畸变特性。     

Semi-sparse deconvolution robust to uncertainties in the impulse responses

The received signal in ultrasonic pulse-echo inspection can be modeled as a convolution between an impulse response, or prototype, and the reflection sequence that is the impulse characteristic of the inspected object. Deconvolution aims at approximately inverting this process to improve the time resolution so that the overlap between echoes from closely spaced reflectors becomes small. For the relatively high contrast reflection sequences often found in non-destructive testing applications, semi-sparse deconvolution algorithms will typically yield better results than the classical Wiener filter solution. However, this requires a prototype that is a good representative for all echo responses found in the signals. Since, in practice, the prototype is often chosen by the operator directly from the inspection data, the prototype may incidentally be a bad representative for modeling the echoes for the remaining part of the object. Because of the sensitivity to deviations in the prototype this can yield deconvolution results with poor reproducibility. This paper presents a new semi-sparse deconvolution algorithm that is robust to deviations in the prototype. The new robust algorithm is based on a modification of an earlier presented non-robust semi-sparse algorithm. The robustness is obtained by including a stochastic model of the variations in the prototypes to the signal model when deriving the algorithm. Experiments performed using simulated data verify that the robust algorithm is less sensitive to deviations in the prototypes compared to the non-robust version of the algorithm and show that the proposed algorithm yields better estimates than its non-robust version and the Wiener filter in scenarios for which the algorithm was derived. Results using real ultrasonic data further show that the algorithm can be useful in practical scenarios where similar deconvolution results are required from slightly different echoes.     

A method of improving overall resolution in ultrasonic array imaging using spatio-temporal deconvolution

In this paper a beamforming method for ultrasonic array imaging is presented that performs both spatial and temporal deconvolution based on a minimum mean square error (MMSE) criteria. The presented MMSE receive mode beamformer performs a regularized inversion of the propagation operator for the ultrasonic array system at hand. The MMSE beamformer accounts for the transmit and receive processes, defined in terms of finite array element sizes, transmit focusing laws and electrical transducer characteristics. The MMSE beamformer is compared to the traditional delay-and-sum (DAS) beamformer with respect to both resolution and signal-to-noise ratio. The two algorithms are compared using both simulated and measured data. The simulated data was obtained using ultrasonic field simulations and the measured data was acquired using a linear phased array imaging wire targets in water. The results show that the MMSE beamformer has superior temporal and lateral resolution compared to DAS. It is also shown that the MMSE beamformer can be expressed as a filter bank, which enables parallel processing at high frame rates.     

三反射式柱面光学系统设计及优化

针对传统单片柱透镜和柱面反射镜成像光束不理想以及视场通常小于1°,提出并设计了一种三反射式柱面结构.对柱面光线追迹及单片柱面镜成像进行了深入分析,分别设计了三反射式圆柱面和二次曲线柱面系统,提出了一种基于抛物柱面镜理想线聚焦的新型像差优化方法,使其在子午面方向各视场调制传递函数得到最佳优化,并达到成像光谱仪等在狭缝方向卜高空间分辨率要求.其子午面总视场均达到了3°,在45 lp/mm分辨率条件下,边缘视场子午面方向的调制传递函数分别优于0.2和0.6.