retrdata.m

改进等高线拟合代码

function retrdata(dummy)

% Part of GUI to retrieve data and draw fits and points
% The retrieved file must be in .mat form containing the variable 
% fit_data consisting of (in this order)
%  	fdatax; fdatay; 
%  	fcolor; fstyle; fmarker; 
%  	xp ;yp ; xsize; ysize
%
% If no arguments are given, a filename is asked for,
% with a dummy argument the file temps.mat is opened.

global FILENAME
global FITHANDLES
global PATHNAME
global FILECHANGED


if nargin ==0,
	seterror('');
	[filename, PATHNAME]= uigetfile('*.mat','Open Data File');
	len = length(filename);
   if (len < 4),      % probally cancel in uigetfile
   	 seterror(' ... Canceled');
		 return;
	end

	filename = lower(filename);
	PATHNAME = lower(PATHNAME);

	lf = length(filename);
	ok =0;
	if lf>4,
   	if strcmp(filename(lf-3:lf),'.mat'),
      	if lf < 13,
         	ok = 1;
	      else
   	      seterror(' Error: Filename to long.');
      	   return;
	      end
   	end
	end 
	if ~ok,  
	   for i =1:lf 
   	   if filename(i) == '.',
      	   seterror(' Error: Incorrect filename.');
         	return;
	      end
   	end
	   if lf<9,
	      filename = [filename,'.mat'];
	   else
	      seterror(' Error: Filename to long.');
	      return
	   end
	end
	if ~exist(filename),
		seterror([' Error: Can''t open file: ''',filename,'''. File doesn''t exist.']);
		return;
	end

	routine('new');
	seterror(' ... reading data.');
	eval(['load ' PATHNAME filename]);
	FILENAME = filename;
   FILECHANGED =0;	
else
	eval(['load temps.mat']);
end

fdatax = fit_data{1};
fdatay = fit_data{2};
fcolor = fit_data{3};
fstyle = fit_data{4}; 
fmarker = fit_data{5}; 
xp = fit_data{6};
yp = fit_data{7}; 
xsize = fit_data{8}; 
ysize = fit_data{9};

nrofpoints = length(xp);
set(FITHANDLES(1,12),'value',nrofpoints);
for i=1:nrofpoints,
   set(FITHANDLES(3,i),'value',xp(i),'string',num2str(xp(i)));
   set(FITHANDLES(4,i),'value',yp(i),'string',num2str(yp(i)));
end

hold on
nroffits = length(fdatax);
for i=1:nroffits,
	h(i) = plot(fdatax{i},fdatay{i});
   set(h(i),'linestyle',fstyle{i},'marker',fmarker{i},'color',fcolor(i,:));
end
hold off

ims=['imagesize',num2str(max(xsize,ysize))];
setting(ims);
setting('showtitle');
setting('drawaxes');

if nargin == 0,
	seterror([' Data read from ''', FILENAME ,''' .']);
end