-
Notifications
You must be signed in to change notification settings - Fork 14
Expand file tree
/
Copy pathusgsMatch.m
More file actions
123 lines (117 loc) · 3.92 KB
/
usgsMatch.m
File metadata and controls
123 lines (117 loc) · 3.92 KB
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
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
%% read binary usgs file
fileHandle= fopen('./data/Cup95RefEm.sli', 'rb');
cupLib = fread(fileHandle, [50, 38], 'single');
cupSpectrum = 1000 * [1.990800, 2.000900, 2.010900, 2.020900, 2.030900, 2.040900, 2.050900, ...
2.060900, 2.071000, 2.081000, 2.091000, 2.101000, 2.111000, 2.121000, ...
2.130900, 2.140900, 2.150900, 2.160900, 2.170900, 2.180900, 2.190800, ...
2.200800, 2.210800, 2.220800, 2.230700, 2.240700, 2.250600, 2.260600, ...
2.270600, 2.280500, 2.290400, 2.300400, 2.310400, 2.320300, 2.330200, ...
2.340200, 2.350100, 2.360000, 2.370000, 2.379900, 2.389800, 2.399700, ...
2.409600, 2.419600, 2.429500, 2.439400, 2.449300, 2.459200, 2.469100, ...
2.479000];
% extract neccesary band
figure(1)
subplot(231)
emMvc = HMvcOrig(:,7)/1000;
emMvcAscl1_2 = HMdcAscl1_2Orig(:,7)/1000;
emNfindr = HNfindrOrig(:,7)/1000;
emRef = cupLib(:,5);
plot(cupSpectrum,emRef, 'k--', 'LineWidth', 2)
ylim([0, 1])
hold on
plot(cupSpectrum, emMvc, 'r', 'LineWidth', 2)
plot(cupSpectrum, emMvcAscl1_2, 'b', 'LineWidth', 2)
plot(cupSpectrum, emNfindr , 'g', 'LineWidth', 2)
ylabel('reflectance(%)')
xlabel('wavelength (nm)')
title('Nontronite', 'FontSize', 20)
legend('Lib', 'MVC', 'MDSC', 'NFINDR')
sad(emRef, emMvc)
sad(emRef, emNfindr)
subplot(232)
emNfindr = HNfindrOrig(:,3)/1200;
emRef = cupLib(:,11);
emMvc = HMvcOrig(:,3)/1400;
emMvcAscl1_2 = HMdcAscl1_2Orig(:,3)/1400;
plot(cupSpectrum, emRef, 'k--', 'LineWidth', 2)
ylim([0, 1])
hold on
plot(cupSpectrum, emMvc, 'r', 'LineWidth', 2)
plot(cupSpectrum, emMvcAscl1_2, 'b', 'LineWidth', 2)
plot(cupSpectrum, emNfindr , 'g', 'LineWidth', 2)
ylabel('reflectance(%)')
xlabel('wavelength (nm)')
title('Chalcedony','FontSize', 20)
legend('Lib', 'MVC', 'MDSC', 'NFINDR')
sad(emRef, emMvc)
sad(emRef, emNfindr)
subplot(233)
emMvc = HMvcOrig(:,4)/500;
emMvcAscl1_2 = HMdcAscl1_2Orig(:,4)/1000;
emNfindr = HNfindrOrig(:,4)/1000;
emRef = cupLib(:,2);
plot(cupSpectrum, emRef, 'k--', 'LineWidth', 2)
ylim([0, 1])
hold on
plot(cupSpectrum, emMvc, 'r', 'LineWidth', 2)
plot(cupSpectrum, emMvcAscl1_2, 'b', 'LineWidth', 2)
plot(cupSpectrum, emNfindr , 'g', 'LineWidth', 2)
ylabel('reflectance(%)')
xlabel('wavelength (nm)')
title('Kaolinite', 'FontSize', 20)
legend('Lib', 'MVC', 'MDSC', 'NFINDR')
sad(emRef, emMvc)
sad(emRef, emNfindr)
subplot(234)
emMvc = HMvcOrig(:,5)/600;
emMvcAscl1_2 = HMdcAscl1_2Orig(:,5)/600;
emNfindr = HNfindrOrig(:,5)/600;
emRef = cupLib(:,1);
plot(cupSpectrum, emRef, 'k--', 'LineWidth', 2)
ylim([0, 1])
hold on
plot(cupSpectrum, emMvc, 'r', 'LineWidth', 2)
plot(cupSpectrum, emMvcAscl1_2, 'b', 'LineWidth', 2)
plot(cupSpectrum, emNfindr, 'g', 'LineWidth', 2)
ylabel('reflectance(%)')
xlabel('wavelength (nm)')
title('Alunite','FontSize', 20)
legend('Lib', 'MVC', 'MDSC', 'NFINDR')
sad(emRef, emMvc)
sad(emRef, emNfindr)
subplot(235)
emMvc = HMvcOrig(:,8)/400;
emMvcAscl1_2 = HMdcAscl1_2Orig(:,8)/400;
emNfindr = HNfindrOrig(:,8)/400;
emRef = cupLib(:,4);
plot(cupSpectrum, emRef, 'k--', 'LineWidth', 2)
ylim([0, 1])
hold on
plot(cupSpectrum, emMvc, 'r', 'LineWidth', 2)
plot(cupSpectrum, emMvcAscl1_2, 'b', 'LineWidth', 2)
plot(cupSpectrum, emNfindr, 'g', 'LineWidth', 2)
ylabel('reflectance(%)')
xlabel('wavelength (nm)')
title('Andradite', 'FontSize', 20)
legend('Lib', 'MVC', 'MDSC', 'NFINDR')
sad(emRef, emMvc)
sad(emRef, emNfindr)
subplot(236)
emMvc = HMvcOrig(:,9)/400;
emMvcAscl1_2 = HMdcAscl1_2Orig(:,9)/400;
emNfindr = HNfindrOrig(:,9)/400;
emRef = cupLib(:,6);
plot(cupSpectrum, emRef, 'k--', 'LineWidth', 2)
ylim([0, 1])
hold on
plot(cupSpectrum, emMvc, 'r', 'LineWidth', 2)
plot(cupSpectrum, emMvcAscl1_2, 'b', 'LineWidth', 2)
plot(cupSpectrum, emNfindr, 'g', 'LineWidth', 2)
ylabel('reflectance(%)')
xlabel('wavelength (nm)')
title('Muscovite', 'FontSize', 20)
legend('Lib', 'MVC', 'MDSC', 'NFINDR')
sad(emRef, emMvc)
sad(emRef, emNfindr)
% plot(HMvcOrig(:,3)/200, '--', 'LineWidth', 2)
fclose(fileHandle);