24 Juni 2014

Convert Coordinate to UTM

Convert Coordinate to Universal Transverse Mercator(UTM)

Saya menyimpulkan bahwa konversi coordinate ke Universal Transverse Mercator (UTM) ke 5 tahapan
  • Konversi Coordinat ke Decimal Derajat
  • Menentukan datum dan nilai nilai constannya
  • Kalkulasi Meridional Arc
  • Menghitung Coefficient  untuk Koordinat UTM
  • Konversi ke Koordinat UTM

HD Video


1.  Konversi Coordinat ke Decimal Derajat

Sebelumnya, kita perlu merubah format coordinate ke format Decimal Derajat dulu, baik dia dalam bentuk degree minute second terpisah ataupun dalam format Coordinate.

Decimal Derajat(dd) = Derajat + menit/60 + detik/3600
Contoh : 30 degrees 15 minutes 22 seconds = 30 + 15/60 + 22/3600 = 30,2561
jadi disini kita perlu mendapatkan nilai Decimal Derajat 30,2561.

Contoh : 3°26'35.40"S = 3 + 26/60 + 35.40/3600 = 3,443166667
karena ada di selatan/south, kita jadikan negatif, jadi hasilnya -3,443166667



Untuk merubah yang memiliki format string coordinate, selengkapnya bisa dilihat di 

2.  Menentukan datum dan nilai nilai constan lainnya

Datum
Equatorial Radius,
meters (a)
Polar Radius, meters (b)
Flattening
(a-b)/a
Use
NAD83/WGS84
6,378,137
6,356,752.3142
1/298.257223563
Global
GRS 80
6,378,137
6,356,752.3141
1/298.257222101
US
WGS72
6,378,135
6,356,750.5
1/298.26
NASA, DOD
Australian 1965
6,378,160
6,356,774.7
1/298.25
Australia
Krasovsky 1940
6,378,245
6,356,863.0
1/298.3
Soviet Union
International (1924) -Hayford (1909)
6,378,388
6,356,911.9
1/297
Global except as listed
Clake 1880
6,378,249.1
6,356,514.9
1/293.46
France, Africa
Clarke 1866
6,378,206.4
6,356,583.8
1/294.98
North America
Airy 1830
6,377,563.4
6,356,256.9
1/299.32
Great Britain
Bessel 1841
6,377,397.2
6,356,079.0
1/299.15
Central Europe, Chile, Indonesia
Everest 1830
6,377,276.3
6,356,075.4
1/300.80
South Asia

lat = latitude/titik lintang
long = longitude/titik bujur
long0 = central meridian of zone
k0 = scale factor = 0,9996
e = SQRT (1-b2/a2) = eccentricity sekitar 0,08, adalah eksentrisitas elips penampang bumi
e'2 = (ea / b) 2 = e2 / (1-e2) = sekitar 0,007
n = (a-b) / (a + b) 
rho = a (1-e2) / (1-e2sin2 (lat)) 3/2. Ini adalah jari-jari kelengkungan bumi pada bidang meridian. 
nu = a / (1-e2sin2 (lat)) 1/2. Ini adalah jari-jari kelengkungan bumi tegak lurus terhadap bidang meridian. Hal ini juga jarak dari titik terhadap sumbu kutub, diukur tegak lurus terhadap permukaan bumi. 
p = (long-long0) dalam radian(delta long)

Disini kita perlu menentukan datum yang digunakan, misalnya WGS84 dengan nilai Equatorial Radius = 6,378,137
Karena menggunakan trigonometri, kita perlu merubah latitude dan longitude ke radian



3. Kalkulasi Meridional Arc (S)

S = A'lat - B'sin(2lat) + C'sin(4lat) - D'sin(6lat) + E'sin(8lat), dimana lat dalam radian.
A' = a[1 - n + (5/4)(n2 - n3) + (81/64)(n4 - n5) ...]
B' = (3 tan/2)[1 - n + (7/8)(n2 - n3) + (55/64)(n4 - n5) ...]
C' = (15 tan2/16)[1 - n + (3/4)(n2 - n3) ...]
D' = (35 tan3/48)[1 - n + (11/16)(n2 - n3) ...]
E' = (315 tan4/51)[1 - n ...]
S = A'lat - B'sin(2lat) + C'sin(4lat) - D'sin(6lat) + E'sin(8lat), dimana lat dalam radians.

4.  Menghitung Coefficient  untuk Koordinat UTM

K1 = Sk0,
K2 = k0 nu sin(lat)cos(lat)/2 = k0 nu sin(2 lat)/4
K3 = [k0 nu sin(lat)cos3(lat)/24][(5 - tan2(lat) + 9e'2cos2(lat) + 4e'4cos4(lat)]
K4 = k0 nu cos(lat)
K5 = (k0 nu cos3(lat)/6)[1 - tan2(lat) + e'2cos2(lat)]

5.  Konversi ke Koordinat UTM

y = northing = K1 + K2p2 + K3p4
x = easting = K4p + K5p3
dan zona


sumber : www.uwgb. edu

1 komentar:

Komputer H Kusuma mengatakan...

Bagaimana Pak Bikin Macronya, saya belumpenah buat maro excel... Mohon pencerahannya