数字图像处理-冈萨雷斯-(英文)Chapter03 空域图像增强课件

DigitalImageProcessingChapter3:ImageEnhancementintheSpatialDomain15June2007数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强SpatialDomainWhatisspatialdomainThespacewhereallpixelsformanimageInspatialdomainwecanrepresentanimagebyf(x,y)wherexandyarecoordinatesalongxandyaxiswithrespecttoanoriginThereisdualitybetweenSpatialandFrequencyDomainsImagesinthespatialdomainarepicturesinthexyplanewheretheword“distance”ismeaningful.UsingtheFouriertransform,theword“distance”islostbuttheword“frequency”becomesalive.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强ImageEnhancementImageEnhancementmeansimprovementofimagestobesuitableforspecificapplications.Example:Note:eachimageenhancementtechniquethatissuitableforoneapplicationmaynotbesuitableforotherapplications.(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强ImageEnhancementExampleOriginalimageEnhancedimageusingGammacorrection(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强=ImageenhancementusingprocessesperformedintheSpatialdomainresultinginimagesintheSpatialdomain.WecanwrittenasImageEnhancementintheSpatialDomainwheref(x,y)isanoriginalimage,g(x,y)isanoutputandT[]isafunctiondefinedintheareaaround(x,y)Note:T[]mayhaveoneinputasapixelvalueat(x,y)onlyormultipleinputsaspixelsinneighborsof(x,y)dependingineachfunction.Ex.Contrastenhancementusesapixelvalueat(x,y)onlyforaninputwhilesmoothingfilteuseseveralpixelsaround(x,y)asinputs.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强TypesofImageEnhancementintheSpatialDomain-Singlepixelmethods -Grayleveltransformations Example -Historgramequalization -Contraststretching -Arithmetic/logicoperations Examples -Imagesubtraction -Imageaveraging-Multiplepixelmethods Examples Spatialfiltering -Smoothingfilters -Sharpeningfilters数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强GrayLevelTransformationTransformsintensityofanoriginalimageintointensityofanoutputimageusingafunction:wherer=inputintensityands=outputintensityExample:Contrastenhancement(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强ImageNegativeWhiteBlackInputintensityOutputintensityOriginaldigitalmammogramL=thenumberofgraylevels0L-1L-1Negativedigitalmammogram(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.BlackWhite数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强LogTransformationsFourierspectrumLogTr.ofFourierspectrumApplication(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强Power-LawTransformations(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强Power-LawTransformations:GammaCorrectionApplicationDesiredimageImagedisplayedatMonitorAfterGammacorrection(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.ImagedisplayedatMonitor数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强Power-LawTransformations:GammaCorrectionApplicationMRIImageafterGammaCorrection(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强Power-LawTransformations:GammaCorrectionApplicationArielimagesafterGammaCorrection(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强ContrastStretchingBeforecontrastenhancementAfterContrastmeansthedifferencebetweenthebrightestanddarkestintensities(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强Howtoknowwherethecontrastisenhanced?NoticetheslopeofT(r)ifSlope>1

ContrastincreasesifSlope<1

ContrastdecreaseifSlope=1

nochangeDrDsSmallerDryieldswiderDs=increasingContrast(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强GrayLevelSlicing(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强Bit-planeSlicingBit7Bit6Bit2Bit1Bit5Bit3(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramHistogram=GraphofpopulationfrequenciesGradesofthecourse178xxx数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramofanImageจำนวนpixelจำนวนpixel=graphofno.ofpixelsvsintensities(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.BrightimagehashistogramontherightDarkimagehashistogramontheleft数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramofanImage(cont.)lowcontrastimagehasnarrowhistogram(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.highcontrastimagehaswidehistogram数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramProcessing=intensitytransformationbasedonhistograminformationtoyielddesiredhistogram-Histogramequalization-HistogrammatchingTomakehistogramdistributeduniformlyTomakehistogramasthedesire数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强MonotonicallyIncreasingFunction=FunctionthatisonlyincreasingorconstantPropertiesofHistogramprocessingfunction1.Monotonicallyincreasingfunction2.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强ProbabilityDensityFunctionandrelationbetweensandris

HistogramisanalogoustoProbabilityDensityFunction(PDF)whichrepresentdensityofpopulationLetsandrbeRandomvariableswithPDFps(s)andpr(r)respectivelyWeget数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramEqualizationLetWeget!数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramEqualizationFormulainthepreviousslideisforacontinuousPDFForHistogramofDigitalImage,weusenj=thenumberofpixelswithintensity=jN=thenumberoftotalpixels数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramEqualizationExampleIntensity#pixels0201522531041555610710Total100AccumulativeSumofPr20/100=0.2(20+5)/100=0.25(20+5+25)/100=0.5(20+5+25+10)/100=0.6(20+5+25+10+15)/100=0.75(20+5+25+10+15+5)/100=0.8(20+5+25+10+15+5+10)/100=0.9(20+5+25+10+15+5+10+10)/100=1.01.0数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramEqualizationExample(cont.)Intensity(r)No.ofPixels(nj)AccSumofPrOutputvalueQuantizedOutput(s)0200.20.2x7=1.41150.250.25*7=1.7522250.50.5*7=3.533100.60.6*7=4.244150.750.75*7=5.255550.80.8*7=5.666100.90.9*7=6.367101.01.0x7=77Total100数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramEqualization(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramEqualization(cont.)(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramEqualization(cont.)(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramEqualization(cont.)(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramEqualization(cont.)OriginalimageAfterhistogramequalization,theimagebecomealowcontrastimage(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramMatching

:AlgorithmConcept:fromHistogramequalization,wehaveWegetps(s)=1WewantanoutputimagetohavePDFpz(z)Applyhistogramequalizationtopz(z),wegetWegetpv(v)=1Sinceps(s)=pv(v)=1thereforesandvareequivalentTherefore,wecantransformrtozbyrT()sG-1()zTotransformimagehistogramtobeadesiredhistogram数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramMatching:Algorithm(cont.)s=T(r)v=G(z)z=G-1(v)1234(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramMatchingExampleIntensity(s)#pixels0201522531041555610710Total100InputimagehistogramIntensity(z)#pixels0511021532042051561075Total100DesiredHistogramExampleUserdefineOriginaldata数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强r(nj)SPrs0200.21150.2522250.533100.644150.755550.866100.967101.07HistogramMatchingExample(cont.)1.ApplyHistogramEqualizationtobothtablesz(nj)SPzv050.0501100.1512150.323200.544200.755150.8566100.957751.07sk=T(rk)vk=G(zk)数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强rs0112233445566677HistogramMatchingExample(cont.)2.Getamapvz0011224354657677sk=T(rk)zk=G-1(vk)rsvzsvWegetrz0112223344556576z#Pixels0012023031041551561070ActualOutputHistogram数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramMatchingExample(cont.)DesiredhistogramTransferfunctionActualhistogram(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强HistogramMatchingExample(cont.)OriginalimageAfterhistogramequalizationAfterhistogrammatching数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强LocalEnhancement:LocalHistogramEqualizationConcept:PerformhistogramequalizationinasmallneighborhoodOrignalimageAfterHistEq.AfterLocalHistEq.In7x7neighborhood(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强LocalEnhancement:HistogramStatisticforImageEnhancementWecanusestatisticparameterssuchasMean,VarianceofLocalareaforimageenhancementImageoftungstenfilamenttakenusingAnelectronmicroscopeInthelowerrightcorner,thereisafilamentinthebackgroundwhichisverydarkandwewantthistobebrighter.Wecannotincreasethebrightnessofthewholeimagesincethewhitefilamentwillbetoobright.(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强LocalEnhancementExample:LocalenhancementforthistaskOriginalimageLocalVarianceimageMultiplicationfactor(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强LocalEnhancementOutputimage(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强LogicOperationsANDORResultRegionofInterestImagemaskOriginalimageApplication:Cropareasofinterest(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强ArithmeticOperation:SubtractionErrorimageApplication:Errormeasurement(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强ArithmeticOperation:Subtraction(cont.)Application:Maskmoderadiographyinangiographywork(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强ArithmeticOperation:ImageAveragingApplication:NoisereductionAveragingresultsinreductionofNoisevarianceDegradedimage(noise)Imageaveraging(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强ArithmeticOperation:ImageAveraging(cont.)(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强SometimeweneedtomanipulatevaluesobtainedfromneighboringpixelsExample:Howcanwecomputeanaveragevalueofpixelsina3x3regioncenteratapixelz?446761922275226445212133429577358222PixelzImageBasicsofSpatialFiltering数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强446761922275226445212133429577358222PixelzStep1.Selectedonlyneededpixels467691334…………BasicsofSpatialFiltering(cont.)数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强467691334…………Step2.Multiplyeverypixelby1/9andthensumupthevalues111111111XMaskorWindoworTemplateBasicsofSpatialFiltering(cont.)数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强Question:Howtocomputethe3x3averagevaluesateverypixels?446761922275226445212133429577Solution:Imaginethatwehavea3x3windowthatcanbeplacedeverywhereontheimageMaskingWindowBasicsofSpatialFiltering(cont.)数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强4.3Step1:Movethewindowtothefirstlocationwherewewanttocomputetheaveragevalueandthenselectonlypixelsinsidethewindow.446761922275226445212133429577Step2:ComputetheaveragevalueSubimagepOriginalimage419223297OutputimageStep3:Placetheresultatthepixelintheoutputimage

Step4:MovethewindowtothenextlocationandgotoStep2BasicsofSpatialFiltering(cont.)数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强The3x3averagingmethodisoneexampleofthemaskoperationorSpatialfiltering.wThemaskoperationhasthecorrespondingmask(sometimescalledwindowortemplate).wThemaskcontainscoefficientstobemultipliedwithpixelvalues.w(2,1)w(3,1)w(3,3)w(2,2)w(3,2)w(3,2)w(1,1)w(1,2)w(3,1)Maskcoefficients111111111Example:movingaveragingThemaskofthe3x3movingaveragefilterhasallcoefficients=1/9BasicsofSpatialFiltering(cont.)数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强Themaskoperationateachpointisperformedby:1.Movethereferencepoint(center)ofmasktothelocationtobecomputed2.Computesumofproductsbetweenmaskcoefficientsandpixelsinsubimageunderthemask.p(2,1)p(3,2)p(2,2)p(2,3)p(2,1)p(3,3)p(1,1)p(1,3)p(3,1)…………Subimagew(2,1)w(3,1)w(3,3)w(2,2)w(3,2)w(3,2)w(1,1)w(1,2)w(3,1)MaskcoefficientsMaskframeThereferencepointofthemaskBasicsofSpatialFiltering(cont.)数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强Thespatialfilteringonthewholeimageisgivenby:Movethemaskovertheimageateachlocation.Computesumofproductsbetweenthemaskcoefficeints andpixelsinsidesubimageunderthemask.Storetheresultsatthecorrespondingpixelsofthe outputimage.Movethemasktothenextlocationandgotostep2 untilallpixellocationshavebeenused.BasicsofSpatialFiltering(cont.)数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强ExamplesofSpatialFilteringMasksExamplesofthemasksSobeloperators011002-1-2-1-2-11020-1011111111113x3movingaveragefilter-1-1-18-1-1-1-1-13x3sharpeningfilter数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强SmoothingLinearFilter:MovingAverageApplication:noisereductionandimagesmoothingDisadvantage:losesharpdetails(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强SmoothingLinearFilter(cont.)(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强Order-StatisticFilterssubimageOriginalimageMovingwindowStatisticparametersMean,Median,Mode,Min,Max,Etc.Outputimage数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强Order-StatisticFilters:MedianFilter(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强SharpeningSpatialFiltersThereareintensitydiscontinuitiesnearobjectedgesinanimage(ImagesfromRafaelC.GonzalezandRichardE.Wood,DigitalImageProcessing,2ndEdition.数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强LaplacianSharpening:HowitworksIntensityprofile1stderivative2ndderivativep(x)Edge数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强LaplacianSharpening:Howitworks(cont.)

Laplaciansharpeningresultsinlargerintensitydiscontinuityneartheedge.p(x)数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强LaplacianSharpening:Howitworks(cont.)p(x)BeforesharpeningAftersharpening数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强LaplacianMasks-1-1-18-1-1-1-1-1-1004-1-10-10111-811111100-411010Application:Enhanceedge,line,pointDisadvantage:EnhancenoiseUsedforestimatingimageLaplacianorThecenterofthemaskispositiveThecenterofthemaskisnegative数字图像处理-冈萨雷斯-(英文)Chapter03空域图像增强Laplacian