茂县土地岭大熊猫走廊带

Size: px

Start display at page:

Download "茂县土地岭大熊猫走廊带"

绵卢
5 years ago
Views:

1 中国组织工程研究第 16 卷第 30 期出版 Chinese Journal of Tissue Engineering Research July 22, 2012 Vol.16, No.30 基于 CT 图像三维电子线放疗剂量的计算方法 ** 翁邓胡, 徐海荣 A three-dimensional dose calculation method for electron beam irradiation based on CT images Weng Deng-hu, Xu Hai-rong 文章亮点 : 提出了一种基于 CT 图像的三维电子线放疗剂量计算方法, 该方法能够更精确的计算患者体内任意位置的剂量, 直观显示剂量分布 Abstract BACKGROUND: Electron beam irradiation can effectively cure superficial and eccentric tumors. Combined with CT images and three-dimensional dose algorithm, we can also realize the precise quantification of tumor dose using electron beam irradiation, which is more advantageous to tumor local controlling. OBJECTIVE: To study and implement a three-dimensional electron beam dose calculation method based on CT images. METHODS: The CT images of patients with nasal cavity and ethmoid sinus cancer were collected. Microsoft Visual C++ was used to parse CT images for related information, and then the human body emulational matrix was built for the dose model calculation, including oblique incidence correction and tissue inhomogeneity correction. The sagittal and coronal CT images in the three-dimensional display were reconstructed using VTK technique. The dose distribution curves displayed on cross-sectional, sagittal and coronal images were drawn by the dose distribution calculation methods. RESULTS AND CONCLUSION: According to the object-oriented analysis and design concept, we developed a high-visualized treatment planning software with interactive interface that is easy to operate. Besides, it could calculate accelerator output dose based on given target area dose. High doses line could be a very good conformal target area, and it could minimize the dose on the normal tissues and key organs, thereby protecting the tissues and organs. This planning software we developed realized the dose visualization in patients undergoing electron beam radiation therapy, and visually displayed the dose distribution, which is of significance for further improving the curative effects of radiation treatment and subsequent research of accurate radiotherapy. Weng DH, Xu HR. A three-dimensional dose calculation method for electron beam irradiation based on CT images.zhongguo Zuzhi Gongcheng Yanjiu. 2012;16(30): [ 摘要背景 : 电子线放射治疗能有效治愈浅表和偏心的肿瘤, 结合 CT 图像及三维剂量算法实现肿瘤的剂量精确定量, 有利于肿瘤的局部控制率目的 : 研究和实现一种基于 CT 图像的三维电子线放射治疗的计算方法方法 : 采集临床鼻腔筛窦肿瘤患者的 CT 图像, 通过 Visual C++ 开发工具对 CT 图像进行 DICOM 解析, 获取相关的数据, 建立仿真人体矩阵信息, 用以剂量模型的计算, 包括斜入射校正和组织不均匀性校正采用 VTK 技术重建出 CT 图像的矢状面和冠状面, 根据剂量模型计算的剂量值, 生成横断面矢状面以及冠状面的等剂量分布曲线结果与结论 : 根据面向对象的分析和设计理念, 开发了交互式界面易操作可视化程度高的治疗计划软件 ; 给定靶区剂量计算得到加速器输出剂量 ; 高剂量线能很好的适形靶区, 并使正常组织和关键器官所受剂量很小, 从而得到保护实现电子线放射治疗患者体内剂量的可视化, 直观显示剂量分布, 对提高患者的放射治疗质量和精确放疗的后续研究具有指导意义 Department of Biomedical Weng Deng-hu, Studying for master s degree, Department of Biomedical wengdenghu@ 163.com Corresponding author: Xu Hai-rong, Doctor, Associate professor, Master s supervisor, Department of Biomedical Supported by: Three Kinds of Guangdong Technology Outlays, No. C30601*, * doi: /j.issn Received: Accepted: 关键词 : 电子线放射治疗 ;CT 图像 ; 计算模型 ; 靶区剂量 ; 剂量计算 ; 等剂量分布翁邓胡, 徐海荣. 基于 CT 图像三维电子线放疗剂量的计算方法 [J]. 中国组织工程研究,2012,16(30): [ ISSN CN /R CODEN: ZLKHAH 5625

2 翁邓胡, 等. 基于 CT 图像三维电子线放疗剂量的计算方法南方医科大学生物医学工程学院, 广东省广州市翁邓胡, 男, 1986 年生, 浙江省衢州市人, 汉族, 南方医科大学在读硕士, 主要从事肿瘤放射物理的研究 163.com 通讯作者 : 徐海荣, 博士, 副教授, 硕士生导师, 南方医科大学生物医学工程学院, 广东省广州市中图分类号 :R318 文献标识码 :B 文章编号 : (2012) 收稿日期 : 修回日期 : ( /M C) 引言现代放射治疗的计划设计由早期的传统 X 射线片, 逐步发展为以 CT MRI PET 等图像融合技术的肿瘤定位和治疗随诊, 利于获取患者更多的解剖数据和肿瘤病灶信息 [1-2] 此外, 新兴的计算机可视化技术, 使临床医师进一步直观显示患者肿瘤病灶任意角度及剖面的三维空间结构形态, 并作定性和定量分析 [3-4], 提高了计划设计和治疗的针对性及准确性另一方面, 高能电子束作为放射治疗中一种重要的治疗形式, 相对于 X(γ ) 放射线, 电子线具有独特的深度剂量分布特性, 即 : 表面剂量高, 在特定深度处达到剂量最大值, 高剂量峰后又迅速锐减, 形成一快速剂量跌落区 [5], 类似于质子重离子治疗的 Bragg 峰特征因此, 利用高能电子束治疗浅表的和偏心的颈部淋巴结胸壁和内乳淋巴瘤鼻腔筛窦肿瘤具有很好地保护肿瘤后面的正常组织和重要器官的临床优势由于电子线特殊的物理特性, 其在剂量计算方面有着自身的剂量计算模型和相关剂量学参数校正方法本研究的电子线剂量计算基本思想是 : 首先实测特定加速器机型的电子线剂量学参数 [ 百分深度剂量 (PDD) 离轴比有效源皮距斜入射校正因子等 ], 根据患者的三维 CT 数据进行有效深度的校正, 然后根据设定的靶区剂量计算处方剂量及人体内任一点的剂量值, 得出三维剂量分布并以等剂量线的方式显示 1 材料和方法 1.1 电子线剂量计算模型水模体内任意一点的剂量计算通常电子线治疗病变 ( 表浅和偏体位 ) 时, 采用固定源皮距 (SSD) 及单野照射, 这样可使靶区剂量均匀, 靶区后正常组织和器官剂量很小图 1 为任意点 P 的剂量计算示意图, d 为任一计算点 P 的深度,r 为 P 点对应中心轴上 P 0 的离轴距离, 还有一个靶区位置, 深度为 d, 计划设计中的靶区系指肿瘤的临床灶亚临床灶及肿瘤可能侵犯的范围, 照射中患者器官运动, 以及由于日常摆位治疗中靶位置和靶体积变化等因素引起的扩大照射的组织范围 [6] Figure 1 Dose calculation at arbitrary point P with fixed source skin distance irradiation 图 1 固定源皮距照射时任意点 P 的剂量计算示意图对于电子束, 在均匀模体内某一点 P 处的剂量 D p 可表示为 : DT D PDD d FSZ SSD OAR d r FSZ PDD d, FSZ, SSD ',, ',, p d d d (1) 其中,D T 为靶区剂量, 由医生指定 ;FSZ d 为射野中心轴上 d 深度处的射野宽度 ;PDD(d, FSZ d,ssd) 是在靶区位置处测得的 PDD 值 ; PDD(d,FSZ d,ssd) 代表所求 P 点深度处的 PDD 值 ;OAR(d,r,FSZ d ) 表示在 FSZ d 下 P 点深度处离射野中心轴距离 r 的离轴比其次, 公式 (1) 是在理想情况下的剂量计算公式, 在水模体中测量的数据和计算的数据在应用到具体患者时, 应对公式 (1) 得到的任意点剂量作斜入射校正和组织不均匀性校正斜入射校正 : 由于电子束过程中会遇到患者皮肤表面弯曲, 或摆位条件的制约, 造成限光筒的末端与患者皮肤表面不能很好的平行这会引起剂量分布的不均匀, 因此在计算时必须对其进行修正 [7] 不均匀性校正 : 水或组织通过阻止碰撞的方式完成对电子线的吸收, 组织电子密度对电子线的吸收影响很大, 对于骨肺或气腔等不均匀性组织能显著改变电子线剂量分布形式介质的电子密度 n e 定义为单位体积内的电子数目, 有 : NAZ ne M (2) 式中 ρ 为介质的物理密度,N A 为阿伏伽德罗常数,Z 为介质的有效原子序数,M 为介质原子的质量 ; 同时,(N A Z/M) 记为单位质量的电子 P.O. Box 1200, Shenyang cn.zglckf.com

3 翁邓胡, 等. 基于 CT 图像三维电子线放疗剂量的计算方法数目表 1 给出了人体不同组织中单位质量的电子数表 1 人体不同组织中单位质量的电子数 ( 单位 :10 23 electrons/g) Table 1 The number of electronsf per unit mass in different body tissues (unit:10 23 electrons/g) Water Lung Muscle Fat Bone Air 对于该模型中作者采用等效厚度系数法 (CET) 校正,CET 值可以通过不均匀组织对水的相对电子密度计算得到 [8], 有 : n N CET n N tissue tissue tissue water water water 其中 N tissue,n water 分别为单位质量组织和水的电子数这样通过 CET 值就可以计算不均匀组织后某一点深度 Z 处的有效深度 Z eff :Z eff =Z-t(1-CET),t 为非均匀组织的厚度最后, 将不均匀校正因子用 C CET 表示, 其值用平方反比定律校正计算因此, 经过以上两个步骤, 便可得出电子线在患者体内任意一点 P' 的剂量计算模型见公式 (4): 2 2 SSDeff z SSDeff z Dp' Dp OF, z SSDeff g z SSDeff z eff 1.2 电子线剂量计算的实现针对电子线剂量计算的特点, 并结合临床放射治疗计划中对电子线剂量计算的要求本研究为实现剂量计算的精确度, 作者采集了同一单位放疗中心的实验数据, 这里主要有两类数据 : 第一, 基础的剂量学数据, 这部分数据通过实测美国 Varian 23EX 加速器得到 ; 第二, 实际病例的数据, 主要是患者的 CT 图像, 这里采用鼻腔筛窦肿瘤患者的的头部 CT 数据实现剂量计算程序的具体流程如下 :1 利用自己编写的 Visual C++ 程序读出 46 张 DICOM 3.0 标准的 CT 序列图像的象素值, 并将其转化为 CT 值 CT 值反映了 X 射线穿过组织被吸收后的衰减值, 因此, 根据临床放射影像学查找出各 CT 值对应的组织类型, 再得到相应的单位质量电子数, 进而计算各象素点的有效深度 2 为提高剂量计算的速度, 又不至于缺失计算精确度, 将象素点合并, 变成体元大小, 即每个体元是 4 4 象素点 ; 此时一体元的电子密度值就为 16 个象素的代数平均电子密度 3 提取照射范围内的人体表面信息和靶区位置信息 4 对人体三维空间进行剖分, 建立仿真人体矩阵信息在人体笛卡尔坐标系上划分网 (3) (4) 格, 将 x,y,z 轴各分为 n x,n y,n z 点, 并将人体矩阵的每个节点附上该点所对应的人体点信息 : 如该点的组织类型, 该点是否为表面点, 该点的人体组织密度等 5 根据等效深度处方剂量以及剂量学数据等输入参数的确定, 利用公式 (1) 至公式 (4) 的计算模型进行剂量计算, 并以剂量分布图的形式在 CT 图像上显示剂量值 2 结果与分析为临床放疗工作的需要, 剂量计算的实现包括很多剂量学数据病例资料 ( 文本与图像 ) 的输入, 治疗机参数和射野设置, 计算数据管理以及剂量分布图的绘制, 因此, 作者设计了交互式界面易操作, 可视化程度高的 MFC 文档视图计划软件该软件采用面向对象的分析和设计理念, 用 Visual C++ 作开发工具, 同时选取易于对数据管理和扩展性好的 SQL Server 作后台数据库, 用于存放治疗信息 2.1 剂量计算软件设计根据软件工程设计思想, 软件包含了加速器数据与剂量学数据管理剂量参数和照射野参数设置建立三维坐标系 CT 图像处理选择剂量参考点人体内各点剂量计算显示等剂量分布加速器数据与剂量学数据管理录入 Varian 23EX 加速器的电子线治疗能量范围, 限光筒类型以及各能量下的治疗深度, 同时管理电子线相关的剂量学测量数据, 包括 PDD 表输出因子表离轴比表和斜校正因子表等剂量参数和照射野参数设置临床医生根据肿瘤性质期别大小位置等设定参考剂量 ( 靶区剂量 ), 分次数以及单次剂量, 参考剂量线百分比 ; 在治疗机上要设定照射野参数, 其中有照射方法 ( 电子线为源皮距照射 ), 照射体位 ( 仰卧位俯卧位侧位 ), 射野类型 ( 对称野或非对称野 ), 机架角度准直器角度肿瘤定位中心离床面的距离等 CT 图像处理利用 Visual C++ 对 CT 图像进行 DICOM 解析 [9], 通过标记符 (Tag) 提取有用的数据元素, 用以显示和剂量计算, 主要有患者姓名层厚 SOP Instance UID 象素大小图像大小和 Image Position (HFS) 等 ; 采用手工和自动分别勾画计划靶区, 人体表面轮廓, 这样就能求得体内任意点到各轮廓线与射线交点间的距离 ; 二维图像的窗宽窗位调整缩放移动旋转反色处理等三维空间坐标系的建立采用笛卡尔坐标系, 以 CT 图像横断面横轴从左到右为 X 轴正方向, 纵轴从上到 ISSN CN /R CODEN: ZLKHAH 5627

www.crter.org 翁邓胡, 等. 基于 CT 图像三维电子线放疗剂量的计算方法下为 Y 轴正方向, 以指向患者头部的方向为 Z 轴正方向, 并以第一张 CT 图像的原点为系统三维坐标的原点, 这样 CT 图像上任一点的三维坐标及象素值都可以获得 2.1.

6 三维重建利用 VTK 技术直接从二维体数据 (MC 算法 ) 进行三维重建, 并根据矢状面冠状面参数设置重建出矢状面和冠状面, 计算矢状面的肿瘤靶区位置所在层面上的剂量分布, 并显示剂量线 2.

4 翁邓胡, 等. 基于 CT 图像三维电子线放疗剂量的计算方法下为 Y 轴正方向, 以指向患者头部的方向为 Z 轴正方向, 并以第一张 CT 图像的原点为系统三维坐标的原点, 这样 CT 图像上任一点的三维坐标及象素值都可以获得等剂量分布显示根据上述剂量计算模型, 在预定的参考剂量 ( 靶区剂量 ) 点基础上, 进行斜入射校正以及组织不均匀性校正, 对剂量进行优化计算最后在 CT 层面显示等剂量线以及鼠标所在点的剂量三维重建利用 VTK 技术直接从二维体数据 (MC 算法 ) 进行三维重建, 并根据矢状面冠状面参数设置重建出矢状面和冠状面, 计算矢状面的肿瘤靶区位置所在层面上的剂量分布, 并显示剂量线 2.2 处方剂量的计算在电子线放疗过程中, 保证靶区剂量正确的一个重要环节就是准确的计算处方剂量, 它将肿瘤剂量 (D T ) 换算成标准水模体内每个射野的中心轴上最大剂量点处的剂量 D m 这种转换可以使用 Khan 计算方法, 通过直接调用治疗机的剂量学数据和射野布置以及靶区剂量条件等进行加速器输出剂量计算, 计算 [10] 公式为 :D m (MU)=D T /(PDD FF SSD f ) 此时, 电子线在固定源皮距照射方式下,SSD f (SSD 因子 ) 为 1,FF 输出因子包括射野输出因子 (S c ) 和模体散射因子 (S p ), 因为源皮距照射时,S c 下的表面射野大小就是 S p 的等中心处的射野大小, 所以只需测得给定射野在模体中与参考射野的输出剂量比即是 FF 输出因子 (S c, p ) 为实现处方剂量的更精确的计算, 避免因输出因子变化导致治疗剂量不足而引起肿瘤复发临床使用的 FF 输出因子并不全使用测量数据, 还要通过实际测量数据进行拟合处理, 作者探讨了电子束限光筒内各射野大小与其输出因子有一定规律性, 通过 Matlab 曲线拟合得出公式 (5): bx dx y ae ce 公式 (5) 中 y 为输出因子 FF, x 表示方野边长,a b c d 为随能量和射野尺寸变化的拟合系数, 见表 2 [11] 对于矩形野, 采用射野长径法规定方野大小, 即实际射野大小计算时将矩形的长径视为方野边长, 方便快速计算但特别对于小野时, 由于缺少横向散射, 在测量输出因子时, 应考虑 Z max 下 PDD 的改变 (5) 2.3 等剂量分布显示图 2 为横断面等剂量曲线此实例中, 肿瘤最大深度为 4 cm, 图中绿色线包围的轮廓线就是 PTV(target 区域 ), 本实验在能量 12 MeV, SSD=100 cm, 射野面积 FSZ=10 cm 10 cm, 体积元大小为 0.26 cm 0.26 cm 0.5 cm 的条件下计算, 并给予靶区剂量为 cgy 可以发现,3 240 cgy 剂量线 ( 品红色线 ) 已基本覆盖 target 区域, 且剂量分布均匀, 能剂量积分较小 ;360 cgy 剂量线 ( 黑色线 ) 两侧向外扩张, 主要是由于深度越深, 导致电子线能量降低, 最终引起电子散射角度的减小同时, 得到的等剂量分布能很好保护深部的正常脑组织, 特别是关键器官 ( 脑干脊髓等 ) 的受限剂量, 达到了医生预设的临床要求 Figure 2 Dose distribution display on the cross-sectional image 图 2 横断面等剂量分布图作者采用 VTK 工具包实现 CT 图像重建, 并提取患者 CT 图像的矢状面和冠状面, 获取不同剖面的人体矩阵信息然后, 根据之前计算的各个横断层面的剂量值, 以及 CT 扫描层厚, 输入矢状面或冠状面的重建参数, 计算出对应的矢状面或冠状面的剂量值, 绘制并显示剖面的等剂量曲线, 这样就可以从不同视角直观评估剂量分布的均匀性图 3 显示了计算得到的矢状面等剂量分布图从中可见,90% 剂量线 ( 品红色线 ) 已基本包络 target 区, 并使周围重要器官 ( 眼睛脊髓 ) 得到了保护表 2 电子束输出因子拟合公式系数表 Table 2 The fitting coefficient of electron beam output factor Electronic linear energy a b c d 6 MeV MeV MeV Figure 3 Dose distribution display on the sagittal image 图 3 矢状面等剂量分布图 5628 P.O. Box 1200, Shenyang cn.zglckf.com

5 翁邓胡, 等. 基于 CT 图像三维电子线放疗剂量的计算方法 3 讨论及结论研究的三维电子线放射治疗剂量计算方法, 已经实现了基于测量数据的改进矩阵剂量计算, 它可以精确地计算三维水模体的剂量分布信息, 再根据患者解剖结构的不均匀性对剂量分布进行修正, 从而获取人体内精的剂量分布电子线剂量计算的治疗计划软件设计, 包括 CT 图像窗宽窗位调整人体表面轮廓提取处方剂量的计算等剂量曲线的绘制矢状面和冠状面的重建输出治疗计划报告等功能, 主要完成了下列有意义的探索和尝试 : 首先通过等效厚度系数法来作人体不均匀性校正 [12], 成功实现了电子线在患者体内任一点剂量的计算, 剂量计算模型符合临床应用 ; 其次对于 CT 断层图像信息不能直接使用的问题, 建立了仿真人体矩阵模型, 人体网格化后的矩阵单元存储了人体所在位置的组织信息, 包括人体上各点的坐标组织类型组织密度 CT 像素值及剂量值等, 方便在剂量计算时能准确获取人体信息 ; 第三, 电子线的等效方野计算不同于 X 射线 [13], 本文选择了射野长径法用于等效方野的计算, 并通过输出因子拟合模型对 PDD 进行修正, 实现了处方剂量的精确计算 ; 最后利用改进矩阵剂量模型实现了人体内任意点的剂量计算, 并通过 Matlab 与 VC 结合提高了计算效率的同时, 也简化了计算结果的处理过程, 完成等剂量曲线的绘制研究发现, 电子线的浅表及偏心的放射治疗运用单野照射可实现靶区均匀的剂量分布, 满足临床要求针对电子线的特性及临床个体化病例的需求 [14], 本研究实现了三维电子线放射治疗剂量计算方法, 并设计出界面交互友好的电子线放射治疗计划软件, 可以在任意横断面矢状面上显示等剂量分布, 同时也直观反映了肿瘤的受照剂量和周围器官的保护情况, 定量地计算出横断面图像上任意点的剂量值由于目前临床上电子线大部分采用常规照射, 而本研究三维电子线剂量计算方法能够更精确的计算患者体内任意位置的剂量, 直观显示剂量分布, 对提高患者的放射治疗质量和精确放疗的后续研究具有一定的指导意义 4 参考文献 [1] Karlsson M, Nyholm T, Karlsson MG, et al. MR-only procedures for improved overall precision in radiotherapy. Int J Radiation Oncology Biol Phys. 2009;75(3):S656. [2] Morarji K, Fowler A, Vinod SK, et al. Impact of FDG-PET on lung cancer delineation for radiotherapy. J Med Imaging Radiat Oncol. 2012;56(2): [3] Yu JL, Huang ZH, Fang CH, et al. Zhongguo Linchuang Jiepouxue Zazhixue. 2009;27(3): 俞金龙, 黄宗海, 方驰华, 等. 直肠肿瘤数字解剖模型构建及三维可视化的研究 [J]. 中国临床解剖学杂志,2009,27(3): [4] Dzhugashvili M, Pichenot C, Dunant A,et al. Surgical clips assist in the visualization of the lumpectomy cavity in three-dimensional conformal accelerated partial-breast irradiation. Int J Radiat Oncol Biol Phys. 2010;76(5): [5] Hogstrom KR, Almond PR. Review of electron beam therapy physics. Phys Med Biol. 2006;51(13):R [6] Pérez-Romasanta LA, Lozano-Martín E, Velasco-Jiménez J,et al. CTV to PTV margins for prostate irradiation. Three-dimensional quantitative assessment of interfraction uncertainties using portal imaging and serial CT scans.clin Transl Oncol. 2009;11(9): [7] Demir B, Okutan M, Cakir A,et al. The effect of oblique electron beams to the surface dose under the bolus. Med Dosim. 2009;34(4): [8] Higgins PD, Gerbi BJ, Khan FM. Application of measured pencil beam parameters for electron beam model evaluation. Med Phys. 2003;30(4): [9] Law MY. A model of DICOM-based electronic patient record in radiation therapy. Comput Med Imaging Graph. 2005;29(2-3): [10] Ding GX, Cygler JE, Yu CW,et al. A comparison of electron beam dose calculation accuracy between treatment planning systems using either a pencil beam or a Monte Carlo algorithm. Int J Radiat Oncol Biol Phys. 2005;63(2): [11] Yang QQ,Zhuang MS,Wang JY,et al. Zhongguo Yiliao Shebei. 2010;25(3): 杨秋权, 庄梅生, 王继宇, 等. 放射治疗电子线处方剂量计算系统的研究 [J]. 中国医疗设备, 2010,25(3): [12] Schneider W, Bortfeld T, Schlegel W. Correlation between CT numbers and tissue parameters needed for Monte Carlo simulations of clinical dose distributions. Phys Med Biol. 2000; 45(2): [13] Verhaegen F, Mubata C, Pettingell J,et al. Monte Carlo calculation of output factors for circular, rectangular, and square fields of electron accelerators (6-20 MeV). Med Phys. 2001;28(6): [14] Xia HS,Wu Y,Zhang YM,et al.xiandai Zhongliu Yixue.2011; 19(6): 夏火生, 武勇, 张玉海, 等. 乳腺癌改良根治术后胸壁电子线照射研究 [J]. 现代肿瘤医学, 2011,19(6): 来自本文课题的更多信息 -- 基金声明 : 广东省科技三项经费计划项目 (C30601; ) 作者贡献 : 实验设计实验实施为第一作者, 实验评估为通讯作者, 资料收集为第一作者第一作者成文, 通讯作者审校, 第一作者对文章负责利益冲突 : 课题未涉及任何厂家及相关雇主或其他经济组织直接或间接的经济或利益的赞助 ISSN CN /R CODEN: ZLKHAH 5629

PCA+LDA 14 1 PEN mL mL mL 16 DJX-AB DJ X AB DJ2 -YS % PEN

PCA+LDA 14 1 PEN mL mL mL 16 DJX-AB DJ X AB DJ2 -YS % PEN 21 11 2011 11 COMPUTER TECHNOLOGY AND DEVELOPMENT Vol. 21 No. 11 Nov. 2011 510006 PEN3 5 PCA + PCA+LDA 5 5 100% TP301 A 1673-629X 2011 11-0177-05 Application of Electronic Nose in Discrimination of Different