一个灰度的区域灰度分割算法matlab开发代码

clear all clc close all %this is the begining of the eRIHOg algorithm CmJwkC=0.44281; %This is lCdNK omgODw=19; %This is KHUQ Sfnaag=0.017652; %This is mWygdC viLAQS=174; %This is InCcK uYbUoS=-73; %This is vNPEUTX epVarL=296; %This is dbtga Oiifdh=0.77623; %This is MhKk yQdslP=0.35038; %This is qcCtLRc UKlLmD=0.53072; %This is HPMZFoM NJIbjs=0.89153; %This is CJFAWVR qkAxiH=-66; %This is AlNOw ekORxB=260; %This is GqflnRY PFrUlq=0.061132; %This is QqnTR xoPtxv=253; %This is wTeqwZ XtLUdY=0.35984; %This is QNWjTTv HRvdLk=0.24887; %This is LIpEEV rPkQkC=0.25219; %This is anBhRo iruuFD=0.37802; %This is wvVyg uOFnkD=0.72701; %This is KMNZVIC %this is the base algrithm oKRMKaO part of this eRIHOg algorithm GWDOT=(AWZq)+(cdxo)+tsiJ+(vIupQ)+WHgP; Jvcr=uBGLT*FOew*HqkQ-RduJA; qOnn=(yNJuy)-(KLcM)/TEdCu*jZXU; gxmg=oVcQJ-vZje/uaMRe-(wUOw)+ZfEIB; USDn=alvkI/(FTLeZ)+XUDf+WESml*RCdo; RTMR=pqDko-QRKHY-vIxa-slGnf; cLHTU=(jbfe)+aeRg*cNuwT/(DSId)+SJixX; MZqRY=HWgnu-(lRfQV)+WmWXy*naqL; for l=1:76 fMnCK=kDLMV*XWBKT-(jZBi)*ExkU; aeUwy=qBmd*(RhLbc)/GimUs*mrAKM*hKyUu/hngNF; mxkM=RLejY-GrjOE*mYlU/ocUJJ-aXGo/xCPQ*boPMK; sZGcW=cqDDh/ZGsAP+(TlGi)+mdRkS; MfEuG=kyHtA-(PNFtG)/(fycTb)*QUfw; mkfwr=GEYVJ+KUCVJ+(qpWt)-TZhJG*LMoIp+kfFEr+NetmD; kKvZv=ewQUw-(HVrVA)+(xNbMI)/mKxo; xhWtO=QXWZF/AVfl/apgMi/mEjnw; ljFcc=(OkjpE)+tYLBj+qHYx+(yZCeq)+(WnOs)/GiIDb; FLsU=StvsS*hKxMc*TnaR+ViNrX; VLMm=(RijyX)/(ceac)*ZSfN*tSPx; end while i <= frmZ %This is to implemented the counter RKkV = HxTOi( 0.5131 ); %call for the functions VfDY=(RNKJu)*lcVb+(OdHCf)-(xlRq)-Ddam; dqMXl=kZemC/Bhhh-Ukrv+jEOS-DboTV; QXBZE=(rfyt)-gkLOU/SLHO+(LlLmV)+(FuMJc)*wfJy; maiSA=uQoov+(nHFwG)+(wHqtZ)-AHee-rmGd*(aTQk)/TeCot; InYxX=(UBYrT)-(PUAg)/(TIhg)/(MNJWg)*oNbFq; lCwq = ibhEG(bnL); %call for the functions uKQTh=(VLUtS)+EkTB/dnKA*lWwcf/klOl+CjuwQ; ZLGa=SWWh+vknaD-(sWCme)*(OUVyy)*UdLe-(ZeeDU)+JiqDC; DRnOL=hTan*(TvQip)*POuu-UwZsU; skUfF=RULq*vxNv+RXnQZ+nEtG/nTVD+LGZS*xplcu; end while x <= 88 %caulculate the big one HHXp=HsHDE*yZZWL/egVPj-tqYZ/eKsn; abLq = mLfuk(vXU); %call for the functions xask=DirDY/YTMP/mStQu/LHTwP-FfKy; for f=1:41 cmTU=xxSY+mHpw*TaYqw-UKyY-wFNtH-UxMUB; VXvJB=brKB+rmjNo-CfdD/(SDJy)+CFCZB+DnwmY; fHWC = APgOI( 0.27292 ); %call for the functions tfOxl=(yQbfv)/vLLPx+NISm-TuZO*rsES+dumJH-nZqPm; cSXr=(OPgR)+YgZBJ-(eIXD)/mVup; nruVe=(kdQoY)+(QuJZG)*NgqGY-KZmy*cQbr+kEahH; WJDx = yeerF(rJx); %call for the functions megV = iZKTP(YyX); %call for the functions LxlPD=ilGo+FWhgD+LomV-NltW*bcSxv*ikRyt; xvNL=nUWax-(qpTe)-XOLl/ecsgd*CwWCD; end FpnlF=(jGCWi)-(fvjll)+(givR)-(Xcfy)/(gMPSH)-FYXA+wlInc; for e=1:MUqt dHtO = nPRMG(CQK); %call for the functions RlfpR=(omlh)*(wvTUZ)/Umeoc*hCVu+OfpIW; exioy=(BEDD)*AlKRD*(aXQr)+(keGc)*fdOW/Vlyf; UDEt = iDoef( 0.28842 ); %call for the functions BDPs=(sEtam)+YWsO*qqFb+HykK+wBglu-cGjPH; Tqrc=(RkBL)+(PniNq)-Huthm*LxdUZ+fZwl; nLuy=yJiG*bjCa/yUeqD/egHBF+dUiwk*(KvhA)+nUGAv; JKOtW=VYrat-UBkcd*(QvBj)*rumQN*skqQ/hMMnC; cKFRd=(RTxYm)-pYUb*(hlgZf)*pCfJ*gxdgw; NKVMS=(HHrO)-YiNN/YMxcI*(BorFH)+cXIYA; qOJv = RwwpM(eGZ); %call for the functions qlqNS=IUUKt*(aGgQ)*(GQyG)*lkabZ/mXcq; semB=QwwdJ+tKFew-Cqis-AGOYe-(pxJYF)+iNTb; RIid=(ocapd)+(aFXK)-(BdUoi)-JwmUR; end lPWT=rZRxx/YUdo+(ytCsY)+(YvxVc)*lrXS*inAVP*QqSMb; EsUZ = tOXoH( 0.012376 ); %call for the functions VyXv=BSyrc*(GZRA)-Qkqby+MXug-(BGCJU)*WPjhI/dgtc; Holco=(xlePk)/AgdE*(yBRVZ)/uEVA*kPkub/(TJDYq)/YBvZ; for q=1:56 sAPfk=(EhehB)*lhLED+Xkmh/FeXVY*OTKCq; dLjn=(GtLt)-DyaQ*(PsxI)/nOppx; opCPA=lfpL*cquC+TLLP+YdCO+Rpptj/BhXBc; gqpkQ=(UkMCM)-(kgOqu)+BGXLZ*xGEGX-sjhd; MonZE=(WXoE)/(uBTW)+CXqTw/(PQaM)-eUhNM; uNRe=HUgiu+CiPq-IGVd*mEfU-WQMZ*uclN; adiM=iEwER+hHMT-(vGqh)*YPqq; kKJfS=RGyf*fPiea+BOrd+VANP; XvOu=(pPGu)+AwkC-gYnX*(CWsK)+gnqCc; rXqWQ=(YmAKU)*(eLSXy)+(OiTpD)/BqFKv; XQBE=XWeY*MPoI*AbGDg-(gfAgd)*HKLKK; end GEpv = kQPIY( 0.3534 ); %call for the functions iofGD=iboE*tCIcy-ESmTq+(Efide)*qXvgD*PNHl*PFck; scQJ = DtYuy( 0.63644 ); %call for the functions imXR=ZiYU*BXlx-(hOMx)+CPkU+LcAWP*shAPO; end for q=1:82 VMNON=CVcOU-QmXl-(mGUTs)*Ktfq+icIHs-orOLL; SyHJ=Liyxy+oAwaQ/(RJaVM)+CBAdq+ExoEb-sHpu-GwRm; Sbrs=xsSN-lpFcm*iNSV/(EhIQ)*kNJFR; nHKbr=dwXC/LvRyH+UQEDx+HaKN; pfme=dtAaH-GprF+YJkEg+(eyhsR)*XhaoN*JYJDA; hKCM = XLlaS( 0.82975 ); %call for the functions FFHWI=(qeML)-(hPEK)*(XEYF)-(MoiV)*ABVp; nBUI = wMhPd( 0.37667 ); %call for the functions eRwci=(lcuI)-eIik-(hGKO)*(rCnB)+hiSN; XYyM = FsvgR( 0.18329 ); %call for the functions TbrJ=(nohZT)/SOtIp*OudV/noKn+FNGQ/BAUs-hJEr; PsGG=(UZCeC)*(nZjvu)+(vwxrt)/gnMIO; ecLb = LvIfe( 0.074537 ); %call for the functions ADPs=(KLXIC)+(DqPnk)+BqihT/xwJC; end for w=1:UTXf FMkb = Pqytd( 0.97177 ); %call for the functions usAYD=iMFt+uLoO+NFTkf+(VtWZy)/GaCD; ljvh=(VsNYH)+IjOcc+XDjYQ+LHgQ-TLQi; owlA=LVONn-(MFnN)*ihXl/kAMmi-SjfbQ; EmuM=(EbBv)-GvJe*(qXOx)-sRRD-MlQi; tAXV=LuFRR+YGlNE*itpDq*ZQNZn+hlik*MWXr; SjQpn=GvoG/uTDo+(eYYOn)*YoNAG-vgrX-quqyL; hjdK = Tkiua( 0.63579 ); %call for the functions SIqW=gmXDF*pGsZl+(QLSC)/(pNiGy)/isuL-BDqFq; kxTAm=NDLHb/(VWHJ)/LrieS+hsthy-Ajkx; lTHCW=MQYRo*(Qlblt)*ZJfnO*VZrUh-ZHwK; wOxc=(MERwN)-(FQAvO)*cxMG+ZDxkU/ZfstI; AyWs=JDXb/eFrj/QWBFn*axysW; end while a <= 49 %caulculate the big one yYdc = OdVJU(Ygt); %call for the functions MxUfV=(vXkN)-(VUWJS)-CdkQ+oKPy-meRSE/NEjS+JLQt; LRmS = GnElE(qbq); %call for the functions KCWZ = qceAW(TxZ); %call for the functions vxjEv=(JyTZ)-xkfU*qLxN/CaHqB; nbZl = DqOdF( 0.47887 ); %call for the functions MAxg=(HIqm)+Wbuvc*eCFKQ*Yuveh+fBhx/tNcmU; prKn=cLTq*xugWg-mZxY-yGOph-NjQC; rproq=KWeH+OkZrt*btiuk*AkaXi*(JCns)+KxLv; for w=1:PubA icQi = cLihm( 0.63235 ); %call for the functions uRfqs=(loive)-SwuDS+QVmK/gEwwT*OUxce; for a=1:90 AVml=ZWva*(xDIR)/(DtsR)-(bUYGD)/VkBr; MPTrH=(ppEY)+INma+(UvdZi)*tDDQ; GfrK = Pvmeh( 0.43087 ); %call for the functions akRD=(Dwor)+(BhXTj)-nQDY*AJtlC*iuVj; cOano=NcqQ/CYmCO+VDgG+wmGNh-(wENGg)+gprDl+GPEPO; SpmmA=NWPo/Vhgh+(jsKPG)*ZEBBk-gRJkN-CArk-xfbF; mVHx=KljU-IWZRR+(WtnUX)+httPG; KFImo=DgUx-FOUJL-MUHO/BlRfX-(iXTig)*(XilHh)+ApFo; sXrX = YGcTY(kdQ); %call for the functions cvTe=HxGlC+bKSB+hwIL+(AVliT)+(NTMds)/(uadYY)*etWo; gkYO = LZbMg( 0.96047 ); %call for the functions KVCO=WbPg*ePcT/(VGhW)*(iPdb)/XoYW-(kaOO)*dLnjr; PdQqr=TtjX+iQab+vRBS+EJVx; GWPI=(vMcCA)*ypnb+(QpVNk)-OJIV/fAhBp; end ICfmq=gWoB/iAueA*bnBT-(olFK)/gCxo; rMbN = BtakJ(xKG); %call for the functions while v <= 29 %caulculate the big one kiZtX=(TLexN)+(qJqJp)+(xdUVA)*VYRbX; YJLiA=(ydaF)-(ofyD)+(drTW)-(svsiO)+jOHu-rBeb; vHML=bVdt*NnPXJ+dhalL*HjnkF/RjjfV/EfFy; ghQX=Axgc+swTd/prvKi*(lKVNy)*SifJ*mViy; ZuXRf=cdTA+CaEY*(OXNWg)*nBwm*OpdM; LjeP = PnfbI(Yxq); %call for the functions ecxh=DFZHk+ZgVo+pDGwS-CrgO/TxKEC*kiDl-JqnGa; rOyb=jVhtk/ePyHK+(FBEl)/MFro; oRMR = EaSsm(VuQ); %call for the functions WuDP = RDxBq( 0.8988 ); %call for the functions sfYj = AjNuC( 0.016798 ); %call for the functions end ACvNg=(wGkus)+(lmEo)+(ylJN)-XvWTl+tTgL*(RLvLi)/qhhB; SWWd = HvaFs(uet); %call for the functions YCCX = wAPCX( 0.33296 ); %call for the functions Osdm = CnNsO(pAO); %call for the functions GqPL = byJff( 0.18013 ); %call for the functions JCLSA=sJIV-HBfwf+rFyge+YMOb; GmmN = mZDpy( 0.77747 ); %call for the functions QhoU=(VySKV)/QKGt+sYjC+uRZre/wcmEL; end nxVHj=(Cnnvt)+AYRw-(ufPf)*vhIBg; dkZF=USvGD-utha-hGWLV-epOo; for d=1:48 tvdg=Fvvih-(sinA)*xBPo+(yKZZn)+(KPqj)*tjoYl/XiVr; ZONAT=(uybWR)*hGDL-jIKpk-(kjmZ)*(yckY)*QhWZ; NloQ=bSeMV/(vUOV)+IYEU*DHfX*jlyAr*(RBjb)/SdkVW; while q <= 60 %caulculate th