PROGRAM DENGAN MATLAB UNTUK VISUALISASI DATA BATIMETRI 2 DIMENSI YANG DIOVERLAY DENGAN ARUS PERMUKAAN

Syntax program matlab tampilan batimetri 2D dan arus permukaan:

% -------------------------------------

% Program Batimetri dan Arus Permukaan

% Oleh : Henky Wibowo

% -------------------------------------

data = load('batiarus.txt');

p = 136.2084:0.001:136.2556; q = -1.9251:0.001:-1.8746;

[XI,YI] = meshgrid(p,q);

ZI = griddata(data(:,1),data(:,2),data(:,3),XI,YI,'v4');

contour(XI,YI,ZI)

v = -199 : 5 : 0;

[c,h] = contourf(XI,YI,ZI,v)

x = (data(:,1));

y = (data(:,2));

u = gradient(data(:,4));

v = gradient(data(:,5));

hold on

h = quiver(x,y,u,v);

hold off

Title('Batimetri 2D dan Arus Permukaan','fontweight','bold','fontsize',12);

xlabel('Bujur','fontweight','bold');

ylabel('Lintang','fontweight','bold');

Penjelasan per baris dari Syntax:

·         Baris 1: memuat data

·         Baris 2: memanggil data dengan batas minimum dan maksimum serta interval tertentu. Untuk p (bujur) minimum 136.2084 dan maksimum 136.2556 dengan interval 0.001. Untuk q (lintang) minimum -1.9251 dan maksimum -1.8746 dengan interval 0.001.

·         Baris 3: XI dan YI adalah p dan q yang telah mengalami meshgrid. Meshgrid berfungsi merubah domain vektor, p dan q menjadi matriks XI dan YI.

·         Baris 4: ZI merupakan data pada kolom 1, kolom 2, kolom 3, XI, dan YI yang telah di grid dengan menggunakan fungsi griddata dengan metode ‘v4’ yang merupakan MATLAB4 griddata method.

·         Baris 5: membentuk kontur dari data XI, YI, dan ZI.

·         Baris 6: data v dipanggil dengan nilai minimum -199 dan maksimum 0 dengan interval 5.

·         Baris 7: memberikan warna pada tiap-tiap interval kedalaman yang berbeda.

·         Baris 8: menjelaskan bahwa x adalah data yang dimuat dari kolom 1.

·         Baris 9: menjelaskan bahwa y adalah data yang dimuat dari kolom 2.

·         Baris 10: menjelaskan bahwa u adalah gradien dari data yang dimuat dari kolom 4. Gradient ini berfungsi untuk menentukan arah dari data.

·         Baris 11: menjelaskan bahwa u adalah gradien dari data yang dimuat dari kolom 5.

·         Baris 12: menahan gambar yang ditampilkan sehingga pada perintah berikutnya gambar dapat ditampilkan bersama dengan gambar sebelumnya.

·         Baris 13: menampilkan kecepatan arus dengan komponen (u,v) pada titik (x,y)

·         Baris 14: me-reset perintah sebelum menampilkan gambar baru atau membuat default.

·         Baris 13: membuat judul gambar.

·         Baris 14: membuat keterangan pada sumbu x.

·         Baris 15: membuat keterangan pada sumbu y.

Oleh: Henky Wibowo