www.pudn.com > TimeFrequency.rar > LogtfCBack.m
function LogtfCBack(action)
persistent fullfile;
persistent line;
persistent depth;
persistent val;
persistent mindep;
persistent maxdep;
persistent wlog;
persistent wavefun;
persistent scales;
persistent func;
persistent scalefun;
persistent cjcode;
persistent cjunit;
switch (action)
case 'Init' % Figure initialization.
func={'db3' 'sym2' 'sym3' 'morl'};
scalefun={'1:32' '1:64' '1:128' '16:32' '1:4:32' '1:4:64' '100:1:400'};
cjcode={'AC/DT' 'CNL' 'DEN' 'GR' 'RT' 'SP'};
cjunit={'us/m' 'us/ft' '%' 'g/cm3' 'API' '¦¸M' 'mV'};
case 'browes'
[filename, pathname] = uigetfile('*.*','Well Log File');
if filename ~= 0
fullfile = strcat(pathname, filename);
hnd = findobj('Tag', 'logdata');
set(hnd,'String',fullfile);
end;
case 'computer'
hnd=findobj('Tag', 'stdep');
mindep=str2num(get(hnd,'String'));
hnd=findobj('Tag', 'eddep');
maxdep=str2num(get(hnd,'String'));
[minval,stindex]=min(abs(depth-mindep));
[minval,edindex]=min(abs(depth-maxdep));
if mod((edindex-stindex),2)==0
edindex=edindex-1
end;
scale=[700,800,900,1000];
hnd=findobj('Tag','wavefun');
tmp=get(hnd,'Value');
wavefun=char(func(tmp));
n=str2num(get(findobj('Tag','wavedecn'),'String'));
if isempty(n)==1
n=wmaxlev(val(stindex:edindex)',wavefun);
end;
new=cwt(val(stindex:edindex)',scale,wavefun);
[tf,L]=wavedec(val(stindex:edindex)',n,wavefun);
[det1,det2,det3,det4]=detcoef(tf,L,[6,7,8,9]);
wlog.curves=horzcat(depth(stindex:edindex),val(stindex:edindex),new');
hnd=findobj('Tag','wellname');
wlog.name=get(hnd,'String');
wlog.curve_info(1,1)={'depth'};
wlog.curve_info(1,2)={'m'};
wlog.curve_info(1,3)={'test depth'};
hnd=findobj('Tag','curvename');
cjnum=get(hnd,'Value');
wlog.curve_info(2,1)=cjcode(cjnum);
hnd=findobj('Tag','unit');
tmp=get(hnd,'value');
wlog.curve_info(2,2)=cjunit(cjnum+tmp);
wlog.curve_info(2,3)={'test log'};
for i=1:length(scale)
wlog.curve_info(i+2,1)={num2str(scale(i),'Wavelet scale %2d')};
wlog.curve_info(i+2,2)={''};
wlog.curve_info(i+2,3)={'test TF'};
end;
wlog.first=depth(stindex);
wlog.last=depth(edindex);
aux=l_plot(wlog);
hold on;
h2 = subplot(1,6,2);
pos=get(h2,'Position');
set(h2,'FontSize',8,'Position',[pos(1),pos(2)-0.04,pos(3),0.75]);
% axis('tight');
plot( ...
det1,...
[1:length(det1)],'-' );
set(h2,'XAxisLocation','top');
xlabel('det1');
h3 = subplot(1,6,3);
pos=get(h3,'Position');
set(h3,'FontSize',8,'Position',[pos(1),pos(2)-0.04,pos(3),0.75]);
% axis('tight');
plot( ...
det2,...
[1:length(det2)],'-' );
set(h3,'XAxisLocation','top');
xlabel('det2');
h4 = subplot(1,6,4);
pos=get(h4,'Position');
set(h4,'FontSize',8,'Position',[pos(1),pos(2)-0.04,pos(3),0.75]);
% axis('tight');
plot( ...
det3,...
[1:length(det3)],'-' );
set(h4,'XAxisLocation','top');
xlabel('det3');
% set(axAct,'Ylim',[ymin ymax],'Tag',tag_axe_dec);
h5 = subplot(1,6,5);
pos=get(h5,'Position');
set(h5,'FontSize',8,'Position',[pos(1),pos(2)-0.04,pos(3),0.75]);
% axis('tight');
plot( ...
det4,...
[1:length(det4)],'-' );
set(h5,'XAxisLocation','top');
xlabel('det4');
h = subplot(1,6,6);
pos=get(h,'Position');
set(h,'FontSize',8,'Position',[pos(1),pos(2)-0.04,pos(3),0.75]);
axis('tight');
hnd=findobj('Tag','scale');
tmp=get(hnd,'Value');
scale=str2num(char(scalefun(tmp)));
tf=cwtnew(val(stindex:edindex)',scale,wavefun,'plot');
set(h,'XAxisLocation','top');
% specgm = Spectrogram(val(stindex:edindex));
case 'readdep'
hnd = findobj('Tag', 'logdata');
datapath=get(hnd,'String');
if ~isempty(datapath)
[depth,val] = textread(fullfile, '%f %f', 'headerlines', 1);
maxdep=max(depth);
hnd=findobj('Tag', 'eddep');
set(hnd,'String',num2str(maxdep));
mindep=min(depth);
hnd=findobj('Tag', 'stdep');
set(hnd,'String',num2str(mindep));
hnd=findobj('Tag', 'computer');
set(hnd,'enable','on');
end;
case 'exit'
delete(gcf);
case 'cjcodesel'
hnd=findobj('Tag','curvename');
i=get(hnd,'Value');
hnd=findobj('Tag','unit');
if i==1
set(hnd,'String',cjunit(1:2));
else
set(hnd,'String',cjunit(i+1));
end;
end;
function aux=l_plot(wlog,varargin)
% Function plots log curves; a curve's properties can be changed by
% right-clicking on the curve and choosing new curve parameters from the
% pop-up menu.
%
% Written by: E. R.: May 6, 2000
% Last update: July 17, 2006: use "mysuptitle" with different location for
% landscape and portrait
%
% aux=l_plot(wlog,varargin)
% INPUT
% wlog log structure
% varargin one or more cell arrays; the first element of each cell array is a
% keyword, the other elements are parameters. Presently, keywords are:
% 'annotation' subplot (curve annotation). Possible values are 'mnemonic'
% and 'description' which refer to columns 1 or 3 of the
% field "curve_info'. Mnemonic is generally much shorter.
% Default: {'annotation','mnemonic'}
% 'figure' Specifies if new figure should be created or if the
% seismic traces should be plotted to an existing figure.
% Possible values are 'new' and any other string.
% Default: {'figure','new'}
% 'color' color (and line style) of curves. Default {'color','r'}
% Colors may be changed interactively by right-clicking on
% a curve and selecting a new color from the pop-up menu.
% 'curves' mnemonics of curves to plot. {'curves',[]} means all
% curves.
% Default: {'curves',[]}
% If S4M.interactive == 1 this brings up a list of curve
% mnemonics for interactive curve selection
% otherwise all curves are plotted
% 'orient' plot orientation; possible values are: 'landscape' and 'portrait'
% Default: for four or fewer curves: {'orient','portrait'
% for more than four curves: {'orient','landscape'}
% 'depths' Depth range (or rather range of values of first column); can be two
% comma-separated numbers or a two-element vector
% Default: the whole depth range
% 'axis_scaling' controls axis scaling (see help axis);
% possible values are: 'auto' and 'tight'
% Default: {'axis_scaling','tight'}
% OUTPUT
% aux optional output argument. Structure with fields "figure_handle" and
% "axis_handles" which contains the handles to the axes of all subplots
%
% EXAMPLES
% l_plot(wlog,{'curves','DT','RHOB'},{'orient','landscape'}) % Plot
% % sonic and density log in landscape orientation
% l_plot(wlog,{'color','b'},{'depths',4000,5000}) % Plot all logs in the depth
% % range from 4000 to 5000 (in terms of log depth units)
% % using line color blue.
% aux=l_plot(wlog) Plot all curves of log structure; output "aux",
% % structure with axis handles; changing, say,
% % the x-axis direction from normal to reverse,
% % of the third curve (e.g. a sonic log) can be
% % achieved by
% set(aux.axis_handles(3),'XDir','reverse') % Reverse the x-axis for
% % the third curve
aux.figure_handle=[];
param.annotation='mnemonic';
param.color='r';
param.curves='';
param.depths=[wlog.first,wlog.last];
param.figure='new';
param.orient=[];
param.axis_scaling='tight';
% Decode and assign input arguments
param=assign_input(param,varargin);
if iscell(param.depths)
param.depths=cat(2,param.depths{:});
end
if isempty(param.curves)
param.curves=wlog.curve_info(2:end,1);
end
if strcmp(param.annotation,'mnemonic')
idescr=1; % Curve mnemonics used as subplot titles
else
idescr=3; % Curve descriptions used as subplot titles
end
if ~iscell(param.curves)
param.curves={param.curves};
end
ncurves=length(param.curves);
if ncurves > 12
alert('The maximum number of curves that can be displayed is 12')
ncurves=12;
end
if ncurves == 1 & strcmp(param.curves,'*')
param.curves=wlog.curve_info(2:end,1)';
ncurves=length(param.curves);
end
aux.figure_handle=[];
if strcmp(param.figure,'new')
if isempty(param.orient)
if ncurves > 4
param.orient='landscape';
else
param.orient='portrait';
end
end
if strcmpi(param.orient,'landscape')
figure_handle=lfigure;
font_size=(60/max([7.5,ncurves]))+1; % Adjust font size to the number of curves to plot
elseif strcmpi(param.orient,'portrait')
figure_handle=pfigure;
font_size=(40/max([5,ncurves]))+1; % Adjust font size to the number of curves to plot
else
alert([' Unknown picture orientation (',param.orient',')'])
end
timeStamp
bgGray
set(figure_handle,'name','²â¾®ÇúÏßÏÔʾ','numbertitle','off');
figure_export_menu(figure_handle);
else
font_size=10;
end
index=find(wlog.curves(:,1) >= param.depths(1) & wlog.curves(:,1) <= param.depths(2));
if isempty(index)
error([' Log has no values in requested depth/time range: ', ...
num2str(param.depths)])
end
ier=0;
hh=zeros(ncurves,1); % Reserve room for axis handles
tracking_button=logical(1);
if ncurves > 1
switch param.orient
case 'landscape'
mysuptitle(mnem2tex(wlog.name),{'yloc',1.09})
case 'portrait'
mysuptitle(mnem2tex(wlog.name))
end
else
mytitle(mnem2tex(wlog.name))
end
for ii=1:ncurves
if ncurves > 1 % Avoid the "subplot" command if there is only one curve
hh(ii)=subplot(1,ncurves+1,ii);
else
hh=gca;
end
pos=get(hh(ii),'Position');
set(hh(ii),'FontSize',font_size,'Position',[pos(1),pos(2)-0.04,pos(3),0.75]);
[idx,ier]=curve_index1(wlog,param.curves{ii});
if isempty(idx)
disp([' Requested curve mnemonic "',param.curves{ii},'" not available'])
ier=1;
elseif length(idx) > 1
error([' More than one curve with mnemonic "',param.curves{ii},'"'])
else
if ~strcmpi(wlog.curve_info(idx(1),2),'logical') % Regular curves
plot(wlog.curves(index,idx(1)),wlog.curves(index,1),param.color)
try
axis(param.axis_scaling)
catch
disp(['Unknown value for "axis_scaling": "',param.axis_acaling,'".'])
end
else % Logical curves
x=wlog.curves(index,idx(1));
y=wlog.curves(index,1);
ya=(y(1:end-1)+y(2:end))*0.5;
lya=length(ya);
iidx=reshape([1:lya;1:lya],1,2*lya);
yy=[y(1);y(1);ya(iidx);y(end);y(end)];
xx=[0;x(iidx);x(end);x(end);0];
dxx=diff(xx);
idx1=find(dxx ~= 0);
idx2=unique([idx1;idx1+1]);
fill([0;xx(idx2);0],[yy(1);yy(idx2);yy(end)],param.color,'EdgeColor','none');
end
% Create a "menu button" if "tracking_button" is true
initiate_2d_tracking(wlog.curve_info(idx(1),:), ...
wlog.curve_info(1,:),tracking_button)
tracking_button=logical(0);
set(hh(ii),'XAxisLocation','top')
title(mnem2tex(wlog.curve_info{idx,idescr}));
units=wlog.curve_info{idx,2};
if ncurves > 1
if ~strcmpi(units,'n/a')
xlabel(units2tex(wlog.curve_info{idx,2}));
end
else
if ~strcmpi(units,'n/a')
xlabel([strrep(wlog.curve_info{idx,idescr},'_','\_'),' (',units2tex(wlog.curve_info{idx,2}),')']);
else
xlabel(strrep(wlog.curve_info{idx,idescr},'_','\_'));
end
end
set(hh(ii),'YDir','reverse')
if strcmpi(wlog.curve_info{1,2},'n/a')
dunits='';
else
dunits=[' (',units2tex(wlog.curve_info{1,2}),')'];
end
if ii == 1
ylabel([wlog.curve_info{1,idescr},dunits])
end
if ii > 1 & ii < ncurves
set(hh(ii),'YtickLabel','')
end
if ii == ncurves & ii ~= 1
set(hh(ii),'YAxisLocation','right')
ylabel([wlog.curve_info{1,idescr},dunits])
end
end
grid
end
if ier == 1
disp([' Available curves: ',cell2str(wlog.curve_info(2:end,1))])
end
% Create linked zoom
if length(hh) > 1
axes(hh(1))
mylinkaxes(hh,'y')
else
zoom on
end
% Create over-all title
linemenu % Allow interactive modification of curves
if nargout > 0
aux.figure_handle=figure_handle;
aux.axis_handles=hh;
else
clear aux
end