Memisahkan data batuan (sifat fisis)berdasarkan gradien dari dua parameter // separate a rocks due to phisical parameters

Memisahkan data batuan berdasarkan gradien dari dua parameter// separate a rocks data due to gradien of 2 parameters

File input.dat sebagai data //  input.dat file as data

File plo.dat sebagai  file untuk parameter garis pembatas //  plo.dat files for gradient  (separator line)

SCRIPT :

————————————————————————————————————–

% plot dan memilah

% load data dulu

clear;

clc;

load input.dat

load plo.dat

x=input(:,3);

y=input(:,2);

z=input(:,1);

a=plo(:,1);

b=plo(:,2);

figure(1)

plot(x,y,’*'),hold on

plot(a,b,’-*r’)

%————–

f1=polyfit(a,b,1)

y1=polyval(f1,x)

%————

figure(3)

for i=1:1:length(x)

if y(i)<=y1(i)

plot(x(i),y(i),’*g’), hold on

depth(i)=z(i);

rock(i)=100; % id batuan

else

plot(x(i),y(i),’*c’),hold on

depth(i)=z(i);

rock(i)=200; % id batuan

end

%plot(x,y1,’o')

end

plot(a,b,’-*r’)

hasil=[depth' rock']

save hasil.dat hasil -ascii

figure(4)

for i=1:1:length(x)

if y(i)<=y1(i)

plot3(x(i),y(i),z(i),’*g’),grid on, hold on

depth(i)=z(i);

rock(i)=100;

else

plot3(x(i),y(i),z(i),’*c’),grid on,hold on

depth(i)=z(i);

rock(i)=200;

end

end

plot(a,b,’-*r’)

———————————————————————————————————-

result :

before separate

before separate

After separate

After separate

Warna biru menunjukkan batuan Dolomit, warna hijau menunjukkan batuan Limestone // blue color show a dolomite, and green to limestone

File hasil pemisahan tersebut disimpan pada file “hasil.dat”. Dolomit dan limestone diidentifikasikan dengan harga yaitu 100 dan 200 (lihat warna kuning) // result  save in ” hasil.dat” . Dolomite and limestone identified with value 100 and 200