水中气泡尺度计算的图像处理方法

361Vo.l36,No.120061LASER&INFRAREDJanuary,2006:1001-5078(2006)01-0075-03葛卫龙,张晓晖,雷选华(,430033):舰船尾流场中气泡幕的存在,使得舰船尾流区域的海水的光特性与周围海水的光特性有很大差异,利用尾流场光学特性的异变,可作为水下航行体跟踪识别和定位目标的重要手段尾流中气泡大小浓度对于研究对光在海水中传输特性有很重要的意义文中给出了在实验室条件下用图像处理方法计算水中微小气泡尺度的方法,实验结果对于研究舰船尾流场的光学特性有实际意义:尾流;气泡;图像处理:O435;TP391.4:AComputingUnderwaterBubbleRadiusUsingImageProcessingGEWe-ilong,ZHANGXiao-hu,iLEIXuan-hua(DepartmentofWeaponryEngineering,NavalUniversityofEngineering,Wuhan430033,China)Abstract:Airbubblescreeninshipwakesmakesseawaterinwakesopticalcharacteristicmoredifferentwithothers,whichcanbeusedtotrack,identifyandorientateobjectsintheseawater.Theradiusandthicknessofairbubblesisveryimportanttolighttransmissioncharacteristicinwakes.Themethodofcomputingtinybubbleradiususingimageprocessingisgiveninthearticle,andtheresultshavepracticalmeaningtostudyopticalcharacteristicinshipwakes.Keywords:wakes;airbubble;imageprocessing1,,,,[1][2],,,,(1),,,5Lm10Lm20Lm25Lm30Lm,LZB-3WB,60~600ml/min,[34],,,,,,,,,,,:973:(1979-),,,,:2005-04-05;:2005-06-301Fig.1bubbleproducinginstrument2,QV-5700Casio,2560@1920dp,i,AF,6cm,,MATLAB,,,,22Fig.2programflowchart33,60ml/min3.1像素尺寸的定标31cm,:(x1y1)(x2y2)(x3y3)(x4y4),1cm:DB1=((x1-x2)2+(y1-y2)2+(x2-x3)2+(y2-y3)2+(x3-x4)2+(y3-y4)2)/3DB=(10000/DB1)Lm,DB1=349.2,DB=28.6Lm3.2选择要研究的气泡区域3,3,,,4,wienerwiener,,3Fig.3originalimage4Fig.4chosenswatchfield3.3气泡边缘提取和改进的填充算法sobel5,,1,0,,,,,:484,;8,76365,sobel,48[5],,(x1,y1)(x2,y2),(x1,y1)(x2,y2),65sobelFig.5bubbleedgepick-upusingsobeloperator6Fig.6bubbleedgefillingin,(),4(4),;1.2,;0.7,(,)784,,920Lm60ml/minLm27Fig.7bubbleafterfilteranditsareadistributinghistogramLmLm8Fig.8bubblediameterhistogram9,n=469Fig.9bubblediameterdistributingrulen=4695,,,,,,:(1)20Lm60200400312Lm386Lm407Lm,,,,(2):[1],,.He-Ne[J].,2001,9(1):20-23.[2],,.[J].,2000,8(1):11-15.[3]JinWu.Bubblepopulationsandspectrainnear-surfaceo-cean:summaryandreviewoffieldmeasurements[J].J.GeophRes,1981,86(C1):457-463.[4]MedvinH,BreitsND.Ambientandtransientbubblespectraldensitiesinquiescentseasandunderspillingbreakers[J].J.GeophRes,1989,94(C9):12751-12759.[5],.MATLAB[M].,2002,112-113.77No.12006