NOTES : A Model of Gas Exchange for Hyperpolarized Xe(129) Magnetic Resonance of the Lung
NOTES : A Model of Gas Exchange for Hyperpolarized Xe(129) Magnetic Resonance of the Lung
背景知识:
Gas exchange is the essential function of the lung. In general, a lung can be viewed as a porous medium(多孔介质) consisting of capillary(毛细管) circuits with blood flowing inside. The blood flows are separated from the air spaces by several layers of tissues, including epithelium(组织上皮细胞), endothelium(内皮),and interstitium(小间隙).These tissues together are called air–blood barrier.(空气血液栅栏).At equilibrium, gas molecules in the alveolar spaces(肺泡空间) are under constant exchange with those dissolved in the barrier and capillary blood.
Remark: 气体由肺泡扩散浓度自高处向低处,平衡状态时,以常数速率和溶解在组织和毛细管血管的气体交换.
Almost all pulmonary diseases can be attributed to deficient gas exchange or delivery in the lung.
气体交换不足或者过量均会导致肺部疾病,联系5个肺部关键参数。
Prob: How measure gas exchange?
Ans: However, none of the established imaging techniques provides direct measurement of gas exchange!!!
没有直接的方法,间接实验。
Exp.1 computerized tomography ( 断层摄影技术 ) measures tissue density
Exp.2 magnetic resonance imaging of hyperpolarized (HP)
He (超极化氦气) images the air spaces in the lung.
Exp.3 MR of HP Xe is capable of providing direct measurements of gas exchange in the lung. As a contrast
agent, not only does xenon yield dramatically enhanced
MR signals in the air spaces, but it also dissolves into lung tissue and blood. As the dissolved xenon in blood follows the same physiological pathways of the normal blood gases (i.e., O2,CO2) (参看下面的图)
量化:
Quantification of dissolved xenon dynamics will lead to a quantified understanding of lung function.
(How? M.R. experiment)
原理:
A particular feature of xenon that permits such a study, is that the xenon dissolved into human lung exhibits two large chemical shifts from the resonance frequency(共振频率) of the free xenon gas—one at 197 ppm(parts per million), for xenon in lung tissue and blood plasma(血浆) (TP xenon), the other at 217 ppm, for the xenon in the red blood cells (RBC xenon)
RBC :Red Blood Cell
TP : Tissue and Plasma(组织和血浆)
Xenon ----> Dissolved Xenon +Free Xenon
Dissolved Xenon = TP Xe + PBC Xe
197ppm: xenon in lung tissue and blood plasma(血浆)
217ppm: the xenon in the red blood cells
核磁共振具体实验技术暂不知。
CSSR has been used frequently, over the past years in various lung diseases, including fibrosis(纤维变性) and emphysema (肺气肿).
Despite the previous applications of CSSR, there has not been a satisfactory theory to interpret xenon uptake dynamics for both dissolved xenon peaks in the lung.
density M_d within the septum(隔膜), perpendicular(垂直的) to the blood flow
M_d is density of dissolved xenon
M_f is the density of free xenon gas (at 0 ppm) in the air spaces
D is the diffusion coefficient of dissolved xenon
\lambda is the Ostwald solubility of xenon in lung parenchyma(软细胞组织)
Prob: How to define the boundary value conditions?
Apply Fourier’s separation of variable obtain series solution:
Where T is the xenon-exchange time constant in the lung
(That is Par.4)
Proof:
。。。。。。。。
初边值问题的数值解(张文生book):
级数解渐进:
Assume:
S_A is the total surface area of air space
V_g is total volume of the air spaces in the lung
The dissolved xenon signal is proportional to M_d *S_A/2
The free xenon gas signal is proportional to M_f * V_g
The normalized signal distribution S_d(x, t) for dissolved xenon can be written as (WHY? Normalized)
和 (4)式进行对比 compare with (4).
S_d= M_d SVR / 2M_f
S_A / V_g is the surface-area-to-volume ratio
(That is Par.1 SVR)
S_d1(t): The Total Signal From The Tissue
(组织中xenon 信号)
It is not difficult to derive (6)--->(7),.calculated as the spatial integral of S_d in Eq.6 over the two regions from (0) to (\delta) and from (d–\delta) to (d).
Where,b is the normalization factor (dimensionless), and \delta/d is what we call barrier-to-septum ratio.
(That is Par.2 BSR)
The dissolved xenon signal from the blood (‘‘blood xenon’’) is more difficult to calculate owing to flow.
S_d2(s): The Total Signal From Blood (No-flow & flow )
(总血液中Xenon 信号)
No flow : integral of S_d at (\delta ,d-\delta)
How is consider blood flow? (Key)
The total xenon signal from the blood is:
Where t_X is pulmonary capillary transit time
(Par.5 运输时间)
To calculate the xenon signal in tissue and blood plasma, e.g., TP xenon signal S_TP(t) at 197 ppm, and signal of RBC xenon S_RBC(t) at 217 ppm, we let h denote the fraction of RBC xenon relative to total xenon in blood
S_PBC = h * S_d2
S_blood plasma =(1-h) * S_d2
S_TP = S_d1 + (1-h) * S_d2(组织和血浆信号)
The hematocrit (HCT) can be calculated from h using
(Par.3 血细胞比容)
Where \lambda_RBC and \lambda_P are Ostwald solubilities of xenon in RBCs and plasma, respectively.
Remark : How calculate h ? MR of HP Xe
上面建立的理论模型,实验信号可观测到(不同的肺部疾病),接下来拟合需要的五个参数与正常的对比下???
文章选取了肺部处于四种不同的状态。
Remark:
COPD (Chronic Obstructive Pulmonary Disease)慢性阻塞肺病
Fibrosis 纤维病变
Anemia 贫血
正常健康人的参数:
\lambda 0.2
\delta 2
d 10
SVR 250
T 40
\eta 0.5
t_X 1.6
NOTES : A Model of Gas Exchange for Hyperpolarized Xe(129) Magnetic Resonance of the Lung的更多相关文章
- Notes(一)
Numerous experimental measurements in spatially complex systems have revealed anomalous diffusion in ...
- RFC3261--sip
本文转载自 http://www.ietf.org/rfc/rfc3261.txt 中文翻译可参考 http://wenku.baidu.com/view/3e59517b1711cc7931b716 ...
- ML—随机森林·1
Introduction to Random forest(Simplified) With increase in computational power, we can now choose al ...
- pm剩余要看的内容
1. thermal profiles tuning, customize thermal daemon for different thermal usage model 2.gas gauge d ...
- .net平台的rabbitmq使用封装
前言 RabbitMq大家再熟悉不过,这篇文章主要整对rabbitmq学习后封装RabbitMQ.Client的一个分享.文章最后,我会把封装组件和demo奉上. Rabbitmq的运作 从下图可以看 ...
- R语言︱XGBoost极端梯度上升以及forecastxgb(预测)+xgboost(回归)双案例解读
XGBoost不仅仅可以用来做分类还可以做时间序列方面的预测,而且已经有人做的很好,可以见最后的案例. 应用一:XGBoost用来做预测 ------------------------------- ...
- Asp.NetCore轻松学-实现一个轻量级高可复用的RabbitMQ客户端
前言 本示例通过对服务订阅的封装.隐藏细节实现.统一配置.自动重连.异常处理等各个方面来打造一个简单易用的 RabbitMQ 工厂:本文适合适合有一定 RabbitMQ 使用经验的读者阅读,如果你还没 ...
- android的USB MTP && USB CDC/USBnet(ECM, NCM, ACM) && USB gardget
MTP的全称是Media Transfer Protocol(媒体传输协议),它是微软公司提出的一套媒体文件传输协议.早在智能手机普及前,数码相机和MP3播放器等都使用了MTP的前身PTP(Pictu ...
- EIT: where is it now and what lies ahead?
EIT: where is it now and what lies ahead? Electrical impedance tomography (EIT) is an emerging clini ...
随机推荐
- XMPP聊天客户端环境搭建
1.服务器选择:ejabberd,具体安装过程请参考:http://blog.csdn.net/linhanmin/article/details/9876819 2.客户端配置: 采用xmppfra ...
- PHP+Mysql无限分类的方法汇总
无限分类是个老话题了,来看看PHP结合Mysql如何实现.第一种方法这种方法是很常见.很传统的一种,先看表结构表:categoryid int 主键,自增name varchar 分类名称pid in ...
- DOM4J介绍与代码示例【转载】
DOM4J是dom4j.org出品的一个开源XML解析包.Dom4j是一个易用的.开源的库,用于XML,XPath和XSLT.它应用于Java平台,采用了Java集合框架并完全支持DOM,SAX和JA ...
- bash 统计文件行数
#假设文件名是:fortest.gtf declare -i fileLinesfileLines=`sed -n '$=' fortest.gtf`echo $fileLines #-------- ...
- poj 2676 Sudoku ( dfs )
dfs 用的还是不行啊,做题还是得看别人的博客!!! 题目:http://poj.org/problem?id=2676 题意:把一个9行9列的网格,再细分为9个3*3的子网格,要求每行.每列.每个子 ...
- Repeater实例应用
在实际开发过程中,涉及到数据绑定,分页,以及一对多展示数据时,遇到这样的需求我们怎么解决呢?下面以帖子展示来逐一说明. 帖子主要由两部分组成,第一部分是发帖人的原创内容部分,第二部分是用户评论部分,这 ...
- [HDU 1695] GCD
GCD Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submiss ...
- C#中哈希表与List的比较
简单概念 在c#中,List是顺序线性表(非链表),用一组地址连续的存储单元依次存储数据元素的线性结构. 哈希表也叫散列表,是一种通过把关键码值映射到表中一个位置来访问记录的数据结构.c#中的哈希表有 ...
- UVALive 3713 Astronauts (2-SAT,变形)
题意: 有A,B,C三种任务,每个人必获得1个任务,大于等于平均年龄的可以选择A和C,小于平均年龄的可以选择B和C.这些人有一些是互相讨厌的,必须不能执行同任务,问能否安排他们工作?若行,输出任意一组 ...
- nginx.conf配置
在此记录下Nginx服务器nginx.conf的配置文件说明, 部分注释收集与网络. #运行用户 user www-data; #启动进程,通常设置成和cpu的数量相等 worker_processe ...