-
Notifications
You must be signed in to change notification settings - Fork 4
/
DepositionToDepthInversion.m
29 lines (29 loc) · 1.23 KB
/
DepositionToDepthInversion.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
%% Deposition to Depth Inversion
%
% Allocation
RadarDepth = RadarDeposition;
% Change to Age Depth Matrix
for ii = 1:nFiles
RadDeposit = RadarDeposition{ii};
RadDepth = zeros(size(RadarDeposition{ii}));
% Axis for Deposition Image
zStak = DepthAxis{ii};
% if isPickAgeHorizons || isLoadIRH
% ageStak = depositionAgeModel{ii};
% DepositionAxis{ii} = (linspace(0,max(ageStak(:)),size(RadarDeposition{ii},1)))';
% DepositAxe = DepositionAxis{ii};
% else
ageStak = AgeModel{ii};
DepositionAxis{ii} = (linspace(0,max(ageStak(:)),size(RadarDeposition{ii},1)))';
DepositAxe = DepositionAxis{ii};
% end
% Grid Search for Depths Occurs in DepthtoDepositionConversion.m
% Convert Deposition Time Image to Depth Image
parfor (kk = 1:size(RadarDeposition{ii},2), nWorkers)
% [Age,Depth] - Array to be resampled
azcurve = [DepositAxe(:),zStak];
[RadDepth(:,kk)] = timeDepthConversion(RadDeposit(:,kk),azcurve,depths(:,kk));
end
RadarDepth{ii} = RadDepth;
end
clear('RadDeposit','RadDepth','DepositAxe','zStak','ageStak');%'RadarDeposition','DepositionAxis','depths'