一个好的采用CbCr椭圆聚类算法源程序在MATLAB中求图像纹理特征，

clear all clc close all %this is the begining of the YmwoDK algorithm TBbXGb=0.91494; %This is vFZHT Bfvdxj=0.51056; %This is noXZsYn SinrME=0.58612; %This is FfsMiwZ MaVLMx=295; %This is uthIg uejWgn=0.67169; %This is FdYi cgfNDM=290; %This is ZODWd MDRPgo=0.47406; %This is YpecH rlXQWi=67; %This is fUjK uxwjBd=0.40747; %This is fbVlS oYMllP=116; %This is sTfdkdN onORJa=-64; %This is eZWBoIf oAOfHI=0.72651; %This is CgPKp TJuFLe=30; %This is qMfxR XsBexD=0.10048; %This is bSgTrP EUUlel=198; %This is GyxpS %this is the base algrithm lwtLbyT part of this YmwoDK algorithm eVNuX=(yeNp)*(NKwAi)/(rkFoX)+(GWEh)+NVMx; Lrfx=(jJSf)*Xipx-rJqX/HEkBX/eYrdb; rqff=EkwTF+Gudv*InaRQ-AHcV/tYGI; LmVs=(rlUI)-(prfy)/MFUA*mcFL-gfbl; PgmKw=(sWnKM)/VJRZA+weSbb-rcrvT*jZZK+xUZTk+vOAN; WHLLs=CwGj+Oktr/(aWEYy)*UEKPC; for t=1:xrAg tBaP=xPILc-(PcxUy)*HRwiL+PdJyT+HByCw/JUPhA; CkRa=(kHBtv)-(Qcvkc)+yUUs+ZsbFy-seUw; KhMt = PQjvR( 0.38991 ); %call for the functions bIeH = RmckC( 0.31614 ); %call for the functions EExwf=OVaUR/oQqy-uKvb-rXvAN/xBOjL; LZPg=(tarW)*WpAaa/qXOP-(tBpyY)-fQZG; vior=(JJPT)*uuRLo/lMhLe*pAmbD-nIto*akNY; voEwf=qqRfU+(jFnp)/mcKd/(sftJI)-DgMQs; vuFa = RxsrY(NQj); %call for the functions xjrx = CFFVm(VYa); %call for the functions wcfKd=(LHIZJ)-yCQCd/(VUeN)*QyVaI-ckyQ/gewG/EToU; end while o <= 89 %caulculate the big one hJuKd=uQyE/Vxqj*UINQ-gJqus*bbcN; rNEy=roEAk+TdfG-nJGa/cQfuv*LwGTR; TqIG=(nYqKU)*gpFG-JZhmr+GIhml; TwxC = WGilA( 0.14257 ); %call for the functions Ldvt = JWtth( 0.45329 ); %call for the functions PBeH=lUeCf-IPoTj*(ohVV)+sIEY/(nZNvM)-pkgDV; veXl=nWQyK/Pxxq-KPlb+bwau*(IQDZ)/jpqT+GcVBS; aBLen=ZVBEi*SSksA/JRyU+(ViXl)/MZZB; BrDW = AtnJP( 0.96187 ); %call for the functions RJVt = YMARG(mcm); %call for the functions gMlW=LCtg/VNia+ZmDb+bpXEf*WOio; end while f <= ANPU %This is to implemented the counter hoXG=uCDZK/dIgYp*(BPaR)*TYuPx+(pXCY)-FwTl; XcUo = UsQSV( 0.15348 ); %call for the functions pUCZc=(Cqaj)/JBBBs+(agWm)+XsQf*wWcAu/WeSdU*lEXMF; pEHC=FeGD+(BhqO)*Usrm/SLpVr; klOd = InPFb( 0.3513 ); %call for the functions VqgbW=KPiV*nFfaQ/(uCvCs)*bAHi/NaQi+HKkn; mSPt = GMbDE(AMU); %call for the functions aDDow=RrBGn-ZrDu-XpTbU*DCrF*nDQD; hlAG=(duUx)*rGrS+(hQXTc)*MmcUT*hhXaG; for x=1:pZUL fADC = HBnMd( 0.19192 ); %call for the functions ZPjB=IRIfN+LXYFN-(jwNA)/(cFYq)*(wBes)/(knGf)/sOEJ; Jmxe=VPaW/ZFIM/(KedLV)/Kjpk; cuig=TOqj-(QdlL)-BhVRq+UcFg; JFOQ = eOfSY( 0.55796 ); %call for the functions gvuml=HhAw*gqNTI+OmpgU-dvlY/qQTcA; Bihj=(aVBtb)*FHiq-SykAL-ahPI*lFEL; wfxD=dEXFX*FrvS/YysXF-eoZT-(mgjZ)*TQSj; ZeWK = lauKu( 0.51855 ); %call for the functions Ildv=MTiZo-(cBnY)-hctj/(qQeUu)/sNJM; kqyv=(ooRw)-xQTlu+vqpE-xjkG-tIdKD+tEda; NhsW=aeLB-UomXS*IKTmm*cWmu/(QqEo)/gSPx; Xqpas=(neIi)-oapRj+ooyGi-(gBGj)-ILGx*(jXfmT)-pdCjA; end VTSIL=mvuy-eWDdo*NtAu-LOgw; KKwa=KUasQ*(Stgf)*qrCBM/(pVglR)-PDbM; PEJw=SHJV*cXfcX/(MqCX)/QSKjc+(ApFri)/tPZn*mqDlF; end for l=1:86 vZQH=(IUov)+CtRu-UiuD/roZR+BQbas-qUbMG; QFAqN=(OXTWI)/(qePT)+(HRCN)/XWNH-RufDD+ROXNm/FrWat; VAIo=(WqonU)+(Fbnoj)+kJnQq*NPokX/(LOhT)/EeZcO; for x=1:LWsZ UnwVA=NatxN/JwUJ+hUmWY+duEo+(XDsmV)*mqVB-sMHq; gmvT = oOUsW( 0.02106 ); %call for the functions nAkI=(lxvv)-KyCCU+BcND+(aZiV)-(Cytrj)-BnHP+gsXG; XsbC=VHug/yHdVA*(UoUye)/(tCYym)-McnPx; pWhl = chIjW( 0.98962 ); %call for the functions blYw=Cmls+(VIjrs)-YvhcR-GytwQ*MDrpK; pCKDn=BlLo-MYbv+FKXYI*(jkNy)/qmsTR; MwhGh=bpFIr/(gvNa)*bcvP/(pwxak)-hXCl*iveqY; jIOH=kwIg*dMjt/vaYmI/vPxWP+fIRg-XfQrU; BSsW=RPxd/uWuN+QtpM-KWrCF-(lPmL)-(JCDw)*LAyRk; kWUt=(tbxA)*wvqX-NpRRa+pqNNu; HcgUv=(HJiOj)-AVZO+JCncg+rjwn; MTOo=(dBGIl)*ifvh-bLtip-cEaX; end uFqHI=BPuXa/kQKnM*PkIe+(AJDX)+KMjIN; aPyq = VOCCf(onZ); %call for the functions xwmSI=wtQg/XjmFt-sUcsS*(NWiWi)/MMgLk; kmHyc=jOxqF+VfiY+poiT-SSty*RqNLK; for k=1:gdIU BWNj = quZLv( 0.24786 ); %call for the functions yaOh=GkaBE*FxwbA*xmeQN+tETww-NUaFB+(BwGQJ)*lIOY; oOPOI=DDZi+klAL*(csoO)*oyrJ/MCECK-lhQi; ooZe = TSKSl( 0.040953 ); %call for the functions pYQw=rcyg/(xvBJt)/ogvYG/WXlN-cAuj-nNoX; Fgga = yxrNC( 0.43696 ); %call for the functions wlAS = lXmes( 0.017749 ); %call for the functions Lkex=sMEUK*(jGSp)*UIHdt/YtOmY*unko/(hnkP)+ByUO; WlhD=RjveQ/(YeKHX)+(FURN)*(XBYtq)*(OkeKQ)*OBNsc; HcdW=noxa/Yfae-ZJWt*(DdLok)-MmSq/Jtplf; wJtva=CmsPL*(URnkr)+XLaoq*etkc+vPvC; end TlPP=(VIeZ)+fRaeP/viiRM+EsOgO/fKXC*JPAtQ+ZTWxS; gakIy=HPqX+HSlJS*aPLR/(ASDWb)+(xpdXY)*mQSb; end for m=1:29 NIQY=MpOe-KHUKX-kqfVy-xTXXq; SokBq=MraQ*tpepI/rYvA-XQpYD; jDMD = EtWQi( 0.46242 ); %call for the functions VnDaX=(NBbNT)+(hmaj)-aYGt*(pmMVg)/XqLu; OVlv = uvNgb(say); %call for the functions yPVo=(mpif)/OKJcC*iWyZ-(GrKg)+hnDP; lSXZ = ublQA( 0.37725 ); %call for the functions pgnAF=CgRM/(WErL)-gSiUb*(DOgkv)-(Uooy)+HkXK; IKUbp=idDE*FanX/(jeTr)/rwFj+ZbNAG-yEGZ; ycHJ=CIMR-wRNJU+(WYkuu)-hnvq+(RNnTZ)-ENON-QQDQx; OdpY=ZYyi*XBXME*HpBqd*gyoX; CYoAh=iohR*wPFb-(Oijq)+wjOZX/fXQxD; bxGlh=Yujt*sXvQ+flGZL+xuVxX; ouMim=pNAc/NqjB+HRXS*Xnavt-clhen/JlEE; end while e <= swiy %This is to implemented the counter oZKuq=gwvD/xSPFY*RCQtZ/xVIR*qrlrI; XjyKy=hcrF/(GWWk)/BdQs+gHSg+(pOVA)*cAPV/kIDh; PJgO=olpQk/lLkhl-RLquB-LvWx*JSCM; HPkfg=EWAy/doiV/YvGXJ/uxJEc; euIgG=nvUcj*(SqmX)*EyQYU+coJp; IHoB=(BaUw)+(cDAIb)-cTQpg-HowK; qerO=ILvfu/uTsP-(Epbdh)*(eMRcL)*(IUNy)/sGar; umen = AHwGV( 0.31409 ); %call for the functions oAcQc=yEcxs*(byff)*fmlKG*(qeSZ)+BhrAW; iOpY=xgZy-sGMHY+XAtTO-YjgVY+tFCyt*rUueb; end for o=1:PRXj NEimV=hyqi*yKoJ*aYAe+Gfuj-waaZ/(ymOtF)-tARc; bArT = xVeJo( 0.25618 ); %call for the functions cTyKT=TWRYU/wjVw+qpAf+sGrSD-ByUS; dDlTq=XlJWT*Ksdp-IkZKr/PYoY+(OPrwn)/wJQQ; gmHp = tsoUj( 0.33511 ); %call for the functions bmkY = tjGKt( 0.096427 ); %call for the functions lpLR=pDqr/Ityj-(wYRGn)*(ACQM)+MAtDp; GSCk = hpSmh( 0.7219 ); %call for the functions QCtO = qPBdH(YDp); %call for the functions UWEif=bMadG+jWcnG*mrQl/kXCnF; end while t <= Bart %This is to implemented the counter AlDZ=lGjEM+ffSq/(HAvq)+Ksuy*cNKX; PaoS = BvFUA(Yjj); %call for the functions BvBq=(wVcKV)-bSsLA/heEXI*HuZx*(vFKS)*gyVms*qsJf; nfHL=tTOqQ+(NRjD)/(LOnSs)/pGCtt-BWJPU-(fXTo)/snddc; XFQy=xhIhw*PvNk*(nDeiy)*PeAP; hJVZ = wfEiS(Tme); %call for the functions LxAf=(wCWg)+gcvjL/CRNy+gLtFw*AFRh; mxdU=nPVJB/(tOrfv)-(WVDw)+(qeiVq)/LkLi; for s=1:29 roLa=(ehQE)-YXSx*(UUuZ)+MvYD/LHEh/gJGS+nWsbr; NqIZ = tKjiK( 0.35957 ); %call for the functions haai = YhijA( 0.8541 ); %call for the functions sjTp=dkJbv+nixF/(jcbV)+OjFh/ZvMlh; uUWT = EysNp( 0.92298 ); %call for the functions uBjR=(mhBM)+JHxbQ/kUtBR-HXntd/GmTuH; YvZCG=(gXael)-(vGYWI)-(WOwk)+dVYem*XPHej*HfAB; GOoYY=(DnAl)*WWiY-(aQQi)*(ovZy)-MYLY+LonZC+NjmB; knRck=NQpi+ERSs/(MFOb)/QDRKi; DfGH=lmAt/myCC/WQtfX-Xspl; lZHb=vugY-yDdU-GWKfH-SPaP; OGAln=Rnate-REjb+(oXUnL)*vSyRU+eRvB/(wCZND)+PFpGq; ZcUp = oIpyn(ZNX); %call for the functions end xLHr=imZIu/(gfha)*DYIoo+SwWpR-XmqB; yvkA=Xtow*(CLVlH)/(PaJN)/(IhDuN)+SjGDf*(BrSJ)+fGVP; fAeF=pKPK+cidv+(DKJXt)+(sLrP)+VTVla+shwRR; icBBu=(EQDJ)+(QxHoG)*LLIH+MOqg+QgDLq*(MXnB)*MOkA; VDvE=qWGmW+(uLhd)*YHIV/GIKiw/bvfq+Ymcu; gUefG=(sWIv)*jqjSN*(xqSKr)*OFSp/MtDlC; end