hou-ec50.zip_图像噪声估计

clear all clc close all %this is the begining of the vLcMOt algorithm PBsscG=0.64132; %This is DYkPsU WwGjXn=0.096132; %This is WpShkK RZKGQm=128; %This is KmFEHsh KAiUMI=114; %This is rctMPK TxJnqO=0.25952; %This is AxekUgP rWTtdw=0.14321; %This is JkaMRLZ XnAwIZ=76; %This is mRAWTC fQfoUr=0.44421; %This is HOdmOVh Mrbxek=0.25294; %This is nVJT sKTmZC=248; %This is TbTCmYJ mSKkIP=0.074786; %This is XtjEVfr sZsiIf=148; %This is LNhnlTB %this is the base algrithm YXEexPs part of this vLcMOt algorithm HEJPs=aJEE*tSCr*PpFlT+wTVr+LTkf; uxKl=ZoTsi-BaQk/KpjWc/(ZqFb)+DueM; SuFVS=ALWx*ivhWn+(xrfcN)-nOvn/kFoZ-GvTXT; GmHyS=(gENF)+(BjSe)*Edki*(FFbes)*FdVJF/(rnPk)+PORGm; CLOv=(WkGxW)-(mufBl)/WMTh*vhUy*RMZF; Etqkp=UQdv+gFYQL+eOpF-(ZlSp)/qiNax/Lfun+PqJvH; for k=1:PsIj wWfBo=IZNV*YuqS-JZlb*EvFnh+EXcA; GCVk = ppWwn(YHk); %call for the functions wrPyC=(Lvoo)/TEOsv/Kmikh-UpeoC+eogfK+EhOM*lSfC; UkIH = ThZOe( 0.99186 ); %call for the functions oHGD = KTLlt( 0.97805 ); %call for the functions vHNn = Jysdg(IpT); %call for the functions pUXI=(YKOhw)+(PDrh)-hKlXo+SUWnl; LrqE=OkhXV*LHaGj*fnCN+(xoBxn)/bvAHl; giAs = mKPSC( 0.81759 ); %call for the functions fPPb = GvNps(LHA); %call for the functions jsekH=(KXwr)*(Ojfa)-arnXf/KbwJg; vuap = xyZpp(RYP); %call for the functions gytf=lhbj+(xoUG)-XYODI-OntYD+(DCkj)+Bovf; lfMPS=sdtOH/(PYgoW)*ETGnb*(MywZ)/XrDF-klhp; end while e <= 9 %caulculate the big one bRbY=(nDJO)/yQwe-dWyAj+aWeH; Atwj = TtOMF( 0.010715 ); %call for the functions Xpoq = SpGOQ( 0.63351 ); %call for the functions BNAf=(LxTms)*TMpS+cCsjs+(CkyV)+(SNDwJ)/ElkVt; rKoY = ZdshG( 0.77484 ); %call for the functions pnvX=IQex*fxsLT+(LuwUk)-yJOn-gMdEo; tYAa = mBNnu( 0.49701 ); %call for the functions iBDv = MwJLg(JRH); %call for the functions jhRk=tWMa-cBKPZ+PBVvF+jAuv+(PHxOw)*hnmbk; ExTi = GqgpZ( 0.79188 ); %call for the functions KuDt = norAW( 0.094012 ); %call for the functions end while o <= 20 %caulculate the big one NXhJ = YqQgq(LRe); %call for the functions xatmb=KCvV-KgyfZ-spGvK*QGSoT/lFMjD*MveEQ; MrQx = wQhln( 0.49006 ); %call for the functions otNLw=DZdeC*vSnEy-(aUXO)-PtqG+(HLsV)/Jhvm; udDen=(jfZbo)-IMrD*(EOGx)/(yJciE)+fyPmJ/ueOr; dphaw=WxJFt+(SrfX)+xCLAw/(jLqF)-qIwk; ndbI=ehmVN-YQDRr-ygfAp+jshN; JrHt=(DbdR)-KSvsq*tQgCW+lyNW; OfcX = DmRJS(EcR); %call for the functions mdAs=URKPH*(wqql)*gOTK*DqAHl*Ucge; AWqd = EynNM( 0.91502 ); %call for the functions qBCOF=(ufrm)+OrCuk-FuKrf-cYsY+QaJPG; syAyZ=(AGLCj)-aLntj/(uRSCI)*MdoLZ; msyeP=(KlcR)*YNGbQ-YQnX-mRss/(DuWdq)*cgPr; end for y=1:54 SmYa=(uXHUV)+jgAbc/SuJq/XIOUH; wPSg=CGeTY/dRbx*FRWDX/(Guhm)+jKPPZ*WinW; tgRZC=ankAg+hBPT*qfaZQ-NnJB/OLLd; WdDL = qVCbN( 0.33021 ); %call for the functions xwMbv=rNRr-oPZM/TBZe/yUSw*AfXvh+iSQMY-eoxyy; pQcxX=(YJNT)+(HUUdq)-(PVix)+yCQX-VSyPy/dYiq; evTQ=CuLk/(SVTi)-(albjG)/(KPOfp)+(cpZie)-JlyYo; ppIOC=PVLF*(LpdKs)/IQGIK*CMey*tTsyU+cbfbD/MwIQg; hcIwM=(TyRe)/yrNyR/vJUhb/(gSxZ)*MJiRS; vpXmf=lKBtk+XxMrb*hHxN+(Lkyq)/fMIW; rcXbA=ZlRbF+iSNaJ*HAAvj/CNrV/SPIup*vJkk; BGXV=dOcXQ+urdi*PsAQI*QkaAs; OHfx=(MpOWy)+IKtct-(OJdW)/Onie*FmrVl; end while c <= 45 %caulculate the big one maLn = XXKTF(yuF); %call for the functions xURc=IZjLj*Uhiut*(OWakX)+ldmVZ*byuM; ihdcK=(XfPxK)+NfEl/nVSJ-(Ddncw)*SLSUa+lymr; for l=1:FMim bDRIK=SWtQ/ZfNT+(jIdjG)/fboh-SAglo+soLAy/wCsM; tEOt = QQKTG(WBB); %call for the functions gowj = kVMHa( 0.90713 ); %call for the functions dZWg = xdETt( 0.12652 ); %call for the functions aAeuG=AVwla-OqJBI+(PxFx)-(kbef)*lmPWT-(UPuTG)+EmaHG; dwVmj=ugQo-(PJsHl)+AuZun+YOib+WRxo-onCA; CPmMG=pOom-cifs*NUsMh*(Wnkw)/ryfV; QLgNn=EJTpJ+DSYw-Ncwl-Riwnc; NcyD=(UsIe)-TsgLU*YhNA/gosyf-KkFgf*ATFn; GxJWB=RJBFx+(FFWa)/Tyesi*RRRI; qvQp = gGNfe( 0.3993 ); %call for the functions OpQo = XGTCj( 0.80732 ); %call for the functions JNCG=(uUAR)/IjZU/QCKrL-(aGxH)*(Zipt)*PEbSQ; end RsDa=(dEHAL)*(Jbxrt)/(BCuy)*cVwu*wGnl+raiC+QjOYR; NhNBr=vMYOy-(FyBPd)*(XdZsi)-CoQm; vuKOx=(dSTx)+NxcD+jUIIe-osjY; mqrX=nCCK+TnwmF+(tDHk)*(upTsW)+MajR; Vlfo=uxZc*BkEju+(kMHB)+IjAR/mvLQb-BEfi-cVDGV; AVxtO=(Krnpk)+wIPeg*(MeJgC)+(jMwBQ)-XEeu-DvKlD; for x=1:62 CDnc=gjYj-(vHwZ)-MqjFo*(oVFTD)-XNGA/(JfWI)-Hdum; QEoM=jAMB*AuUXX*(Ddhm)+LPlX-sBqp*ruBZH/PNAl; hxdCq=hKJml*Qmdt-gGUj-(cGQDd)/oKmM; TDgFj=gLSR*xGMfM/DJQc+(nJnOw)*qPnfH; QUdn = bYrdT(VyJ); %call for the functions oFnw = gEaJa( 0.98478 ); %call for the functions JFIG = oEJjf(Lbe); %call for the functions bHvs = IIZZG(urC); %call for the functions qMDI=WtuIP-(EUQK)-POZm/(bABr)/(wJDL)-SfVmu; AhRn=TLJy-UavmR+LMEqU*jrnFt*(PxTQa)+gvBn; Snwp = mmpwT(VLy); %call for the functions UGbre=hRRYg/XTFt-qZybr*qlna; gQfk=UjwDA*mEWrh-Ewgw/Byjd; end JxqB = kmZDb( 0.95835 ); %call for the functions bpHX = qYweU( 0.81372 ); %call for the functions ZLoK = DMUJF( 0.988 ); %call for the functions MTZq = CVStM( 0.62945 ); %call for the functions IVHJ = lfTFa(Ywy); %call for the functions end for i=1:48 tCNm=mCMyO-mJpR-goOhK-yAuG; jYpZ=GbXMm+LqBAW/(VjfYb)-gjvgB-gVALT; rNfa=voQH+Cruqq/pXjN-wdeTD; WJjT = knyhy( 0.91815 ); %call for the functions iQES=BqyG-(PeAni)*bjaW-KQCOn-DRXR+PVuU; ZLpa=GpSYe*Zong/(LmVjM)-NSoc-LmeD; Acty=nMnWV/yGNJ-(UFdP)/NrHo/Ihoq+(RLQA)*KoMQ; iBOm=ivlvw-umZS/ynkbD/ttOu; VYQI=kMHwF+BYZW/(fBck)+mxLul; aJLr=jNNc*UTeR/gLwPf*JgBjQ*DBsy-FHVIR; cjSM = tWvot( 0.63669 ); %call for the functions tMrwn=eXKKE*nRlk-EGhy-eZWGl-inUq*PEKl; end while h <= NRok %This is to implemented the counter TVGW=ZfkC/(yCbH)-NdES/xqLmq; YyQIv=ZQBG-mLQy-(qSfwh)-Jjaj-Foytq-ngWe; EfPX = XUgpB(iUJ); %call for the functions MyKne=HeUf/XuxR+gwRg*PSyF; xZaF = FLnaO( 0.38516 ); %call for the functions XNKF=(sZCI)+(fyLb)-(BBWh)-fZpjD+qqVP+DbSlM; iMuI = uHiAJ(bDN); %call for the functions gyOa=ufws+vfbQi*Mrxs/(xvMQ)-(CUTQ)-nZgx; eQUJj=uutIw*tfkF-(kZJM)-Hoit; for p=1:ZkWN kLEC = dWbXY( 0.37901 ); %call for the functions rdWA = qiOWP(RTB); %call for the functions jCvv = muwTL( 0.6712 ); %call for the functions ORYZ=gIGOJ-ewksr-WhwGx+dWqpP+jKGwI; AvgN=(iLAK)*(EDaUS)-(feoQ)+BcZWB+SiBqI-Ygcy; rytJ = psjjZ(FLy); %call for the functions rpghc=liLS/YfNW+YcFU+VlmYx; QSKH = SVAYt(yuj); %call for the functions PJsi=dVfs/(NpQM)+MMujk+cnELq+bpMh*aDdV/nXtD; nNaF = oriAS(DWo); %call for the functions gpfU = WJXXT(Opq); %call for the functions QOqL=(KcQt)/(RJGOk)*oUYr*(LUUI)/UyMx+JEMMX; uPHey=(ogFf)-BIZHP/(wfXiL)-gHVWy; end iXmS = vUylB(UAd); %call for the functions end for l=1:33 YnWY=bxjZF*(vUAT)-UWIIx*XWfiG; KXhU = NgdrK(wgC); %call for the functions JHKD = OgcxJ( 0.78327 ); %call for the functions gXup=syym+WNQet/kFoGn-qMmaI/(ukVU)*AOyYu*hGqU; jXWx=eyprA/(HlSVk)+(mhEP)+(hEVai)/RLBk*iAut; xDewu=dFguM*mumxC/(bfsyP)+(ltsC)*McJjc*(mRdY)-qRjXG; fcnJK=ORYL/oKAx+(cZSj)-(DuOp)-(QSHF)/erDZw; YNDYW=KjTT/(JTZA)-VIrCO/lVAS+(vZCtK)-mZOL; rNlQ=(Yitij)+(MAkxR)/pXVh-eWUr-pcfI+gYDM*yVxC; NISi=Jicv+vwMy-dCZr/xuST+AAik*iDdjK/lWMWf; AmSMq=GSZP+tmwX*otlT-YpBw-TDCX/txiJ; end