ideal=avsread('phtmesh.fld'); ii=ideal(:);
%new=avsread('fbplcmine3.fld');
%old=avsread('fbplchuber.fld');
new=avsread('fbphcmine3.fld');
old=avsread('fbphchuber.fld');
N=50; x=zeros(128*128,N); y=zeros(128*128,N);
zimg=zeros(128,128);

for i=1:N, tt=new(:,:,i); x(:,i)=tt(:); tt=old(:,:,i); y(:,i)=tt(:); end
a=std(x'); b=std(y'); z=a./b;
%zimg(:)=z'; image(zimg*128);
zz=z(find(ii>0));
hist(zz,30);
%print histhc.ps
