Bagian Kedua

GD2211 IHG 2
GEOMETRI
ELLIPSOID
Dosen : Kosasih Prijatna
Wedyanto Kuntjoro

Versi 2006
 ELLIPS dan ELLIPSOID
z

b b

a a y

x2 z2 x2 + y2 z2
+ =1 + =1
2 2 2 2
a b a b
 Menggambar ellips secara grafis

lingkaran
berjari-jari a
(sb. panjang)

ellips

lingkaran
berjari-jari b
(sb. pendek)
 ELLIPSOID
Kutub Utara

b
Rotational Biaxial
a
Ellipsoid a

ekuator

Parameter-parameter bentuk dan dimensi ellipsoid :

Sumbu pendek : b
Sumbu panjang : a
Pegepengan : f = (a-b)/a
Kosasih Prijatna, 2005
 ELLIPS dan ELLIPSOID
F dan F’ masing-masing adalah titik fokus ellips.

Q
Eksentrisitas pertama :
an

nia
i

2 2
id

2 a −b
rid
er

e =
m

me

a2

Eksentrisitas kedua :

2 a 2 − b2
e' =
b2

FQ + F ′Q = konstan
Ellips adalah tempat kedudukan titik-titik yang mempunyai jumlah
jarak yang tetap ke kedua titik fokusnya.
Hubungan antar parameter ellipsoid
pegepengan eksentrisitas pertama eksentrisitas kedua

a −b a2 − b2 a 2 − b2
f = e2 = e'2 =
a a2 b2

e2
(1 − e )(1 + e′2 ) = 1
2 e′2 =
1 − e2

2 e′2
e = =2f − f 2 f = 1 − 1 − e2
1 + e′2

(1 − e 2 ) = 1 − 2 f + f 2
 Beberapa Ellipsoid Referensi
Thn. Nama a (m) b (m) 1/f

1830 Airy 6377563 6356257 299.325

1830 Everest 6377276 6356075 300.802
1841 Bessel 6377397 6356079 299.153
1866 Clarke 6378206 6356584 294.978
1907 Helmert 6379200 6356818 298.300
1909 Hayford 6378388 6356912 297.000
1927 NAD-27 6378206.4 6356912 294.9786982
1948 Krassovsky 6378245 6356863 298.300
1960 WGS-60 6378165.0 6356783.3 298.3
1966 WGS-66 6378145 6356760 298.25
1967 GRS-67 6378160.0 6356774.5 298.247167427
1972 WGS-72 6378135.0 6356751 298.26
1980 GRS-80 6378137.0 6356752 298.257222101
1984 WGS-84 6378137.0 6356752 298.257223563

Kosasih Prijatna, 2005

BENTUK dan UKURAN BUMI
PENAMPANG EKUATORIAL
dari bumi (geoid global).
Pada gambar ini, perbedaan
dengan ellipsoid diperbesar
sekitar 10000 kali;
a adalah sumbu panjang
ellipsoid referensi,
Sekitar 6378 km.

Ref. Vanicek & Krakiwsky, 1986

Hasanuddin Z. Abidin, 2001

BENTUK dan UKURAN BUMI
PENAMPANG MERIDIAN NOL
dari bumi (geoid global).
Pada gambar ini, perbedaan
dengan ellipsoid diperbesar
sekitar 10000 kali;
a adalah sumbu panjang
ellipsoid referensi,
Sekitar 6357 km.

Ref. Vanicek & Krakiwsky, 1986

Hasanuddin Z. Abidin, 2001

 SISTEM KOORDINAT GEODETIK
Z Parameter sistem koordinat :
Kutub Q
• Lokasi titik nol sistem koordinat
Pusat hQ
• Orientasi sumbu-sumbu koordinat
ellipsoid
ich

ZQ • Besaran yang digunakan untuk

enw

mendefinisikan posisi titik pada

gre

sistem koordinat tersebut

Y
ϕQ Koordinat kartesian :
λQ XQ
X (XQ , YQ , ZQ )
YQ
el
lip Koordinat geodetik :
s oi
d (ϕQ , λ Q ,hQ )

XQ = (N + hQ ) cos ϕQ cos λQ
ϕ = lintang geodetik
YQ = (N + hQ ) cos ϕQ sin λ Q
λ = bujur geodetik
ZQ = [N(1 − e 2 ) + hQ ] sin ϕQ h = tinggi geodetik
N = radius lengkung vertikal utama Hasanuddin Z. Abidin, 1997
 LINTANG GEODETIK
no
z rm
al
ϕ = lintang geodetik
an xQ Q
i
id
er
m
zQ
ϕ 90o+ϕ
x

a cos ϕ a (1 − e 2 ) sin ϕ
x= z=
2 2
1 − e sin ϕ 1 − e 2 sin 2 ϕ

dz
= tan(90o + ϕ) = − cot ϕ
dx
 LINTANG REDUKSI
z
Q1 θ = lintang reduksi

Q
a QQ2 b
=
Q1Q2 a
θ ϕ
a x
O Q2

elli b x = OQ2 = a cos θ

pso
id
z = Q1Q2 = a sin θ
bo
la
 LINTANG GEOSENTRIK
z

Q ψ = lintang geosentrik

a ψ
x
O

elli b
pso
id x = r cos ψ
z = r sin ψ
r = x2 + z 2
 Hubungan antar Lintang
Berdasarkan hubungan sebagai berikut :
z
tan ψ =
x
dapat diturunkan :

( )
tan ψ = 1 − e 2 tan θ = 1 − e 2 tan ϕ

tan θ = 1 − e 2 tan ϕ ϕ : lintang geodetik

ψ : lintang reduksi
θ : lintang geosentrik
tan ϕ = 1 + e'2 tan θ
 RADIUS KELENGKUNGAN
Dari kalkulus, kelengkungan sebuah kurva y = f (x) :

y′′ dy d2y
κ= y′ = y′′ =
[
1 + ( y′)2 ]
3
2 dx dx 2

• kurva dengan kelengkungan tinggi Î κ besar

• kurva dengan kelengkungan kecil Î κ kecil

1
Radius kurva di satu titik : R=
κ

Soal :
Perlihatkan bahwa lingkaran yang berjari-jari R
mempunyai kelengkungan κ = 1/R !
 Irisan Normal pada Ellipsoid
Bidang normal adalah sebuah bidang yang berimpit dengan garis
normal ellipsoid di satu titik dan memotong permukaan ellipsoid.

P
irisan Irisan normal adalah kurva yang
normal dibentuk oleh perpotongan antara
bidang normal dengan permukaan
garis normal
ellipsoid.

Umumnya, radius kelengkungan

irisan normal di satu titik pada
permukaan ellipsoid tidak sama,
ell
ips tergantung orientasi dari bidang
oid
normalnya.

Radius kelengkungan irisan normal di setiap titik pada permukaan

bola adalah sama, tak tergantung dari orientasi bidang normalnya.
 Irisan Normal pada Ellipsoid
Untuk mengetahui kelengkungan kurva irisan normal berdasarkan formula
kalkulus, perlu diketahui terlebih dahulu model fungsi kurva tersebut.

Fungsi kurva di setiap titik di permukaan ellipsoid dan di setiap orientasi

adalah berbeda-beda.
Untuk menentukan kelengkungan kurva di setiap titik dan pada berbagai
orientasi dengan menggunakan formula kalkulus, terlebih dahulu fungsi
kurva irisan normalnya harus diketahui Î tidak praktis !

Pada setiap titik ada nilai kelengkungan kurva minimum dan maksimum.

Kelengkungan maksimum Î pada meridian

Radius kelengkungan Î minimum

Kelengkungan minimum Î pada bidang normal yang tegak lurus

terhadap meridian (vertikal utama)
Radius kelengkungan Î maksimum
 Irisan Normal pada Ellipsoid

Radius lengkung meridian (M) : Minimum di ekuator

Maksimum di kutub

Radius lengkung vertikal utama (N) : Minimum di ekuator

Maksimum di kutub

Radius kelengkungan kurva irisan normal yang orientasinya di

antara arah meridian dan irisan vertikal utama dapat ditentukan
melalui formula Euler sebagai berikut :

MN
R=
M sin 2 α + N cos 2 α
α = asimut
 Radius Lengkung Meridian
dz
= − cot ϕ
dx

−1
d 2z dϕ
1 1  dx 
= =  
2 2 dx 2  dϕ 
dx sin ϕ sin ϕ

( )
3
ell 2 2 2
ips d z − 1 − e sin ϕ 2
oid =
dx 2 a sin 3 ϕ (1 − e 2 )

( )
  dz  2  2
1 +    a 1 − e2
  dx   M=
M=
(1 − e2 sin 2 ϕ) 2
3
d 2z
dx 2
 Radius Lengkung Vertikal Utama

elli elli
pso pso
id id

Menurut teorema Meusnier :

a
( o
)
x = N sin 90 − ϕ = N cos ϕ
N=
(1 − e 2 2
sin ϕ 2 )
1
Perbandingan antara M dengan N

Pada umumnya :
N >M
elli
pso Kecuali di kutub :
id
N =M
 Perbandingan antara M dengan N

Di ekuator (ϕ = 0o) : (
M 0 = a 1 − e2 )
N0 = a

a2
Di kutub (ϕ = 90o atau ϕ = -90o) : M 90 =
b
a2
N 90 =
b
M 90 > M 0 dan N 90 > N 0

Soal : Hitung besaran-besaran di atas dengan menggunakan

ellipsoid Bessel, GRS67, dan WGS84 !
 Radius−radius lainnya
Radius rata-rata Gauss :
2π
1
RG =
2π ∫ R dα = NM R : radius Euler α : asimut
0
Radius rata-rata sumbu ellipsoid :
a+a+b
Rm =
3
Radius bola (luas bola = luas ellipsoid) :
 e 2 17 4 67 6  2
luas bola : 4πRL
L
RL = = a 1 − − e − e − ........ 
4π  6 360 3024  L : luas ellipsoid
 
Radius bola (volume bola = volume ellipsoid) :
3 2 4 3 4 2
RV = a b Volume bola : πRV Volume ellipsoid : πa b
3 3
 Panjang Busur Meridian
KU panjang busur dSϕ :
ϕ +dϕ
dSϕ ϕ dSϕ = M dϕ

M panjang busur Sϕ :
M ϕ2 ϕ2
2 dϕ
O dϕ Sϕ = ∫ M dϕ = a(1 − e ) ∫ W 3
ϕ1 ϕ1
me
rid Bentuk integral eliptik !
ia n
3
dengan W =  1 − e 2 sin 2 ϕ 
3
 
Integral di atas tidak dapat langsung diintegrasikan secara elementer.
Salah satu solusinya adalah dengan terlebih dahulu menguraikan W-3
dengan menggunakan deret MacLaurin.
 Panjang Busur Meridian
Deret Taylor :
( x − xo ) 2 ( x − xo )3
f ( x) = f ( xo ) + ( x − xo ) f ′( xo ) + f ′′( xo ) + f ′′′( xo ) + .........
2! 3!
Deret MacLaurin adalah bentuk khusus dari Deret Taylor, yaitu
untuk xo= 0 :
x2 x3
f ( x) = f (0) + x f ′(0) + f ′′(0) + f ′′′(0) + .........
2! 3!

Contoh : Uraian deret MacLaurin untuk f(x) = sin(x)

x3 x5 x 7 x9
sin( x) = x − + − + − .......... (x dalam radian)
3! 5! 7! 9!

Soal : Uraikan W-3 dengan menggunakan deret MacLaurin !

 Panjang Busur Meridian
Multiple angle formulas :
1 1
sin 2 x = − cos 2 x
2 2
3 1
sin 3 x = sin x − sin 3 x
4 4
3 1 1
sin 4 x = − cos 2 x + cos 4 x
8 2 8
5 5 5 1
sin x = sin x − sin 3 x + sin 5 x
8 16 16
5 15 3 1
sin 6 x = − cos 2 x + cos 4 x − cos 6 x
16 32 16 32

7
Soal : Tentukan pula sin x sin 8 x sin 9 x ……. !
 Panjang Busur Meridian
Apabila W -3 diuraikan dengan deret MacLaurin, diperoleh :
1 3 15 35 315 8 8
= 1 + e 2 sin 2 ϕ + e 4 sin 4 ϕ + e 6 sin 6 ϕ + e sin ϕ + ...........
3 2 8 16 128
W

Untuk mempermudah integrasi, gunakan multiple angle formulas :

1
= A − B cos 2ϕ + C cos 4ϕ − D cos 6ϕ + E cos 8ϕ − F cos 10ϕ + .............
3
W

Sehingga panjang busur meridian antara ϕ1 dan ϕ2 adalah :

 B C
S ϕ = a (1 − e 2 )  A(ϕ 2 − ϕ1 ) − (sin 2ϕ 2 − sin 2ϕ1 ) + (sin 4ϕ 2 − sin 4ϕ1 )
 2 4
D E F
− sin 6ϕ 2 − sin 6ϕ1 + sin 8ϕ2 − sin 8ϕ1 − (sin 10ϕ 2 − sin 10ϕ1 ) ] + ......
( ) ( )
6 8 10
 Panjang Busur Meridian
3 45 175 6 11025 8 43659 10
A = 1 + e2 + e4 + e + e + e + ..........
4 64 256 16384 65536

3 15 525 6 2205 8 72765 10

B = e2 + e4 + e + e + e + ...........
4 16 512 2048 65536

15 4 105 6 2205 8 10395 10

C= e + e + e + e + ............
64 256 2048 16384
35 6 315 8 31185 10
D= e + e + e + ............
512 2048 131072
315 8 3465 10
E= e + e + ............
16384 65536
693 10
F= e + ............
131072
 Panjang Busur Paralel
KU λ2
lingkaran paralel

p
p Sλ

ϕ KU
λ2−λ1
O N p
λ1
me
rid pa
ia n ral
el

Radius lingkaran paralel : p = N cos ϕ

Panjang busur paralel : Sλ = (λ 2 − λ1 ) p = (λ 2 − λ1 ) N cos ϕ

(λ2−λ1) dalam radian
 Luas Permukaan Ellipsoid
Luas elemen permukaan
paralel
dL = MN cos ϕ dϕ dλ
λ 2 ϕ2 ϕ2
merid

M dϕ dL L= ∫ ∫ MN cos ϕ dϕ dλ = (λ 2 − λ1 ) ∫ MN cos ϕ dϕ
i an

N cos ϕ dλ λ1 ϕ1 ϕ1
ϕ2
b 2  sin ϕ 1 1 + e sin ϕ  ϕ2
∫ MN cos ϕ dϕ = 
2 2
+ ln 
2 1 − e sin ϕ 2e 1 − e sin ϕ  ϕ1
ϕ1
Luas setengah permukaan ellipsoid ( λ 2 − λ1 = 2π ϕ1 = 0 ϕ2 = π 2 ) :
1 1 + e
2 1
L0o −90o = πb  + ln 
1 − e2 2e 1 − e
Luas seluruh permukaan ellipsoid : Lellipsoid = 2 L0o −90o
 IRISAN NORMAL
z
∆ ′′ Umumnya, irisan normal dari arah P1 ke P2
tidak berimpit dengan irisan normal dari
arah kebalikannya (P2 ke P1).

Bidang normal di P1 : P1 – n1 – P2
Bidang normal di P2 : P2 – n2 – P1

ae 2 sin ϕ p
zn =
(1 − e 2 2
sin ϕ p )
1/ 2

Bila kedua titik tidak terletak pada bujur dan lintang yang sama, maka :
z n1 < z n2 untuk ϕ p1 < ϕ p2

z n1 > z n2 untuk ϕ p1 > ϕ p2

 IRISAN NORMAL
Sudut perbedaan antara dua bidang normal (direct & inverse) dapat
dihitung melalui persamaan sebagai berikut :

1  ϕ p1 + ϕ p2 2 × jarak
∆ ′′ = ρ′′ e 2 σ 2 cos 2 ϕ m sin 2α p1 p2  ϕm = σ=
4  2 N p1 + N p2

Contoh :
o o
Kondisi maksimum (ϕ m = 0 dan α p1 p2 = 45 ) , jarak 200 km : ∆′′ = 0.36′′
C
Arah pengukuran sudut-sudut segitiga maupun
asimut di permukaan ellipsoid dari dua arah
yang berbeda akan tidak konsisten ! Kenapa ?

Pada praktisnya (poligon dsb), keadaan tersebut

dapat diabaikan.
A B
 Geodesik atau Garis Geodetik
Geodesik adalah garis hubung terpendek
antara dua titik di permukaan ellipsoid.

Di setiap titik sepanjang geodesik, arah

vektor radius berimpit dengan arah
normal ellipsoid.

Perbedaan antara jarak sepanjang

irisan normal dengan jarak sepanjang
geodesik (∆s) dapat dihitung melalui :

ae 4
∆s = sin 2 2α12 cos 4 ϕ m σ 5
360

Untuk jarak = 600 km, ∆s adalah sekitar 9x10-6 meter.

(dalam praktis dapat diabaikan !)