Fits multiple transects at a single location with a single line. Inputs are number of files (z) and data matrix (Z)(see ReadFiles.m). Outputs are the intercept (b) and slope (b) of the best fit line. (adapted from code by J. Czuba) P.R. Jackson, USGS, 12-9-08
0001 function [m,b] = VMT_MeanXSLine(z,A) 0002 % Fits multiple transects at a single location with a single line. Inputs 0003 % are number of files (z) and data matrix (Z)(see ReadFiles.m). Outputs are 0004 % the intercept (b) and slope (b) of the best fit line. 0005 % 0006 % (adapted from code by J. Czuba) 0007 % 0008 % P.R. Jackson, USGS, 12-9-08 0009 0010 0011 0012 %% Determine the best fit mean cross-section line from multiple transects 0013 % initialize vectors for concatenation 0014 0015 x = []; 0016 y = []; 0017 0018 for zi = 1 : z 0019 0020 % concatenate long and lat into a single column vector for regression 0021 x=cat(1,x,A(zi).Comp.xUTM); 0022 y=cat(1,y,A(zi).Comp.yUTM); 0023 0024 plot(A(zi).Comp.xUTM,A(zi).Comp.yUTM,'r') 0025 plot(A(zi).Comp.xUTMraw,A(zi).Comp.yUTMraw,'b') 0026 0027 end 0028 0029 % find the equation of the best fit line 0030 whichstats = {'tstat','yhat'}; 0031 stats = regstats(y,x,'linear', whichstats); 0032 0033 % mean cross-section line slope and intercept 0034 V.m = stats.tstat.beta(2); 0035 V.b = stats.tstat.beta(1); 0036 0037 clear x y stats whichstats zi