clear;clc;
IM=imread（‘4040.JPG‘）;
IM=rgb2gray（IM）;
IM=double（IM）;
figure（1）
imshow（uint8（IM））
[maxXmaxY]=size（IM）;
IMM=cat（4IMIMIMIM）;
%%%%%%%%%%%%%%%%
cc1=8;
cc2=90;
cc3=180;
cc4=250;

ttFcm=0;
while（ttFcm<15）
????ttFcm=ttFcm+1
????
????c1=repmat（cc1maxXmaxY）;
????c2=repmat（cc2maxXmaxY）;
????c3=repmat（cc3maxXmaxY）;
????c4=repmat（cc4maxXmaxY）;
????c=cat（4c1c2c3c4）;
????
????ree=repmat（0.000001maxXmaxY）;
????ree1=cat（4reereereeree）;
????
????distance=IMM-c;
????distance=distance.*distance+ree1;
????
????daoShu=1./distance;
????
????daoShu2=daoShu（::1）+daoShu（::2）+daoShu（::3）+daoShu（::4）;
????distance1=distance（::1）.*daoShu2;
????u1=1./distance1;
????distance2=distance（::2）.*daoShu2;
????u2=1./distance2;
????distance3=di