function   [cpaxyz,cpixy,pictureim]=getpicturedata(cpfname)   
  fid=fopen(cpfname,'r+t');   
  [cpnum,count]=fscanf(fid,'%4d',1);   
  cpnum   
  cpaxyz=zeros(cpnum,3);   
  for   i=1:cpnum   
        for   j=1:3   
              [cpaxyz(i,j),count]=fscanf(fid,'%e',1);   
        end   
  end   
  cpaxyz   
  [picturenum,count]=fscanf(fid,'%4d',1);   
  picturenum   
  for   k=1:picturenum   
        [picturename,count]=fscanf(fid,'%s',1);   
        picturename   
        pictureim(:,:,k)=getapicture(picturename);   
  %       figure;   
  %       imshow(pictureim(:,:,k));   
        for   i=1:cpnum     
        for   j=1:2   
          [cpixy(i,j,k),count]=fscanf(fid,'%e',1);   
  end     
          end   
        cpixy(:,:,k)   
  end   
  fclose(fid);   
  函數getapicture.m取得圖片的灰度矩陣   
  function   im=getapicture(bmpfile)   
  [x,map]=imread(bmpfile);   
  im=ind2gray(x,map);