Calviewplot.m
Jump to navigation
Jump to search
% NOTES:
% Enters new times, temp into qcaled window
%
function out=calviewplot(action)
switch action
case 0
%initialize figure
fig=figure;
out=fig;
zoom on % leave here, erases buttondownfcn if moved afterward
set(findall(0,'Tag','calviewBut'),'Userdata',fig);
%put the pointer into the go button
set(findall(0,'Tag','goBut'),'Userdata',...
findall(0,'Tag','calviewBut'));
set(0,'CurrentFig',fig);
set(fig,'pos',[130 30 600 520]);
set(fig,'Tag','calviewFig');
%initialize the plots
subplot(221); %calplot
set(gca,'Tag','calplot');
set(gca,'ButtonDownFcn','calviewplot(3)')
subplot(425); %tempplot
set(gca,'Tag','tempplot');
subplot(414); %specplot
set(gca,'Tag','specplot');
h=subplot(222); %dummy for table
get(gca,'Pos');
wavepos=[ans(1:2) ans(3)/3 ans(4)];
meanpos=[wavepos(1)+wavepos(3) wavepos(2) wavepos(3:4)];
stdpos=[meanpos(1)+wavepos(3) wavepos(2) wavepos(3:4)];
table(fig,wavepos,meanpos,stdpos);
delete(h);
%
rbutpos=[meanpos(1) wavepos(2) wavepos(3) wavepos(4)/10];
rbut(fig,rbutpos)
subplot(426); %stdplot
set(gca,'Tag','stdplot');
signature('calview.m');
%if this window gets closed, calviewBut should be cleared
%set(fig,'Close','disp(''Use Go for next file, or Quit to close all'')');
case 1
[ax_han,fig]=get_ax;
axes(ax_han(1));
udata=get(fig,'Userdata'); %global for all files
prop='a';
if udata.par.prop,
prop='c';
end
ylabel([prop '-' prop '_w (m^{-1})']);
xlabel('time (min)');
axes(ax_han(2));
ylabel('instr. temp. (C)');
axes(ax_han(3));
ylabel([prop '-' prop '_w (m^{-1})']);
xlabel('\lambda (nm)');
axes(ax_han(5));
ylabel('Std. Dev.');
xlabel('\lambda (nm)');
case 2
ax_han=get_ax;
udata=get(ax_han(1),'Userdata');
%initialize calplot
calviewplot1(udata);
axis auto
t_lim=udata.t_lim;
calviewmask(t_lim)
zoom reset
calviewplot2(udata);
%calculate the stats
calviewplot(6)
case 3
%calplot callback
ax_han=get_ax;
udata=get(ax_han(1),'Userdata');
newdata=get(ax_han(3),'Userdata');
%calculate the stats
calviewstats(udata,newdata);
%plot new spectra
newdata=get(ax_han(3),'Userdata');
calviewplot3(udata,newdata);
%update the stats table
newdata=get(ax_han(3),'Userdata');
calviewplot4(udata,newdata);
%update stdplot
calviewplot5(udata,newdata);
case 4
%AllChan callback
ax_han=get_ax;
axes(ax_han(1));
udata=get(gca,'Userdata');
t_lim=get(gca,'Xlim');
zoom out
calviewplot1(udata);
calviewmask(t_lim);
case 5
%AllTime callback
%calviewplot1;
calviewmask([0 0]);
calviewplot(3)
case 6
[ax_han,fig]=get_ax;
%ResetStats callback
l_han=findobj(ax_han(5),'Type','line')
delete(l_han)
set(ax_han(5),'userdata',[])
l_han=findobj(ax_han(3),'Type','line')
delete(l_han)
set(ax_han(3),'userdata',[])
calviewplot(3)
case 7
%ResetTime callback
ax_han=get_ax;
udata=get(ax_han(1),'Userdata');
t_lim=udata.t_lim;
calviewmask(t_lim)
calviewplot(3)
case 8
[ax_han,fig]=get_ax()
axes(ax_han(1));
ax=axis
minmin=ax(1);
maxmin=ax(2);
t1=sprintf('%5.2f',minmin)
t2=sprintf('%5.2f',maxmin)
set(findobj('Tag','t1Box'),'String',t1);
set(findobj('Tag','t2Box'),'String',t2);
mtemp=sprintf('%6.1f',get(ax_han(2),'Userdata'))
set(findobj('Tag','InttempBox'),'String',mtemp);
qcaled_cb('update');
return
end
return
function [ax_han,fig]=get_ax()
fig=get(findall(0,'Tag','calviewBut'),'Userdata');
ax_han(1)=findobj(fig,'Tag','calplot');
ax_han(2)=findobj(fig,'Tag','tempplot');
ax_han(3)=findobj(fig,'Tag','specplot');
ax_han(5)=findobj(fig,'Tag','stdplot');
return
function calviewstats(udata,newdata)
ax_han=get_ax;
axes(ax_han(1));
ax=axis
minmin=ax(1);
maxmin=ax(2);
%new values from calplot axes
t1=sprintf('%5.2f',minmin)
t2=sprintf('%5.2f',maxmin)
mtemp=sprintf('%6.1f',get(ax_han(2),'Userdata'))
%pass array to qcaled_cb to update text fields
isishandle = ishandle(findobj('Tag','QcaledFig1'))
if ishandle(findobj('Tag','QcaledFig1')),
%file text boxes with current file parameters
qcaled_cb('gettextfields');
%update with new values from calplot axes
findobj('Tag','t1Box')
set(findobj('Tag','t1Box'),'String',t1);
set(findobj('Tag','t2Box'),'String',t2);
set(findobj('Tag','InttempBox'),'String',mtemp);
qcaled_cb('update');
else
display('qcal.txt file has not been updated, it is up to you');
end
% $$$ udata=get(ax_han(1),'Userdata');
par=udata.par;
time2=udata.time2;
chan2=udata.chan2;
itemp2=udata.itemp2;
I=udata.I;
ichannels=udata.ichannels;
range=find(time2>minmin & time2<maxmin);
% $$$ newdata=get(ax_han(3),'Userdata');
newdata.avespec=mean(chan2(range,I(ichannels)));
newdata.npts=length(range)*ones(1,9);
newdata.stdspec=std(chan2(range,I(ichannels)));
set(ax_han(3),'Userdata',newdata);
return
function calviewmask(t_lim)
ax_han=get_ax;
axes(ax_han(1));
h=findobj(gca,'Type','line');
%check for t_lim and adjust
if exist('t_lim') & ~isempty(t_lim),
minmin=t_lim(1);
maxmin=t_lim(2);
if maxmin==0,
t_data=get(h(1),'xdata');
maxmin=ceil(t_data(end));
end
set(ax_han(1),'XLim',[minmin maxmin])
end
return
function [ax]=calviewplot1(udata)
%this makes the initial plot of the data using times from file
ax_han=get_ax;
axes(ax_han(1));
h_old=findobj(gca,'Type','line');
if ~isempty(h_old),
delete(h_old)
end
% $$$ udata=get(gca,'Userdata');
title(strrep(udata.file,'_','\_'));
%change plot to line to avoid reinitializing
h=line(udata.time2,udata.chan2(:,udata.I(udata.ichannels)));
mylines(h,'c');
% $$$ %check for t_lim and adjust
% $$$ if exist('t_lim') & ~isempty(t_lim),
% $$$ minmin=t_lim(1);
% $$$ maxmin=t_lim(2);
% $$$ if maxmin==0,
% $$$ t_data=get(h(1),'xdata');
% $$$ maxmin=ceil(t_data(end));
% $$$ end
% $$$ set(ax_han(1),'XLim',[minmin maxmin])
% $$$ end
%zoom reset
l_han=legend(udata.channels(udata.ichannels,:),-1);
l_pos=get(l_han,'pos');
a_pos=get(ax_han(1),'pos');
set(l_han,'pos',[l_pos(1) a_pos(2)+a_pos(4)-l_pos(4) l_pos(3: ...
4)]);
return
function calviewplot2(udata)
%returns max. int. temp, only called once per file
%change plot to line to avoid reinitializing
ax_han=get_ax;
% $$$ udata=get(ax_han(1),'Userdata');
axes(ax_han(2));
%remove previous lines first
h=get(gca,'chil');
delete(h);
line(udata.time2,udata.itemp2);
mtemp=max(get(findobj(gca,'Type','line'),'ydata'));
%store mtemp in axis rather than returning it
set(gca,'Userdata',mtemp);
return
function calviewplot3(udata,newdata)
%plot/replot spec data and store in axis
ax_han=get_ax
% $$$ udata=get(ax_han(1),'Userdata')
axes(ax_han(3));
% $$$ newdata=get(ax_han(3),'Userdata');
l_han=errorbar(udata.par.wave,newdata.avespec,newdata.stdspec);
hold on
if isfield(newdata,'specdata'),
plot(udata.par.wave,newdata.specdata,'k')
else
newdata.specdata=[];
end
hold off
newdata.specdata=[newdata.specdata; newdata.avespec];
set(gca,'Tag','specplot');
set(gca,'Userdata',newdata);
return
function calviewplot4(udata,newdata)
ax_han=get_ax;
% $$$ udata=get(ax_han(1),'Userdata');
% $$$ newdata=get(ax_han(3),'Userdata');
strvcat(sprintf('%3.0f\n',udata.par.wave));
set(findobj('Tag','WaveBox'),'String',ans);
sprintf('%7.4f\n',newdata.avespec);
set(findobj('Tag','MeanBox'),'String',ans);
sprintf('%7.4f\n',newdata.stdspec);
set(findobj('Tag','StdBox'),'String',ans);
return
function calviewplot5(udata,newdata)
%plot/replot std of calplot spectra
ax_han=get_ax;
% $$$ udata=get(ax_han(1),'Userdata');
axes(ax_han(5));
%stddata=get(gca,'Userdata');
% $$$ newdata=get(ax_han(3),'Userdata');
hold on
plot(udata.par.wave,newdata.stdspec);
if isfield(newdata,'stddata'),
plot(udata.par.wave,newdata.stddata,'k')
else
%first time through, initialize these fields
newdata.stddata=[];
newdata.ndata=[];
%problem is that if the plot rescales, the markers don't move
%with the line always, huh?
for ip=1:9,
h(ip)=line(udata.par.wave(ip),newdata.stdspec(ip));
end
mylines(h,'mc');
end
hold off
newdata.stddata=[newdata.stddata; newdata.stdspec];
newdata.ndata=[newdata.ndata; newdata.npts];
title(['npts = ' num2str(mean(newdata.npts))]);
set(ax_han(5),'Tag','stdplot');
set(ax_han(3),'Userdata',newdata);
return
function []=rbut(fig,rbutpos)
uicontrol('Parent',fig, ...
'Callback','calviewplot(2)',...
'Units','norm', ...
'BackgroundColor',[.8 .8 .8], ...
'Position',rbutpos, ...
'String','Reset time', ...
'Style','Push', ...
'Tag','ResetTimeBut');
uicontrol('Parent',fig, ...
'Callback','calviewplot(4)', ...
'Units','norm', ...
'BackgroundColor',[.8 .8 .8], ...
'Position',rbutpos, ...
'String','All chan', ...
'Style','Push', ...
'Tag','AllChanBut');
uicontrol('Parent',fig, ...
'Callback','calviewplot(5)',...
'Units','norm', ...
'BackgroundColor',[.8 .8 .8], ...
'Position',[rbutpos(1) rbutpos(2)+1.1*rbutpos(4) rbutpos(3:4)], ...
'String','All time', ...
'Style','Push', ...
'Tag','AllTimeBut');
uicontrol('Parent',fig, ...
'Callback','calviewplot(6)', ...
'Units','norm', ...
'BackgroundColor',[.8 .8 .8], ...
'Position',[rbutpos(1)+1.1*rbutpos(3) rbutpos(2) rbutpos(3:4)], ...
'String','Reset stats', ...
'Style','Push', ...
'Tag','ResetStatsBut');
uicontrol('Parent',fig, ...
'Callback','calviewplot(7)', ...
'Units','norm', ...
'BackgroundColor',[.8 .8 .8], ...
'Position',[rbutpos(1)+1.1*rbutpos(3) rbutpos(2)+1.1*rbutpos(4)...
rbutpos(3:4)], ...
'String','Reset time', ...
'Style','Push', ...
'Tag','ResetTimeBut');
return
function []=table(fig,wavepos,meanpos,stdpos)
uicontrol('Parent',fig, ...
'Units','norm', ...
'BackgroundColor',[.8 .8 .8], ...
'Position',wavepos, ...
'hor','right', ...
'String','', ...
'Style','text', ...
'Tag','WaveBox');
uicontrol('Parent',fig, ...
'Units','norm', ...
'BackgroundColor',[.8 .8 .8], ...
'Position',meanpos, ...
'hor','right', ...
'String','', ...
'Style','text', ...
'Tag','MeanBox');
uicontrol('Parent',fig, ...
'Units','norm', ...
'BackgroundColor',[.8 .8 .8], ...
'Position',stdpos, ...
'hor','right', ...
'String','', ...
'Style','text', ...
'Tag','StdBox');
return
% TODO
% Add colored symbols to std and ts plots to match the colors on calplot
% Make buttondown on legend names plot that/those channels for all
% time with the greatest resolution.
% Unclick names to get all the lines back
%should newdata be saved in the figure udata rather than the
%axes. Make the figure udata the place to pass data to other figures
