%   read the calibration data
%   function needs the input calibration grid data file '.byt'
%   Npoints number of calibration points 
function [coods] = read_calib_data(infname,Npoints)
%
fp = fopen(infname,'rb');
calib_data = fread(fp,'float32');
fclose(fp);
%read the coordinate information
[m,n] = size(calib_data);
coods = double(calib_data(2:m));
coods=coods';
coods = reshape(coods,7,Npoints)';
%coods = reshape(coods,8,Npoints)';

%