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
| Lon=load([rep 'Longitude.txt']);
Lat=load([rep 'Latitude.txt']);
CN=load([rep 'sur_refl_b01_1.txt']);
ReflB1=CN/10000/pi; % reflectance
%Grid data
nlon=1;
for t=2:size(Lat)
if(Lat(t)~=Lat(t-1)), break; end;
nlon=nlon+1;
end
nlat=length(Lat)/nlon;
Xp=zeros(nlat,nlon);
Yp=zeros(nlat,nlon);
RB1=zeros(nlat,nlon);
for j=1:nlat
Xp(j,1:nlon)=Lon((j-1)*nlon+1:j*nlon);
Yp(j,1:nlon)=Lat((j-1)*nlon+1:j*nlon);
for i=1:nlon
ReflB1t=ReflB1((j-1)*nlon+i);
RB1(j,i) =ReflB1t; % reflectance à cartographier
end
end
fig=figure(1);clf %% Polynomiale - Turb
set(fig, 'visible', 'off')
pcolor(Xp,Yp,RB1); shading interp;
caxis([0 0.05]) ; axis equal tight ; colormap(jet2);colorbar('FontSize',20,'Fontweight','bold');
saveas(gcf,[rep 'RB1_Insitu_Prof.png'],'png'); |
Partager