医学图像处理实验报告(要)hidx.docx

医学图像处理实验报告医学图像像处理实实验报告告Lab Repportt off Meediccal Imaage Proocesssinng系 部：学 号：姓 名：指导教师师：实验1 MATTLABB基本操操作及函函数用法法1.1 实验原原理介绍绍1.1.1读取取函数函数immreaad可以以将图像像读入MMATLLAB环环境，语语法为：imreead('fiilenamme')其中fiilennamee是一个个含有图图像文件件全名的的字符串串（包括括任何可可用的扩扩展名）。例例如命令令行>>f = iimreead ( ''cheestxxrayy.jppg');将JPEEG图像像cheestxxrayy读入图图像数组组f中。函数siize可可以给出出一副图图像的行行数和列列数：>>siize(f)ans-1024 10244函数whhos可可以显示示一个数数组的基基本附加加信息。例例如，语语句>>whhos fwhoss行结尾尾处的分分号对结结果没有有影响。1.1.2 显示示图像一般使用用函数iimshhow显显示，其其基本语语法为：imshhow (f, G)其中，ff是一个个图像数数组，GG是显示示该图像像的灰度度级数。若若G省略略，默认认为2556。语语法：imshhow f, llow hiigh会将所有有小于或或等于llow的的值都显显示为黑黑，所有有大于或或等于hhighh的值都都显示为为白色。介介于loow和hhighh之间的的值将以以默认的的级数显显示为中中等亮度度值。当用immshoow显示示另外一一幅图像像时，MMATLLAB会会在屏幕幕上用新新图像替替换旧图图像，为为保持第第一幅图图像并同同时显示示第二副副图像，可可以使用用如下ffiguure函函数：>>fiigurre, imsshoww(g)要了解图图像文件件的其他他详细信信息，可可以使用用imffinffo函数数，其语语法结构构为：imfiinfoo fiilennamee其中，ffileenamme是存存储在磁磁盘中的的图像全全名。1.1.3 保保存图像像使用函数数imwwritte可以以将图像像写在磁磁盘上，该该函数语语法为：imwrritee (ff, fillenaame)下面的命命令可以以把图像像f写为为TIFFF格式式，在默默认情况况下使用用语句>>immwriite (f,'paatieent110_rrun11', 'tiif')或>>immwriite (f, 'ppatiientt10_runn1.ttif'')会将文件件保存到到当前的的工作目目录中。另外，如如果是JJPEGG图像，其其语法可可以为：imwrritee (ff, ''fillenaame.jpgg', 'quualiity'', qq)其中q是是一个在在0到1100之之间的整整数（qq越小，图图像退化化越严重重）。函数immwriite另另一种常常用但只只用于ttif图图像的语语法为：imwrritee (gg, fillenaame.tiff, commpreessiion, parrameeterr,ressoluutioon, collress rrowrres)其中，parrameeterr有三种种情况：为nonne表示无无压缩，packbits表示比特压缩包压缩（非二值图像默认参数），ccitt表示ccitt压缩（二值图像默认参数）。colres为列分辨率。1.1.4 数数据类型型见下表：表格 Error! No text of specified style in document.1 数据据类型。前前8项数数值类，第第9项字字符类，第第10项项逻辑类类。名称描述doubble双精度浮浮点数，范范围【-103308,103308】（88比特每每像素）unitt8无符号88比特整整数，范范围【00,2555】（11比特每每像素）unitt16无符号116比特特整数，范范围【00,6555355】（22比特每每像素）unitt32无符号332比特特整数，范范围【00,4229499672295】（44比特每每像素）int88有符号88比特整整数，范范围【-1288,1227】（11比特每每像素）int116有符号116比特特整数，范范围【-327768,327767】（22比特每每像素）int332有符号332比特特整数，范范围【-2144748836448,22147748336477】（44比特每每像素）singgle单精度浮浮点数，范范围【-10338，110388】（44比特每每像素）charr字符（22比特每每像素）logiicall值为0或或者1（11比特每每像素）1.1.5 图图像类型型主要分为为4种：l 亮度图像像（Inntennsitty iimagge）l 二值图像像（Biinarry iimagge）l 索引图像像（Inndexxed imaage）l RGB图图像（RRGB imaage）注意：一一副图像像的特征征是由数数据类和和图像类类这两者者来表征征的。例例如，“uniit8亮亮度图像像”表示一一副像素素都是uunitt8数据据类的亮亮度图像像。1.1.6 数数据类与与图像类类型之间间的转换换（1）数数据类间间的转换换通用语法法为：B=daata_claass_namme(AA)其中，ddataa_cllasss_naame可可以为表表1第一一列中的的任意一一项。（2）图图像类和和数据类类型间的的转换见下表：表格 Error! No text of specified style in document.2 IPPT中用用于进行行图像类类和数据据类转换换的函数数名称将输入转转换为有效的图图像数据据类im2uuninnt8unitt8logiicall,uiint88,uiint116,ddoubbleim2uuintt16unitt16logiicall,uiint88,uiint116,ddoubblemat22graaydoubble,范围【00,1】doubbleim2ddoubbledoubblelogiicall,uiint88,uiint116,ddoubbleim2bbwlogiicalluintt8,uuintt16,douublee要把一个个douublee类的任任意数组组转换成成取值范范围为0,11的归归一化ddoubble类类数组，可可以通过过函数mmat22graay完成成，其基基本语法法为：g=maat2ggrayy(A, AAminn, AAmaxx)其中，图图像g的的取值范范围在00（黑色色）到11（白色色）之间间。转换换时，小小于Ammin的的值在gg中变为为0，大大于Ammax的的值变为为1。1.1.7 数数组索引引在MATTLABB中，vv=11 3 5 77 9代表一一维行向向量，使使用转置置运算符符（.）可以以把行向向量转为为列向量量。要存取元元素数据据块，可可以使用用冒号（:）。例例如，下下列语句句代表存存取第三三个到最最后一个个元素：>>v(3:eend)ans= 5 77 9注意：符符号1:2：eend表表示索引引从1开开始计数数，步长长为2，直直到结束束。步长长可以为为负。例如，求求一个矩矩阵所有有元素之之和：>>s=summ(A(:)1.1.8 运运算符MATLLAB运运算符一一般分为为以下三三种主要要类别：l 执行数值值计算的的算术运运算符l 在数量上上比较操操作数的的关系运运算符l 执行函数数ANDD、ORR和NOOT的逻逻辑运算算符常用的运运算符和和算术函函数如下下面几个个表所示示：表格 Error! No text of specified style in document.3 数数组和矩矩阵运算算符运算符名称MATLLAB函函数注释和示示例+数组和矩矩阵加pluss(A,B)a+b,A+BB或a+A-数组和矩矩阵减minuus(AA,B)a-b,A-BB或a-A.*数组乘timees(AA,B)C=A.*B,C(II,J)=A(I,JJ)*BB(I*J)*矩阵乘mtimmes(A,BB)标准矩阵阵乘./数组右除除rdivvidee(A,B)C=A./B.数组左除除ldivvidee(A,B)C=A.B/矩阵右除除mrdiividde(AA,B)矩阵左除除mldiividde(AA,B).数组求幂幂poweer(AA,B)C=A.B,有C(I,JJ)=AA(I,J)B(II,J)矩阵求幂幂mpowwer(A,BB).|向量和矩矩阵转置置trannspoose(A)标准的向向量和矩矩阵转置置|向量和矩矩阵复共共轭转置置ctannspoose(A)标准的向向量和矩矩阵共轭轭转置+一元加upluus(AA)+A=00+A-一元减uminnus(A)-A=00-A:冒号表格 Error! No text of specified style in document.4 IPPT支持持的图像像算术函函数函数描述imaddd两幅图像像相加或或把常数数加到图图像imsuubtrractt两幅图像像相减或或从图像像中减去去常数immuultiiplyy两幅图像像相乘，在在相应的的像素对对间进行行；或图图像乘以以一个常常数imdiividde两幅图像像相除，在在相应的的像素对对间进行行；或图图像除以以一个常常数imaddsdiiff计算两幅幅图像间间的绝对对差imcoompllemeent对图像求求补imliicommb计算两幅幅或多幅幅的线性性组合表格 Error! No text of specified style in document.5 关关系运算算符运算符名称<小于<=小于等于于>大于>=大于等于于=等于=不等于表格 Error! No text of specified style in document.6 逻逻辑运算算符运算符名称&AND（与与）|OR（或或）NOT（非非）表格 Error! No text of specified style in document.7 逻逻辑函数数函数注释xor若两个操操作数逻逻辑上不不同，则则函数xxor返返回1；否则，返返回0all若一个向向量中的的所有元元素都非非零，返返回1；否则，返返回0any若一个向向量中的的任何元元素都非非零，返返回1；否则，返返回0。该该函数在在矩阵中中按列操操作1.2 实验内内容1.2.1 在在“素材”文件夹夹下，读读入“Figg2.119(aa).jjpg”并且提提取该图图像的基基本信息息，写出出实验步步骤和实实验结果果。图 1 图像像的读取取与显示示>> ff=immshoow(''Figg2.119(aa).jjpg'');>> wwhoss f Naame SSizee Byttes Cllasss AAttrribuutess f 11x1 8 dooublle 1.2.2读取取图像“Figg1.114(dd).jjpg”，使用用函数iimwrritee和suubpllot函函数分别别显示质质量参数数q=550,225,55,0的的结果如如图所示示。写出出实验步步骤。>> ff=immreaad(''Figg1.114(dd).jjpg'');>> ssubpplott(3,2,11),iimshhow(f);>> ssubpplott(3,2,22),iimshhow('buubblles550.jjpg'');>> iimwrritee(f,'buubblles225.jjpg'','qquallityy',225)>> ssubpplott(3,2,33),iimshhow('buubblles225.jjpg'');>> iimwrritee(f,'buubblles115.jjpg'','qquallityy',115)>> ssubpplott(3,2,44),iimshhow('buubblles115.jjpg'');>> iimwrritee(f,'buubblles55.jppg','quualiity'',5)>> ssubpplott(3,2,55),iimshhow('buubblles55.jppg');>> iimwrritee(f,'buubblles00.jppg','quualiity'',0)>> ssubpplott(3,2,66),iimshhow('buubblles00.jppg');观察哪副副图像压压缩最大大后没有有“伪轮廓廓”。此时键入入命令：>> iimfiinfoo buubblles225.jjpg会得到如如下信息息：ans = Fiilennamee: ''bubbblees255.jppg' FiileMModDDatee: ''15-十月-20112 115:449:444' FiileSSizee: 1138775 Forrmatt: ''jpgg' FFormmatVVerssionn: ''' Wiidthh: 7714 HHeigght: 6882 BiitDeepthh: 88 CollorTTypee: ''graaysccalee' ForrmattSiggnatturee: ''' NummberrOfSSamppless: 11 CoddinggMetthodd: ''Hufffmaan' CCodiingPProccesss: ''Seqquenntiaal' CCommmentt: 其中，使使用wiidthh乘以hheigght乘乘以BiitDeepthh再除以以8，得得到：7714*6822*8/8=44869948字字节。把把这个结结果除以以FilleSiize就就可以的的得到压压缩比：48669488/1338755=355.100。此时时压缩比比是在保保持图像像质量与与应用要要求一致致的前提提下得到到的，它它可以使使得单位位时间内内传输的的数据量量大约是是压缩前前的355倍。利用上述述原理，请请写出在在MATTLABB中如何何计算得得到压缩缩比。>> kk=immfinnfo('buubblles225.jjpg'');>> iimagge_bbytees=kk.Wiidthh*k.Heiightt*k.BittDeppth/8;>> ccomppresssedd_byytess=k.FilleSiize;>> ccomppresssioon_rratiio=iimagge_bbytees/ccomppresssedd_byytesscomppresssioon_rratiio = 335.0095441.2.3 如如图（Figg5.007(aa).jjpg）所示，把把原始图图像以ttifff格式存存储为无无压缩的的名为“sf.ttif”的图像像，且像像素不变变，但尺尺寸减小小为1.5*11.5英英寸。写写出实验验步骤。1.2.4 将将douublee类图像像“ >>>f=1 22;3 4”转换为为二值图图像，其其中1、22变为00,3、44变为11。请写写出实验验步骤。1.2.5 利利用数组组索引的的概念，写写出实现现下列图图像的实实验步骤骤。1.2.5 设设A=1 22 3;4 55 6; 7 8 99和BB=00 2 4; 3 55 6; 3 4 99计算算A=B,AA>=BB的结果果。1.2.6 设设A=1 22 0; 0 4 55,BB=11 -22 3; 0 1 11,计计算A&&B, xorr(A,B), alll(AA), anyy(A), aall(B), anny(BB)的结结果。实验2 亮度变变换与空空间滤波波2.1 实验原原理2.1.1 亮亮度变换换函数函数immadjjustt是对灰灰度图像像进行亮亮度变换换的基本本IPTT工具。语语法为：g = imaadjuust (f, llow_in hiigh_in, loww_ouut higgh_oout, ggammma)次函数将将图像ff中的亮亮度值映映射到gg中的新新值，即即将loow_iin至hhighh_inn之间的的值映射射到loow_oout至至higgh_oout之之间的值值。参数数gammma指指定了曲曲线的形形状，该该曲线用用来映射射f的亮亮度值，以以便生成成图像gg。2.1.2 直直方图处处理与函函数绘图图（1）函函数immhisst是处处理图像像直方图图的核心心函数，其其语法为为：h = imhhistt (ff, bb)其中，ff为输入入图像，hh为其直直方图hh(rkk)，bb是用于于形成直直方图的的“收集箱箱“的个数数（即灰灰度级的的个数）。bb默认为为2566。（2）绘绘制直方方图还有有另外三三个函数数barr、sttem和和ploot分别别以条形形图、杆杆状图和和线图的的形式体体现直方方图，其其语法为为：bar (hoorz, v, wiidthh)stemm (hhorzz, vv, collor_linnesttylee_maarkeer, filll)plott (hhorzz, vv, collor_linnesttylee_maarkeer)其中，hhorzz是一个个与v有有着相同同维数的的向量，它它包含水水平标度度值的增增量；wwidtth是一一个在00和1之之间的数数，为00则竖条条是简单单垂直线线，为11则竖条条明显；collor_linnesttylee_maarkeer来自自表格 008的三个个值，默默认颜色色为bllackk，默认认线条为为sollid（实线线），默默认标记记点形状状为ciirclle。表格 Error! No text of specified style in document.8 函数数steem和函函数poolt属属性。nnonee属性只只适用于于函数pplott，且必必须单独独指定。符号颜色符号线型符号标记k黑-实线+加号w白-虚线圆形r红:点线*星号g绿-.点划线.点b蓝nonee无线x叉形c青s方形y黄d菱形m品红nonee无标记（3）直直方图均均衡化由由工具箱箱中的函函数hiisteeq实现现，该函函数语法法为：g = hissteqq (ff, nnlevv)其中，ff为输入入图像；nleev是为为输出图图像指定定的灰度度级数，其其默认值值为644，但一一般设定定为2556。（4）直直方图匹匹配（也也称直方方图规定定化）使使用hiisteeq实现现，语法法为：g = hissteqq (ff, hhspeec)其中，ff为输入入图像，hhspeec为指指定的直直方图（一一个由指指定值构构成的行行向量），gg为输出出图像，其其直方图图类似指指定的直直方图hhspeec。（5）线线性空间间滤波是是利用工工具箱中中的函数数imffiltter来来实现，通用语法为：g = imffiltter (f, w, fiilteerinng_mmodee, bbounndarry_ooptiionss, ssizee_opptioons)其中，ff是输入入图像，ww为滤波波掩模，gg为掩模模结果，其其他参数数参见下下表。表格 Error! No text of specified style in document.9 函数数imffiltter的的选项选项描述滤波类型型corrr滤波通过过使用相相关来完完成。该该值默认认。connv边界类型型P输入图像像的边界界通过用用值P（无无引号）来来填充扩扩展。默默认为00。reppliccatee图像大小小通过复复制外边边界的值值来扩展展。symmmettricc图像大小小通过镜镜像反射射外边界界来扩展展。cirrcullar图像大小小通过将将图像看看成一个个二维周周期函数数的一个个周期来来扩展。大小选项项fulll输出图像像大小与与被扩展展图像大大小相同同。samme输出图像像大小与与输入图图像大小小相同。这这可将通通过滤波波掩模的的中心点点的偏移移限制到到原始图图像中包包含的点点来实现现。该值值默认。（6）线线性空间间滤波器器可通过过函数ffspeeciaal来实实现，语语法为：w = fsppeciial (typpe, pparaametterss)其中，typpe表示滤滤波器类类型，pparaametterss进一步步定义了了指定的的滤波器器，详见见下表。表格 Error! No text of specified style in document.10 函数数fsppeciial支支持的空空间滤波波器类型函数和参参数aveeraggefspeeciaal(aveeragge,rr c)即大大小为rr*c的的一个矩矩阵平均均滤波器器。默认认为3*3。若若由一个个数来代代替rr c，则表表示方形形滤波器器。disskfspeeciaal(dissk,r)一个圆圆形平均均滤波器器（包含含在2rr+1大大小的正正方形内内），半半径为rr。默认认半径为为5。gauussiianfspeeciaal(lapplacciann,rr c,siig)一一个大小小为r*c的高高斯低通通滤波器器，标准准偏差为为sigg（正）。默默认值为为3*33和0.5。若若由一个个数来代代替rr c，则表表示方形形滤波器器。lapplacciannfspeeciaal(lapplacciann,allphaa) 一一个大小小为3*3的拉拉普拉斯斯滤波器器，其形形状由aalphha决定定，allphaa是范围围内00,1的一个个数.aalphha默认认为0.5。loggfspeeciaal(logg,rr c,siig) 一个大大小为rr*c的的高斯-拉普拉拉斯（LLoG）滤滤波器，标标准偏差差为siig（正正）。默默认值为为5*55和0.5。若若用一个个数代替替r c，则则表示一一个方形形滤波器器。mottionnfspeeciaal(mottionn,leen,tthetta)围围绕一副副有着llen像像素的图图像线性性运动时时（类似似照相机机与景物物的关系系），输输出一个个滤波器器。运动动方向为为theeta，单单位为度度，即从从水平方方向逆时时针转动动的角度度，默认认值为99和0，表表示沿水水平方向向9个像像素的运运动。preewitttfspeeciaal(preewittt)输出出一个大大小为33*3的的preewittt掩模模wv，它它近似于于垂直梯梯度。水水平梯度度掩模可可以通过过置换结结果whh=wvv获得。sobbelfspeeciaal(sobbel)输出一一个大小小为3*3的掩掩模wvv，它近近似于垂垂直梯度度。水平平梯度掩掩模可以以通过置置换结果果sh=sv获得。unssharrpfspeeciaal(unssharrp,allphaa)输出出一个大大小为33*3的的非锐化化滤波器器。参数数alppha控控制形状状；它必必须是一一个大于于等于11.0的的数，默默认值为为0.22。（6）非非线性空空间滤波波，主要要介绍二二维中值值滤波函函数：g = meddfillt2 (f, mm n, ppadoopt)其中，数数组mm n定义了了一个大大小为mm*n的的邻域，中中值就在在该邻域域上计算算，而ppadoopt指指定了三三个可能能的边界界填充选选项之一一：zerro（默认认值）；symmmettricc，此地地f按照照镜像反反射方式式对称地地沿着边边界扩展展；inddexeed，若ff是dooublle类图图像，则则以1来来填充图图像，否否则以00来填充充图像，该该函数默默认形式式为：g = meddfillt2 (f)它使用一一个大小小为3*3的邻邻域来计计算中值值，并用用0来填填充图像像的边界界。2.2 实验内内容2.2.1 如如图是数数字乳房房X射线线原始图像像、负片图图像、亮亮度范围围为00.5 0.775的的图像和和令gaammaa=2的的图像。写写出实验验步骤来来显示这这4副图图像。（提提示：使使用immadjjustt函数）>> ff= iimreead('Fiig3.04(a).jpgg');>> ssubpplott(2,2,11),iimshhow(f)>> gg1=iimaddjusst(ff, 0 11,1 00);>> ssubpplott(2,2,22),iimshhow(g1)>> gg2=iimaddjusst(ff,00.5 0.775, 00 1);>> ssubpplott(2,2,33),iimshhow(g2)>> gg3=iimaddjusst(ff, , , 2);>> ssubpplott(2,2,44),iimshhow(g3)2.2.2读取取图像''Figg1.006(bb).jjpg''，并绘制制4种类类型直方方图。效效果如下下图所示示。>> ff=immreaad(''Figg1.006(bb).jjpg'');>> iimhiist(f);>> hh1=hh(1:10:2566);>> hhorzze=11:100:2556;>> bbar(horrze,h1)>> aaxiss(00 2555 00 1550000)>> sset(gcaa, ''xtiick'', 00:500:2555)>> sset(gcaa, ''ytiick'', 00:20000:150000)>> sstemm(hoorzee, hh1, 'fiill'')>> aaxiss(00 2555 00 1550000)>> sset(gcaa,'xxticck',0:50:2555)>> sset(gcaa,'yyticck',0:20000:1150000)>> pplott(h)>> aaxiss(00 2555 00 1550000)>> sset(gcaa, ''xtiick'',00:500:2555)>> sset(gcaa, ''ytiick'',00:20000:150000)2.2.3 如如图（Figg3.115(aa)1ttop.jpgg）是电子子显微镜镜下放大大7000倍的花花粉图像像，为得得到如图图所示结结果，请请写出实实验步骤骤。>> ff=immreaad(''Figg3.115(aa)1ttop.jpgg');>> iimshhow(f)>> ffiguure,imhhistt(f)>> yylimm('aautoo')>> gg=hiisteeq(ff,2556);>> ffiguure,imsshoww(g)>> ffiguure,imhhistt(g)>> yylimm('aautoo')>> hhnorrm=iimhiist(f)./nuumell(f);>> ccdf=cummsumm(hnnormm);>> pplott(x,cdff)>> aaxiss(00 1 0 11)>> sset(gcaa,'xxticck',0:.2:11)>> sset(gcaa,'yyticck',0:.2:11)>> xxlabbel('innputt inntennsitty vvaluues'','ffonttsizze',9)>> yylabbel('输出出亮度值值',''fonntsiize'',9)>> ttextt(0.2,00.5,'变换换函数'','ffonttsizze',9)实验3 图像复复原复原的目目的是在在预定义义的意义义上改善善给定的的图像。图图像增强强主要是是主观的的处理，图图像复原原主要是是客观的的处理，趋趋向于将将退化模模型化，并并用相反反的处理理来恢复复原图像像。3.1 实验原原理3.1.1噪声声模型（1）工工具箱使使用immnoiise函函数来污污染一副副图像，基基本语法法为;g = imnnoisse (f, typpe, parrameeterrs)其中，ff是输入入图像，该函数常见语法形式如下：l g = imnnoisse (f, gauussiian, mm, vvar)将均值值为m，方方差为vvar的的高斯噪噪声加到到图像ff上。默默认值为为0,、方方差为00.011的噪声声。l g = imnnoisse (f, loccalvvar, VV)将均均值为00、局部部方差为为V的高高斯噪声声添加到到图像ff上，其其中V是是与f大大小相同同的一个个数组，它它包含每每个点的的理想方方差值。l g = imnnoisse (f, loccatllvarr, iimagge_iinteensiity, vaar)将均值值为0的的高斯噪噪声添加加到图像像f上，其其中噪声声的局部部方差vvar是是图像ff的亮度度值的函函数。参参量immagee_inntennsitty必须须包含在在范围0,11内的的归一化化亮度值值，它与与varr是大小小相同的的向量。本本函数绘绘制出噪噪声方差差和图像像亮度的的函数关关系。l g = imnnoisse (f, sallt && peeppeer, dd)用椒椒盐噪声声污染图图像f，其其中d是是噪声密密度（即即包含噪噪声值的的图像区区域的百百分比）。因因此，大大约有dd*nuumell(f)个像素素受到影影响，默默认噪声声密度为为0.005。l g = imnnoisse (f, speecklle, vvar)用方程程g=ff+n*f将乘乘法性质质噪声添添加到图图像f上上，其中中n是均均值为00、方差差为vaar的均均匀分布布的随机机噪声。vvar默默认为00.044。l g = imnnoisse (f, poiissiion)从数据据中生成成泊松噪噪声，而而不是将将人工的的噪声添添加到数数据中。为为了遵守守泊松统统计，uunitt8类和和uinnt166类图像像的亮度度必须和和光子（或或其他任任何的量量子信息息）的数数量相符符合。当当每个像像素的光光子数量量大于6655335（但但必须小小于10012）时时，就要要使用双双精度图图像。亮亮度值在在0和11之间变变化，并并且对应应于光子子的数量量除以110122。（2）使使用指定定的分布布产生空空间随机机噪声。空间噪声值是随机数，以概率密度函数（PDF）或相应的累积分布函数（CDF）表征。利用已有M文件“imnoise2.m”。（3）周周期噪声声，一般般产生于于图像采采集过程程中的电电气和/或电机机干扰。利利用已有有M文件件“imnnoisse3.m”完成实实验。（4）估估计噪声声参数。函数statmoments计算均值和n阶中心矩，并返回向量v。由于0阶矩阵总为1,1阶矩总为0，所以statmoments忽略这两个矩，语法形式如下：v, unvv = sttatmmomeentss (pp, nn)其中，pp为直方方图向量量，n是是计算的的矩的数数量。噪声参数数通常直直接由带带噪声的的图像或或一组带带噪声的的图像来来估计。在MATTLABB中，使使用函数数roiipolly来选选择一个个感兴趣趣区域（RROI），语法为：B = roiipolly （ff, cc, rr）其中，ff是我们们感兴趣趣的图像像，c和和r是多多边形顶顶点的相相应列坐坐标和行行坐标，输输出的BB是一个个二值图图像，它它是一个个与f大大小相同同且在感感兴趣区区域外为为0、在在感兴趣趣去野外外为1的的二值图图像，它它是一个