function [result] = mat5write(matfile,imageC,imageI,xw,yw,zw);
% open the world coordinate files
fnamex = strcat(matfile,'X.byt');fnamey = strcat(matfile,'Y.byt');
fnamez = strcat(matfile,'Z.byt');fnameC = strcat(matfile,'C.bmp');
fnameI = strcat(matfile,'I.bmp');
%get the dimensions
[my,nx,pz] = size(imageI);
%reshape the 1D world coordinate vectors
x = reshape(xw',1,my*nx)';y = reshape(yw',1,nx*my)';z = reshape(zw',1,nx*my)';
%xw
fpx = fopen(fnamex,'wb');fwrite(fpx,x,'float32');fclose(fpx);
%yw
fpy = fopen(fnamey,'wb');fwrite(fpy,y,'float32');fclose(fpy);
%zw
fpz = fopen(fnamez,'wb');fwrite(fpz,z,'float32');fclose(fpz);
%C
imwrite(imageC,fnameC,'bmp');imwrite(imageI,fnameI,'bmp');
%
result = 1;