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 

convert dms format decimal

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