DAFTAR PUSTAKA
Djunarsjah, E. 2007. Aspek Teknik Hukum Laut. Diktat Kuliah. Penerbit ITB. Bandung.
Djunarsjah, E. 2007. Konsep Penetapan Batas Laut. Slide Kuliah. KK Sains dan Rekayasa Hidrografi ITB. Bandung.
IHO. 2006. A Manual on Technical Aspect of the United Nations Convention on the Law of The Sea, Special Publication No. 51 (4th edition), International Hydrographic Bureau, Monaco.
Krakiwsky, E. J. 1978. Geodetic Position Computation. Lecture Notes Geodesy dan Geomatics Engineering UNB. Kanada.
Lesmana, H. 2005. Penetapan Batas Laut Teritorial Antara Indonesia dengan Malaysia di Sekitar Pulau Sipadan dan Ligitan. Tugas Akhir. ITB. Bandung.
Pemerintah RI. 2002. Peraturan Pemerintah No. 38 Tahun 2002 tentang Daftar Koordinat Geografis Titik-titik Garis Pangkal Kepulauan Indonesia. Jakarta.
Rais, J. 2003. Pedoman Penentuan Batas Wilayah Laut Kewenangan Daerah Menurut UU No.22/1999. Koleksi Dokumen Proyek Pesisir 1997 – 2003. Jakarta.
Umaryono. 2001. Hitung dan Proyeksi Geodesi. Diktat Kuliah. Teknik Geodesi ITB. Bandung.
DP - 1
Vincenty. 1975. Direct and Inverse Solutions of Geodesics on the Ellipsoid with Application of Nested Equations. Survey Review Vol. XXIII. Directorate of Overseas Surveys of the Ministry of Overseas Development. Surrey.
Veness, C. 2006. Vincenty formula for distance between two Latitude/Longitude points. http://www.movable-type.co.uk/scripts/LatLongVincentyDirect.html.
_______. Direct and Inverse Geodetic Problem. Surveying Engineering Department. Ferris State University.
_______. Geocentric Datum of Australia Technical Manual Version 2.2. Intergovernmental Committee on Surveying and Mapping. Australia.
DP - 2
LAMPIRAN
Lampiran A
Program Matlab Rumus Vincenty
Lampiran B
Daftar Koordinat Titik-titik Dasar Di Sebelah Barat Pulau Sumatra Sesuai PP No. 38 tahun 2002
LP‐1
LAMPIRAN A Program Matlab Rumus Vincenty
LP‐2
Program Matlab Untuk Rumus Vincenty :
format long g
%masukkan nilai lintang
derajatli =input('Masukkan Nilai derajat lintang =');%Dalam koma menitli =input('Masukkan Nilai menit lintang =');%Dalam koma sekonli=input('Masukkan Nilai sekon lintang =');%Dalam koma
%penentuan lintang secara koordinat
lat1=(derajatli+(menitli/60)+(sekonli/3600))
%masukkan nilai bujur
derajatbu =input('Masukkan Nilai derajat bujur =');%Dalam koma menitbu =input('Masukkan Nilai menit bujur =');%Dalam koma sekonbu=input('Masukkan Nilai sekon bujur =');%Dalam koma
%penentuan bujur secara koordinat
lon1=(derajatbu+(menitbu/60)+(sekonbu/3600))
if abs(lat1)>90 error('nilai yang dimasukkan harus antara -90 sampai +90 derajat.') end a = 6378137; % semimajor axis b = 6356752.31424518; % semiminor axis f = 1/298.257223563; % flattening coefficient LP‐3
lat1 = lat1 * .1745329251994329577e-1; % intial latitude in radians lon1 = lon1 * .1745329251994329577e-1; % intial longitude in radians % correct for errors at exact poles by adjusting 0.6 millimeters: kidx = abs(pi/2-abs(lat1)) < 1e-10;
if any(kidx); lat1(kidx) = sign(lat1(kidx))*(pi/2-(1e-10)); end
%masukkan sudut jurusan
derajatjrsn =input('Masukkan Nilai derajat sudut jurusan =');%Dalam koma menitjrsn =input('Masukkan Nilai menit sudut jurusan =');%Dalam koma sekonjrsn=input('Masukkan Nilai sekon sudut jurusan =');%Dalam koma
alpha=(derajatjrsn+(menitjrsn/60)+(sekonjrsn/3600))
alpha1 = alpha * .1745329251994329577e-1; % inital azimuth in radians sinAlpha1 = sin(alpha1); cosAlpha1 = cos(alpha1); tanU1 = (1-f) * tan(lat1); cosU1 = 1 / sqrt(1 + tanU1*tanU1); sinU1 = tanU1*cosU1; sigma1 = atan2(tanU1, cosAlpha1); sinAlpha = cosU1 * sinAlpha1; cosSqAlpha = 1 - sinAlpha*sinAlpha; uSq = cosSqAlpha * (a*a - b*b) / (b*b); A = 1 + uSq/16384*(4096+uSq*(-768+uSq*(320-175*uSq))); B = uSq/1024 * (256+uSq*(-128+uSq*(74-47*uSq)));
LP‐4
%masukkan nilai jarak s=input('Masukkan Nilai Jarak='); %dalam meter
sigma = s / (b*A); sigmaP = 2*pi; while (abs(sigma-sigmaP) > 1e-12) cos2SigmaM = cos(2*sigma1 + sigma); sinSigma = sin(sigma); cosSigma = cos(sigma); deltaSigma = B*sinSigma*(cos2SigmaM+B/4*(cosSigma*(-1+... 2*cos2SigmaM*cos2SigmaM)-... B/6*cos2SigmaM*(-3+4*sinSigma*sinSigma)*(-3+... 4*cos2SigmaM*cos2SigmaM))); sigmaP = sigma; sigma = s / (b*A) + deltaSigma; end tmp = sinU1*sinSigma - cosU1*cosSigma*cosAlpha1; lat2 = atan2(sinU1*cosSigma + cosU1*sinSigma*cosAlpha1,... (1-f)*sqrt(sinAlpha*sinAlpha + tmp*tmp)); lambda = atan2(sinSigma*sinAlpha1, cosU1*cosSigma - ... sinU1*sinSigma*cosAlpha1); C = f/16*cosSqAlpha*(4+f*(4-3*cosSqAlpha)); L = lambda - (1-C) * f * sinAlpha * (sigma + C*sinSigma*(cos2SigmaM+... C*cosSigma*(-1+2*cos2SigmaM*cos2SigmaM))); lon2 = lon1 + L; % output degrees lat2 = lat2 * 57.295779513082322865 lon2 = lon2 * 57.295779513082322865; lon2 = mod(lon2,360) % follow [0,360] convention if nargout > 2 a21 = atan2(sinAlpha, -tmp); LP‐5
a21 = 180 + a21 * 57.295779513082322865; % note direction reversal a21=mod(a21,360); end return
LP‐6
LAMPIRAN B Daftar Koordinat Titik-titik Dasar Di Sebelah Barat Pulau Sumatra Sesuai PP No. 38 tahun 2002
LP‐7
LP‐8
LP‐9
LP‐10
