-
Notifications
You must be signed in to change notification settings - Fork 4
/
DepthToDepositionConversion.m
44 lines (43 loc) · 1.99 KB
/
DepthToDepositionConversion.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
%% Depth to Deposition
%
% Allocation
RadarDeposition = RadarDepth;
depositionAgeModel{ii} = cell(1,nFiles);
DepositionAxis = cell(1,nFiles);
for ii = 1:nFiles
RadDeposition = zeros(size(RadarDeposition{ii}));
% newAge = zeros(size(AgeModel{ii}));
RadDepth = RadarDeposition{ii};
% Axis for Deposition Image
zStak = zStack{ii};%DepthAxis{ii};
ageStak = AgeModel{ii};
DepositionAxis{ii} = (linspace(min(ageStak(:)),max(ageStak(:)),size(RadarDeposition{ii},1)))';
DepositAxe = DepositionAxis{ii};
n = length(DepositAxe);
% depths = zeros(size(RadarDeposition{ii}));
% % Find the Max Depths
% % Grid Search for Depths of Isochrones at an evenly sampled rate.
% DepositAxeMat = repmat(DepositAxe,1,length(DepositAxe));
% % for (kk = 1:size(RadarDeposition{ii},2))
% parfor (kk = 1:size(RadarDeposition{ii},2),nWorkers)
% [~,tmpIx] = min(abs(DepositAxeMat - ageStak(:,kk)'));
% maxDepth = max(zStak(tmpIx));
% depths(:,kk) = linspace(0,maxDepth,n);
% end
% Do the Converison
parfor (kk = 1:size(RadarDeposition{ii},2), nWorkers)
% for kk = 1:size(RadarDeposition{ii},2)
% [Depth, Age] - Array to be resampled
% zacurve = [depths(:,kk),ageStak(:,kk)];
zacurve = [zStak(:,kk),ageStak(:,kk)];
% Perform Depth Conversion
[RadDeposition(:,kk)] = timeDepthConversionNonLin(RadDepth(:,kk),zacurve,DepositAxe);
% Convert Age-Depth Model
% Yes this Works, but it is simpler to compute this directly
% [newAge(:,kk)] = timeDepthConversion(ageStak(:,kk),zacurve,DepositAxe);
end
RadarDeposition{ii} = RadDeposition;
% depositionAgeModel{ii} = newAge;
depositionAgeModel{ii} = ones(1,size(RadarDeposition{ii},2)).*DepositionAxis{ii};
end
clear('RadDeposition','RadDepth','DepositAxe','zStak','ageStak');