SpectroscopyAnalysis/Spectroscopy_Data_Analysis.m at master · abhaykunayak/SpectroscopyAnalysis · GitHub

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
% Reference 'LineSpec.m' by Jonathan. \MATLAB\STM Data Acquisition\LineSpec.m
% Read spectroscopy data from *.3ds file using Jonathan's modified code
% 'load3dsV.m'

%% Load .3ds file
[FILENAME, PATH] = uigetfile('*.3ds');
[X, Y, V, Z, DATA, header, pm] = load3dsV([PATH FILENAME]);

%% Auto load relevant channel
chs = 1:numel(header.channels);
% Spectroscopy channel
LS_ch_no = chs(cellfun(@(x) strcmp(x, 'LIX 1 omega (A)')...
    |strcmp(x, 'LIX 1 omega [AVG] (A)')...
    |strcmp(x, 'Input 2 [AVG] (V)'),header.channels));
LS = squeeze(DATA(:,:,:,LS_ch_no));

LS_op_ch_no = chs(cellfun(@(x) strcmp(x, 'LIY 1 omega (A)')...
    |strcmp(x, 'LIY 1 omega [AVG] (A)')...
    |strcmp(x, 'Input 3 [AVG] (V)'),header.channels));
LS_Y = squeeze(DATA(:,:,:,LS_op_ch_no));

% Current channel
I_ch_no = chs(cellfun(@(x) strcmp(x, 'Current (A)')...
    |strcmp(x, 'Current [AVG] (A)'),header.channels));
I = squeeze(DATA(:,:,:,I_ch_no));

%% Remove bad sweeps;
[LS, I] = remove_bad_sweeps(DATA, 2);

%% Fix data scale
LS = LS.*1e12;
I = I.*1e9;
X = X*1e9;
Y = Y*1e9;
Z = Z*1e9;

% offs = 0;                   % Offset for clarity in nS
% sensitivity = 20e-3;        % LI7265 sensitivity
% scale = 0.1*sensitivity;
% dV = 3e-3;                  % OSC amplitude
% LS = LS*scale/dV;

%% Select lines
lines_selected = 1:7;
LS = LS(lines_selected,:,:);
I = I(lines_selected,:,:);
Z = Z(lines_selected,:);
Y = Y(lines_selected);

%% Subtract Linear fit
[Z, slope] = subtract_linear_fit(Z);

%% Correct Slope
slope_avg = mean(slope);
theta = atand(slope_avg);
X = X./cosd(theta);

%%
X = X./1.1715;
%% Drift Correction
[Z, LS] = drift_correct(Z, LS);

%% Plot Topography
% plot_topography(X,Y,diff(Z,1,2));
plot_topography(X,Y,Z);

%% Data Analysis Each Energy SLice
browse_fft_slider(X,Y,V,LS,I,Z);

%% Data Analysis Energy Dispersion
data_analysis(X,Y,Z,V,LS,I);