3317


GEODEZJA WYŻSZA

WYKŁAD DLA 3. ROKU GEODEZJI BY DR MAREK PLEWAKO

POWIERZCHNIE ODNIESIENIA

Geoida rzeczywista - fizyczna powierzchnia Ziemi, nie da się przedstawić jej modelem matematycznym, więc nie nadaje się do wykorzystania jako powierzchnia odniesienia.

Powierzchnia odniesienia powinna spełniać następujące warunki:

  1. geoida zerowa - przedłużenie średniego poziomu mórz i oceanów pod lądami; jest to powierzchnia ekwipotemcjalna (o stałym potencjale), wyznaczana mechanicznie (średni poziom morza wyznacza się za pomocą mareografów); geoida zerowa jest bardzo skomplikowana matematycznie, nie spełnia więc 2. warunku, służy jednak jako poziom odniesienia w niwelacji precyzyjnej.

b) elipsoida trójosiowa

0x08 graphic
0x08 graphic

c) elipsoida dwuosiowa (obrotowa) - powstaje wskutek obrotu elipsy wokół małej osi (bb); przekroje południkowe to elipsy, równoleżnikowe to koła; jest to główna elipsoida w geodezji wyższej.

d) kula - wystarcza w niektórych zagadnieniach - zależy od kryteriów dokładnościowych np.:skali (kartografia małoskalowa)

e) płaszczyzna - styczna w środku obszaru - wystarcza dla małych obszarów

[a, b, c, d, e - kolejne przybliżenia powierzchni Ziemi]

GEOMETRIA ELIPSOIDY OBROTOWEJ

0x08 graphic

0x08 graphic

Historia:

- elipsoida Bessela 1841, ulubiona w Niemczech

- elipsoida międzynarodowa Hayforda 1909

- elipsoida Krasowskiego 1940, narzucona dla bloku wschodniego, nadal wykorzystywana w Polsce ze względu na dużą ilość istniejących materiałów opartych na elipsoidzie Krasowskiego

- WGS 84, elipsoida międzynarodowa 1984r.

WSPÓŁRZĘDNE GEODEZYJNE I PROSTOKĄTNE

0x08 graphic

0x08 graphic

Równanie elipsy:

0x08 graphic
0x01 graphic

0x08 graphic
0x01 graphic

po podniesieniu do kwadratu i podstawieniu:

0x01 graphic

otrzymujemy związek między współrzędnymi XY a BL:

0x08 graphic

0x01 graphic

SZEROKOŚĆ GEODEZYJNA, GEOCENTRYCZNA, ZREDUKOWANA.

0x08 graphic

0x08 graphic

Związki:

0x01 graphic

wyprowadza się:

0x01 graphic

w przybliżeniu:

0x01 graphic

0x01 graphic

KRZYWIZNA LINII NA POWIERZCHNI, TWIERDZENIA MEUSNIERA I EULERA.

0x08 graphic
0x08 graphic

0x08 graphic

0x08 graphic

PROMIENIE KRZYWIZNY GŁÓWNYCH PRZEKROJÓW NORMALNYCH ELIPSOIDY ZIEMSKIEJ (OBROTOWEJ)

Na powierzchni elipsoidy obrotowej głównymi przekrojami normalnymi w dowolnym punkcie są:

- przekrój południkowy (ponieważ zawiera normalną w danym punkcie)

- przekrój poprzeczny, zwany też przekrojem pierwszego wertykału,

0x08 graphic

0x08 graphic

ŚREDNI PROMIEŃ KRZYWIZNY W DANYM PUNKCIE

0x01 graphic

PODSTAWOWE FUNKCJE DLA ELIPSOIDY OBROTOWEJ, W, V.

Z definicji:

0x01 graphic

Oznaczenie:

0x01 graphic

Sprawdzenie, który przekrój jest większy:

0x01 graphic

Więc:

0x01 graphic

0x01 graphic
gdy 0x01 graphic
, czyli na biegunie

WYKŁAD 2 17.10.02

OBLICZENIE DŁUGOŚCI ŁUKU POŁUDNIKA

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

Wynik rozwinięcia w szereg i całkowania:

0x01 graphic

Współczynniki 0x01 graphic
obliczamy dla każdej elipsoidy.

Dla krótkich łuków 0x01 graphic
można stosować wzór uproszczony:

0x01 graphic

0x08 graphic
0x01 graphic
- średni promień krzywizny dla 0x01 graphic
. Dokładność tego wzoru 0x01 graphic
(na naszej szerokości 0x01 graphic
).

0x08 graphic
0x08 graphic

0x08 graphic

Promień krzywizny przekroju

normalnego sięga do osi figury.

0x01 graphic
- azymut

0x01 graphic
- krzywa przekroju normalnego

0x01 graphic
i 0x01 graphic
są skośne więc łączymy punkty 0x01 graphic
i 0x01 graphic

0x01 graphic
- kąt między 0x01 graphic
i tym połączeniem

0x01 graphic

z tego wzoru wyliczamy przez odwrócenie wzoru

0x08 graphic
0x08 graphic

0x01 graphic

0x08 graphic

ORIENTACJA ELIPSOIDY WZGLĘDEM GEOIDY. ODCHYLENIA PIONU.

Orientacja elipsoidy względem geoidy na całej jej powierzchni wymaga spełnienia trzech warunków:

Tak zorientowane elipsoidy nazywamy ogólną elipsoidą ziemską.

WGS 84 - odstępy geoidy od elipsoidy nie przekraczają na całej ziemi 100m (tzw. anomalia geoidy),

w Krakowie 40m

W lokalnym aspekcie zagadnienia ważne będą inne 3 warunki:

0x08 graphic
0x08 graphic

WZAJEMNE PRZEKROJE NORMALNE

0x08 graphic

Szczególne przypadki:

0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic

TRÓJKĄT SFEROIDALNY (NA POWIERZCHNI ELIPSOIDY)

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

Dla elipsoidy o średnim promieniu obliczonym dla:

0x01 graphic

0x08 graphic

0x08 graphic
A na kuli:

Istnieją wzory:

0x01 graphic

i podstawiamy: 0x01 graphic

otrzymujemy: 0x01 graphic

Więc można małe trójkąty sferoidalne obliczać jak sferyczne, a wynikające z tego błędy będą znacznie mniejsze od dokładności pomiarów.

ROZWIĄZYWANIE MAŁYCH TRÓJKĄTÓW SFERYCZNYCH

  1. Metoda Legendre'a

  2. Metoda Soldnera (additamentów)

Metoda Legendre'a opiera się na twierdzeniu Legendre'a:

„Mały trójkąt sferyczny można rozwiązać z wielkim przybliżeniem jako trójkąt płaski o tych samych bokach i o kątach zmniejszonych o 1/3 nadmiaru sferycznegio.”

0x08 graphic
0x08 graphic

0x01 graphic

- tutaj przybliżenie i tracimy wyrazy powyżej czwartego rzędu.

Naprawdę powinno być:

0x01 graphic

Rozszerzone twierdzenie Legendre'a (uwzględnia wyraz czwartego rzędu):

0x01 graphic

Metoda Soldnera:

Kąt pozostają sferyczne a operacje wykonywane są na bokach:

Wzór sinusowy dla trójkąta sferycznego: 0x01 graphic

I jego rozwinięcie w szereg Taylora: 0x01 graphic

0x08 graphic
0x01 graphic

oznaczenie:

0x08 graphic
0x01 graphic

po zlogarytmowaniu równania *:

0x01 graphic

jeśli 0x01 graphic
to: 0x01 graphic

więc:

0x08 graphic
0x08 graphic
0x01 graphic

ogólnie:

0x08 graphic
0x01 graphic
- additament do logarytmu boku

Praktyczne postępowanie:

  1. Zmniejszamy bok wyjściowy o jego additament

0x01 graphic

  1. Obliczamy pozostałe logarytmy boków z wzorów sinusowych

0x01 graphic

  1. Do tak wyliczonych logarytmów boków dodajemy additamenty

0x01 graphic

WYKŁAD 3 24.10.02.

LINIA GEODEZYJNA (ORTODROMA)

Taka krzywa na powierzchni, która:

Na płaszczyźnie ortodroma to linia prosta, na kuli koło wielkie, na walcu każda linia śrubowa (rys.) o dowolnym skoku, na elipsoidzie obrotowej południki tak, równoleżniki nie, inne linie tak, pod warunkiem, że spełniają równanie.

0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

Równanie uwikłane powierzchni:

0x01 graphic

w teorii powierzchni udowadnia się, że:

(I) 0x01 graphic

(II) 0x01 graphic

Dla elipsoidy obrotowej:

0x08 graphic
0x01 graphic

liczymy pochodne cząstkowe:

0x01 graphic

Proporcja z równań (I) i (II):

0x08 graphic
0x01 graphic

Całkujemy

0x08 graphic
0x01 graphic

0x08 graphic

0x08 graphic

0x08 graphic

0x01 graphic

0x01 graphic

0x08 graphic
0x01 graphic

0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic
Iloczyn promienia równoleżnika i sinusa azymutu linii geodezyjnej jest stały na całej jej długości.

0x08 graphic
0x01 graphic

KSZTAŁT LINII GEODEZYJNEJ NA POWIERZCHNI ELIPSOIDY OBROTOWEJ

0x08 graphic

PRZEBIEG LINII GEODEZYJNEJ WZGLĘDEM WZAJEMNYCH PRZEKROJÓW NORMALNYCH

0x08 graphic

Dowodzi się, że

0x01 graphic

0x01 graphic

PRZENOSZENIE WSPÓŁRZĘDNYCH

0x08 graphic

Zadanie wprost:

0x01 graphic

Zadanie odwrotne:

0x01 graphic

Zadanie przeniesienia współrzędnych dzieli się na 3 działy:

Lub inny podział:

Dokładność:

0x01 graphic

0x01 graphic

na ziemi 0x01 graphic

PRZENIESIENIE WSPÓŁRZĘDNYCH METODĄ SZEREGÓW POTĘGOWYCH

0x08 graphic

metoda Legendre'a:

0x01 graphic

0x08 graphic
0x01 graphic

0x01 graphic

itd.

METODA ŚREDNIEJ SZEROKOŚCI GAUSSA

0x08 graphic

Obieramy punkt pomocniczy 0x01 graphic
, gdzie: 0x01 graphic

0x01 graphic

Jeśli oba te równania raz zsumujemy, a raz odejmiemy, to:

z sumy:

0x01 graphic

z różnicy:

0x01 graphic

0x08 graphic
0x01 graphic

Ostateczne wzory w podręcznikach.

Rozwiązanie zadania odwrotnego metodą średniej szerokości Gaussa uważa się za najlepsze z istniejących.

METODA SCHREIBERA

(do s= 150km)

0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic

METODA KRUGERA

Wymaga obliczeń w dwóch etapach:

na kuli:

0x01 graphic

promień kuli: 0x01 graphic

Kruger korzysta z szeregu Legendre'a dla kuli znajduje:

0x08 graphic
0x08 graphic
0x01 graphic

METODA CLARKE'A

(dla małych 0x01 graphic
mniejsze od 30 - 50km)

0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

Łuk południka:

0x01 graphic

0x01 graphic
jest nieznane, dlatego traktujemy tę wartość jako przybliżenie

0x08 graphic
0x08 graphic
0x01 graphic

Znając 0x01 graphic
obliczamy 0x01 graphic

0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic
0x01 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x01 graphic

Metodę tę można stosować w okolicach podbiegunowych, można też stosowa zadania odwrotne, nie potrzeba korzystać z tablic.

PRZENOSZENIE WSPÓŁRZĘDNYCH NA ŚREDNIE I DALEKIE ODLEGŁOŚCI

Na średnie odległości: (800 - 1000km)

Metoda Jordana (na szeregach Legendre'a)

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x01 graphic

inne metody:

  1. met. Bessela

  2. H - met. Helmerta

  3. L - D - met Levallois

0x01 graphic
- całka Wallisa

POMIARY STOPNIA - WYZNACZANIE ELEMENTÓW ELIPSOIDY ZIEMSKIEJ METODĄ ASTRONOMICZNO - GEODEZYJNĄ

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic

Redukcja na południk:

0x08 graphic

Podobnie można wyznaczyć elementy elipsoidy ziemskiej na podstawie równości równoleżników.

ODWZOROWANIA UŻYWANE W POLSKIEJ GEODEZJI

0x08 graphic
WSPÓŁRZĘDNE IZOMETRYCZNE

0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x01 graphic

Jeżeli dla danej powierzchni będą:

0x01 graphic

takie, że:

0x01 graphic

to takie współrzędne 0x01 graphic
nazywamy izometrycznymi.

Jeżeli

0x01 graphic
w 0x01 graphic
to 0x01 graphic
w obu przypadkach dostanie ten sam przyrost

a) na płaszczyźnie:

0x01 graphic
to współrzędne izometryczne

b) na kuli

0x01 graphic

okazuje się, że 0x01 graphic
nie są współrzędnymi izometrycznymi, zmiana o 0x01 graphic
jest mnożona przez 0x01 graphic
a zmiana o 0x01 graphic
przez 0x01 graphic
.

Przekształcamy wyrażenie na 0x01 graphic
:

0x01 graphic

0x08 graphic

0x08 graphic

0x01 graphic

0x01 graphic
- współrzędne izometryczne dla kuli

0x01 graphic
- szerokość izometryczna

0x01 graphic
- długość izometryczna

dla 0x01 graphic

c) dla elipsoidy obrotowej

0x01 graphic
- 0x01 graphic
nie są współrzędnymi izometrycznymi

0x01 graphic

0x08 graphic

0x08 graphic

0x01 graphic

0x01 graphic

0x01 graphic
- szerokość izometryczna dla elipsoidy obrotowej

Na przykład:

Kula i płaszczyzna - odwzorowanie wiernokątne:

0x08 graphic
0x08 graphic
0x01 graphic

Elipsoida i płaszczyzna:

0x08 graphic
0x08 graphic
0x01 graphic

ODWZOROWANIE GAUSSA - KRUGERA

(zwane w krajach anglosaskich poprzecznym odwzorowaniem Merkatora)

I. 0x01 graphic
wiernokątnie

0x08 graphic
0x01 graphic

II. 0x01 graphic
wiernokątnie

działamy w zmiennych zespolonych

założenia Gaussa:

1. 0x01 graphic
- obrazem południka środkowego 0x01 graphic
jest linia prosta

2. 0x01 graphic
- 0x01 graphic
przedstawia wierny obraz południka osiowego, tzn. skala na nim równa się 1

0x08 graphic

0x08 graphic
0x01 graphic

warunki odwzorowawcze:

0x08 graphic
0x01 graphic

Funkcję 0x01 graphic
można rozbijać w szereg Taylora w okolicy części rzeczywistej, przy czym:

0x01 graphic

0x01 graphic

Oddzielamy część rzeczywstą od urojonej:

0x08 graphic
0x01 graphic

Pierwszy warunek Gaussa 0x01 graphic
jest spełniony

dla 0x01 graphic
mamy 0x01 graphic
, czyli obraz osi 0x01 graphic
jest osią 0x01 graphic

Drugi warunek realizuje się poprzez obliczenie pochodnej:

0x01 graphic

0x01 graphic

obliczamy 0x01 graphic

0x01 graphic

0x01 graphic

Po podstawieniu pochodnych dostajemy ostateczne wzory Gaussa - Krugera:

0x08 graphic
0x01 graphic

długość łuku południka

0x01 graphic

0x08 graphic

Zadanie odwrotne w odwzorowaniu Gaussa - Krugera

Dane: 0x01 graphic

Szukane: 0x01 graphic

0x01 graphic

ZBIEŻNOŚĆ POŁUDNIKÓW NA ELIPSOIDZIE

0x08 graphic
0x01 graphic

Def.

Kąt zawarty między płaszczyzną południka punktu P a płaszczyzną równoległą do płaszczyzny południka początkowego oznaczamy 0x01 graphic
i nazywamy zbieżnością południków na elipsoidzie w punkcie P.

0x01 graphic

Na płaszczyźnie:

0x08 graphic

0x01 graphic
- zbieżność południków na płaszczyźnie Gaussa - Krugera

0x01 graphic

Porównując oba wzory widzimy, że:

0x08 graphic
0x08 graphic

SKALA LINIOWA W PUNKCIE P

0x01 graphic

po przekształceniu

0x01 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

Widać, że:

- 0x01 graphic
dla wszystkich punktów leżących poza południkiem osiowym

- 0x01 graphic
dla punktów na południku osiowym (bo 0x01 graphic
, a więc 0x01 graphic
0x01 graphic
II założenie Gaussa)

W odwzorowaniu wiernokątnym skala jest funkcją położenia punktów.

W odwzorowaniu Gaussa - Krugera skala rośnie z kwadratem y. Im dalej od południka osiowego tym jest większa. Inaczej jest blisko równika, inaczej w naszych szerokościach.

0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic

0x08 graphic

REDUKCJE ODWZOROWAWCZE

0x08 graphic
Redukcje liniowe:

0x08 graphic
0x01 graphic

Łuk praktycznie identyfikuje się z cięciwą (linią prostą).

0x08 graphic
0x08 graphic
0x01 graphic

Z wielkości pomierzonych na elipsoidzie uzyskujemy wielkości w odwzorowaniu.

Redukcja dla kierunków i kątów.

0x08 graphic

0x01 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x01 graphic

0x01 graphic

0x08 graphic
W trójkącie jest 6 takich poprawek 0x01 graphic
, więc:

0x01 graphic

Odwzorowanie Gaussa - Krugera w strefach 3 stopniowych:

0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x01 graphic

Np.: 0x01 graphic

ODWZOROWANIE SOLDNERA

W tym odwzorowaniu wykonano wszystkie mapy kopalń z południkiem osiowym Sucha Góra, w okresie międzywojennym wykonano też mapę Warszawy. Odwzorowanie Soldnera nadaje się tylko dla małych obszarów.

0x08 graphic

0x01 graphic

0x01 graphic
jeżeli 0x01 graphic
leży na wschód od południka centralnego

0x01 graphic
jeżeli 0x01 graphic
leży na zachód od południka centralnego

Zadanie wprost 0x01 graphic
są wzory.

Zadanie odwrotne 0x01 graphic
.

Odwzorowanie to jest częściowo wiernoodległościowe (tylko południk osiowy i wszystkie igreki), dlatego są zniekształcenia kątów, pól i długości.

0x01 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

różnica kątów (redukcja):

0x01 graphic

różnica długości (redukcja)

0x01 graphic

różnica pól

0x01 graphic

0x08 graphic

ODWZOROWANIE WIG - QUASI STEREOGRAFICZNE (ROUSSILHE'A)

0x08 graphic
0x08 graphic

TRANSFORMACJA WSPÓŁRZĘDNYCH

  1. Transformacja afiniczna jest stosowana dla małych obszarów (nie jest wiernokątna) - min. 4 punkty dostosowania.

0x08 graphic

Jeżeli są 3 punkty wspólne to można przetransformować każdy punkt leżący wewnątrz tzw. trójkąta afinicznego.

0x08 graphic
0x08 graphic
0x01 graphic

Jeżeli obszar jest duży, dzielimy go na trójkąty afiniczne.

Transformacja afiniczna przestrzenna:

0x01 graphic

0x08 graphic
0x08 graphic
0x01 graphic

  1. Przekształcenie rzutowe (używane w fotogrametrii, np. przetwornik, również niewiernokątne), dla małych obszarów o małych deniwelacjach terenu. Jest to przekształcenie płaskie tzn. płaszczyzna w płaszczyznę.

0x01 graphic

0x08 graphic
0x08 graphic
0x01 graphic

  1. Transformacja wiernokątne współrzędnych płaskich.

0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x01 graphic

Po przesunięciu obu układów do pokrycia się można napisać funkcję analityczną zmiennej zespolonej.

0x01 graphic

Dla małych obszarów można rozwijać funkcję w szeregi:

0x01 graphic

gdzie:

0x01 graphic

0x01 graphic

0x01 graphic
- wyrażają przesunięcie dlatego można opuścić 0x01 graphic
przy założeniu nasunięcia obu układów

0x01 graphic
- wyrażają obrót, skręt i zmianę skali

wobec czego:

0x01 graphic

Po wymnożeniu i rozdzieleniu części rzeczywistej i urojonej i koniec.

Trzy główne przyczyny wpływają na zmianę wartości g na powierzchni Ziemi:

1. przyśpieszenie odśrodkowe,

2. spłaszczenie Ziemi (odległość od środka masy M),

3. rozkład gęstości mas w skorupie ziemskiej.

Inne przyczyny to czynniki zmienne (przyciąganie Słońca, Księżyca, planet i skład grawimetryczny planet).

Przesunięcie jednorodnej masy kulistej m o dr wbrew sile przyciągania w polu działania siły o masie M wymaga elementarnej pracy dL:

0x08 graphic

0x08 graphic


0x01 graphic
- potencjał siły przyciągania

Potencjał siły odśrodkowej:

Praca wykonana w płaszczyźnie równoleżnika przez siłę odśrodkową f wynosi:

0x01 graphic

0x01 graphic
- potencjał siły odśrodkowej

Potencjał siły ciężkości nazywamy funkcją W

0x01 graphic

0x01 graphic
( porównaj z równaniem *)

Powierzchnie ekwipotencjalne (jednakowego potencjału):

0x01 graphic

0x01 graphic
- geoida, zawiera w sobie powierzchnię mórz i oceanów, przebieg powierzchni ekwipotencjalnych bliska Ziemi jest podobny do elipsoid obrotowych, dwie takie powierzchnie nie stykają się i nie przecinają, przez jeden punkt przechodzi tylko jedna powierzchnia ekwipotencjalna.

przebieg powierzchni ekwipotencjalnych

0x08 graphic
0x01 graphic

0x01 graphic

linie sił ciężkości = linie pionowe

0x08 graphic

0x08 graphic

Własność powierzchni ekwipotencjalnej:

- są to powierzchnie obrotowe, gdzie osią obrotu jest oś Ziemi;

- w pobliżu Ziemi są to powierzchnie zamknięte, ciągłe zbliżone do elipsoid obrotowych lub dalej od Ziemi rozciągające się do ∞ powierzchnie walcowe;

- przyspieszenie g jest na powierzchni ekwipotencjalnej najmniejsze na równiku, największe na biegunie;

- gdy g≠0 istnieje tylko jedna normalna linia pionowa do powierzchni ekwipotencjalnej;

- są powierzchniami nierównoległymi, odstępy są większe w pobliżu równika;

- w odległości 0x01 graphic
zachodzi g'= -c czyli g=0.

0x08 graphic

Przyśpieszenie siły ciężkości - wzory, pomiary, redukcje.

Wzór Clairaulta.

0x08 graphic
0x08 graphic
0x01 graphic
lub 0x01 graphic

0x08 graphic

0x08 graphic

0x08 graphic

Wyznaczenie ga i β:

w n punktach o znanym ϕ wyznaczamy g:

0x08 graphic
0x01 graphic

g - wynik pomiaru przyśpieszenia siły ciężkości

γ - oznaczenie przyśpieszenia normalnego, teoretyczne, obliczone ze wzorów typu Clairaulta dla ziemi zregularyzowanej (przyjmując pewien rozkład mas i pewną figurę Ziemi)

np.

- wzór Helmerta (1909r.) 0x01 graphic

- wzór Cassinis (1936r.) 0x01 graphic

* wynik pomiaru (g) nie musi zgadzać się z modelem teoretycznym (γ)

Pomiary przyśpieszenia siły ciężkości:

Dzielą się na:

Bezwzględne od dawna były wykonywane wahadłami, w ostatnich trzydziestu latach doszły trzy metody:

1. swobodnego spadku ciała;

2. balistycznego ruchu ciała w aparacie pomiarowym;

3. metoda z wykorzystaniem długiego giętego wahadła.

ad. 1

wahadło fizyczne

0x08 graphic

Okres:

0x01 graphic

Okres - to odstęp czasu w jakim wahadło przechodzi przez ustalony punkt toru w tej samej fazie.

Ad. 2

Największe znaczenie ma metoda druga, przeprowadził ja profesor Sakuma w Serve

0x08 graphic
1970r. ± 0,003 mgala (koło Paryża)

1973r. ± 0,001 mgala

t - czujniki laserowe mierzące momenty czasu

0x01 graphic
- podlega bezpośredniemu pomiarowi (interferometr laserowy)

0x01 graphic
0x01 graphic
0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic

0x01 graphic

W pomiarach względnych potrzebne nam jest nawiązanie do punktów, które wcześniej były wyznaczone metodami bezwzględnymi. Najczęściej wykonuje się te pomiary przy pomocy grawimetrów (przyrządów do pomiaru siły ciężkości). Najdokładniejszy grawimetr to LaCaste - Romberg (od 1930r.) typ G (geodezyjny). W Polsce tworzyło się ciągi grawimetryczne i punkty grawimetryczne (1956 - 1959r. dokładność ±0,2 - 0,3 mgala).

0x08 graphic

Redukcja pomiarów grawimetrycznych.

1. Redukcja wodnopowietrzna (Redukcja Faye'a)

0x08 graphic

0x08 graphic
0x01 graphic

„+” - tzn. przyspieszenie dla punktów nad

geoidą jest za małe

2. Redukcja Bouguara (Bugie) - uwzględnia masy

`-` - wpływ masy nad geoidą, jest więc za duże przyspieszenie

0x08 graphic

0x01 graphic

0x01 graphic
- gęstość Ziemi pod punktem P

0x01 graphic
- gęstość średnia Ziemi

3. Redukcja topograficzna - uwzględnia zarówno wzniesienie mas ponad poziom punktu P jak i ich ubytek (zasięg 30 - 40 km)

„t” topografia powoduje, że przyspieszenie jest za małe.

0x08 graphic
Mapa - sektory

HIPOTEZY IZOSTAZJI

1. Pratt (Hayforda) (1855) „różne gęstości dla różnych wysokości słupów”

0x08 graphic

2. Airy Heisanen (1856)

0x08 graphic

Skorupa o mniejszej gęstości SiAl jest zanurzona w podłożu SiMa tym głębiej im jej wzniesienie nad poziom morza jest większe, stąd nazwa „hipoteza korzeni gór”. Gęstość w obrębie słupów jest stała, natomiast zmienne są wysokość słupów, ich masa, głębokość zanurzenia, grubość skorupy.

Anomalie 0x01 graphic
przyspieszenia siły ciężkości i II twierdzenie Stokes'a

g - wynik pomiaru i redukcji (na geoidę)

0x01 graphic
- przyspieszenie normalne (teoretyczne) siły ciężkości

0x01 graphic
- anomalie przyspieszenia siły ciężkości.

0x08 graphic
wzory wiernokątnego

odwzorowania

płaskiego

0x08 graphic

transformacja Helmerta ( dla małych obszarów) - wystarczą 2 punkty wspólne -> 4 niewiadome (go, ho, g1, h1)

Współrzędne bieguna: 0x01 graphic
.

Δx, Δy - duże to transformacja Helmerta jest za mało dokładna, dlatego na większych obszarach trzeba mieć więcej niż dwa punkty dostosowania.

4. Transformacja z pasa na pas w odwzorowaniu Gaussa - Kr*gera.

0x08 graphic
0x01 graphic

0x01 graphic

0x01 graphic

h2=... g3, h3 - odpowiednie

wzory

0x01 graphic
0x01 graphic
0x01 graphic
0x01 graphic

x0, y0 można obliczyć dla danej elipsy.

Tablice Tarczy - Hornocha i Hristowa - są w nich przeliczenia między pasami. Wszystkie tablice zabezpieczają dokładności centymetrowe.

5. Transformacja współrzędnych geodezyjnych (B, L) z jednej elipsy na drugą.

Elipsoida E'

Elipsoida E

0x01 graphic

0x08 graphic
0x01 graphic


„e“ jest inne niż przy z'1.

0x08 graphic
0x08 graphic
0x08 graphic
0x01 graphic

przesunięcie współczynniki kierunkowe skale

Transformacja między odwzorowaniami:

0x01 graphic

0x01 graphic

GEODEZJA FIZYCZNA ELEMENTÓW

Elementy geodezji dynamicznej.

Prawo Newtona (o wzajemnym przyciąganiu): 0x01 graphic

Przyspieszenie siły przyciągania, siły ciężkości i odśrodkowej:

0x08 graphic
0x01 graphic

- siła przyciągania: 0x01 graphic

0x01 graphic
- przyspieszenie siły przyciągania; 0x01 graphic

- siła odśrodkowa: 0x01 graphic

0x08 graphic

F - siła przyciągania ziemskiego 0x01 graphic

f - siła odśrodkowa 0x01 graphic

P - siła ciężkości 0x01 graphic
ρ - promień równoleżnika

υ - prędkość liniowa

ω - prędkość kątowa Ziemi

c - przyspieszenie odśrodkowe.

- siła ciężkości (ciężar)

0x01 graphic
(*)

0x01 graphic

0x01 graphic

0x08 graphic
Mapy anomalii 0x01 graphic

II twierdzenie Stokes'a (1849 r.)

0x01 graphic
(całka po całej Ziemi (0;Π) gdzie: a - promień równika, G - średnia wartość ze wszystkich pomierzonych i zredukowanych wartości 0x01 graphic
, 0x01 graphic
- odległość sferyczna między punktem o znanej anomalii 0x01 graphic
a punktem gdzie wyznaczamy odstęp elipsoidy od geoidy (N), N - odstęp geoidy od elipsoidy, 0x01 graphic
- anomalie przyspieszenia siły ciężkości, 0x01 graphic
- funkcja Stokes'a

0x01 graphic

10°

30°

60°

90°

120°

180°

0x01 graphic

+1,00

+1,21

0,47

-0,90

-1,08

+0,08

0,00

0x08 graphic
0x01 graphic

Odchylenie pionu.

- kąt między kierunkiem linii pionu i kierunkiem normalnej do elipsoidy.

0x08 graphic

zg - zenit geodezyjny

za - zenit astronomiczny

B - biegun geoidy

θ - odchylenie pionu

0x01 graphic

0x01 graphic
- składowa południkowa odchylenia pionu

η - składowa w I wertykale odchylenia pionu

0x01 graphic

0x01 graphic
ma znak „+” gdy linia pionu odchyla się na północ od normalnej do elipsoidy

0x01 graphic
ma znak „+” gdy za odchyla się na wschód od zg.

Dla małych rozmiarów trójkąta za, zg, z'

0x08 graphic
0x01 graphic
0x01 graphic

0x08 graphic
0x01 graphic
0x01 graphic

0x08 graphic
0x01 graphic
- równanie Laplace'a

Wzory Vehig - Meinesza

0x01 graphic
Znając anomalie 0x01 graphic
otrzymamy odchylenia pionu (ich składowe)

NIWELACJA ASTRONOMICZNA

Pozwala wyznaczyć odstępy geoidy od elipsoidy gdy znane są względne odchylenia pionu. W tym celu obserwuje się 0x01 graphic
na punktach, gdzie B i L znane są z pomiarów geodezyjnych. Mając parami współrzędne na tych samych punktach (astronomiczne i geodezyjne) można obliczyć składowe odchyleń pionów.

0x08 graphic
0x01 graphic

0x01 graphic
0x01 graphic
0x01 graphic

0x01 graphic

0x01 graphic

0x08 graphic

0x08 graphic

0x01 graphic

0x01 graphic
- dla skończonych wartości s; na równinach si=10 - 20km, w górach si=1 - 3km.

0x01 graphic

Można również wykreślić mapę odstępów n

0x08 graphic

Można wykreślić też mapę ξ i mapę η.

REDUKCJE NA GEOIDĘ POMIARÓW ASTRONOMICZNO - GEODEZYJNYCH.

1. Redukcje bazy na geoidę.

Pomiary geodezyjne na fizycznej powierzchni Ziemi, gdzie są różne wzniesienia nad geoidą, a obliczeni są na elipsoidzie odniesienia, stąd redukcje w dwóch etapach:

I. na geoidę po liniach pionowych;

II. na elipsoidę po normalnych dla elipsoidy odniesienia.

0x08 graphic

0x01 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x01 graphic

0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic

2. Redukcja szerokości astronomicznych na geoidę.

0x08 graphic

Przez punkt P na wysokości H nad geoidą poprowadzona jest linia siły ciężkości (wklęsła ku biegunowy) i przecina ona geoidę w punkcie G. Kierunek prostej pionowej w pkt. P różni się od kierunku prostej pionowej w pkt. G. Zatem kąty jakie tworzą te proste z płaszczyzną równika, czyli szerokości ϕP , ϕG są różne. Wyprowadza się wzór:

0x01 graphic
dla R=6371km, β=0,0053

0x01 graphic

Redukcja równa się 0 dla ϕ=90° i ϕ=0° (biegun i równik), ponieważ linie sił są tam proste. Dla ϕ=45° redukcja jest największa, dla H=5830m redukcja wyniosła by 1”.

1” - 30m 0,1” - 3m

3. Redukcje azymutu astronomicznego.

3.1 Redukcja azymutu ze względu na wzniesienie celu ponad geoidę.

0x08 graphic

S - stanowisko instrumentu (na geoidzie elipsoidy)

P - cel na wysokości H nad geoidą

P' - rzut normalny pkt. P na geoidę

BS, BP' - łuki południków punktów P i S

Ns,NP - przecięcia normalnych z S i P z osią figury,

normalne s i p są wichrowate.

Prowadząc przez Ns i P płaszczyznę normalną otrzymamy przekrój normalny SP”, którego azymut A'=BSP” jest azymutem zmierzonym. Wyprowadza się :

0x01 graphic

np. dla H=1km max redukcji wynosi 0,11”.

Ponieważ azymut obarczony jest większymi błędami obserwacyjnymi, to jeśli cel nie był bardzo wysoko często redukcję można zaniedbać.

3.2 Redukcja azymutu ze względu na wzniesienie stanowiska.

0x08 graphic

P - stanowisko na wysokości H nad geoidą;

P' - rzut punktu P po linii pionu;

P” - rzut punktu P po stycznej do linii pionu;

α - azymut zredukowany;

α' - azymut zmierzony punktu B (kąt pionowy z P)

wtedy:

0x01 graphic
.

Przeważnie h w sieciach poziomych jest małe i redukcję można zaniedbać. Dodatkowo jeśli celowa znajduje się w płaszczyźnie południka to redukcja równa się 0.

4. Redukcja kątów poziomych na geoidę.

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic

0x08 graphic

Kąt poziomy jest różnicą azymutów:

0x01 graphic
, z pominięciem wzniesienia stanowiska.

HL, HK - wzniesienia punktów celowych nad geoidą.

NIWELACJA PRECYZYJNA

Systemy wysokościowe precyzyjnej niwelacji geometrycznej.

0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

Zmiana potencjału W o dW wiąże się ze zmianą wysokości H o dh.

0x01 graphic
czyli 0x01 graphic

gdzie: W - potencjał siły ciężkości

h - kierunek pionowy zwrócony do zenitu

g - przyspieszenie siły ciężkości

Ciąg niwelacji od punktu O (stacja mareograficzna) do punktu B:

0x01 graphic
(1)

dla zaobserwowanych dh i g możemy obliczyć całkę (1) tzn. różnicę potencjałów dla zaobserwowanych ciągów

1) wysokość dynamiczna i poprawka dynamiczna zaobserwowanego wzniesienia:

„Wysokością dynamiczną punktu B nazywamy ujemną różnicę potencjału siły ciężkości w rozważanym punkcie i w poziomie morza, podzielona przez normalną wartość przyspieszenia siły ciężkości w poziomie morza dla

0x01 graphic
- 0x01 graphic

0x01 graphic
(2)

dla wzniesienia dynam. Z punktu B nad punktem A otrzymujemy

0x01 graphic

przechodząc od różniczki dh do przyrostów skończonych otrzymamy:

0x08 graphic
0x01 graphic

Jak widać niezbędne są oprócz pomiarów wysokości pomiary g.

Własności wysokości dynamicznych:

1. Wysokości dynamiczne nie zależą od drogi niwelowania.

2. Wysokości dynamiczne punktów tej samej powierzchni poziomej są jednakowe.

2) Wysokość ortometryczna i poprawka ortometryczna zaobserwowanego wzniesienia.

„Wysokością ortometryczną punktu nazywamy jego odległość od geoidy (poziomu morza) mierzona po linii pionowej.”

0x08 graphic

0x08 graphic
0x08 graphic

0x01 graphic
(4)

0x01 graphic
- zmienna wartość przyśpieszenia siły ciężkości we wnętrzu Ziemi na odcinku linii pionowej B0 B

0x01 graphic
- przeciętna wartość tego przyspieszenia na tej samej linii B0 B

Różnicę potencjału na geoidzie i punktu B możemy przedstawić jako całkę po linii pionowej od punktu B0 do B.

(4)=(1)

0x01 graphic

takie przyjęcie 0x01 graphic
uniemożliwia dokładne wyznaczenie wysokości ponieważ trzeba znać wartość 0x01 graphic
w każdym punkcie pionu B0 B.

Własności wysokości ortometrycznych:

1. Wysokość ortometryczna nie zależy od drogi niwelacji do punktu B i będzie różna dla

różnych punktów tej samej powierzchni ekwipotencjalnej przechodzącej przez punkt B.

0x01 graphic
nie są równe (WA i WB - nie są to dwie równoległe powierzchnie)

Wzniesienie ortometryczne punktu B nad punktem A: 0x01 graphic

Wzniesienie ortometryczne punktu A nad punktem B: 0x01 graphic

Ponieważ powierzchnie nie są równoległe 0x01 graphic
ciąg od A do B

0x01 graphic
- suma wzniesień ortometrycznych w zamkniętej pętli niwelacyjnej 0x01 graphic
0

0x01 graphic
- „rozwartość ciągu”

i=1, 2, ... n - poszczególne stanowiska w ciągu niwelacyjnym między punktami A i B

Poprawka ortometryczna:

  1. wzniesienia: 0x01 graphic

0x01 graphic
0x01 graphic

poprawka ortometryczna wzniesienia zależy od średniej wysokości ciągu, od średniej szerokości ciągu oraz od ukierunkowania (ϕB - ϕA)

b) wysokość punktu B: 0x01 graphic

0x01 graphic

3) Wysokość normalna (geopotencjalna) i poprawka normalna zaobserwowanego wzniesienia

0x08 graphic

Wprowadził ją Rosjanin M. S. Mołodieński w latach 40-tych.

Biorąc pod uwagę trudności określenia poprawnych wysokości ortometrycznych i ścisłego wyznaczenia geoidy wynikające z nieznajomości rozkładu mas zakłócających w skorupie ziemskiej M. S. Mołodieński opracował metodę wyznaczenia figury fizycznej powierzchni Ziemi na podstawie obserwowanych tylko na powierzchni Ziemi „wartości potencjału siły ciężkości i jego gradientu „g”.

W metodzie tej rezygnuje się z wyznaczenia geoidy, otrzymuje się natomiast ściśle określona powierzchnię, bliską geoidy tzw. quasigeoidę, od której liczone odległości fizycznej powierzchni Ziemi nazywamy wysokościami normalnymi lub potencjalnymi.

Wysokość normalna w stosunku do wysokości ortometrycznej:

0x01 graphic

gdzie: 0x01 graphic
- przeciętna wartość przyśpieszenia normalnego na odcinku B0'B i może być obliczona ze wzoru:

0x01 graphic

Wysokości normalne punktów położonych w jednakowych szerokościach geograficznych są równe, ponieważ obliczone 0x01 graphic
jest funkcją ϕ. 0x01 graphic
to stałe. Quasigeoida niewiele odbiega od geoidy (maksymalnie w górach 1-2 m), obie powierzchnie pokrywają się na obszarach oceanów i mórz oraz wszędzie tam gdzie 0x01 graphic
. Na równinach różnice wynoszą kilka centrymetrów.

Jak wyznaczyć geoidę?

dh - przewyższenie zaobserwowane na łatach

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x01 graphic

gdzie: 0x01 graphic
- są zredukowane wolnopowietrznie na quasigeoide

0x01 graphic
- jest anomalią wolnopowietrzną na quasigeoidzie

0x01 graphic

0x01 graphic
- wysokość punktu bieżącego linii całkowania (ciągu niwelacyjnego)

Wpływ krzywizny powierzchni ekwipotencjalnej na wyniki niwelacji

- w przypadku niwelowania ze środka wpływ ten jest zaniedbywalny, ale w geodezji nie zawsze niweluje się ze środka.

0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x01 graphic
0x01 graphic
dla 0x01 graphic

0x08 graphic
Refrakcja w niwelacji precyzyjnej (Lallemand)

0x08 graphic
0x01 graphic
wg Guy Rusacka

0x01 graphic
- współczynnik rozszerzalności gazu 0x01 graphic

0x01 graphic
- temperatura w oC

0x01 graphic
- gęstość powietrza wyjściowa

0x01 graphic
- ciśnienie

logarytmuje a następnie różniczkuje ⇒

0x01 graphic

największy wpływ na zmianę gęstości powietrza ma temperatura

0x08 graphic
0x01 graphic
zmiana gęstości

I prawo Lallemanda

„Warstwy powietrza o jednakowej gęstości są równoległe nie do powierzchni ekwipotencjalnej, ale do terenu.”

0x08 graphic

II prawo Lallemanda

0x08 graphic

0x08 graphic
0x01 graphic

gdzie: a, b, c - stałe

t - temperatura

w - wysokość od Ziemi

„Jeżeli wysokość nad Ziemią zmienia się w postępie geometrycznym to temperatura zmienia się w postępie arytmetycznym.”

0x08 graphic

0x01 graphic
- łączy wpływ refrakcji na stanowisku

Ep - wpływ refrakcji na łacie w przód

Et - wpływ refrakcji na łacie wstecz

0x01 graphic

gdzie:

s - spadek terenu

tt - temperatura na łacie wstecz

tp - temperatura na łacie w przód

Niwelowanie jednego i tego samego odcinka w kierunku przeciwnych nie usunie z wyniku refrakcji.

Wpływ refrakcji wg T. J. Kukkamäki

0x01 graphic
gdzie: a, b, c - stałe

t - temperatura oC

z - wysokość od Ziemi

R - refrakcja

0x01 graphic

gdzie:

(t2-t1) - różnica temperatur w oC między wysokościami 0,5m a 2,5m na Lacie 3 metrowej

d - długość celowej w m

h - różnica wysokości wyrażona w pół centymetrach (pomierzona)

γ - stała zależna od stałej c, skomplikowany wzór

Refrakcja ma czasem charakter systematyczny, aby zmniejszyć jej wpływ trzeba skracać długość celowych. W nocy refrakcja jest przeciwnego znaku, ale wartości tylko połowy dziennej. Najkorzystniejsze dla pomiaru są okresy: ranny (1 - 1,5 godziny po wschodzie Słońca) i wieczorny (1 godzina przed zachodem Słońca), wtedy refrakcja ma najmniejszy wpływ.

Różnica między niwelatorem precyzyjnym i niwelatorem technicznym.

1. powiększenie lunety: w precyzyjnym jest większe 40x, 50x, w technicznym 20x, 25x.

2. czułość libeli

3. krzyż nitek: w precyzyjnym

0x08 graphic

w technicznym

0x08 graphic

4. układ mechaniczno - optyczny z płytką płaskorównoległą przed obiektywem połączoną z wyskalowaną śrubą mikrometryczną.

Wielkość przesunięcia pionowego można odczytać na bębenku śruby mikrometrycznej. Odczyt w niwelacji precyzyjnej jest sześciocyfrowy: 3 cyfry z łaty i 3 cyfry z bębenka śruby.

Łata do niwelacji precyzyjnej:

0x08 graphic

Znaki w niwelacji precyzyjnej;

W niwelacji precyzyjnej stosuje się znaki wiekowe osadzone na litych skałach (znaki I i II rzędu), są to tzw. „Pola Staussa”. I dalej są jeszcze 4 rodzaje znaków wysokościowych stabilizuje się je głęboko poniżej strefy zamazania gruntu.

Zasady wykonywania niwelacji precyzyjnej:

Niwelację wykonuję się 2 razy „tam” i dwa razy „z powrotem” na kilku klinach i ze zmianą wysokości osi celowej. Kolejność odczytów powinna być t p p t - program „A” (t- wstecz, p- w przód) t1 t2 p1 p2 - program „C”

W niwelacji mierzy się temperaturę i ciśnienie oraz rozmierza się dokładnie ruletką odległość od łat do instrumentu z dokładnością do 1 dm

0x08 graphic

Dawny Zabór Pruski ma punkt zerowy, poziom odniesienia odniesiony do mareografu w Amsterdamie (+11 cm)

W Zaborze Austriackim to średnia z poziomu trzech mórz: Bałtyku, Czarnego i Azowskiego.

Współcześnie poziomem odniesienia jest mareograf w Kronsztadzie.

Błędy niwelacji precyzyjnej:

  1. systematyczne

  2. przypadkowe lub

  1. błędy instrumentalne

-bł. łaty

-bł. instrumentu

  1. błędy na stanowisku

- bł. obserwacji

-bł. Obserwatora

Bł. łaty:

-wygięcie pod wpływem wilgotności

- niepionowość spowodowana wadliwą libelą

-zmiana długości łaty pod wpływem temperatury

Bł. instrumentu:

-niepoziomośc nitki poziomej

-nierównoległość osi libeli do osi celowej

-martwy ruch śruby mikrometrycznej

-błędy samej śruby mikrometrycznej

-błędy samej śruby mikrometrycznej

-nierównoległość ścian płytki płaskorównoległej

-bezwładność pęcherzyka i gęstości cieczy

Bł. obserwatora:

-celowanie i nastawianie na ostrość

-doprowadzanie libeli do poziomu

-odczytywanie na łacie lub nastawianie na środek kreski

-błąd w pomiarze odległości instrument - łata

Bł: inne ( ze strony środowiska):

-grunt w którym mamy wbite kliny i osadzamy legistratywy

    1. osiadanie instrumentu i łaty

    2. wypychanie instrumentu i łat przez sprężysty grunt

    3. wpływ wiatru

    4. wpływ zmian temperatury na instrument i łaty

    5. wibracja czyli drganie powietrza

    6. wpływ refrakcji

Ilość stanowisk powinna być parzysta, aby wyeliminować błąd zera łaty.

0x08 graphic

Wyrównanie sieci niwelacyjnej zazwyczaj metoda pośredniczącą.

Podstawowa sieć pionowa w Polsce.

Początki niwelacji precyzyjnej to druga połowa XIX w.

1867r. Szwajcaria (CH)

1871r. Prusy, Rosja

1874r. Francja

Pierwsze ciągi były wzdłuż linii kolejowej.

1875 - 1894 Austria ok. 3 tys. Punktów z poziomem odniesionym do Triestu.

Były trzy rożne poziomy odniesienia, więc i trzy różne niwelacje po 1918r. Czyli po odzyskaniu niepodległości.

1926 - 1938 - polska sieć niwelacji precyzyjnej z poziomem odniesienia do Amsterdamu, było siedem punktów fundamentalnych 36 oczek, niwelatorem Zeiss-a, z nasadką zawierającą płytkę płaskorównoległą , łaty były inwarowe, ilość linii: 121, średnia odległość między znakami w sieci I rzędu to 1,7 km. Wyrównanie metoda warunkową, długości oczek od czterech do kilkanaście km.

Punkt główny był na ratuszu w Toruniu

1931r. - mareograf w Gdyni (polski);

uzyskiwany błąd to ±1mm/1km podwójnej niwelacji przy błędzie systematycznym ±0,2mm

1945 - 1947r. - inwentaryzacja

do 1952r. - niwelacja o dokładności 1,5mm/1km - 4 400km tej niwelacji z zach i płn

1953 - Konferencja w Sofii odnośnie standardów niwelacji we wszystkich państwach w Europie, zarządziła ona warunki identyczne dla wszystkich:

a - pomiar ze środka z celowa ograniczona do 40m (nie krótsza niż 8m)

b - na czterech klinach ma się odbywać niwelacja

c - obowiązuje podwójna niwelacja tam i z powrotem

d - wysokość celowej nad terenem do 0,5m nie niżej

e - wspólne wyrównanie, duże oczka

f - identyczne typy znaków zalecone do osadzenia

g - system wysokości normalnych i poziom odniesienia „Kronsztad”

Takie podejście korzystne jeśli chodzi o badanie współrzędnych ruchów skorupy ziemskiej.

1953 - 1956r. - robiono w Polsce pomiar sieci niwelacyjnej precyzyjnej uzupełnionej o pomiary niwelacją grawimetryczną. Na terenie kraju realizowano osnowy 1 rzędu dokładność - 1mm/km

1955 - 1959r. - dalsze rzędy osnowy: 2 rzędu (dokładność 2mm/km), 3 rzedu (7mm/km), 4 rzędu (20mm/km).

1971r. - II Konferencja w Sofii - zażądano ponownego pomiaru sieci KDL

1974 - 1979r. - powtórny pomiar 1 rzędu (10 400km linii), około 170 linii o sredniej długości 50 - 60 km, odcinek niwelacji około 1km; stosowano: tppt pttp; średni błąd wyznaczony z odchyłek zamknięć poligonów to ±0,94mm/km; stosowano niwelatory Zeiss Ni 002 i opton Ni 1.

Niwelacja trygonometryczna.

d ± 1mm α ± 3cc

Inne techniki:

ML - Motorized Levelling

MTL - Motorized Trigonometring Levelling

FL - Food Levelling

BL - Bice Levelling

W ciągu sezonu pomiarowego można wykonać 1 400km ML, a 304km FL, 504km BL.

0x01 graphic

0x08 graphic

41

Przekroje: południkowy i równoleżnikowy mają kształt elips. Elipsoida trójosiowa jest dobrym przybliżeniem, jednak jest bardzo skomplikowana matematycznie, więc jest używana tylko w pracach naukowo-badawczych "najwyższej próby".

a

b

c

Spłaszczenie p: 0x01 graphic
, dla Ziemi 0x01 graphic

I mimośród e: 0x01 graphic
, 0x01 graphic

II mimośród e': 0x01 graphic

Aby określić elipsoidę obrotową wystarczy znać dwie z następujących wielkości: a, b, p, e2, e'2

Oficjalnie 0x01 graphic
ale właściwie wartość ta znana jest dokładniej (elipsoida WGS 84 - wyznaczona na podstawie badań satelitarnych)

a

b

a

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
- szerokość geodezyjna 0x01 graphic

- tangens nachylenia stycznej do osi X

- z równania elipsy

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
- szerokość geodezyjna

0x01 graphic
- szerokość geocentryczna

promieniem 0x01 graphic
zakreślam łuk kołowy i odrzutowuję punkt 0x01 graphic
prostopadle do równika

0x01 graphic
- szerokość zredukowana

0x01 graphic

Krzywa przechodzi przez punkt0x01 graphic
i leży na powierzchni, więc można powiedzieć, że powstała przez przecięcie powierzchni z płaszczyzną (a przynajmniej każdy nieskończenie mały element krzywej). Wśród krzywych leżących na powierzchni i przechodzących przez punkt 0x01 graphic
można wyróżnić:

1° - krzywe powstałe przez przecięcie płaszczyznami zawierającymi normalną do powierzchni w tym punkcie (tzw. przekroje normalne)

2° - wszystkie inne krzywe to przekroje skośne

0x01 graphic

Twierdzenie Meusniera: promień krzywizny przekroju skośnego równa się rzutowi na płaszczyznę ściśle styczną tego przekroju promienia krzywizny przekroju normalnego przechodzącego przez styczną do przekroju skośnego.

0x01 graphic
0x01 graphic
- kąt, pod którym się rzuca

Twierdzenie Eulera. Euler zauważył, że dwa przekroje normalne mają ekstremalną krzywiznę. Nazwał je przekrojami minimalnym (0x01 graphic
) i maksymalnym (0x01 graphic
). Są one wzajemnie prostopadłe, a wszystkie inne mają promień krzywizny z tego przedziału. W dodatku każdy z pozostałych przekrojów:

0x01 graphic

0x01 graphic
- kąt między obliczanym przekrojem normalnym 0x01 graphic
a promieniem krzywizny o minimalnej wartości 0x01 graphic
.

0x01 graphic
- promień krzywizny przekroju normalnego

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

Równoleżnik jest przekrojem skośnym (oprócz równika).

Przekrój poprzeczny leży w płaszczyźnie prostopadłej do kartki, a przekrój południkowy w płaszczyźnie kartki.

0x01 graphic
- promień krzywizny przekroju południkowego (jest przekrojem minimalnym)

0x01 graphic
- promień pierwszego wertykału

0x01 graphic

Rzucamy na płaszczyznę równoleżnika:

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic

. 0x01 graphic
- elipsoida lokalna

0x01 graphic
- geoida

0x01 graphic
- elipsoida obrotowa

0x01 graphic
- normalna do elipsoidy lokalnej

0x01 graphic
- normalna do elipsoidy obrotowej

0x01 graphic
- względne odchylenie pionu

0x01 graphic
- bezwzględne odchylenie pionu

Odchylenie pionu w Krakowie: kilka sekund (w ogóle kilkadziesiąt sekund, a w najwyższych górach nawet minuta).

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

OBLICZENIE DŁUGOŚCI ŁUKU DOWOLNEGO PRZEKROJU NORMALNEGO

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
- łuk równoleżnika - łuk kołowy

0x01 graphic

OBLICZENIE DŁUGOŚCI ŁUKU RÓWNOLEŻNIKA

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

Jest to wzór ścisły umożliwiający obliczenie długości łuku południka o dowolnej długości. Ponieważ jednak jest to całka eliptyczna niemożliwa do przedstawienia funkcjami elementarnymi, dlatego obliczamy ją rozwijając funkcję podcałkową w szereg i całkując wyraz po wyrazie. Rozwiązanie jest tym dokładniejsze im więcej wyrazów całkujemy.

0x01 graphic

Dwoistość przekrojów normalnych:

0x01 graphic
- przekrój wprost 0x01 graphic

0x01 graphic
- przekrój odwrotny 0x01 graphic

0x01 graphic
- skośne

0x01 graphic
nie leżą na tym samym południku ani

równoleżniku

N

S

równoleżnik

przekrój normalny

przekrój normalny

równoleżnik

Dowodzi się, że 0x01 graphic

0x01 graphic
- największy rozstęp liniowy

0x01 graphic

0x01 graphic
- średnia azymutów

0x01 graphic
- długość przekroju normalnego

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
-łączymy wzajemnymi przekrojami normalnymi

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

A

B

C

c a

b

B'

A' C'

c a

b

0x01 graphic

*

- additamenty liniowe, które oznaczamy: 0x01 graphic

additamenty boków logarytmiczne

ortodroma

normalna do krzywej

0x01 graphic

równanie parametryczne krzywej na powierzchni

0x01 graphic

normalna do powierzchni

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

Kawałek krzywej 0x01 graphic

0x01 graphic
rzutujemy na płaszczyznę 0x01 graphic
i otrzymujemy 0x01 graphic

0x01 graphic

Element pola:

0x01 graphic

0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic

0x01 graphic

równanie Clairaut

0x01 graphic
- szerokość zredukowana

0x01 graphic
- duża półoś elipsoidy inna postać równania

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
- południki linie

0x01 graphic
- równik geodezyjne

równoleżniki nie spełniają 1., 3. i 4. definicji...

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

szukane:

0x01 graphic
, 0x01 graphic

+ wyrazy małe drugiego rzędu

- wyrazy małe drugiego rzędu

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

Dane: 0x01 graphic

Szukane: 0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
- obliczamy z twierdzenia sinusów w tym trójkącie

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

.

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

ortodroma

równoleżnik

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

obliczone 0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

.

0x01 graphic

0x01 graphic

0x01 graphic
bo jest to mały kąt

0x01 graphic

mały kąt

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

.0x01 graphic

dł. łuku w Peru

dł. łuku w Laponii

2 równania z 2 niewiadomymi: 0x01 graphic

to 0x01 graphic

nowa wsp. izom. ozn. 0x01 graphic

ozn. 0x01 graphic

odwzorowanie Merkatora dla kuli

odwzorowanie Merkatora dla elipsoidy zwane też odwzorowaniem Gaussa

0x01 graphic
- długość pewnego południka środkowego

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

te wzory różnią się wyrazami trzeciego rzędu

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

- redukcja długości

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

szczeciński

bydgoski

18°

warszawski

21°

białostocki

24°

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
- soldnerowska zbieżność południków

0x01 graphic

azymut

geodezyjny

azymut

soldnerowski

G' B' D'

F' A' C'

0x01 graphic

G B D

F A C

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

Wiedeń Bielsko - Biała Lwów

Lwów

WCI 4000 OCII OCIII X

Y

mapy katastralne:

1:2880

w małych miastach 1:1440

1:2000

skala 0,9995

184km

22°E

52°N

koło bez zniekształceń

układ pierwotny (stary)

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

układ wtórny (nowy)

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

6 niewiadomych - 3 punkty wspólne

12 niewiadomych - 4 punkty wspólne

8 niewiadomych - 4 punkty wspólne nie leżące na jednej prostej

układ pierwotny

układ wtórny

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
- łuk od 0x01 graphic

0x01 graphic

Kąty zachowują się między krzywymi na elipsoidzie i kuli, a nas interesują kąty między prostymi na płaszczyźnie.

r

m

dr

M

0x01 graphic

100m

ΔHbiegun

ΔHrównik

45o

Ziemia

P

prosta pionowa (pion)

P

linia pionowa

prosta pionowa jest tylko styczną

do zakrzywionej linii pionowej

0x01 graphic

g=0

R

przyśpieszenie

na Ziemi

w punkcie

pod szerokością

geograficzną

przyśpieszenie

siły

ciężkości

na równiku

spłaszczenie grawimetryczne

na biegunie

0x01 graphic

n równań dla 2 niewiadomych (ga i β)

ϕ

g

C

B

A

l

C

B

A

klosz próżniowy

t1

t2

t3

t4

S1

S2

pow. Ziemi

R

geoida

g

g0

H

P

H

P'

0x01 graphic

100 - 110 km

powierzchnia izostatyczna (równowaga hydrostatyczna)

0x01 graphic

100 - 110 km

powierzchnia izostatyczna (równowaga hydrostatyczna)

SiAl

SiMa

0x01 graphic

0x01 graphic

y0

x0

x0'

P0 (x0', y0')

γ0'

γ0

m

M

R

Ziemia

F

f

P

ρ

oś obrotu

N0x01 graphic

0x01 graphic

ν2

ν3

θ

A

S

N

90°-B

z'

η

zo

A

θ

90°-ϕ

za

λ-L

B

ξ

Θ'

geoida

elipsoida

N1

N2

normalna do elipsoidy

normalna do geoidy

ds

-10km

P

b'

K

H

HP

HK

Powierzchnia fizyczna

geoida - kula

elipsoida obrotowa

N

P0

pion

K0

K

P

b

b'

pion

P'

K'

Θ1

Θ2

redukcja ze względu na odchylenie pionu

redukcja bazy z geoidy na elipsoidę

redukcja bazy na geoidę

geoida

ϕG

ϕP

G

H

P

linia pionowa lekko zakrzywiona

styczna do zakrzywionej linii pionu

P

A

A'

Ns

NP

P'

P”

H

S

B””

normalna w pkt. P

linia pionu

(zakrzywiona)

prosta styczna

geoida

pow. ekwipotencjalna

H

α'

α

B

P'

P”

P

P'L

P'K

K'

S'

K

PL

PK

S

Kąt poziomy (wszystkie punkty nad geoidą)

Rzut na geoidę

W=Wi

W=WB

B

g

dh

W=Wi+dW

W=WA

W=W0

potencjał

A

poziom geoidy

zaobserwowane przewyższenie

poprawka dynamiczna

do wielkości pomierzonych

W=WB

A'

B

dh

W=Wi+dW

g

HBortom

W=Wi

B'

A

W=WA

A0

geoida

B0

W=WA

W=W0

0

quasigeoida

W=WB

B

A

B'

B'0

geoida

A'

B''0

C

D

B

E

A

HAnorm

HBnorm

HCnorm

HDnorm

HEnorm

normalne

quasigeoida

geoida

poprawka normalna

d

warstwy gęstości pow.

równoległe do terenu

pow.

ekwipotencjalna

2

1

łata wstecz

łata w przód

Et

Ep

służy do naprowadzania kreski podziału miedzy te dwa wąsy

kreska podziału łaty

taśma inwarowa z naniesionym podziałem

5mm - kreska

5mm - odstęp

3m

Zakres mikrometru pozwala na odczyt tylko jednej kreski podziału, podział może być główny i kontrolny przysunięty o pewną stałą wielkość.

d1≈d2 ± 0,1m

d < 35m

d1

d2

nogi statywu na stanowiskach

2-3m



Wyszukiwarka