医学影像系统原理课后作业
一医学影像系统原理课后作业 Part 1 1、影响X射线模糊度的因素有哪些? 主要因素有几何学模糊、运动型模糊、增感屏-胶片系统产生的模糊和散射线性模糊。 2、如何获得能量高的硬X射线? 提高管电压。 3、普通X射线摄影技术与CR, DR有什么不同? CR是一种X线摄影的数字采集技术,使用常规X线摄影的采集结构,利用成像板取代 传统的屏-片体系,它在光激励荧光体中记录X线影像,并使其影像信息以电信号方式提取 出来;CR的应用突破了常规X线摄影技术的固有局限性,实现了常规X线摄影信息数字化。 然而,CR也存在不足之处,如时间分辩力较差,不能满足动态器官形态的显示;间接转换 可出现信息的丢失,其空间分辨率仍有不足,难于显示细微的组织结构;成像板(IP)为易 耗品,图像质量会随使用次数过多而下降;不能改变传统X线摄影检查的工作流程。 DR较之CR具有更高的空间分辩率,更高的动态范围和DQE,更低的X线照射量,图 像层次更丰富,在曝光后几秒内即可显示图像,大大改善了工作流程,提高了工作效率。 DR系统最重要的部件是平板探测器,直接数字化X线成像的平板探测器利用了非晶硒的光 电导性,将X线直接转换成电信号,形成全数字化影像。 Part 2 能否采用比体元数少的投影方程数重建CT图像?为什么?如果可能,请给出你的实现方案与 步骤。 可以。 Part 3 什么是纵向驰豫时间?什么是横向驰豫时间?二者有什么区别和联系。 纵向驰豫时间T1是指纵向磁化矢量从最小值恢复至平衡态的63%所经历的弛豫时间。其 物理学意义相当于一个“弛豫周期”,每经过一个T1时间则纵向磁化恢复其剩余值的63%。 T1是不同组织的弛豫特征值,反应不同组织的纵向弛豫率的快慢差别。山于纵向弛豫是高 能态自旋释放能量恢复低能态的过程,所以高能态自旋必须通过有效的途径将能量传递至周 围环境(晶格)中去,故又称其为自旋一晶格弛豫。 横向驰豫时间T2是射频脉冲停止后,横向磁化矢量衰减至其最大值的37%时所经历的时 间。T2也是不同组织的弛豫特征值,反应不同组织横向磁化弛豫率的快慢差别,其物理意 义与T1相似,只是T2代表横向磁化的“衰减周期”,每过一个T2时间,横向磁化减少至 其剩余值的37%,与放射性元素的半衰期意义相近。 Part 4 功能成像是一种还是一类磁共振成像方法? 是一类磁共振成像,MRI是相对于MR形态学而言的,具有较广泛的含义,包括扩散、 灌注加权成像、皮层活动功能定义及MR波谱成像等。 Part 5 SPECT与PET也是一种CT成像方式吗? SPECT与PET不是CT成像方式,SPECT与PET都是ECT的一种。发射式计算机断层 (emission computed tomography, ECT)是以引入体内的放射性药物为放射源,应用可以旋转 的或环状闪烁探测器在体外从不同角度采集脏器中放射性药物分布的信息,并通过计算机进 行数据处理和断层重建,获得脏器或组织的横断面、矢状面及冠状面三维图像。SPECT是 单光子发射式计算机断层成像,PET是正电子发射式计算机断层成像。 医学图像处理 Part 1 1. In averaging de-nosing , how much improvement has been made for SNR of averaged image? Answer: If there are N items of data for the average, the SNR increased N times. 2. Derive the ula of averaging using discrete convolution Answer: 3. Program a median filter and filter a noise image Answer: load an image, then median filtering, and then show dear: de; [im,map]=imread( 4 lena.gif1 ): 11=ind2gray(im,map); imshow(ll); I2=medfilt2(tl): % 中值建次 figure; imshow(l2); Part 2 1. What is the mathematical relationship of the mean filter and Sobel operator? Answer: They are all first-order differential operator, the er general average, which is the weighted average filtering 2. Describe how to approximate of LoG by DoG Answer: Approximately described by Gradient Algorithm: A Gni - G„2 Part 3 1. Program the run-length encoding and calculate the compression ratio using different types of image. Answer: matlab RLE code : function data = rle(x) if iscell(x) % decoding i = cumsum([ 1 x{2}]): j = zerojfl. i(end)-l); J(i(l:end-1)) = 1; data = x{l}(cumsum(j)); else % encoding if size(xj) > size(x,2)« x = x ; end % if x is a column vector, tronspose i = [ find(x(l:end-l) ~= x(2:end)) length(x)]; data{2} = diff([0i]); data{l} = x(I): end To different type of image, the compression ratio of RLE is different. dear: [im,map]=imreadfcameraman.gif): 11=ind2gray(im,map); compresjedDatal = rlefll); [im.map] = imreadfhead.gif): 12=ind2gray (im.map): compres$edData2 = He (12); [im.map] = imreadfbob.gif); !3=ind2grayGm,map): compre$5edData3 = rle(l3); The compression ratio of the three images: ratiol =