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
|
%slab chilien
clear all
close all
catalogues={'ISC.txt' 'NEIC.txt' 'CSEM.txt' 'GUC.txt' 'SJA.txt'};
Donnee=cell(numel(catalogues),1);
for n=1:numel(catalogues)
Donnee{n}= load(catalogues{n});%dlmread(catalogues{n},',');
end
Latitude=[Donnee{1}(:,1);Donnee{4}(:,1);Donnee{5}(:,1)];
Longitude=[Donnee{1}(:,2);Donnee{4}(:,2);Donnee{5}(:,2)];
Profondeur=[-Donnee{1}(:,3);-Donnee{4}(:,3);-Donnee{5}(:,3)];
Magnitude=[Donnee{1}(:,4);Donnee{4}(:,4);Donnee{5}(:,4)];
for n=1:numel(Latitude)
Latitude(n)=round((10*Latitude(n)))/10;
Longitude(n)=round(10*Longitude(n))/10;
end
C=[Latitude,Longitude,Profondeur,Magnitude];
figure
plot(Magnitude,Profondeur,'r.')
xlabel('Magnitude')
ylabel('Profondeur')
title('Magnitude fonction de la profondeur')
C=sortrows(C,[3]);
x=Latitude;
y=Longitude;
z=Profondeur;
xi=linspace(min(x),max(x),30);
yi=linspace(min(y),max(y),30);
[XI YI]=meshgrid(xi,yi);
ZI = griddata(x,y,z,XI,YI);
contour3(XI,YI,ZI);
colormap hsv
colorbar
figure
%subplot(121)
%h=surf(XI,YI,ZI)
h=plot3(XI,YI,ZI)
%colorbar(blue)
%subplot(122)
hold on
fittedmodel = fit([x, y], z, 'cubicinterp');
evaluation=fittedmodel(x,y);
plot(fittedmodel)
%Comparaison |
