Remote Sensing (Pengolahan Data Angin Negara Panama September 1999)

clear all;

% rubah nama file sesuai kebutuhan !!!

namafile   = 'usept99.dat'; % file yg dibuka (format binary)

namahasil  = 'u1999.tmp'; % file hasil temporary (ascii)

namahasil1 = 'u1999_hsl'; % file hasil (ascii)

% buka dan baca file data yang mau diolah...

fid = fopen(namafile,'rb')

% buka di sini utk thn 1988 --> tahun kabisat jadi ada 366 hari x 4 = 1464 waktu

%[p,count]= fread(fid,[475,1464],'real*4');  %26barisx34kolomx1464waktu utk th.1988

% buka di sini utk thn 1989 --> bukan tahun kabisat jadi ada 365 hari x 4 = 1460 waktu

[p,count]= fread(fid,[12,120],'real*4'); %26barisx34kolomx1460waktu utk th.1989

fclose(fid)

% proses tahap 1

%for i=1:1464;   %tutup di sini kalau mau buka yg tahun 1989

% buka di sini utk thn 1989

for i=1:120;

fid=fopen(strcat(namahasil,num2str(i),'.dat'),'wb');

var1=p(:,i);

var1=var1';

var1=fliplr(var1);

% jumlah kolom data sesuai gambar

% (34 kolom)

fprintf(fid,'%6.2f %6.2f %6.2f %6.2f  n',var1);

fclose(fid);

end;

% proses tahap 2

%for i=1:1464;  %tutup di sini kalau mau buka yg tahun 1989

% buka di sini utk thn 1989

for i=1:120;

fid = fopen(strcat(namahasil1,num2str(i),'.dat'),'wb');

namafile1=strcat(namahasil,num2str(i),'.dat');

s=load(namafile1);

s=fliplr(s);

s2=s';

% (34 kolom)

fprintf(fid,'%6.2f %6.2f %6.2f %6.2f n',s2);

fclose(fid);

end;

% hapus file yg gak perlu yg dihaslkan dari proses tahap 1

%for i=1:1460;   %tutup di sini kalau mau buka yg tahun 1989

% buka di sini utk thn 1989

for i=1:120;

namafile1=strcat(namahasil,num2str(i),'.dat');

delete (namafile1);

end;

whos

Sedangkan untuk vector V prosesnya hampir sama kita tinggal mengganti u dengan v

clear all;

% rubah nama file sesuai kebutuhan !!!

namafile   = 'vsept99.dat'; % file yg dibuka (format binary)

namahasil  = 'v1999.tmp'; % file hasil temporary (ascii)

namahasil1 = 'v1999_hsl'; % file hasil (ascii)

% buka dan baca file data yang mau diolah...

fid = fopen(namafile,'rb')

% buka di sini utk thn 1988 --> tahun kabisat jadi ada 366 hari x 4 = 1464 waktu

%[p,count]= fread(fid,[475,1464],'real*4');  %26barisx34kolomx1464waktu utk th.1988

% buka di sini utk thn 1989 --> bukan tahun kabisat jadi ada 365 hari x 4 = 1460 waktu

[p,count]= fread(fid,[12,120],'real*4'); %26barisx34kolomx1460waktu utk th.1989

fclose(fid)

% proses tahap 1

%for i=1:1464;   %tutup di sini kalau mau buka yg tahun 1989

% buka di sini utk thn 1989

for i=1:120;

fid=fopen(strcat(namahasil,num2str(i),'.dat'),'wb');

var1=p(:,i);

var1=var1';

var1=fliplr(var1);

% jumlah kolom data sesuai gambar

% (34 kolom)

fprintf(fid,'%6.2f %6.2f %6.2f %6.2f  n',var1);

fclose(fid);

end;

% proses tahap 2

%for i=1:1464;  %tutup di sini kalau mau buka yg tahun 1989

% buka di sini utk thn 1989

for i=1:120;

fid = fopen(strcat(namahasil1,num2str(i),'.dat'),'wb');

namafile1=strcat(namahasil,num2str(i),'.dat');

s=load(namafile1);

s=fliplr(s);

s2=s';

% (34 kolom)

fprintf(fid,'%6.2f %6.2f %6.2f %6.2f n',s2);

fclose(fid);

end;

% hapus file yg gak perlu yg dihaslkan dari proses tahap 1

%for i=1:1460;   %tutup di sini kalau mau buka yg tahun 1989

% buka di sini utk thn 1989

for i=1:120;

namafile1=strcat(namahasil,num2str(i),'.dat');

delete (namafile1);

end;

whos

Setelah kita mengekstrak File u2009.dat dan file v2009.dat, kita lanjutkan dengan merata-ratakan semua, dengan listing program sebagai berikut :
Untuk u :

clear

jrt=zeros(4,3);

jrt=0;

for i=1:120;

namafile1=strcat('u1999_hsl',num2str(i),'.dat');

s=load(namafile1);

jrt=jrt+s;

end;

rt=(jrt)'/120;

fid = fopen(strcat('urata99.dat'),'wb');

fprintf(fid,'%6.2f %6.2f %6.2f %6.2f n',rt);

fclose(fid);

Untuk v :

clear

jrt=zeros(4,3);

jrt=0;

for i=1:120;

namafile1=strcat('v1999_hsl',num2str(i),'.dat');

s=load(namafile1);

jrt=jrt+s;

end;

rt=(jrt)'/120;

fid = fopen(strcat('vrata99.dat'),'wb');

fprintf(fid,'%6.2f %6.2f %6.2f %6.2f n',rt);

fclose(fid);

Setelah program tersebut di jalankan maka akan tampil dua file yang akan kita olah di transform yaitu urata99.dat dan vrata99.dat

Kedua file tersebut kita buka di transform. Pilih file open. Pilih teks matriks kemudian OK.

1383193709796273028

Setelah itu ubah generate scales dengan menukar first value dengan last value, diperoleh :

1383193730558889838

Lakukan dengan cara yang sama untuk vrata99.dat, maka diperoleh :