56
Lampiran 1. Alat dan Bahan yang digunakan di Lapangan
Scientific Echosounder Simrad EY 60
Kapal Survei
Pipa Paralon berdiameter 7,6 cm (3 inch) dan Sekop
Dongle Echoview
57
Lampiran 2. Foto Tipe Substrat Dasar Perairan di Lokasi Penelitian
Substrat Pasir
Substrat Pasir berlumpur
Lampiran 3. Alat Pengukur Parameter Fisik Sedimen
Shaker (ASTM E – 11. USA Standard)
Timbangan Digital
58
Lampiran 4. Listing Program Matlab Rick Towler untuk menampilkan Grafik Echogram, SV dan SS (Purnawan, 2009) % readEKRaw_EY60.m %----------------------------------------------------------% % Rick Towler % National Oceanic and Atmospheric Administration % Alaska Fisheries Science Center % Midwater Assesment and Conservation Engineering Group %
[email protected] %----------------------------------------------------------% % Dimodifikasi oleh Manik, H.M Dosen Akustik, P.S. Ilmu Kelautan, Institut Pertanian Bogor % readEKRaw_ChunkExample.m % define paths to example raw and bot files rawFile = 'nama_file.raw'; botFile = 'nama_file.bot'; awal=input('masukkan ping awal = '); akhir=input('masukkan ping akhir = '); % membaca file .raw - hanya pada frekuensi 120 kHz disp('membaca .raw file...'); [header, rawData] = readEKRaw(rawFile, 'SampleRange', [1 500],... 'PingRange', [awal akhir]); calParms = readEKRaw_GetCalParms(header, rawData); % membaca file .bot - data yang kembali sebagai range disp('membaca .bot file...'); [header, botData] = readEKBot(botFile, calParms, rawData, ... 'ReturnRange', true); % konversi power ke Sv data = readEKRaw_Power2Sv(rawData, calParms); % konversi sudut electrical ke sudut physical data = readEKRaw_ConvertAngles(data, calParms); % mensortir kembali data yang digunakan % sehingga mempermudah dalam pengolahan data dasar perairan c=1546.35;%kecepatan suara tau=0.000128;%panjang gelombang x=data.pings.number; y=data.pings.range; Z=data.pings.Sv;% Z= Sv logaritma z=10.^(Z/10); ss=z*(c*tau/2);%ini untuk cari ss SS=10*log10(ss);%in untuk cari SS log along=data.pings.alongship; %sudut alongship athw=data.pings.athwartship; % sudut athwartship Svbottom=Z; along1=along; bd=botData.pings.bottomdepth; [k l]=size(Z); % data tbd pada 1 ping terakhir memberikan nilai yang tidak akurat % sehingga perlu dihilangkan l=l-1; for ll=1:l; m=0; for kk=1:k; % mengambil data dasar perairan, dari permukaan hingga 1/2 meter % data yang lainnya diberikan pada kedalaman lain adalah nol if y(kk,1)<(bd(1,ll)+0.05); Svbottom(kk,ll)=-1000; %svbottom(kk+1,ll)=0;
59
along1(kk,ll)=0; elseif y(kk,1)>(bd(1,ll)+0.5); Svbottom(kk,ll)=-1000; along1(kk,ll)=0; else svbottom(kk,ll)=Z(kk,ll); along1(kk,ll)=along(kk,ll); % mengambil data hanya pada dasar perairan hingga setengah meter, svbonly m=m+1; Svbottomonly(m,ll)=Z(kk,ll); along2(m,ll)=along(kk,ll); end;end; end; % agar jumlah data tiap kolom sama % ditentukan ketebalan lapisan yang digunakan, hlyr hlyr=0.1; for ll=1:l; for i=1:m; if y(i,1)<=hlyr; Svbonly(i,ll)=Svbottomonly(i,ll); along3(i,ll)=along2(i,ll); end; end; end Svbottommean=mean(mean(Svbonly)); [i l]=size(Svbonly); for ll=1:l;Zmax(ll)=-999; for ii=1:i; if Svbonly(ii,ll) > Zmax(ll) ; Zmax(ll) = Svbonly(ii,ll); alongmax(ll)=along3(ii,ll); end end end zmax=10.^(Zmax/10); % linier func ratazmax=mean(zmax); ssmax=zmax*(c*tau/2); ssmean=mean(ssmax); SSmax=10*log10(ssmean) stdsv=std(zmax); rataZmax=10*log10(ratazmax) stdSv=10*log10(stdsv); % membuat gambar echogram dan anglogram disp('Plotting...'); nFreqs = length(data.pings); for n=1:nFreqs % plot echogram readEKRaw_SimpleEchogram(SS,x,y, 'Threshold', [-50,0]);% disini ngerubahnya!!!!! % plot the bottom hold on plot(data.pings(n).number, botData.pings.bottomdepth(n,:), 'c'); hold off % plot anglogram readEKRaw_SimpleAnglogram(data.pings(n).alongship, ... data.pings(n).athwartship, data.pings(n).number, ... data.pings(n).range, 'Title', ... ['Angles ' num2str(calParms(n).frequency)]); % plot bottom
60
hold on plot(data.pings(n).number, botData.pings.bottomdepth(n,:), 'c'); hold off end akhir=akhir-1; %mengembalikan nilai dari 'akhir' di atas %Zmax1=0; alongmax1=0;% untuk merubah kembali pingnya sp=akhir-l; %selisih ping yang dimasukkan dengan untuk looping for ll=awal:akhir; Zmax1(ll)=Zmax(ll-sp); alongmax1(ll)=alongmax(ll-sp); end figure subplot(2,1,1); plot(Zmax1); axis([awal akhir -30 0]) xlabel('ping','fontsize',16); ylabel('Scattering volume (dB)','fontsize',16); legend('Sv max (dB)') subplot(2,1,2); plot(alongmax1); axis([awal akhir -10 10]) xlabel('ping','fontsize',16); ylabel('sudut (derajat)','fontsize',16); legend('sudut alongship (derajat)') figure plot (SS(:,1),y,'r') hold on plot (Z(:,1),y,'b') legend('SS','Sv') xlabel ('Intensitas acoustic backscattering strength (dB)') ylabel ('Depth (m)')
61
Lampiran 5. Listing Program Matlab untuk menampilkan Echo Envelope %readEKRaw_EY60.m %----------------------------------------------------------% % Rick Towler % National Oceanic and Atmospheric Administration % Alaska Fisheries Science Center % Midwater Assesment and Conservation Engineering Group %
[email protected] %----------------------------------------------------------% % Dimodifikasi oleh Manik, H.M Dosen Akustik, P.S. Ilmu Kelautan, Institut Pertanian Bogor %readEKRaw_ChunkExample.m % define paths to example raw and bot files rawFile = 'nama_file.raw'; botFile = 'nama_file.bot'; ping_awal = input('masukkan ping awal = '); ping_akhir = input('masukkan ping akhir = '); disp('Reading .raw file...'); % read in the first chunk of the file using PingRange to define chunk size. % Note that we specify the optional 3rd return argument "rstat" that will % contain the reader state when the function exits. % % also note that we do not read in angle data [header, firstRaw, rstat] = readEKRaw(rawFile, 'Frequencies', 120000, ... 'SampleRange',[1 800],'PingRange',[ping_awal ping_akhir],'Angles',false); % extract calibration parameters from the first raw data structure calParms = readEKRaw_GetCalParms(header, firstRaw); disp('Reading .bot file...'); % read in the .bot file - by passing the optional 3rd argument we force % readEKBot to only return data for pings contained in the firstRaw structure. % again, we set the rstat return argument. [header, firstBot, rstat] = readEKBot(botFile, calParms, firstRaw, ... 'ReturnRange', true); % convert power to Sv firstRaw = readEKRaw_Power2Sv(firstRaw, calParms); % plot up the two blocks of data disp('Plotting...'); % plot the first chunk echogram readEKRaw_SimpleEchogram(firstRaw.pings(1).Sv, firstRaw.pings(1).number, ... firstRaw.pings(1).range, 'Threshold', [-70,0], 'Title', ... ['Sv']); hold on % plot the bottom plot(firstRaw.pings(1).number, firstBot.pings.bottomdepth(1,:), 'c'); hold off
62
% colorbar; colorbar('YTickLabel',{'-70 dB','-58 dB','-47 dB','-35 dB','-23 dB','-12 dB'}) xlabel ('Ping') ylabel ('Depth (m)') % plot echo envelope (digunakan setelah selesai menampilkan echogram) Sv1=firstRaw.pings.Sv; Sv1mean=mean(Sv1'); plot(Sv1mean); xlabel ('Time (ms)') ylabel ('Intensitas energi backscattering strength (dB)')
63
Lampiran 6. Gambar Grafik Echogram (Lanjutan…)
Stasiun 1
Stasiun 2
Stasiun 5
Stasiun 6
Stasiun 7
64
Lampiran 7. Gambar Grafik Pola SV dan SS (Lanjutan…)
Stasiun 1
Stasiun 2
Stasiun 5
Stasiun 6
Stasiun 7
65
Lampiran 8. Gambar Grafik Intensitas Energi Backscattering (Lanjutan…)
Stasiun 1
Stasiun 2
Stasiun 5
Stasiun 6
Stasiun 7
66
Lampiran 9. Hasil Olahan Fraksi Sedimen di Balai Penelitian Tanah Laboratorium Fisika Tanah IPB Bogor No.
Tekstur (%) Liat
Urut Seri/ Substrat
Lanau
Lanau
Lanau
Pasir
Pasir
Pasir
Pasir
Pasir
0 – 2 μm 2 –10 μm 10 – 20 μm 20 – 50 μm 50 – 100 μm 100 – 200 μm 200 – 500 μm 50 – 1000 μm 1000 – 2000 μm
1.
11 F/ STA 1
0.9
3.27
5.67
12.98
7.58
11.27
23.57
29.87
4.89
2.
11 F/ STA 2
0.82
4.31
6.83
15.67
4.57
13.61
20.15
31.26
2.78
3.
11 F/ STA 3
1.15
2.63
5.61
8.25
12.78
15.65
29.33
17.85
6.75
4.
11 F/ STA 4
0.89
2.39
6.62
11.74
4.56
7.88
53.57
7.84
4.51
5.
11 F/ STA 5
1.08
2.13
4.55
9.84
6.41
5.58
23.65
31.28
15.48
6.
11 F/ STA 6
1.28
3.67
6.88
15.31
14.64
19.84
21.38
12.36
4.64
7.
11 F/ STA 7
0.24
0.95
4.25
7.58
4.43
18.47
34.12
17.58
12.38
8.
11 F/ STA 8
0.73
1.25
2.13
5.63
20.34
16.94
31.26
18.45
3.27
9.
11 F/ STA 9
0.38
1.78
3.6
9.35
5.16
13.62
35.41
22.14
8.56
