forked from Coastal-Imaging-Research-Network/cBathy-Toolbox
-
Notifications
You must be signed in to change notification settings - Fork 0
/
showcBathyTimeSeriesMovieNotCIL.m
76 lines (66 loc) · 2.5 KB
/
showcBathyTimeSeriesMovieNotCIL.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
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
function map = showcBathyTimeSeriesMovieNotCIL(sn,pa, map)
%
% map = showcBathyTimeSeriesMovienNotCIL(stackName,pixArraySpecs,map)
%
% load a full cBathy data collection and show as a movie. This is to
% allow debugging, i.e. letting you see if actual waves are there.
% The pixels are converted from discrete locations to a mapping in a
% regular array under the assumption that they were designed as an array.
% The array is described by user input information
%
% if pixArraySpecs is missing or empty, create a default array covering
% the entire stack collection area at a spacing of 20m. If is is a 2
% element array, it is [dx dy]. Otherwise
% pixArraySpecs = [xMin dx xMax yMin dy yMax].
%
% The mapping between stack
% columns and this array is returned as map and can presumably be used many
% times if the geometry is unchanging.
%
% If nargin is three, the third argument is assumed to be the mapping and
% a new map is not calculated.
load(sn)
% if there was no map passed in, create an empty one
if nargin < 3
map = [];
end
% if nargin < 2, create a default pa
if nargin < 2
pa = [ 20 20 ];
end
if( length(pa) == 2 )
pa = [min(xyz(:,1)) pa(1) max(xyz(:,1)) ...
min(xyz(:,2)) pa(2) max(xyz(:,2)) ];
end
% and then call this function. If it has a valid map, it will
% create only x and y. This is to keep from duplicating code
% to do that.
[x,y,map] = findInterpMap( xyz, pa, map );
Nx = length(x);
Ny = length(y);
figure(2); clf; colormap gray
data = double(data);
I = nan(Ny,Nx);
for ii = 1: length(t)
I = reshape(data(ii,map),Ny,Nx);
imagesc(x,y,log(I)); axis xy; axis equal; axis tight
pause(0.2)
end
% Copyright (C) 2017 Coastal Imaging Research Network
% and Oregon State University
% This program is free software: you can redistribute it and/or
% modify it under the terms of the GNU General Public License as
% published by the Free Software Foundation, version 3 of the
% License.
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
% You should have received a copy of the GNU General Public License
% along with this program. If not, see
% <http://www.gnu.org/licenses/>.
% CIRN: https://coastal-imaging-research-network.github.io/
% CIL: http://cil-www.coas.oregonstate.edu
%
%key cBathy
%