生物组织光声粘弹显微成像

提出一种反演生物组织粘弹信息的新型无损光声粘弹显微成像方法,它是以强度调制激光作为激发源,通过检测光声（Photoacoustic,PA）信号的相位重建组织粘弹特性分布的成像方法.实验利用不同浓度的琼脂样品来验证光声粘弹显微测量中相位随浓度变化的依赖关系.利用埋有头发丝的琼脂样品来测试这种显微方法的成像分辨率.利用具有不同粘弹性的离体生物组织来验证系统的成像能力.实验结果表明,这种新方法能够高分辨率和高对比度地重建出具有不同粘弹性的生物组织的光声粘弹显微图像,有望实现组织结晶类病变水平的显微在体检测.     

光声结构与功能成像技术研究进展

光声成像技术利用短脉冲激光激发产生光声信号,可重建出组织的光吸收分布图像,它结合了纯光学成像的高对比度和纯声学成像的高分辨率特性.光声成像技术不仅能够有效的刻画生物组织结构,还能够精确实现无损功能成像,为研究生物组织的形态结构,生理、病理特征,代谢功能等提供了全新手段.本文简要分析了光声信号产生的机理,总结报道了目前实验室几套典型的成像系统及其最新应用进展,指出光声成像作为一种新型的生物医学成像方法,可望引发生物医学影像领域的一次革新.     

基于光谱编码的血管内光声组分成像

光声成像突破了传统的光学成像和超声成像在生物组织成像领域的困境,该技术基于光声(Photoacoustic,PA)效应,脉冲激光激励下的生物组织产生超声信号,超声信号被接收后,通过反投影算法将其携带的时间信息和强度信息转化为能够反映生物组织吸收结构和分布的可视化图像。基于不同生物组织的光吸收差异,当激发光强度均匀且稳定时,光声成像反映的就是该物质对于该波长光的吸收特性。本文中,我们基于导管式的血管内光声断层扫描平台结合多波长激发的光声成像算法开发了基于光谱编码的血管内光声组分成像系统,实现了在离体血管斑块中脂质组分的定量成像,高分辨获得了脂质核心的大小形态和边界信息,表征了斑块内的脂质相对含量。     

生物组织弹性成像方法及其新技术进展

弹性是一种描述物质物理意义的重要参数,在描述物质在热力学和动力学的变化过程中有着重要的意义。在医学上,弹性的变化往往和病变联系在一起。然而,绝大多数生物组织在他们的力学特性上所表现出的复杂性并不是弹性模量一项参数就可以完全表述的,在对于他们的粘弹性表征和流变学行为的描述中,粘滞性往往和弹性一样的重要。现在被广泛用来对生物组织机械特性表征的成像技术是弹性成像,其基本原理是给组织施加一个激励,组织会产生一个响应,而该响应的分布结合技术的处理方法,可以反映出其弹性模量等力学属性的差异。本文介绍了生物组织常见的弹性成像方法:超声弹性成像,磁共振弹性成像以及光学相干弹性成像;详细阐述了新发展起来的技术-光声弹性成像和光声粘弹成像,并讨论分析其应用前景。     

用光声光谱技术测量鼻咽癌细胞的光谱吸收特性

为了解人鼻咽癌细胞的光谱吸收特性,用双光束归一化光声光谱技术测量三株人鼻咽癌细胞(SUNE-1,CNE-2,CNE-1)的归一化吸收光谱和光谱吸收系数.结果表明(1)在波长4000～6000A范围内,这三株人鼻咽癌细胞的光谱吸收系数大小为SUNE-1＞CNE-2＞CNE-1;(2)在波长4200A处平均相对光声信号值为SUNE-11.36,CNE-21.20,CNE-11.10,三者之间的差异具有显著性(P值＝0.012).人鼻咽癌细胞株SUNE-1,CNE-2和CNE-1具有不同的光声光谱特性.     

叶绿素的脉冲激光光声研究

用光声法研究了叶绿素的光声信号与入射激光能量及激光照射时间的关系。用液体流动光声池与高效液相色谱联机,在610nm处得到叶绿素a、b的光声色谱图,研究了叶绿素a/b比值与激光辐照时间的函数关系。     

基于环形阵列探测器的快速光声成像

构建一套基于环形阵探测器的快速光声成像系统用于生物组织的结构成像。该系统以环形阵探测器探测光声信号,采用八通道的采集系统采集光声信号,再利用有限场滤波反投影算法重建光声图像。利用埋有铅笔芯的琼脂样品来测试该系统的分辨率,利用离体猪眼和在体老鼠头部血管成像来验证系统的成像能力。实验结果表明,该系统能方便快速地实现生物组织的结构成像,有望实现早期乳腺癌的临床检测应用。     

光声成像及其在生物医学中的应用

光声成像是一种新近迅速发展起来、基于生物组织内部光学吸收差异、以超声作媒介的无损生物光子成像方法,它结合了纯光学成像的高对比度特性和纯超声成像的高穿透深度特性的优点,以超声探测器探测光声波代替光学成像中的光子检测,从原理上避开了光学散射的影响,可以提供高对比度和高分辨率的组织影像,为研究生物组织的结构形态、生理特征、代谢功能、病理特征等提供了重要手段,在生物医学临床诊断以及在体组织结构和功能成像领域具有广泛的应用前景.对光声成像技术的机理、光声成像技术和方法、光声图像重建算法以及光声成像在生物医学上的应用情况作一个简单介绍,希望有助于推动我国在该领域的科研和开发应用工作的迅速发展.     

一体化的光声-荧光显微镜用于活体肿瘤成像

报道了一种利用单一波长激发的同时产生光声和荧光信号的显微成像系统,本成像系统具有超高的成像分辨率(<6μm)。借助外源的造影剂在近红外的吸收特性,利用光声-荧光显微成像系统对活体肿瘤进行光声/荧光成像。实验结果表明,光声-荧光显微镜在早期肿瘤的成像和检测等方面具有潜在的应用价值。因此,通过研究和选择适当的双模态造影剂,该系统在不同病理模型中可以提供更准确的组织信息及生理参数。     

基于非线性热扩散效应的亚衍射极限光声二次谐波显微成像技术

本文提出了一种基于非线性热扩散效应的光声二次谐波显微SH-PAM成像技术,用于实现亚衍射极限光声成像。生物组织受到强度调制的高斯激光束辐射时,组织吸收光子形成高斯分布的温度场,由于热扩散系数非线性热效应引起的非线性光声PA效应,从而产生光声二次谐波信号。模拟和试验结果均表明,重建后的光声二次谐波成像的横向分辨率超过了传统光学成像分辨率。本文通过仿体样品验证了该方法的可行性,并且对人表层皮肤细胞进行了成像,以证明其对生物样品的成像能力。该方法扩展了传统光声成像的范围,为超分辨成像开辟了新的可能性,为生物医学成像和材料检测提供了新的方法。     

Constitutive equations for the lung tissue

The mechanical behavior of the lung tissue (expressed by its constitutive equations) has considerable influence on the normal and pathological function of the lung. It determines the stress field in the tissue, thus affecting the impedence and energy consumption during breathing as well as the localization of certain lung diseases. The lung tissue has a complex mechanical response. It arises from the tissue's structure--a cluster of a very large number of closely packed airsacks (alveoli) and air ducts. Each of the alveoli has a shape of irregular polyhedron. It is bounded by the alveolar wall membrane. In the present study, a stochastic approach to the tissue's structure will be employed. The density distribution function of the membrane's orientation in space is considered as the predominant structural parameter. Based on this model the present theory relates the behavior of both the alveolar membrane and that of its liquid interface to the tissue's general constitutive properties. The resulting equations allow for anisotropic and visco-elastic effects. A protocol for material characterization along the present model is proposed as well. The methodology of the present theory is quite general and can be similarly used with other structural models of the lung tissue (e.g., models in which the effect of the alveolar ducts is included).     

A transient heating technique for the measurement of thermal properties of perfused biological tissue

Knowledge of tissue thermal transport properties is imperative for any therapeutic medical tool which employs the localized application of heat to perfused biological tissue. In this study, several techniques are proposed to measure local tissue thermal diffusion by heating with a focused ultrasound field. Transient as well as near steady-state heat inputs are discussed and examined for their suitability as a measurement technique for either tissue thermal diffusivity or perfusion rate. It is shown that steady-state methods are better suited for the measurement of perfusion; however the uncertainty in the perfusion measurement is directly related to knowledge of the tissue's intrinsic thermal diffusivity. Results are presented for a transient thermal pulse technique for the measurement of the thermal diffusivity of perfused and nonperfused tissues, in vitro and in vivo. Measurements conducted in plexiglas, animal muscle, kidney and brain concur with tabulated values and show a scatter from 5-15 percent from the mean; measurements made in perfused muscle and brain compare well with the nonperfused values. An estimate of the error introduced by the effect of perfusion shows that except for highly perfused kidney tissue the effect of perfusion is less than the experimental scatter. This validation of the tissue heat transfer model will allow its eventual extension to the simultaneous measurement of local tissue thermal diffusivity and perfusion.     

Quantitative measurements of porphyrin pigments in tissues via photoacoustic spectroscopy

We investigated the possibility of using photoacoustic spectroscopy as an analytical technique for the quantitative measurement of injected porphyrins in tissues. Samples of liver excised from treated (i.e., injected 24 h in advance with 100 mg/kg Photofrin II) mice and control mice were lyophilized and reduced to powder; then, about 16-mg powder samples were compacted to equal volumes inside the photoacoustic cell. Amplitude and phase spectra were measured in the range 280-760 nm. From these data we computed the photoacoustic absorbance spectra; they were suitably normalized in order to account for differences in hemoglobin concentration among the livers of different mice. The absorbance difference spectra (treated minus control) were computed in the region of wavelengths above 450 nm, where porphyrins and hemoglobin exhibit major differences. Finally, by estimating the value of the thermal diffusion length of the powdered sample and those of the extinction coefficients of the most relevant Photofrin II components in the spectral region considered, we were able to evaluate the local drug concentration and determined a value (260 micrograms/g of wet tissue) that is in the range expected for the dose of Photofrin II injected.     

Biorheology and fluid flux in swelling tissues. I. Bicomponent theory for small deformations, including concentration effects

A theory for the rheological behavior and fluid flux in swelling tissues under small deformations is presented. Tissues are considered as bicomponent solid-fluid mixtures. Concentration effects are included. The driving forces (body, surface and interactive), are discussed and their constitutive relationships to the tissue's deformation are specified. Mass and momentum balance equations are developed for each component and for the tissue as a whole. The concept of swelling stress emerges from the theory as an anisotropic generalization of the commonly used swelling pressure. It is shown to be a measure of the total chemical potential combining both mechanical and concentration effects. The theory shows that concentration effects modify the tissue's bulk stiffness in a manner consistent with experimental observations.     

Statistical finite element method for real-time tissue mechanics analysis

The finite element (FE) method can accurately calculate tissue deformation. However, its low speed renders it ineffective for many biomedical applications involving real-time data processing. To accelerate FE analysis, we introduce a novel tissue mechanics simulation technique. This technique is suitable for real-time estimation of tissue deformation of specific organs, which is required in computer-aided diagnostic or therapeutic procedures. In this method, principal component analysis is used to describe each organ shape and its corresponding FE field for a pool of patients by a small number of weight factors. A mapping function is developed to relate the parameters of organ shape to their FE field counterpart. We show that irrespective of the complexity of the tissue's constitutive law or its loading conditions, the proposed technique is highly accurate and fast in estimating the FE field. Average deformation errors of less than 2% demonstrate the accuracy of the proposed technique.     

Stochastic Dynamic Model for Estimation of Rate Constants and Their Variances from Noisy and Heterogeneous PET Measurements

Tissue heterogeneity, radioactive decay and measurement noise are the main error sources in compartmental modeling used to estimate the physiologic rate constants of various radiopharmaceuticals from a dynamic PET study. We introduce a new approach to this problem by modeling the tissue heterogeneity with random rate constants in compartment models. In addition, the Poisson nature of the radioactive decay is included as a Poisson random variable in the measurement equations. The estimation problem will be carried out using the maximum likelihood estimation. With this approach, we do not only get accurate mean estimates for the rate constants, but also estimates for tissue heterogeneity within the region of interest and other possibly unknown model parameters, e.g. instrument noise variance, as well. We also avoid the problem of the optimal weighting of the data related to the conventionally used weighted least-squares method. The new approach was tested with simulated time–activity curves from the conventional three compartment – three rate constants model with normally distributed rate constants and with a noise mixture of Poisson and normally distributed random variables. Our simulation results showed that this new model gave accurate estimates for the mean of the rate constants, the measurement noise parameter and also for the tissue heterogeneity, i.e. for the variance of the rate constants within the region of interest.     

Axial speed of sound is related to tendon's nonlinear elasticity

Axial speed of sound (SOS) measurements have been successfully applied to noninvasively evaluate tendon load, while preliminary studies showed that this technique also has a potential clinical interest in the follow up of tendon injuries. The ultrasound propagation theory predicts that the SOS is determined by the effective stiffness, mass density and Poisson's ratio of the propagating medium. Tendon stiffness characterizes the tissue's mechanical quality, but it is often measured in quasi-static condition and for entire tendon segments, so it might not be the same as the effective stiffness which determines the SOS. The objectives of the present study were to investigate the relationship between axial SOS and tendon's nonlinear elasticity, measured in standard laboratory conditions, and to evaluate if tendon's mass density and cross-sectional area (CSA) affect the SOS level. Axial SOS was measured during in vitro cycling of 9 equine superficial digital tendons. Each tendon's stiffness was characterized with a tangent modulus (the continuous derivative of the true stress/true strain curve) and an elastic modulus (the slope of this curve's linear region). Tendon's SOS was found to linearly vary with the square root of the tangent modulus during loading; tendon's SOS level was found correlated to the elastic modulus's square root and inversely correlated to the tendon's CSA, but it was not affected by tendon's mass density. These results confirm that tendon's tangent and elastic moduli, measured in laboratory conditions, are related to axial SOS and they represent one of its primary determinants.     

Constitutive behavior of ocular tissues over a range of strain rates

The constitutive behavior of bovine scleral and corneal tissues is measured in tension and compression, at quasi-static and moderate strain rates. Experiments are conducted at strain rates up to about 50 strain per second by a pneumatic testing system developed to overcome noise and measurement difficulties associated with the time dependent test of low impedance materials. Results for the tissues at room and the natural bovine body temperatures are similar and indicate that ocular tissue exhibits nonlinear stiffening for increasing strain rates, a phenomena termed rate hardening. For example, at a tensile strain rate of 29/s, corneal tissue is found to develop 10 times the stress that it does quasi-statically at the same strain. Thus, conventional constitutive models will grossly underpredict stresses occurring in the corneo-scleral shell due to moderate dynamic events. This has implication to the accuracy of ocular injury models, the study of the stress field in the corneo-scleral shell for glaucoma research and tonometry measurements. The measured data at various strain rates is represented using the general framework of a constitutive model that has been used to represent biological tissue mechanical data. Here it is extended to represent the measured data of the ocular tissues over the range of tested strain rates. Its form allows for straightforward incorporation in various numerical codes. The experimental and analytical methods developed here are felt to be applicable to the test of human ocular tissue.     

Tracer measurements of non-equilibrium volumes of distribution

According to conventional theory the product of the transport flowrate and the mean transit time of a tracer through a system yields the equilibrium volume of distribution for the tracer, regardless of tracer kinetics or space geometry. Experimental results do not support this notion. The influence of measurement time on the volume measured with a bolus technique is addressed using systems theory to analyze a tissue-impedance form of the Sangren-Sheppard model. Assymptotic solutions show that the volume estimates are governed by a time constant, tau, related to diffusion in the tissue, to tissue capacity, and to wall permeability, and by a dimensionless ratio, f, describing a relation of tau to vascular transport time. A third parameter, g, describing the relative contributions to overall resistance to diffusion of effective permeability and of limited diffusivity in the tissue, is shown to be of less importance. The derived tau is similar to but not equivalent to the often cited "characteristic time". The "equilibrium" volume of distribution is defined as that which would be measured if equilibrium were allowed to establish. The "non-equilibrium" volume of distribution is defined as that which would be measured given finite times and is shown to approach the "equilibrium" volume as such times increase. Tracer equilibration is not required to accurately measure the "equilibrium" volume. When there is no flow limitation (f much less than 1) a measurement time of tau (plus vascular transit time) would yield a "non-equilibrium" volume only 33% of the "equilibrium" volume; a time of 2 tau would yield 55%; a time of 10 tau would yield effectively the total equilibrium volume. Finite diffusivity in tissue and permeability restrictions can have significant effects on the proportion of the volume measured.     

Quantitation of time- and frequency-resolved optical spectra for the determination of tissue oxygenation.

The recent development of near-infrared time- and frequency-resolved tissue spectroscopy techniques to probe tissue oxygenation and tissue oxygenation kinetics has led to the need for further quantitation of spectroscopic signals. In this paper, we briefly review the theory of light transport in strongly scattering media as monitored in the time and frequency domains, and use this theory to develop algorithms for quantitation of hemoglobin saturation from the photon decay rate (delta log R/delta t) obtained using time-resolved spectroscopy, and from the phase-shift (theta) obtained from frequency-resolved, phase-modulated spectroscopy. To test the relationship of these optical parameters, we studied the behavior of delta log R/delta t and theta as a function of oxygenation in model systems which mimicked the optical properties of tissue. Our results show that deoxygenation at varying hemoglobin concentrations can be monitored with the change in the photon decay kinetics, delta delta log R/delta t in the time-resolved measurements, and with the change in phase-shift, delta theta, in the frequency-resolved technique. Optical spectra of the adult human brain obtained with these two techniques show similar characteristics identified from the model systems.