正常角膜和圆锥角膜的特征提取.pdf

第 2 3 卷第 1 O期 2 0 1 5 年 1 O月 光学 精密工程 Op t i c s a n d P r e c i s i o n En g i n e e r i n g Vo 1 2 3 No 1 0 文章 编号1 0 0 4 9 2 4 X(2 0 1 5)1 0 2 9 1 9 0 8 正常角膜 和 圆锥角膜 的特征提取 余锦华 ，季春红，李添捷。，田 磊 ，黄一飞 ，汪源源 ，郑永平。(1 复旦大学 电子工程系，上海 2 0 0 4 3 3；2 医学影像计算及计算机协助介入重点实验室 上海 2 0 0 4 3 3；3 香港理工大学 跨学科生物医学工程部门，香港；4 首都医科大学 附属北京同仁医院 北京市 眼科研究所，北京 1 0 0 0 0 5；5 中国人民解放军总医院 眼科系，北京 1 0 0 8 5 3)摘要：基 于角膜测量仪器 C o r v i s S T采集的图像视频，提出提取新特征参数 以便 准确 区分 正常角膜 和圆锥角膜。首 先对 图像进行滤波、分割等预处理，检测 角膜上下边界，并计算前 角膜 曲率 值；用 小波变换 分析角膜 曲率变 化，获 取与角 膜运 动趋势相关 的特征，包括角膜运动 的整体趋势和角膜振动 的范数 和标准差。然后，基 于均方误差 最小化 法，提取特 征参 数，构建最优参数。最后，用支持 向量 机(S VM)对正 常角膜 和圆锥 角膜进行 分类。从 频率 的角度实施 的实 验显示 角膜 在基本运动趋势上存在 着振 动过 程。此外，提 出的参数优 于形变 幅度(DA)、峰值 距离(P D)等 传统参数，使准 确度、灵敏 度和特异性分别提 高了 1 O 2 ，5 7 和 6 9 。受试者工作特征 曲线(RO C)下 面积 为 0 9 4 8，接近 于 1。结果 显示 本文 方法 自动提取 的特征参 数可提高正常角膜和圆锥角膜 区分 的准确性，对临床诊断有辅助作用。关键词：角膜；圆锥 角膜；特征提取；最小均方误差算 法；C o r v i s S T；小波变换 中图分类 号：R 7 2 2 2 1；TP 3 9 1 文献标识码：A d o i：1 0 3 7 8 8 O P E 2 O 1 5 2 3 1 0 2 9 1 9 Ex t r a c t i ng f e a t u r e s f r o m n o r ma l c o r ne a s a n d k e r a t o c o nu s b a s e d o n wa v e l e t a n a l y s i s YU J i n h u a ，J I Ch u n h o n g ，LI Ti a n j i e。，TI AN L e i ，H UANG Yi f e i ，W ANG Yua n y u a n ZH ENG Yo ng p i ng。(1 De p a r t me n t o f El e c t r o n i c En gi n e e r i n g，Fu da n Un i v e r s i t y，S h a n gh a i 2 0 0 4 3 3，Ch i n a：2 Ke y L a b o r a t o r y o f Me d i c a l I ma g i n g C o mp u t i n g a n d C o mp u t e r As s i s t e d I n t e r v e n t i o n o f Sh a n g h a i，S h a n g h a i 2 0 0 4 3 3，Ch i n a；3 I n t e r d i s c i pl i n a r y Di v i s i o n o f Bi o me d i c a l En g i n e e r i n g，Ho n g Ko n g Po l y t e c h n i c，Ho n g Ko n g，C h i n a；4 B e Oi n g I n s t i t u t e o f Op h t h a l m o l o g y，T o n g r e n Ho s p i t a l，C a p i t a l Me d i c a l Un i v e r s i t y，Be i j i n g 1 0 0 0 0 5，C h i na；5 De p a r t me n t o f Op h t h a t mo t o g y，Ch i n e s e PL A Ge n e r a l Ho s p i t a l，Be i j i n g 1 0 0 9 5 3，Ch i n a)*CD r P 户 0 g a u t h o r，E-ma i l：j h y u f u d a n e d u c n Ab s t r a c t：On t he ba s i s o f v i d e o i ma ge c a p t u r e d b y t he c o r ne a m e a s u r i ng i ns t r ume nt Co r v i s ST，t h i s p a pe r pr o po s e s a n i de a t o i m p r ov e t h e a c c ur a c y o f d i s t i ng ui s hi n g n or ma l c o r ne a s f r o m k e r a t oc on i c c o t n e a s by e xt r a c t i ng n e w f e a t u r e pa r a me t e r s，Fi r s t l y，t he o r i gi n a l i m a ge s we r e p r e pr o c e s s e d b y f i l t e r i n g a n d s e g me n t i n g t o d e t e c t t h e u p p e r a n d l o we r b o u n d a r i e s o f t h e c o r n e a a n d c a l c u l a t e t h e c u r v a t u r e o f 收稿 日期：2 0 1 5 0 5 2 1；修订 日期：2 0 1 5 0 7 2 O 基金项 目：国家 9 7 3重 点 基 础研 究 发展 计 划 资 助项 目(No 2 0 1 5 C B 7 5 5 5 0 0)；国家 自然科 学 基 金 资 助项 目(N o 6 1 4 7 1 1 2 5，N o 8 1 2 7 1 0 5 2)；香港理工大学 内地联合监督计 划资助项 目(No G-UB 5 8)2 9 2 0 光学精 密工程 第2 3 卷 a nt e r i or c or n e a The n，t he c ha n ge of c o r n e a l c u r v a t u r e wa s a n a l y z e d by wa v e l e t t r a n s f or ma t i on me t h o d t O o b t a i n f e a t u r e s r e l a t e d t o t h e t r e n d o f c o r n e a l mo v e me n t，i n c l u d i n g t h e t r e n d o f t h e wh o l e c o r n e a l m o t i o n a s we l l t he n or m a nd t he s t a n da r d d e v i a t i o n o f c o r ne a l v i br a t i o n Fur t he r mo r e，t h e f e a t u r e pa r a m e t e r s we r e e xt r a c t e d i n s u c c e s s i on a nd t he o p t i ma l pa r a m e t e r wa s o bt a i ne d by t h e mi ni m um m e a n s q ua r e e r r o r a l go r i t h m The Sup po r t Ve c t o r M a c h i ne(S V M)wa s f i n a l l y a pp l i e d t o d i s t i nc t i o n of n or ma l c or ne a s f r o m ke r a t o c o n i c c o r ne a s The e x pe r i m e nt r e s u l t s o n t he f r e q ue nc y i n di c a t e t h a t t he r e a r e c or n e a l v i br a t i o ns a l o n g wi t h t he b a s i c m o ve m e nt p r oc e s s Be s i de s，t he p r o po s e d p a r a me t e r s a r e b e t t e r t h a n t r a d i t i o n a l p a r a me t e r s s u c h a s De f o r ma t i o n Amp l i t u d e(DA)，Pe a k Di s t a n c e(P D)a t t he hi g he s t c o nc a v i t y，wh i c h i mpr o v e s t he a c c u r a c y，s e ns i t i v i t y a n d s p e c i f i c i t y b y 1 0 2 ，5 7 a nd 6 9 ，r e s p e c t i v e l y Mo r e o v e r，t h e a r e a u n d e r t h e r e c e i v e r o p e r a t i n g c h a r a c t e r i s t i c c u r v e(ROC)i s 0 9 4 8，c l o s e t o u n i t y Th e a u t o ma t i c e x t r a c t e d f e a t u r e p a r a me t e r s i n t h i s p a p e r a r e a b l e t O i mp r o v e t he a c c u r a c y o f c l a s s i f i c a t i o n be t we e n no r ma l a nd k e r a t o c on i c c or n e a s a nd c o nt r i b ut e t o t h e c l i n i c a l d i agnose s Ke y wo r d s：c o r ne a；ke r a t o c o n i c c o r ne a；f e a t u r e e x t r a c t i o n；mi ni mum m e a n s q ua r e e r r o r a l go r i t hm；Cor v i s S T；wa ve l e t t r a ns f o r m 1 引 口 圆锥角 膜(Ke r a t o c o n i c C o r n e a s，KC)是 角膜 中央变薄并扩张，呈圆锥形特征的眼类疾病，多为 双侧性发病，会 导致不规 则散光 和视 力下 降口 。它是 一 种逐 步 扩展 的角 膜 疾 病_ 2 ，治 疗 的 方 法 主 要取决于严重程度，较轻时可进行屈光治疗，严重 时则需进行角膜移植手术。早期发现有助于更好 地选择治疗方式，也可以最大程度地避免并发症 的产生和发展，有效保护视觉功能。近年来，许多新技术用于测量角膜。主要包 括光学相干层析成像(Op t i c a l C o h e r e n c e To p o g r a p h y，OC T)_ 3“、眼前节 分析诊 断 系统 P e n t a cam_ 5 I 6l、眼反应 分 析 仪(O c u l a r Re s p o n s e An a l y z e r，0R A)【7 等。最近，德 国 Oc u l u s 公 司结合 诱导位移的空气脉冲和高速 S h e i mp f l u g照相机，推 出新 的角膜测量仪器 C o r v i s S Tl g 。利用其实 施测量时，先以恒定压力的空气脉冲对着角膜吹，使角膜逐渐 向内凹陷到达最低点，然后移 除空气 脉 冲，角膜 逐 渐反 弹恢 复 到 正 常 状 态。在此 过 程 中，S h e i mp f l u g照相机以 4 3 0 0 f r a me s 的速度记 录角膜 的移动过程，并在 内置控制面板上 实时显 示 l】。整 个形 变过 程 可 以分为 吹气 前、吹气 中和 吹气后 3个 状 态。除 了记 录 角 膜 的 形 变 进 程，C o r v i s S T也 用 于测 量 一 些有 利 于分 析 KC 的生 物力 学参 数。本 文 的研究 就 是基 于 C o r v i s S T 仪 器 的。C o r v i s S T仪器 可计算 出角膜 中心 厚度、眼 内压、第 一扁 平 时 间、峰值 距 离、最 大 凹 陷时 间 等 参数l】。已有研究报道指出，最大 凹陷时刻角膜 中心的形变幅度(DA)E 1 2 、峰值距离(P D)和 中心 曲率(C C c)l_】在 区分正常角膜和 KC时最具 准 确 性。近 期研 究表 明，前 角 膜 曲率 的测 量 在 屈 光 手 术_】、角膜硬度和眼内压的评估 中 1 6 具有重 要意义。此外，前角膜 曲率还反映了角膜形变度，在区分正常角膜和 KC时具 有关键作用L 1。然 而，C o r v i s S T没有直接提供角膜曲率变化参数，故对有关角膜 曲率变化参 数在 区分正常 角膜和 KC 中的作用 也未进 一 步展 开研究。因此，本文通过计算整个形变过程 中的前角 膜曲率半径和离散小 波分析。来获取角膜 中心 曲率半径信号 的高低频 成分并将其作 为特征参 数；然 后，基 于 均方 误 差 最 小化 方 法 提取 特 征，构 建最优 参数。最 后，用支 持 向量机(S u p p o r t Ve c t o r Ma c h i n e，S VM)c 1 9-z 1 对 角 膜 进 行 分 类，测 试 提 出参 数 的有效 性。2 前 角膜 曲率半径的计算 2 1 角膜 边界 的分 割 C o r v i s S T仪 器采 集 的角 膜 图像 中不 仅 存 在 斑点噪声，而且有眼睫毛、眼泪等引起 的杂质噪 第 1 O 期 余锦华，等：正常角膜和圆锥角膜的特征提取 声。针 对这 一 特 征，本 文 主要 利 用 杂 质 和 角 膜边 界间的位置和大小关 系进行去噪，然后采用基于 相位的分割方法进行边界的 自动检测。，算法 流程如图 1 所示。首先，运用线性结构元素对角 膜 图像进行膨胀 学闭运算，将 由于光照不均匀引 起的角膜断裂连接起来。接着，考虑杂质的面积 远小于角膜，故用面积阈值将其滤除。然后，生成 相位对称和反对称 图像 2 。最后，利用相位对称 图像拟合 中心线，基于 中心线 和相位反对称 图像 获得 角膜 的上 下边 界。图 1 角膜边界分割流程 图 Fi g 1 F l o w c h a r t o f c o r n e a l b o u n d a r y s e g me n t a t i o n 2 2曲率 半径 的估计 由于前 角膜 边 界 点是 离 散 的，如 果 采 用 定 义 的方法计算 曲率，则小 的波动将导致较大 误差。因此，本文采用 T h o ma s 提 出的基于面积误差最 小的方法估计 圆弧 中心和半径汹。假设 属于一 段 圆弧上 点 的坐标 为(，Y )(z ，Y )(z ，Y )；圆弧 中心 为(z ，Y )；半 径 为 R，则 目标 函数 为：N J=：=E n R 一7 c (z lz )+(一Y c)。-=1 (1)通过 函数 对 R，Y 分别求偏微分，获 得 R，z ，Y 的表 达式。这种 方 法 不 仅不 需 要 设 定初始值，而且不用迭代算法，很大程度地节省了 计算时间。3 最优参数提 取 3 1离 散小波 分析 小波滤波器是一种多尺度、多分辨率 的信号 分析工具，能够提取角膜 时变运动的动态特征。为 了观察角膜运 动，对中心点 曲率半径随时间变 化的信号进行小波分析，来检测角膜振动的快速 变 化。首先通过一组高通 和低通滤波器(H，L)将 原始信号分解成高频和低频成分。再将低频成分 通过高通和低通滤波器分解成高低成分，以此类 推，经过 4组高低通滤波器后，原始信号将被分解 成 5个 重要 成分 S ，S ，S。，S ，S 。然 后，对 各 个 节 点进行重构，获得重构信号，s r ，s r 。再将信号进行重组，构成 S h i g h 一 +s r 2+s r 3+；Sl o w一 ，如 图 2所 示。从 图 中可 以看 出，低频成分 S l o w基本反映了角膜受空气脉冲 作用先往下凹陷，到达最低点后再往上弹起恢复 到 自然状态的整体运动趋势；高频成分则反应了 角膜受空气脉冲、弹力、粘黏性等相互作用的振动 过 程。这一 发 现与 Z h a o l o n g Ha n等研 究 者 的 研 究结果相一致 2 ，后文从频率的角度对其进行了 解 释。2 5 0 宣2 0 0 鳖 1 5 0 譬 1 0 0 6 5 O 0 (a)原始信号(a)Or i gi na l s i gna l t m s (b)低频成分 l 善 砉 (c)高频成分(b)L o w-f r e q u e n c y c o mp o n e n t(c)Hi g h-f r e q u e n c y c o mpo n e n t 图 2 信 号重组 Fi g 2 Si gna 1 r e c o mbi na t i on 3 2 基 于 最小均 方误 差 的特征 选择 为了减少特征维数且使特征更具鲁棒性，特 征 选择 和 提 取 是 必 需 的 步 骤 o 。本 文 基 于 最 小 均 方误差，提 出新 的 特 征选 择 方 法。设 样 本 个数 为 N；统计的特征个数为 M；第 J个特征第 i 个值 为 叫。由于不 同特征具有不 同的单位和范 围，需要 各 自进行 归一 化，即：一：二 兰 d L i,j ma x(x J)-mi n(x J)其 中 一1，2，N；j 一1，2，M。(2)2 9 2 4 光学精密工程 第2 3 卷 面积(Ae r a)、截点(C u t o f f p o i n t)等参数。从表中 可以看出，C o m2 对应的曲线下面积达到最 大，为 0 9 4 8，此 时对 应 的 截点 为 一0 5 5 1。截 点 对 应 的 灵敏度和特异性分别为 7 8 6 和 9 5 7 。6 结 论 本文首先利用杂质和角膜边界的位置、大小 等先验关系，基于相位分割方法检测 出整个运动 过程 中的角膜边界，并计算 出了前角膜曲率半径 值；然后根据角膜受空气脉冲和弹力等作用的运 动特征，运用小波变换提取 出高低频成分。从频 参 考文 献：1 2 3 4 E 5 E 6-DUDAK0VA L，L I S K0VA P，j I RS 0VA K。I s c o p p e r i mb a l a n c e a n e n v i r o n me n t a l f a c t o r i n f l u e n c i n g k e r a t o c o n u s d e v e l o p me n t?J Me d i c a l Hy p o t h e s e s，2 0 1 5，8 4(5)：5 1 8 5 2 4 KRACHMER J H，F EDER R S，BE LI N M W Ke r a t o c o n u s a n d r e l a t e d n o n-i n f l a mma t o r y c o r n e a l t h i n n i n g d i s o r d e r s J S u r v O p t h a l mo l，1 9 9 4，28：2 93 3 2 2 S UL L A C，VU0NG L N，P RI CE L L，e t a 1 Co mp a r i s o n o f s p e c t r a l F o u r i e r d o ma i n o p t i c a l C O h e r e n c e t o mo g r a p h y i n s t r u me n t s f o r a s s e s s me n t o f n o r ma l ma c u l a r t h i c k n e s s I-J Re t i n a-T h e J o u r n a l 0 fRe t i n a l a n d Vi t r e o u s Di s e a s e s，2 0 1 0，3 0(2)：2 35 2 45 王立科，张佳 莹，田磊，等基 于光学 相干层 析气 冲 印压技术研究角膜生物力学特性 J 光学 精密工 程，2 0 1 5，2 3(2)：3 2 5 3 3 3 W ANG L K，Z HANG J Y，TI AN L，n Z 0CT b a s e d a i r j e t i n d e n t a t i o n f o r c o r n e a l b i o me c h a n i c a l a s s e s s me n t E J O p t P r e c i s i o n E n g ，2 0 1 5，2 3(2)：3 2 5 3 3 3 (i n Ch i n e s e)ANAY0L M A，GULER E，YAGCI R，e t a 1 Comp a r i s o n of c e nt r a l c o r ne a l t h i c kne s s，t hi nne s t c or ne a l t hi c kn e s s，a nt e r i o r c h a mbe r de pt h，a nd s i mu l a t e d k e r a t o me t r y u s i n g g a l i l e i，p e n t a c a m，a n d s i r i u s d e v i c e s I-J C o r n e a，2 0 1 4，3 3(6)：5 8 2 5 8 6 AM BROS I O R，CAI AD0 A L C，GUERRA F P，e t a 1 No v e l p a c h y me t r i e p a r a me t e r s b a s e d o n c o t 一 率的角度，论证 了角膜基本运动趋势存在着振动 过程。此外，还提取 了其他 的特征参数，包括形变 幅度、峰值距离、中心曲率半径等；接着，以最小均 方误差为 目标函数，利用基于共轭方向的 P o we l l 算法求解系数，提取出最优特征参数；最后，基于 S VM 算法对正常角膜和圆锥角膜进行分类。与已提出的经验证过的传统参数相 比，本文 所提 出参 数的特征分类准确率、灵敏度、特异性、阳性预测值和阴性预测值分别提高 了1 0 2 、5 7 、6 9 、3 5 和 2 1 7 。R OC 曲线 下 的 面 积为 0 9 4 8。实验 结果 验证 了本 文 参 数对 于 圆锥 角膜和正常角膜分类的有效性。7 8 9 1 O 1 1 1 2 n e a l t o mo g r a p h y f o r d i a g n o s i n g k e r a t o c o n u s J J o u r n a l o f Re fr a c t i v e Su r g e r y，2 0 1 1，2 7(1 O)：75 3 75 8 L UCE D A De t e r mi n i n g i n v i v o b i o me c h a n i c a l pr op e r t i e s o f t he c o r n e a wi t h a n o c ul a r r e s p on s e a n a l y z e r r J 1 J o u r n a l o f C a t a r a c t&Re f r a c t i v e S u r g e r y，2 0 0 5，3 1(1)：1 5 6 1 6 2 LUCE D，TAYLOR D Re i c he r t oc u l a r r e s p on s e a n a l y z e r m e a s u r e s c o r n e a l b i o me c h a n i c a l p r o p e r t i e s a n d I OP J R e i c h e r t Op h t h a l mi c I n s t r u me n t s，2 0 06 AM BROSI O R，RAM OS I，LUZ A，a1 Dy n a mi c ul t r a hi gh s pe e d Sc he i mpf l ug i ma gi ng f o r a s s e s s i n g c o r n e a l b i o me c h a n i c a l p r o p e r t i e s J Re v i s t a Br a s i l e i r a d e Of t al mo l o g i a，2 0 1 3，7 2(2)：99 1 0 2 HONG J X，XU J J，WEI A J，e t a 1 A n e w t o n o me t e r-t h e Co r v i s S T t o n o me t e r：c l i n i c a l c o mp a r i s o n wi t h n onc o nt a c t a n d Go l d ma nn ap pl a na t i on t o n o me t e r s J I n v e s t i g a t i v e Op h t h a l mo l o g y&s u a l S c i e n c e，2 0 1 3，5 4(1)：6 5 9 6 6 5 AMB R OS I O J R R，NO GUE I R A L P，C AL DA S D L nZ Eva l ua t i on o f c or ne a l s ha pe a n d bi ome c h a n i c s b e f o r e L A S I K J I n t e r n a t i o n a l Op h t h a l mo l o g y Cl i n i c s，2 0 1 1，5 1(2)：1 1 3 8 TI AN L，HUANG Y F，WANG L Q，e t a 1 Co r n e a l b i o me c h a n i c a l a s s e s s me n t u s i n g c o r n e a l v i s ua l i z a t i on s c he i mpf l u g t e c h no l o gy i n ke r a t oc oni c a n d n o r ma l e y e s F J J o u r n a l o f O p h t h a l mo l o g y，2 0 14，2 01 4：1 475 1 6 第 1 O 期 余锦 华，等：正常角 膜 和 圆锥 角膜 的特 征提 取 2 9 2 5 1 3 1 4 1 5 1 6 1 7 1 8 1 9 2 O AMBROS 1 0 J R R，CAL DAS D，RAM0S I，e t 口 Z Co r n e a 1 b i o me c h a n i c a l a s s e s s me n t u s i n g d y n a m i c u l t r a h i g h-s p e e d S c h e i mp f l u g t e c h n o l o g y n o n c o n t a c t t o n o me t r y(U HS-S T NC T)：p r e l i mi n a r y r e s u l t s C Ame r i c a n S o c i e t y o f C a t a r a c t a n d Re fr a c t i v e S u r g e r y Ame r i c a n S o c i e t y o f Op h t h a l mi c Admi n i s t r a t o r s (AS CRS-AS OA)Sympo s i um a nd Co n gr e s s，San Di e go，CA，2 01 1 ROBERTS C Bi o me c ha ni c a l c us t omi z a t i on：t he n e x t g e n e r a t i o n o f l a s e r r e f r a c t i v e s u r g e r y J J o u r n a l o f Ca t a r a c t&Re fr a c t i v e S u r g e r y，2 0 0 5，3 1(1)：2 5 TANG M L，LI Y，AVI LA M，e t a 1 M e a s u r i n g t o t a l c o r n e a l p o we r b e f o r e a n d a f t e r l a s e r i n s i t u k e r a t o mi l e u s i s wi t h h i g h s p e e d o p t i c a l c o h e r e n c e t o mo g r a p h y J J o u r n a l o f C a t a r a c t Re f r a c t i r e S u r g e r y，2 0 0 6，3 2(1 1)：1 8 4 3 1 8 5 0 LANZ A M，CENNAM 0 M，LAC CARI N0 S，e t n Z Ev a l u a t i o n o f c o r n e a l d e f o r ma t i o n a n a l y z e d wi t h s c h e i mp f l u g b a s e d d e v i c e i n h e a l t h y e y e s a n d d i s e a s e d o n e s J B i o Me d R e s e a r c h I n t e r n a t i o n al，2 01 4，20 14：7 48 67 1 H 0N Y，LAM A K CCo r n e a 1 d e f o r ma t i o n me a s u r e me n t u s i n g S c h e i mp f l u g n o n c o n t a c t t o n o me t r y J O p t o me t r y&V i s i o n S c i e n c e，2 0 1 3，9 0(1)：e 1 一 e 8 张浩，袁怡宝，张峰巴特沃思小波在表面形貌信 号分离 中的应 用 J 光 学 精 密 工程，2 0 1 0，1 8 (7)：1 6 6 1-I 6 6 7 ZHANG H，YUAN Y B，ZHANG F Ap p l i c a t i o n o f But t e r wo r t h wa ve l e t t o s u r f a c e t op o gr ap hy s i g n a l s e p a r a t i o n J O p t Pr e c i s i o n E n g ，2 0 1 0，1 8(7)：1 6 6 1-1 6 6 7 (i n Ch i n e s e)高恒振，万建伟，粘永健，等组合核 函数 支持 向量 机高光谱 图像融合分类 J 光学 精 密工程，2 0 1 l，1 9(4)：8 7 8 8 8 3 GA0 H Z H，W AN J W，NI AN Y J，口 Z F u s i o n c l a s s i f i c a t i o n o f h y p e r s p e c t r a l i ma g e b y c o m p o s i t e k e r n e l s s u p p o r t v e c t o r ma c h i n e J O p t Pr e c i s i o n En g，2 0 1 1，1 9(4)：8 7 8 8 8 3 (i n Ch i n e s e)M ELGANI F，BRUZZ0NE L Cl a s s i f i c at i o n o f hyp e r s pe c t r a l r e mot e s e ns i ng i ma g es wi t h s up po r t v e c t o r m a c h i n e s J I E E E T r a n s G e o s c i R e mo t e S e n s，2 0 0 4，4 2(8)：I 7 7 8 1 7 9 0 2 1 谭 熊，余旭初，张鹏强，等基 于多核支 持向量机 的 高光谱影像 非线 性混合像 元分 解 J 光 学 精 密 工 程，2 0 1 4，2 2(7)：1 9 I 2-I 9 2 0 TAN X，YU X CH，ZHANG P Q，e t a 1 No n l i n e a r mi x e d p i x e l d e c o mp o s i t i o n o f h y p e r s p e c t r a l 2 2 2 3 2 4 2 5 3 2 6 2 7 2 8 2 9 3 O i ma g e r y b a s e d o n mu l t i p l e k e r n e l S VM J Op t Pr e c i s i o n En g。2 0 1 4，2 2(7)：1 9 1 2 1 9 2 0 (i n Ch i n e s e)J I C H，YU J H，W ANG Y Y，e t a 1 Au t o ma t i c p h a s e-b a s e d e d g e d e t e c t i o n o f c o r n e a l S h e i