www.pudn.com > logtools.rar > log_datumadj.m, change:2004-08-13,size:2988b


function ylogout=log_datumadj(ylog,xlog,xd,yd,xdold,ydold)  % ylogout=log_datumadj(ylog,xlog,xd,yd,xdold,ydold) % ylogout=log_datumadj(ylog,xlog,xd,yd) % % LOG_DATUMADJ performs datum adjustment on a log's vertical coordinates. % The datum shift is accomplished % by removing the old datum and installing the new one. The % installation of a datum is done % by subtracting the y coordinate of the datum at the logs x coordinate % from the log's y coordinate vector. % Datum removal is the opposite of datum installation. % Datum shifts are rounded to the nearest whole sample (sample rate is  %	taken to be ylog(2)-ylog(1) ) % % ylog ... vector of y (vertical) coordinates for log % xlog ... scalar giving the logs x coordinate % (xd,yd) ... piecewise linear specification of the new datum % (xdold,ydold) ... piecewise linear specification of the old %		datum % % Default for the old datum is all zeros %  % NOTE: It is illegal for you to use this software for a purpose other % than non-profit education or research UNLESS you are employed by a CREWES % Project sponsor. By using this software, you are agreeing to the terms % detailed in this software's Matlab source file.   % BEGIN TERMS OF USE LICENSE % % This SOFTWARE is maintained by the CREWES Project at the Department % of Geology and Geophysics of the University of Calgary, Calgary, % Alberta, Canada.  The copyright and ownership is jointly held by  % its author (identified above) and the CREWES Project.  The CREWES  % project may be contacted via email at:  crewesinfo@crewes.org %  % The term 'SOFTWARE' refers to the Matlab source code, translations to % any other computer language, or object code % % Terms of use of this SOFTWARE % % 1) Use of this SOFTWARE by any for-profit commercial organization is %    expressly forbidden unless said organization is a CREWES Project %    Sponsor. % % 2) A CREWES Project sponsor may use this SOFTWARE under the terms of the  %    CREWES Project Sponsorship agreement. % % 3) A student or employee of a non-profit educational institution may  %    use this SOFTWARE subject to the following terms and conditions: %    - this SOFTWARE is for teaching or research purposes only. %    - this SOFTWARE may be distributed to other students or researchers  %      provided that these license terms are included. %    - reselling the SOFTWARE, or including it or any portion of it, in any %      software that will be resold is expressly forbidden. %    - transfering the SOFTWARE in any form to a commercial firm or any  %      other for-profit organization is expressly forbidden. % % END TERMS OF USE LICENSE  if(nargin< 5 ) 	xdold=xd; 	ydold=zeros(size(xd)); end  	%compute the shifts 	%first the shifts to remove the old datum 	dy=ylog(2)-ylog(1); 	delyold=interpextrap(xdold,ydold,xlog,0); 	delyold=dy*round(delyold/dy); 	%now to install the new one 	dely=interpextrap(xd,yd,xlog,0); 	dely= dy*round( dely/dy );  	%combine the shifts 	dely=dely-delyold;  	ylogout=ylog-dely;