-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathLK_to_NDVI_Multibanda.m
209 lines (162 loc) · 5.88 KB
/
LK_to_NDVI_Multibanda.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
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
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
function LK_to_NDVI_Multibanda (chkImagenes,chkProceso,val_Master,val_B1,Master_db,B1_db,calib_Master,calib_B1,x1_Master,x1_B1,reflectancias_Master,reflectancias_B1,tipoCmap,boolHist,status_cuad,cuad_div,rgb_g_limits,auxiliar_limits,status_suelo,check_aux,status_NDVI,path_Master,path_B1)
%---------------------------------------------------------------%
% SW de Obtencion de NDVI a partir de irradiancia (Lk)
%---------------------------------------------------------------%
%
%
% LK: Nivel de radiancia
% NDVI: Indice de vigor de planta
%---------------------------------------------------------------%
%% Obtener valores de irradiancia de cada banda
ND_IR = ND_IR_to_LK (chkImagenes,chkProceso,val_Master,Master_db,calib_Master,x1_Master,reflectancias_Master,path_Master);
ND_RED = ND_R_to_LK (chkImagenes,chkProceso,val_B1,B1_db,calib_B1,x1_B1,reflectancias_B1,path_B1);
%% Calcular valores NDVI
switch status_NDVI
case 1
ima=(ND_IR-ND_RED)./(ND_IR+ND_RED);
case 2
ima=((ND_IR-ND_RED)./(ND_IR+ND_RED+0.5))*1.5; %SAVI
case 3
ima=(((ND_IR-ND_RED)./(ND_IR+ND_RED))+1)/2; %%NDVI infragram
case 4
ima=(ND_IR)./(ND_RED); %PaD
case 5
ima=(((ND_IR)./(ND_RED))+1)/2; %PaD Infragram
end
%Eliminación NAN=-1
s_ima=size(ima);
for i=1:s_ima(1)
for j=1:s_ima(2)
if isnan(ima(i,j))==1
ima(i,j)=-1;
end
end
end
minimo = min(min(ima));
maximo = max(max(ima));
if chkProceso==1
ima2=figure;
imshow(ima);
title('Índices de Vegetación');
colorbar;
end
s_cmap=size(tipoCmap);
for i=1:s_cmap(2)
if tipoCmap(i)==1
selection=i;
end
end
switch selection
case 1 %RGB/Gris
g_div=64;
cmap=zeros(g_div,3);
for i=1:g_div %Preparar colormap escala grises
cmap(i,1)=(i-1)/(g_div-1); cmap(i,2)=cmap(i,1); cmap(i,3)=cmap(i,1);
end
if status_suelo==1 %Opción suelo blanco
if rgb_g_limits(1)*g_div>=2 %Comprobar donde empieza el primer límite
for i=1:round(rgb_g_limits(1)*g_div) %Asignar color blanco
cmap(i,1)=1; cmap(i,2)=cmap(i,1); cmap(i,3)=cmap(i,1);
end
end
end
if rgb_g_limits(1)*g_div<0.5 %Comprobar el primer límite para índice diferente 0
for i=1:round(rgb_g_limits(2)*g_div)
cmap(i,1)=0; cmap(i,2)=0; cmap(i,3)=1;
end
else
for i=round(rgb_g_limits(1)*g_div):round(rgb_g_limits(2)*g_div)
cmap(i,1)=0; cmap(i,2)=0; cmap(i,3)=1;
end
end
for i=round(rgb_g_limits(3)*g_div)+1:round(rgb_g_limits(4)*g_div)
cmap(i,1)=0; cmap(i,2)=1; cmap(i,3)=0;
end
for i=round(rgb_g_limits(5)*g_div)+1:round(rgb_g_limits(6)*g_div)
cmap(i,1)=1; cmap(i,2)=0; cmap(i,3)=0;
end
if check_aux(1)==1
for i=round(auxiliar_limits(1)*g_div):round(auxiliar_limits(2)*g_div)
cmap(i,1)=1; cmap(i,2)=cmap(i,1); cmap(i,3)=cmap(i,1);
end
end
if check_aux(2)==1
for i=round(auxiliar_limits(3)*g_div):round(auxiliar_limits(4)*g_div)
cmap(i,1)=1; cmap(i,2)=cmap(i,1); cmap(i,3)=cmap(i,1);
end
end
if check_aux(3)==1
for i=round(auxiliar_limits(5)*g_div):round(auxiliar_limits(6)*g_div)
cmap(i,1)=1; cmap(i,2)=cmap(i,1); cmap(i,3)=cmap(i,1);
end
end
colormap(cmap);
case 2 %Jet
cmap=jet;
colormap(jet);
case 3 %Rainbow
rainbow=zeros(10,3);
rainbow(1,1)=1; rainbow(1,2)=0; rainbow(1,3)=0;
rainbow(2,1)=1; rainbow(2,2)=127/255; rainbow(2,3)=0;
rainbow(3,1)=1; rainbow(3,2)=1; rainbow(3,3)=0;
rainbow(4,1)=0; rainbow(4,2)=1; rainbow(4,3)=0;
rainbow(5,1)=0; rainbow(5,2)=1; rainbow(5,3)=127/255;
rainbow(6,1)=0; rainbow(6,2)=1; rainbow(6,3)=1;
rainbow(7,1)=0; rainbow(7,2)=127/255; rainbow(7,3)=1;
rainbow(8,1)=0; rainbow(8,2)=0; rainbow(8,3)=1;
rainbow(9,1)=127/255; rainbow(9,2)=0; rainbow(9,3)=1;
rainbow(10,1)=1; rainbow(10,2)=0; rainbow(10,3)=1;
cmap=rainbow;
colormap(rainbow);
end
if status_cuad==1
% s_ima=size(ima);
for i=1:cuad_div(1)-1
hold on
x_cuad(1)=i*s_ima(2)/cuad_div(1);
x_cuad(2)=x_cuad(1);
y_cuad(1)=0;
y_cuad(2)=s_ima(1);
plot(x_cuad,y_cuad,'Color','k','LineStyle',':');
end
for j=1:cuad_div(2)-1
hold on
x_cuad(1)=0;
x_cuad(2)=s_ima(2);
y_cuad(1)=j*s_ima(1)/cuad_div(2);
y_cuad(2)=y_cuad(1);
plot(x_cuad,y_cuad,'Color','k','LineStyle',':');
end
end
% current_process=waitbar(1,'Calculando Índice de Vegetación');
% pause on; pause(1.5); pause off;
% close(current_process);
if boolHist==1
ima_aux=ima;
for i=1:s_ima(1)
for j=1:s_ima(2)
if ima_aux(i,j)<0
ima_aux(i,j)=0;
end
end
end
figure;
imhist(ima_aux);
end
if chkProceso==1
max_min=msgbox(strcat('Máximo: ',num2str(maximo),' Mínimo: ',num2str(minimo)),'Valor Máximo y Mínimo');
waitfor(max_min);
choice=questdlg('Recalibrar índices?','Recalibración Índices de Vegetación','Sí','No','Sí');
if strcmp(choice,'Sí')==1
NDRVI_Recalibration_IV(0,ima,cmap,maximo,minimo,status_cuad,cuad_div);
end
choice=questdlg('Guardar imagen?','Guardar imagen','Sí','No','Sí');
if strcmp(choice,'Sí')==1
% Guardar imagen procesada con colormap
NDRVI_im_save(ima,cmap);
end
else
[~, name_im, ~] = fileparts(path_Master);
%imwrite(ima,strcat('Procesadas/',name_im,'_NDVI','.jpg'));
NDRVI_im_save(name_im,chkProceso,ima,cmap);
end