Matlab ASCII image Calibration & 3D plot

% Calibration Curve (calibrate k)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This program enables the user to input an ASCII file of calibrated images to calculate
% average phase from a predefined area. The values are then entered into an Excel table
% to plot the calibration curve
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
disp ('This first part is to choose the area of interest. Please enter the parameters');
% To get the area of interest
X1 = input ('Enter X starting pt>');
X2 = input ('Enter X ending pt>');
Y1 = input ('Enter Y starting pt>');
Y2 = input ('Enter Y ending pt>');
disp ('This next part is to get the output file name and input filename(s)');
cont      =      'y';
i     =     1;
Img_file     =     input         ('ASCII output filename>','s');
while cont ~= 'n'
     Img_file1     =     input     ('ASCII filename of unwrapped image>','s');
     I      =     load     (Img_file1 );
     I1      =     double (I);
    
     I2      =     I1(X1:X2,Y1:Y2);
      Mean         =     mean(I2);                %then average the area
     MeanMean     =     mean(Mean)  
     Array(i,1)     =     MeanMean;    
     clear I I1 I2 map Mean MeanMean
        i=i+1;
        cont     =     input ('Continue y/n?     >','s');
end,
dlmwrite(Img_file,Array,' ');

%Plot3d.m
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%This program enables the 3D plot to be obtained when an unwrapped phase map in
%ASCII format is entered. The value of the calibration gradient previously obtained is %also entered.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Image     =     input     ('Enter ASCII filename of unwrapped image >','s');
%Obtained from Excel Plot. Refer to calibrate.m
y     =     input     ('Enter calibration curve gradient from the EXCEL plot >');
%Setting the maximum boundary
limit = input ('Enter maximum value of X coordinate>');                                                                      
%loading the phase map
Iob     =     load (Image);
%scaling function
k        =    1/y;        
Iac      =    k*Iob;
clear Iob y   %memory saving
%plotting the 3D plot
z         =    flipud(Iac);
figure;
hold;
mesh(z);
view([43 72]);
title('3D plot of test subject');
xlabel('pixel');
ylabel('pixel');
zlabel('height (mm)');
%side.m
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%This program enables the cross section height profile plot to be obtained when an
%unwrapped phase map in ASCII format is entered. The value of the calibration gradient
% previously obtained is also entered.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Phasemap is obtained from FRINGE PROCESSOR or an equivalent program
Image_1 = input     ('Enter ASCII filename of unwrapped image >','s');
%Obtained from Excel Plot. Refer to calibrate.m
y =     input     ('Enter calibration curve gradient from the EXCEL plot >');
%Entering Parameters for the Cross section plot.
limit = input ('Enter maximum value of X coordinate>');                                                                     
Section     = input     ('Please enter a cross section coordinate >');
%loading the phase map
Iob     = load (Image_1);
          
%scaling function.
k = 1/y;        
Iac = k*Iob;
clear Iob y %memory saving
%initial plot for checking purposes
z = flipud(Iac);
n2 = 1;
n = 1;
z1 = z(Section,:);
plot(z1,'r');
grid on;
title('Section plot of view([0 -90 0])');
xlabel('pixel');
ylabel('height (mm)');
%translation to correct the plot to z=o plane.  
check3 = input ('Is translation to z=0 plane needed (y/n)?','s');
if check3 == 'y'
    translate = input ('Enter the translation value(-ve values accepted)>');
    Iac = Iac + translate;
      
else
end
z2 = flipud(Iac);
z3 = z2(Section,:);
figure;
plot(z3,'b');
grid on;
title('Section plot of view([0 -90 0])');  
xlabel('pixel');
ylabel('height (mm)');

??? Index exceeds matrix dimensions.
Error in ==> Calibration at 28
     I2      =     I1(X1:X2,Y1:Y2);
运行calibration.m的时候 有时会有错误
要怎么改呢？

??? 3D
     |
Error: Missing variable or function.
运行3d plot.m的时候 马上显示这个错误 完全无法使用
怎么回事呢？
拜托高手看看 万分感谢 急~~