40
LAMPIRAN
41
Lampiran 1 Dokumentasi survei, peralatan instrumentasi akustik yang digunakan dalam penelitian.
Contoh tampilan echogram pada saat perekaman
Komputer dan stabilizer tranduser
Kapal yang digunakan untuk penelitian
Jenis tranduser split beam yang digunakan untuk sounding data akustik
GPS digunakan untuk pencatatan Posisi
42
Lampiran 2 Jenis lifeform terumbu karang yang ditemukan pada lokasi penelitian. Titik
Bentuk lifeform
Kode
Kedalaman
1
Foliose
(CF)
2,5 meter
Faviidae
2
Masive
(CM)
2 meter
Faviidae
3
Branching
(CB)
2,5 meter
Acropidae
4
Acropora Tabular
(ACT)
2 meter
Poritidae
5
Sponge
(SP)
2,5 meter
Agaricidae
6
Acropora Digitale
(ACD)
2,5 meter
Acropidae
Mushroom
(CMR)
2,5 meter
Fungiidae
7
8
Acropora Branching
(ATB)
Photo/dokumentasi
family
2 meter Acropidae
43
Lampiran 3 Uraian Listening program Rick Towler (Purnawan 2009) %readEKRaw_EY60.m % % % %
Richard Towler NOAA Alaska Fisheries Science Center Midwater Assesment and Conservation Engineering Group
[email protected]
% %
dimodifikasi oleh Jafry Manuhutu, c552070011 Mahasiswa
site='lifeform\'; file = input ( 'file yang akan diolah? ( site - site 8)','s'); rawFile = strcat ( '../jefri-1//D20090714-T120017.raw'); botFile =strcat ( '../jefri-1//D20090714-T120017.bot'); awal = input ('masukkan nilai ping awal'); akhir = input ('masukkan nilai ping akhir'); %memberikan +1 ping, agar data yang diproses adalah = jumlah akhir akhir=akhir+1; % read in raw file - only reading 120 khz disp('Reading .raw file...'); [header, rawData] = readEYRaw(rawFile, 'SampleRange',[1 600],'PingRange',[1 1800]); calParms = readEYRaw_GetCalParms(header,rawData); % read in the .bot file - return data as range disp('Reading .bot file...'); [header, botData] = readEYBot(botFile, calParms, rawData, ... 'ReturnRange', true); % convert power to Sv data = readEYRaw_Power2Sv(rawData,calParms); % convert electrical to physical angles data = readEYRaw_ConvertAngles(data, calParms); % mensortir kembali data yang digunakan %sehingga mempermudah dalam pengolahan data dasar perairan c=1.543;% kecepatan suara tau=0.1280000;%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); SS=10*log10(ss); 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
44
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+l,ll)=0; alongl (kk,ll)=0; elseif y(kk,1)>(bd(1,ll)+0.5); Svbottom(kk,ll)=-1000; alongl(kk,ll)=0; else svbottom(kk,ll)=Z(kk,ll); along (kk,ll)=along (kk,ll); % mengambil data hanya pada dasar perairan hingga setengah meter, svbonly m=m+1; Svbottomly(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.3; for ll=1:l; for i=1:m; if y(i,1)<=hlyr; Svbonly(i,ll)=Svbottomly (i,ll); along3(i,ll)=along2(i,ll); end;end;end Svbottommean=mean(mean(Svbonly)); Svbperping=mean(Svbonly) [i l]=size(Svbonly); for ll=1:ll;Zmax(ll)=-999; for ii=1:i; if Svbonly(ii,11) >Zmax(ll); Zmax(ll)=Svbonly(ii,ll); alongmax(ll)=along3(ii,ll); end end end zmax=10.^(Zmax/10); ratazmax=mean(zmax); ssmax=zmax*(c*tau/2); SSmax=10*log10(ss); 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 %plote chogram readEKRaw_SimpleEchogram(Z,x,y,'Threshold',[60,0])%colormap('jet'); %colorbar;
45
colorbar('YTickLabel',{'-120','-80','-60','-40','-20','0'}) %plot the bottom hold on plot(data.pings(n).number,botData.pings.bottomdepth(n,:),'c'); hold off end akhir=akhir-1;%mengembalikan nilai dari 'akhir'diatas %Zmax1=0; alongmax1=0;%untuk merubah kembali pingnya sp=akhir-1;%selisih ping yang dimasukkan dengan untuk looping for ll=awal:akhir; Zmax1(ll)=Zmax(ll); end %plot(Svmax) %plot(Svmax(:,1)) figure plot(Zmax1); axis([awal akhir -50 0]); xlabel('pings','fontsize',8); ylabel('Sv Maksimum(db)','fontsize',8); legend('Sv max (dB)') %merubah nilai Sv maksimum menjadi SS maksimum zmax1=10.^(Zmax1/10); ssmakping=zmax1*(c*tau/2); SSmakping=10*log10(ssmakping) %plot(SSmax) %plot(SSmax(:,1)) figure plot(SSmakping); axis([awal akhir -65 0]); xlabel('pings','fontsize',8); ylabel('SS Maksimum(db)','fontsize',8); legend('SS max (dB)') %plot(rata2 sv bottom tiap ping) figure plot(Svbperping); axis([awal akhir -60 0]); xlabel('pings','fontsize',8); ylabel('Sv rata-rata(db)','fontsize',8); legend('Sv rata2 per ping (dB)') %merubah nilai Sv rata2 menjadi SS rata2 svrata1=10.^(Svbperping/10); ssrataping=svrata1*(c*tau/2); SSrataping=10*log10(ssrataping) %plot(SS rata2) %plot(SSrata(:,1)) figure plot(SSrataping); axis([awal akhir -40 0]);
46
xlabel('pings','fontsize',8); ylabel('SS rata2(db)','fontsize',8); legend('SS rata-rata per ping (dB)') %axis([awal akhir -10 10]) %xlabel('pings','fontsize',16); %ylabel('sudut (derajat)','fontsize',16); %legend('Sudut alongship (derajat)') bottom1=data.pings.Sv(100:200,:); bottom2=data.pings.Sv(200:300,:); mean1=mean(bottom1'); minE1=min(mean1) maxE1=max(mean1) mE1=mean(mean1) svbottom2=z; bd2=bd*2; [k l]=size(z); for ll=1:l; for kk=1:k; if y(kk,1)
(bd2(1,ll)+100); svbottom2(kk,ll)=0; else svbottom2(kk,ll)=z(kk,ll); svbottomonly2(kk,ll)=z(kk,ll); end;end;end; mean2=mean(bottom2'); minE2=min(mean2) maxE2=max(mean2) mE2=mean(mean2) %subplot(SS(:,1),y,'r') %hold on %plot (z(:,1),y,'b') %legend ('SS','SV')
