1、医学图像处理绪论2022-10-151主讲教师：计算机教研室 王 静办公室电话：Email:办公室地址：逸夫楼819医学图像处理绪论2022-10-152医学图像处理绪论3上课安排6：30-7：30 上课7：30-9：30 上机医学图像处理绪论2022-10-154课程考核平时作业（60%）平时考勤及课堂表现（15%）期末作业（25%）医学图像处理绪论2022-10-155教学目的理解医学图像处理和分析的主要方法了解matlab图像处理工具 在今后的学习工作中能更好的理解和使用图像处理软件或硬件 在今后的学习工作中能根据自己的科研或工作需要，编写一些小的图像处理程序医学图像处理绪论2022-1
2、0-156主要参考书田捷 医学影像处理与分析 电子工业出版社罗述谦 医学图象处理与分析 科学出版社冈萨雷斯 数字图像处理(第二版)冈萨雷斯 数字图像处理（Matlab版）医学图像处理绪论2022-10-157讲义及资料.54 医学图像处理文件夹下大家可以在这里下载讲义及其他资料每节课的作业也上传至此医学图像处理绪论8数字图像处理的概念及应用范畴数字图像处理的概念及应用范畴医学图像处理和分析的意义医学图像处理和分析的意义医学图像处理和分析的内容医学图像处理和分析的内容数字图像处理在医学中的应用数字图像处理在医学中的应用常见的医学图像常见的医学图像医学图像处理和分析的内容及实例医学图像处理和分析的
3、内容及实例Matlab基础基础医学图像处理绪论2022-10-159数数 字字 图图 像像 处处 理理I Imma ag ge e P Pr ro oc ce es ss si in ng g医学图像处理绪论2022-10-1510百闻不如一见One picture is worth more than ten thousand words.Anonymous医学图像处理绪论2022-10-1511 人类通过眼、耳、鼻、舌、身接受信息，感知世界。约有75%的信息是通过视觉系统获取的。数字图像处理是用数字计算机处理所获取视觉信息的技术。医学图像处理绪论2022-10-1512数字图像目的：便于计
4、算机处理。连续图象数字图象离散化RC医学图像处理绪论2022-10-1513数字图像的表示数字图像一般以矩阵形式表示：)1,1()1,1()0,1()1,1()1,1()0,1()1,0()1,0()0,0(),(NMfMfMfNfffNfffyxf医学图像处理绪论2022-10-1514数字图像离散化：空间采样：512X512 象素 Pixel分辨率 Resolution幅值量化：256级,28级,8bit,灰度级 Grey level医学图像处理绪论2022-10-1515一幅图像可定义为一个二维函数 ，这里x和y是空间坐标，而在任何一对空间坐标（x,y）上的幅值f称为该点图像的强度或灰度
5、。当x,y和幅值f为有限的、离散的数值时，称该图像为数字图像数字图像处理的概念),(yxf数字图像处理是指借用数字计算机处理数字图像。数字图像是由有限的元素组成的，每个元素都有特定的位置和幅值，这些元素称为图像元素或像素医学图像处理绪论2022-10-1516 分辨率不同 的图象比较分辨率 640 x480分辨率 320 x240医学图像处理绪论19分辨率 160 x120医学图像处理绪论20分辨率 80 x60医学图像处理绪论21医学图像处理绪论22分辨率 640 x480医学图像处理绪论23分辨率 320 x240医学图像处理绪论24分辨率 160 x120医学图像处理绪论25分辨率 80
6、 x60医学图像处理绪论26从粗到细地观察医学图像处理绪论2022-10-1527查视力=检测分辨率？图像数据格式图像在计算机中是以图像文件的形式存储的。存储的格式有多种，较常用的有BMP,GIF,JPEG,TIFF,PCX，dicom等BMP图像文件的结构分三部分：文件头、位图信息数据块以及图像数据 文件头：文件类型、大小和打印信息 位图信息数据：图像数据医学图像处理绪论2022-10-1531数字图像处理的范畴数字图像处理的范畴图像处理（低级处理）：便于处理对图像信息进行改进，比如降低噪声，对比度的增强和图像锐化等。特点：输入图像，输出图像。图像分析（中级处理）涉及分割，及特征提取。输入是
7、图像，输出为图像的特征计算机视觉（高级处理）：涉及在图像分析中被识别物体的总体理解，以及执行与视觉相关的行为医学图像处理绪论2022-10-1532医学图像处理绪论2022-10-1533医学图像处理绪论2022-10-1534医学图像处理绪论2022-10-1535医学图像处理绪论2022-10-1536医学图像处理绪论2022-10-1537数字图像处理的发展概况和应用发展概况：（1）二十世纪二十年代：图像远距离传输。（2）二十世纪五十年代：数字计算机发展到一定水平，数字图象处理引起巨大关注。（3）1964年：美国喷气推进实验室用计算机对“徘徊者七号”太空船发回的大批月球照片进行处理。数字
8、图象处理应用从空间研究计划扩展到生物医学，工业生产，军事侦察等领域。医学图像处理绪论2022-10-1538图像处理最先应用于空间探索医学图像处理绪论2022-10-1539图像处理最先应用于空间探索医学图像处理绪论2022-10-1540（4）二十世纪六十年代末：数字图象处理较完整的理论体系已形成，成为一门新兴的学科。（5）二十世纪六十年代至八十年代：随着离散数学理论的创立和完善，数字图象处理理论和方法进一步完善，应用范围更加广泛。（6）二十世纪八十年代以来：数字图象处理向更高级的方向发展：实时性，智能化，普及化，网络化，低成本。医学图像处理绪论2022-10-1541应用领域：通信：图象传
9、输，电视电话，HDTV等。宇宙探测：星体图片处理。遥感：地形、地质、矿藏探查，森林、水利、海洋、农业等资源调查，自然灾害预测，环境污染的监测，气象云图。生物医学：CT，NMR，PET，SPECT，DSA，X射线成象，B超，红外图象，显微图象。工业生产：产品质量检测，生产过程控制，CAD，CAM。医学图像处理绪论2022-10-1542应用领域：军事：军事目标侦察，制导系统，警戒系统，自动火器控制，反伪装等。公安：现场照片，指纹，手迹，印章，人像等处理和鉴别。档案：过期的文字、图片档案的修复和处理。机器人视觉娱乐：电影特技，动画，广告，MTV等。医学图像处理绪论2022-10-1543应用举例非
10、生物医学领域医学图像处理绪论2022-10-1544医学图像处理绪论45医学图像处理绪论46医学图像处理绪论47医学图像处理绪论2022-10-1548医学图像处理绪论2022-10-1549医学图像处理绪论50医学图像处理绪论51 医学图像处理绪论2022-10-1552 医学图像处理绪论53 医学图像处理绪论54医学图像处理绪论55医学图像处理绪论2022-10-1556 医学图像处理绪论2022-10-1557医学图像处理绪论2022-10-1558医学图像处理绪论2022-10-1559医学图像处理绪论60医学图像处理绪论2022-10-1561医学图像处理绪论62医学图像处理绪论63
11、医学图像处理绪论2022-10-1564医学图像处理绪论2022-10-1565医学图像处理绪论2022-10-1566图像处理在医学中的应用医学图像处理和分析是近几年兴起的新兴交叉学科，正方兴未艾。1。借助图像处理技术的有力手段，医学影像的质量和显示方式得到了极大的改善，从而借助于医学图像处理与分析手段使得诊疗水平大大提高。2。为医学培训、医学研究与教学、计算机辅助临床外科手术等提供数字实现手段，为医学的研究与发展提供坚实的基础，具有不可估量的价值医学图像处理绪论2022-10-1567医学图像处理研究的内容医学图像成像技术（Medical Imaging）医学图像处理和分析（医学图像后处理
12、post-processing）图象增强技术 图象分割技术 图像配准技术 图象融合技术 图像三维可视化技术 对象识别技术医学图像处理绪论2022-10-1568医学图像成像技术（Medical Imaging）主要研究如何将实际人体信息转换为 计算机的二进制数据的问题。经过一定的重建算法把这些数据转换成与二维图像像素或三维体数据对应的原始图像数据集医学图像处理绪论2022-10-1569医学影像设备目前常见的医学影像 B超扫描图像、彩色多普勒超声图像、核磁共振图像、CT图像、PET图像、SPECT图像、数字X光机（DX）图像、X射线透视图像、各种电子内窥镜图像，显微镜下病理切片图像等。医学图像
13、处理绪论2022-10-1570医学图像处理绪论2022-10-1571伦琴开创了人体图像的先河X射线成像基于待成像物体各部分的密度不同，对x射线的吸收不同，透射X射线强度不同，从而在胶片上成像的。X光图片是X射线在通路上物体对射线吸收的积分效果。X光图片不能反映组织或病灶的三维空间位置医学图像处理绪论2022-10-1572X射线应用-血管造影术血管造影术是对比增强辐射成像领域中的另一主要应用，这过程用于得到血管图像（血管造影照片）。一个导管（小的柔软中空的管子）插入动脉或静脉，导管穿过血管并被引导到要研究的区域。当导管到达探查部位时，x射线对比介质通过导管注入，这就增强了血管的对比并可以看
14、到任何病变或阻塞医学图像处理绪论2022-10-1573医学图像处理绪论2022-10-1574医学图像处理绪论2022-10-1575CT计算机体层摄影（computed tomography,CT）扫描仪利用X线对人体某一范围进行逐层的横断扫描，取得信息，经计算机处理后获得重建的图像。获得的图像为人体的横断解剖图，并可通过计算机处理得到三维的重建图像。（X射线断层扫描影像装置）医学图像处理绪论2022-10-1576CT工作原理医学图像处理绪论2022-10-1577CT工作原理人体各种组织（包括正常和异常组织）对X线的吸收不等。CT即利用这一特性，将人体某一选定层面分成许多立方体小块，这
15、些立方体小块称为体素。它们排列成行列方阵，形成图像矩阵。当X线球管从一方向发出X线束穿过选定层面时，沿该方向排列的各体素均在一定程度上吸收一部分X线，使X线衰减。当该X线束穿透组织层面（包括许多体素）为对面探测器接收时，X线量已衰减很多，为该方向所有体素X线衰减值的总和。然后X线球管转动一定角度，再沿另一方向发出X线束，则在其对面的探测器可测得沿第2次照射方向所有体素X线衰减值的总和；以同样方法反复多次在不同方向对组织的选定层面进行X线扫描，即可得到若干个X线衰减值总和。在上述过程中，每扫描一次，即可得一方程。该方程中X线衰减总量为已知值，而形成该总量的各体素X线衰减值是未知值。经过若干次扫描
16、，即可得一联立方程，经过计算机运算可解出这一联立方程，而求出每一体素的X线衰减值，再经数/模转换，使各体素不同的衰减值形成相应各像素的不同灰度，各像素所形成的矩阵图像即为该层面不同密度组织的黑白图像。医学图像处理绪论2022-10-1578某物质的CT值等于该物质的衰减系数与水的吸收系数之差再与水的衰减系数相比之后乘以1000。其单位名称为HU。值水 0值空气 1000值骨 10003）常见人体组织的值（）组织 值 组织 值骨组织 400 肝脏 5070钙值 80300 脾脏 3560血块 6484 胰腺 3055脑白质 2534 肾脏 2550脑灰质 2844 肌肉 4055脑脊液 38 胆
17、囊 1030血液 1332 甲状腺 5090血浆 314 脂肪 20100渗出液 15 水 0医学图像处理绪论2022-10-1579CTCT的密度的密度分析CT图像，一方面是观察解剖结构，另一方面是了解密度改变。后者可通过测定CT值而知，亦可与周围组织的密度对比观察。人体内肿瘤组织因部位、代谢、生长及伴随情况不同，其密度变化各异。CT对组织的密度分辨率较高，且为横断面扫描，提高了肿瘤诊断的准确率。医学图像处理绪论2022-10-1580造影剂的应用造影剂的应用虽然CT较普通X线摄影有更高的密度分辨率，但有些病变与正常组织间的密度差异很小，需要利用造影剂使上述密度差异加大，以帮助诊断。CT扫描
18、用造影剂可分两大类：一类为用于空腔脏器的造影剂，另一类为静脉注射用造影剂（偶尔也通过动脉注射造影剂）。医学图像处理绪论2022-10-1581近20年来，由于CT装置、成像软件及扫描技术的不断更新和改进，出现了CT电影，CT血管造影，超高速CT，高分辨率CT，螺旋CT等技术，使得CT技术在临床的许多领域得到应用，提高了诊断的准确性和可信度。CT是当今影像技术中不可缺少的手段骨骼在CT图像中可以显示得非常清楚医学图像处理绪论2022-10-1582CT检查的优点：（1）CT为无创性检查，检查方便、迅速，易为患者接受。（2）有很高的密度分辨力，密度相差5-6H的不同组织能被区分。能测出各种组织的C
19、T值。（3）CT图象清晰，解剖关系明确。（4）CT能提供没有组织重叠的横断面图象，并可进行冠状和矢状面图象的重建。（5）用造影剂进行增强扫描，不仅提高了病变的发现率，而且有的能做定性诊断。医学图像处理绪论2022-10-1583CT检查的限制CT扫描虽有广泛的适应范围，但仍有限度。虽然发现病变的敏感性极高，但在定性诊断上仍有很大的限制。由于CT机测定的是物理参数，即人体组织对X线的衰减值或物理密度，医生就是根据正常组织和异常组织呈现的衰减值差异作为诊断的依据，如果衰减值无差异，再大的肿瘤也无法鉴别。可见CT扫描尽管有许多优越性，但也有其局限性，只有与其他设备，其他诊断手段相配合，才能充分发挥其
20、作用。医学图像处理绪论2022-10-1584MRI核磁共振成像把病人放在强磁场中，并使无线电短脉冲通过病人身体，每个脉冲将导致一个有病人组织发射的无线电响应脉冲，这些信号发生的位置和强度由计算机确定，从而产生一个病人的横截面图像MRI可以在任何平面产生图像医学图像处理绪论2022-10-1585医学图像处理绪论2022-10-1586医学图像处理绪论2022-10-15871.无电离辐射危害，是一种安全的检查的方法。2.多参数成像，可以提供丰富的诊断信息。3.高对比度成像，在所有医学影像技术中，MRI的软组织对比分辨力最高。4.具有任意方向断层的能力，能够从不同角度直视地观察分析组织结构及其
21、病变。5.与传统血管造影相比，它最大的优点是无创性，无需使用对比剂，可直接显示心脏和血管结构，即磁共振血管成像，从图像质量来看，目前它至少可以显示大血管及各主要脏器的一、二级分支血管。6.无骨伪影干扰，后颅凹病变清晰可辨。7.可进行功能、组织化学和生物化学方面的研究。MRI的优点医学图像处理绪论2022-10-1588临床应用中枢神经系统 MR已成为颅颈交界区，颅底、后颅窝及椎管内病变的最佳检查方法。该技术对脑肿瘤、脑部炎症、脑血管病、脑发育不全、脑变性病与白质脑病有良好的诊断价值 对颅脑损伤与癫痫的诊断起辅助诊断作用对脊柱与椎管内病变有重要诊断意义 如脊管肿瘤、脊柱炎症病变、脊柱先天畸形、退
22、行性病变、脊柱外伤手术与手术后改变等应用肌肉关节系统 MRI已成为肌肉、肌腱、韧带、软骨病变影像检查的主要手段之一。对于关节周围病变、股骨头缺血性坏死，半月板损伤、松质骨细微结构的破坏，骨小梁骨折（骨挫伤）以及骨髓腔内病变均具有重要的诊断价值。电影MRI技术还可进行关节功能检查。fMRI功能成像研究脑功能对肺、密质骨，如四肢长骨、胃肠道疾病不如CT医学图像处理绪论2022-10-15891 Hz2 Hz3 Hz%100755025Schlaug,et al,1995,Harvard Medical School and Beth Israel Hospital医学图像处理绪论2022-10-1
23、590Schlaug,et al,Science 267:699(1995)医学图像处理绪论2022-10-1591MRI的缺点成像速度慢，尤其是对于低场强者，这是MRI的主要缺点。不适合运动性器官，如肺部疾病患者和危重病人的检查，对于躁动或丧失自制能力的人，如不使用镇静剂，也是难以成像的。对钙化灶和骨皮质病灶不够敏感。禁忌症相对较多。MRI设备的强磁场和RF(射频)场有可能使心脏起搏器失灵，也容易使各种金属性植入物移位。在激励电磁波作用下，体内的金属还会因为发热而对病人造成伤害。因此，安置心脏起搏器、假肢、人工髋关节的病人、疑有眼球异物的病人以及动脉瘤银夹结扎术后的病人都禁做MRI检查。装假
24、牙的病人不能进行颌面水平MRI检查。放置宫内节育环者，如检查中出现不适应立刻停止检查。4、图像易受多种伪影干扰。MRI伪影主要来自患者的运动，金属异物以及设备三个方面。医学图像处理绪论2022-10-1592超声图像（sonography）医用超声诊断仪是将声纳原理和雷达技术相结合生产的为临床应用的医疗仪器。其基本原理是高频超声脉冲波辐射到生物作内，由生物体内不同界面反射出不同波形并形成图像从而判断生物体内是否有病变。超声诊断仪由起初的一维超声扫描显示，发展为二维甚置三维、四维的超声扫描和显示，大大增加了回波信息量，使生物体内的病灶清晰，易辨，因此，它将被越来越广泛地应用医用超声诊断仪 医学图
25、像处理绪论2022-10-1593超声优点1.有高的软组织分辨力。组织只要有1的声阻抗差异，仪器就能检测出，并显示其反射回波。2具有高度的安全性 当严格控制声强低于安全阂值时，超声可能成为一种无损伤的诊断技术，对医务人员更是十分安全。3实时成像 它能高速实时成像，可以观察运动的器官，并节省检查时间。4使用简便，费用较低，用途广泛。医学图像处理绪论2022-10-1594B超缺点比如B超在清晰度、分辨率等方面，较弱。B超对肠道等空腔器官病变易漏诊。气体对超声影响很大，患者容易受到患者肠气干扰等多方面因素影响检查结果 B超检查需要改变体位屏气等，对于骨折和不能配合病人不适用。检查结果也易受医师临床
26、技能水平的影响。孕妇滥查B超可能易致胎儿畸形医学图像处理绪论2022-10-1595医学图像处理绪论2022-10-1596医学图像处理绪论2022-10-1597Photo courtesy Philips Research 医学图像处理绪论2022-10-1598Photo courtesy Philips Research 医学图像处理绪论2022-10-1599彩色多普勒超声图像应用案例 借助多普勒频移原理，可探测到人体内血流动态。经过伪彩色编码，以不同的颜色标识血流流向，色彩的灰度表示血流速度，色块面积表示流量。通常肝脏病变，例如肝癌或肝硬化，都会在肝脏供血上反映出来。注射造影剂后，
27、血流动态显像还会有更多血流信息供临床诊断应用医学图像处理绪论2022-10-15100医学图像处理绪论2022-10-15101伽马射线成像在核医学中，将放射性同位素注射到病人体内，当同位素衰变时放射出伽马射线，用伽马射线监测器收集放射物产生图像。如图显示了一幅利用伽马射线成像得到的骨骼扫描图像，这种图像用于骨骼病理（例如感染或肿瘤）定位医学图像处理绪论2022-10-15102正电子断层扫描仪图像 正电子发射断层扫描仪PET(Positron Emission Tomography)是当今最高层次的核医学技术，它也是当前医学界公认的最先进的大型医疗诊断成像设备之一。电子发射断层成像已成为肿瘤
28、、心、脑疾病诊断与病理研究中的不可缺少的重要方法该技术主要研究人体内部组织与器官功能的重要手段医学图像处理绪论2022-10-15103PET显像主要是在分子水平上提供有关脏器及其病变的功能信息。PET所用的示踪药物主要是18F-FDG短寿命正电子核素，它们是组成有机体组织的基本成分并能参与代谢过程。大多数疾病的生化变化先于解剖学的变化，并且PET 对于示踪剂浓度的灵敏度非常高，能高精度地定量地检测出代谢过程的非正常增加并给出清晰的图像（所谓热源成像），因此 PET 能提供很多疾病在发展过程中的早期信息，可以进行超前诊断，尤其适合于肿瘤的早期诊断医学图像处理绪论2022-10-15104如图是
29、构成三维再现图像序列的一幅样品，图像显示脑部和肺部各有一个肿瘤，如图是构成三维再现图像序列的一幅样品，图像显示脑部和肺部各有一个肿瘤，很容易看到的小白块很容易看到的小白块医学图像处理绪论2022-10-15105SPECT技术，即单光子发射计算机断层成像术（，），是核医学的一种技术。由于它是对从病人体内发射的射线成像，故亦称发射型计算机断层成像术 ，。与相比，二者不同的是射线的入射方式，接受的射线是由受检者体内发出来的，而是由射线从体外穿透人体到达接受器的。医学图像处理绪论2022-10-15106SPECT临床应用在神经系统它可以提供局部脑血流的供应数据，诊断、癫痫等，可诊断和鉴别诊断交通性
30、脑积水和梗阻性脑积水；在心血管方面可显示心脏并测定心脏的功能提供、等值，特别在冠心病诊断方面有其独特的意义，与冠状动脉造影相比，的灵敏度可达到，特异性可达到，能直观地提供心肌缺血的部位、范围及严重程度，也能提示冠状动脉病变的部位；还可以在心肌梗死发生后行心肌显像以判断心肌梗死区内的心肌是否存活，这对于下一步的治疗如、放置支架、冠脉搭桥术等冠状动脉血运重建非常重要 医学图像处理绪论2022-10-15107在肿瘤诊疗方面，利用亲肿瘤的放射性药物进行肿瘤显像，可以发现原发肿瘤及转移的肿瘤，特别是转移到骨骼上的肿瘤，它可以在病人还未感到不适或疼痛时发现病变的存在，可比线提早月发现病灶。所以当肿瘤病人
31、感到全身不适或疼痛时对其及时进行检查是合适和必要的。医学图像处理绪论2022-10-15108在泌尿系统疾病方面，可以动态和静态观察两个肾脏的形态、位置、肾脏的功能并计算出两肾、输尿管的通畅情况；在呼吸系统方面，可进行肺灌注显像和肺通气显像诊断肺梗塞及下肢静脉显像寻找栓子存在；医学图像处理绪论2022-10-15109在新生儿黄胆时可行肝胆显像鉴别诊断新生儿肝炎和先天性胆道闭塞；消化道出血时可定位出血部位和诊断美克尔憩室（），后者的诊断率为。当然还可以在器官移植、妇科病、血液病方面以及科研等多学科、多领域进行广泛应用。医学图像处理绪论2022-10-15110仪器是医院必备的常规检查设备，以检
32、查无创、安全、简便、经济、实用等优点而被广大病人所接受，同时也为临床、科研、教学提供有力帮助而被医务人员所接受。医学图像处理绪论2022-10-15111成像技术CTMRIPETSPECT优点1.成像速度快2.对骨组织敏感3.空间分辨率较高4.造价相对较低1.软组织成像效果好2.空间分辨率高3.对人体无危害4.扫描角度灵活5.无骨伪影1.对比度较高2.受体成像3.属于化学成像1.对比度较高2.价格较低缺点1.X射线对人体有危害2.对软组织成像较差3.骨的边缘在成像中易产生条状伪影4.成像角度不灵活1.对骨组织成像效果较差2.成像时间较CT长1.价格昂贵，必须配备回旋加速器2.分辨率较差1.分辨
33、率低2.成像速度慢医学图像处理绪论2022-10-15112DICOM标准DICOM（Digital Imaging and Communication in Medical）医疗设备的国际标准通讯协议以解决不同厂商的各种医疗设备的互连问题医学图像处理绪论2022-10-15113医学图像处理和分析的应用辅助医生诊断 突出感兴趣区域，消除噪声，提供更多的信息 三维建模 图像融合 虚拟内寇镜治疗规划辅助外科手术 制定手术规划 手术导航和术中监护医学教学教学 数字解剖模型 手术教学训练科学研究 脑功能和结构的研究医学图像处理绪论2022-10-15114医学图像处理绪论1152022-10-151
34、16医学图像处理绪论116HeadS/I ProjectionMRI Center,University of Rochester Medical School医学图像处理绪论2022-10-15117Chen,Dougherty,Totterman,Hornak,Magn.Reson.Med.29:358(1993)NormalCysticScleroticOsteoporosis医学图像处理绪论2022-10-15118授课内容数字图像处理绪论（一节课）医学图像增强（3,4节课）医学图像分割、配准及融合（一节课）基于医学影像的计算机辅助诊断（一节课）医学图像三维显示（一节课）虚拟人体计划研
35、究概述（一节课）分子影像学中图像处理与分析技术（一节课）医学图像处理绪论2022-10-15119医学图像处理临床应用现状和前景 目前，临床上医学图像后处理主要以应用随机软件进行直接处理为主，由于计算机技术的进步，各种影像学设备随机安装大量后处理软件，使图像后处理的临床应用越来越广泛。这使图像的显示、其诊断和鉴别诊断能力均不断提高，受到影像学科和其他临床学科医生的欢迎和好评。但是近年来脱机应用工作站进行图像处理有逐渐增多的趋势。以图像融合为例：对同一患者而言，将密度分辨率最高、显示骨质结构和钙化最佳的CT与软组织对比分辨率最高的MRI，或者将显示解剖结构清楚的CT或MRI与显示功能和代谢改变的
36、PET或SPECT图像进行融合，获取一种新型图像，可增加诊断信息，或者使病灶的定位更准确，其形态结构显示得更直观。由于医学影像学设备联网和PACS的建立和临床应用，数字化图像包含大量有用信息，所以，图像后处理的方法有待于进一步开发。例如：计算机辅助诊断(computed aidded diagnosis,CAD)技术，利用计算机进行图像的识别和辅助诊断，在乳腺癌和周围型肺癌的诊断上，有肯定的临床应用价值，并且正在向其他疾病的诊断方面扩展。有足够迹象表明：医学图像后处理的新方法和新技术有广阔的发展前景，应该引起我们的关注。医学图像处理绪论2022-10-15120matlabMATLAB是由美国
37、mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。医学图像处理绪论2022-10-15121Matlab环境医学图像处理绪论2022-10-15122图像的输入输出读取图像读入Tiff,bmp,jpeg,gif,png图像X=imread();显
38、示读入的图像imshow(X)同时显示多幅图像 figure,imshow(g)保存图像imwrite(f,)读入dicom图像X,map=dicomread()显示读入的图像montage(X,map);保存dicom图像dicomwrite(X,)医学图像处理绪论2022-10-15123图像类型亮度图像 像素值表示亮度 uint8 uint16 double二值图像索引图像RGB图像医学图像处理绪论2022-10-15124Matlab中的数据类名称描述doubleUint8Uint16Uint32Int8Int16Int32SingleCharlogical双精度浮点数，范围：-103
39、0810308（8字节/像素）无符号整数 0 255 1字节/像素 0 65535 2字节/像素 0 4294967295 4字节/像素有符号整数 -128 127 1字节/像素 -32768，32767 2字节/像素 -2147483648，2147483647 4字节/像素单精度浮点数，范围：-10381038（4字节/像素）字符（2字节/像素）值为0或1 1字节/像素医学图像处理绪论2022-10-15125数据类和图像类型之间的转换数据类间的转换 B=data_class_name(A)图像类和类型间的转换医学图像处理绪论2022-10-15126数组索引向量（行向量）创建 向量的元素
40、用方括号括起，并由空格或逗号分开。例如：v=1 3 5 7 9W=1，2，3，4,5 单个像素的访问:向量名(下标)如：v(3)访问连续元素：例如要存取w的中间三个元素 w(2:4)存取w中的最后三个元素w(3:end)访问不连续数据w(1:2:end)从1开始计数步长为2(步长可以为负数)end结束 一个向量可以做另一个向量的索引 W(1 4 5)从行向量产生列向量向量转置（.）x=w.w(:)注意与w(1:end)的区别医学图像处理绪论2022-10-15127矩阵索引矩阵的创建 用一列被方括号括起来并用分号隔开的行向量表示例如：A=1 2 3；4 5 6；7 8 9矩阵元素的访问 访问单
41、个元素A(2,3)访问二维元素块 A(:,3)A(2,:)A(1:2,2:3)A(2:end,end:-2:1)使用向量作为矩阵的索引 E=A(1 3,2 3)使用逻辑数组 D=logical(1 0 0;0 0 1;0 0 0);A(D)A(:)逐列的形式将元素组成列向量sum(A(:)sum(A)区别医学图像处理绪论2022-10-15128一些重要的标准数组zeros(M,N)生成一个大小为M*N的double类矩阵，其元素均为0ones(M,N)ture(M,N)false(M,N)rand(M,N)医学图像处理绪论2022-10-15129M函数编程简介M文件可以是简单的执行一系列M
42、atlab语句的脚本，也可以是接受变量并产生一个或多个输出的函数M文件扩展名.mM文件组成：函数定义行 H1行 帮助文本 函数体 命令函数定义行的形式为：函数定义行的形式为：functionoutputs=name(input)例：例：function s,p=sumprod(f,g)调用s,p=sumprod(f,g)医学图像处理绪论2022-10-15130运算符算术运算符关系运算符逻辑运算符医学图像处理绪论2022-10-15131算术运算符运算）数组运算符（逐个元素数规则运算）矩阵运算符（按线性代算术运算符+-两种运算是一致的，其余用圆点（.）来区分两种运算例如：C=A*B表示传统意义
43、上的矩阵成积D=A.*B表示数组乘法，要求A B两个数组大小 相同，逐像素相乘 即：D(I,J)=A(I,J)*B(I,J)医学图像处理绪论2022-10-15132运算符 名称Matlab函数示例+数组和矩阵加减plus(A,B)-minus(A,B).*数组乘times(A,B)*矩阵乘mtimes(A,B)./数组右除rdivide(A,B)C=A./B c(i,j)=A(i,j)/B(i,j)/矩阵右除mrdivide(A,B).数组左除ldivide(A,B)C=A.B c(i,j)=B(i,j)/A(i,j)矩阵左除mldivide(A,B).数组求幂power(A,B)C=A.B
44、 C(I,j)=A(I,j)B(i,j)矩阵求幂mpower(A,B).向量和矩阵转置transpose(A)Ctranspose(A)所有运算的操作数可以是所有运算的操作数可以是矩阵、数组或标量。数据矩阵、数组或标量。数据类型可以实数或复数类型可以实数或复数（double类型的数据）类型的数据）医学图像处理绪论2022-10-15133函数描述imaddimsubtractimmultiplyimdivideimabsdiffimcomplementimlincomb两幅图像相加或把常数加到图像两幅图像相减或从图像减去常数两幅图像相乘，相应像素相乘，或图像乘以常数两幅图像相除，相应像素相除，
45、或图像除以常数计算两幅图像间的绝对差对图像求补计算两幅或多幅图像的线性组合IPT主持的图像算术函数优点是支持整数数据类优点是支持整数数据类医学图像处理绪论2022-10-15134关系运算符运算符 =-=对相同维数的数组中的相应元素进行比较A=1 2 3 4 5 6 7 8 9 B=0 2 4;3 5 6;3 4 9B=0 2 4 3 5 6 3 4 9 A=Bans=0 1 0 0 1 1 0 0 1A=1 2 3 4 5 6 7 8 9 B=0 2 4;3 5 6;3 4 9B=0 2 4 3 5 6 3 4 9 A=Bans=1 1 0 1 1 1 1 1 1医学图像处理绪论2022-1
46、0-15135逻辑运算符与函数逻辑运算符：&、|、运算符既能作用于逻辑数据也能作用于数值数据。逻辑1和所有非零数据作为TRUE 来处理将逻辑0和数值0做false来处理A=1 2 0;0 4 5;B=1-2 3;0 1 1;A&Bans=1 1 0 0 1 1医学图像处理绪论2022-10-15136流控制if expression statemants end If if expression statemants1else statemants2 end if expression1 statemants1else if expression2 statemants2else statem
47、ants2 end 医学图像处理绪论2022-10-15137For语句for index1=start:increment:end statementend 例：count=0for i=0:0.1:1 count=count+1;end 医学图像处理绪论2022-10-15138whileWhile express statement;endwhile express1 statement1;while express2 statement2;endendBreak:跳出循环Continue :跳出本次循环医学图像处理绪论2022-10-15139switchSwitch switch_expressCase case_express1 statement1;Case case_express2 statement2;Case case_express3 statement3;Otherwise statementn;end医学图像处理绪论2022-10-15140作业一结合自己所学专业，谈谈图像处理在本专业的应用熟悉matlab环境医学图像处理绪论2022-10-15141
