07.07.2014 Views

Byte av referenssystem i plan i Oskarshamns kommun - Kungliga ...

Byte av referenssystem i plan i Oskarshamns kommun - Kungliga ...

Byte av referenssystem i plan i Oskarshamns kommun - Kungliga ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>Byte</strong> <strong>av</strong> <strong>referenssystem</strong> i<br />

<strong>plan</strong> i <strong>Oskarshamns</strong><br />

<strong>kommun</strong><br />

Henrik Nilsson<br />

Examensarbete i geodesi nr. 3112<br />

TRITA-GIT EX 09-06<br />

Avdelningen för Geodesi<br />

<strong>Kungliga</strong> Tekniska Högskolan (KTH)<br />

100 44 Stockholm<br />

Maj 2009


Förord<br />

När jag påbörjade en anställning som mätningsingenjör hos <strong>Oskarshamns</strong> <strong>kommun</strong> 2006,<br />

hade <strong>kommun</strong>en precis tagit beslut om att genomföra byte <strong>av</strong> <strong>referenssystem</strong> för <strong>kommun</strong>ens<br />

lägesbundna information. För mig blev detta <strong>referenssystem</strong>sbyte ett mycket spännande och<br />

intressant examensarbete, där jag i högsta grad blev involverad i såväl mätningarna i fält som<br />

det under hela projekttiden fortlöpande analysarbetet.<br />

Vid det här laget är det ett par år sedan <strong>referenssystem</strong>sbytet i <strong>Oskarshamns</strong> <strong>kommun</strong><br />

genomfördes. Rapportskrivandet har <strong>av</strong> olika anledningar dragit ut på tiden. Det positiva är<br />

dock att jag nu med lite mera perspektiv kan se tillbaka på hur det var före övergången till<br />

SWEREF 99, men också på hur det blev och är efter. Min förhoppning är att beskrivningen i<br />

denna rapport <strong>av</strong> hur <strong>referenssystem</strong>sbytet gjordes, kan vara till nytta för de <strong>kommun</strong>er som<br />

<strong>av</strong> olika anledningar ännu inte kommit igång med motsvarande arbete.<br />

Ett stort tack ska riktas till mina handledare Dr Huaan Fan på KTH i Stockholm och Leif<br />

Carlsson, chef på Kart- och <strong>plan</strong><strong>av</strong>delningen i Oskarshamn, som båda på olika sätt styrt mig i<br />

rätt riktning när det gäller såväl teoretiska som praktiska funderingar kring det analys- och<br />

mätningsarbete som examensarbetet omfattade. Tack säger jag också till Bengt Andersson,<br />

geodet vid Lantmäteriets geodesienhet i Gävle, som under hela <strong>referenssystem</strong>sbytets<br />

genomförande villigt ställde upp och delade med sig <strong>av</strong> kunskap, tips och råd.<br />

Maj 2009, Silverdalen<br />

Henrik Nilsson<br />

I


Sammanfattning<br />

Flera olika lokala koordinatsystem hade använts i <strong>Oskarshamns</strong> <strong>kommun</strong> under 1900-talet.<br />

Fram till år 2000 hade dock flera mindre system kunnat anslutas till två större, vilket var bra<br />

men inte tillräckligt. De två större koordinatsystemen hade rikets trianguleringsnät som grund,<br />

men var behäftade med geometriska deformationer som behövde rätas upp. I samband med en<br />

sådan geometrisk renovering var det även naturligt och lätt att se fördelarna med att gå över<br />

till SWEREF 99 – ett <strong>referenssystem</strong> som är anpassat till global satellitn<strong>av</strong>igeringsteknik.<br />

Beslut om en övergång togs och arbetet påbörjades.<br />

Eftersom arbetet med ett <strong>referenssystem</strong>sbyte i en <strong>kommun</strong> sker i samråd med Lantmäteriet,<br />

användes de metoder som Lantmäteriet tagit fram och ger rekommendationer och råd kring.<br />

Rapporten beskriver och ger ett verkligt exempel på hur metoderna helt konkret kan fungera. I<br />

en teoridel beskrivs Lantmäteriets Rix 95-projekt, möjligheten att med hjälp <strong>av</strong> en<br />

direktprojektionsmetod projicera ett lokalt system direkt till ett geodetiskt globalt system,<br />

samt begreppet restfel och olika sätt att redovisa och presentera restfel i vektorkartor och<br />

variationsbilder.<br />

Utgående från resultatet <strong>av</strong> Rix 95-projektet påbörjades en restfelsanalys som blev ett iterativt<br />

förfarande enligt modellen: analys – resultat – mätningsåtgärder – ny analys – nytt resultat –<br />

nya mätningsåtgärder – slutlig analys – slutligt resultat. Mättekniken som användes var<br />

nätverks-RTK. Hjälpmedel i analysarbetet var restfelsvektorkartor och restfelsvariationsbilder<br />

tillsammans med historisk kunskap om och erfarenhet <strong>av</strong> de två <strong>kommun</strong>ala stomnäten.<br />

Några olika tester och kontroller <strong>av</strong> de slutliga transformationssambanden gjordes innan<br />

primärkartan transformerades. Dessa tester inkluderade bland annat kontrollmätningar på ett<br />

antal kartdetaljer.<br />

Det kunde konstateras att transformationsparametrarna tillsammans med restfelsmodellen<br />

hade överfört kartdetaljerna på ett tillförlitligt sätt så att primärkartan uppnått en god kvalitet<br />

som även var fullt jämförbar med andra <strong>kommun</strong>ers.<br />

III


Abstract<br />

A number of different local co-ordinate systems had been used in Oskarshamn municipality<br />

during the 20 th century. However, before the year of 2000 it had been possible to connect a<br />

number of smaller systems to two bigger ones, which was good but not enough. The two<br />

bigger co-ordinate systems were dialects of national systems but only weakly connected to<br />

them. Furthermore, because of the way that the systems had been developed, the two local coordinate<br />

systems were distorted and the distortions needed to be reduced. As the geometry of<br />

the co-ordinate system was renovated, the advantages of a changeover to SWEREF 99 – a<br />

reference system that is adapted to global satellite n<strong>av</strong>igation technology – were natural and<br />

could easily be seen. A decision on the changeover was made and the project could begin.<br />

Since the work on the changeover of reference system in a municipality is done after<br />

consultation with Lantmäteriet (National Land Survey of Sweden), methods prepared and<br />

recommended by Lantmäteriet were used. The thesis describes and gives a real-life example<br />

of how the methods work quite concretely. In a theoretical part the Rix 95 project of<br />

Lantmäteriet is described together with the possibility to project a local system directly to a<br />

global one, by means of the direct projection method, the concept of residual and different<br />

means of presenting residuals in vector maps and variation images.<br />

With the result of the Rix 95 project as a starting point, an analysis of the residuals started<br />

according to an iterative procedure: analysis – result – surveys – new analysis – new result –<br />

new surveys – final analysis – final result. The surveying technique used was network RTK<br />

survey. Aids used for the analysis were residual vector maps and variations of deformations<br />

images or thematic maps, together with historic knowledge about and experience of the two<br />

local geodetic control networks.<br />

Some tests and checks were done of the final transformation parameters before the municipal<br />

large scale map was transformed. These tests included among other things surveys of a<br />

number of control points.<br />

It could be established that the transformation parameters together with the correction model<br />

had transferred the map details in a reliable way so that the municipal large scale map had<br />

obtained a good quality completely comparable with those in other municipalities.<br />

IV


Innehållsförteckning<br />

FÖRORD.................................................................................................................................................................I<br />

SAMMANFATTNING....................................................................................................................................... III<br />

ABSTRACT .........................................................................................................................................................IV<br />

INNEHÅLLSFÖRTECKNING .......................................................................................................................... V<br />

1. GEODETISKT UTGÅNGSLÄGE OCH FÖRBEREDANDE ÅTGÄRDER......................................... 1<br />

1.1 BAKGRUND ........................................................................................................................................... 1<br />

1.2 PROBLEM MED TVÅ KOORDINATSYSTEM............................................................................................... 3<br />

1.3 VARFÖR GÅ ÖVER TILL SWEREF 99?................................................................................................... 3<br />

1.4 VARFÖR BYTA REFERENSSYSTEM I PLAN MEN INTE I HÖJD? .................................................................. 3<br />

1.5 MANDAT ATT GENOMFÖRA ÖVERGÅNGEN ............................................................................................ 4<br />

2. TEORI .......................................................................................................................................................... 5<br />

2.1 RIX 95................................................................................................................................................... 5<br />

2.2 DIREKTPROJEKTION .............................................................................................................................. 6<br />

2.2.1 Bakgrund ......................................................................................................................................... 6<br />

2.2.2 Direktprojektionsmetoden................................................................................................................ 6<br />

2.3 RESTFEL.............................................................................................................................................. 11<br />

2.3.1 Restfelsmodell................................................................................................................................ 12<br />

2.3.2 Redovisning och presentation <strong>av</strong> restfel ........................................................................................ 12<br />

3. RESTFELSANALYS OCH TRANSFORMATION ............................................................................... 15<br />

3.1 RESULTATET AV RIX 95-MÄTNINGARNA I SMÅLAND OCH OSKARSHAMN........................................... 15<br />

3.2 PLANERING AV KOMPLETTERINGSMÄTNINGAR ................................................................................... 17<br />

3.3 VAL AV MÄTMETOD OCH MÄTTEKNIK................................................................................................. 18<br />

3.3.1 Statisk mätning .............................................................................................................................. 18<br />

3.3.2 RTK-mätning ................................................................................................................................. 19<br />

3.3.3 Förutsättningarna i Oskarshamn .................................................................................................. 21<br />

3.4 PERSONAL, UTRUSTNING OCH UPPLÄGG.............................................................................................. 21<br />

3.5 MÄTNINGARNAS FÖRSTA ETAPP.......................................................................................................... 22<br />

3.6 FORTSATT ANALYS.............................................................................................................................. 23<br />

3.6.1 Regionsystemet .............................................................................................................................. 25<br />

3.6.2 <strong>Oskarshamns</strong> lokala koordinatsystem ........................................................................................... 29<br />

3.7 PLANERING AV ANDRA MÄTETAPPENS MÄTNINGAR ............................................................................ 32<br />

3.8 MÄTNINGARNAS ANDRA ETAPP........................................................................................................... 32<br />

3.9 ANALYSENS SLUTFAS.......................................................................................................................... 32<br />

3.9.1 Regionsystemet .............................................................................................................................. 35<br />

3.9.2 <strong>Oskarshamns</strong> lokala koordinatsystem ........................................................................................... 38<br />

3.10 UPPRÄTTANDE AV RESTFELSMODELLER.............................................................................................. 39<br />

3.11 FRAMTAGANDE AV TRANSFORMATIONSSAMBAND.............................................................................. 44<br />

3.12 TEST OCH KONTROLL AV TRANSFORMATIONSSAMBAND ..................................................................... 44<br />

3.13 TRANSFORMATION AV OSKARSHAMNS KOMMUNS PRIMÄRKARTA...................................................... 49<br />

4. DISKUSSION OCH SLUTSATSER........................................................................................................ 50<br />

4.1 VALET AV MÄTTEKNIK........................................................................................................................ 50<br />

4.2 ANALYSARBETE OCH HISTORISK STOMNÄTSKUNSKAP ........................................................................ 50<br />

4.3 HJÄLPMEDEL I ANALYSPROCESSEN ..................................................................................................... 51<br />

4.4 TRANSFORMATIONENS KVALITET ....................................................................................................... 52<br />

LITTERATURFÖRTECKNING....................................................................................................................... 54<br />

V


1. Geodetiskt utgångsläge och förberedande åtgärder<br />

1.1 Bakgrund<br />

Liksom i många andra <strong>kommun</strong>er hade lägesbunden information i <strong>Oskarshamns</strong> <strong>kommun</strong> i<br />

Kalmar län kommit att hanteras i flera olika lokala koordinatsystem parallellt. Genom att<br />

ansluta mindre lokala koordinatsystem till större hade dock antalet koordinatsystem fram till<br />

2000-talet kunnat reduceras till två. Detta var visserligen ett steg i rätt riktning och en<br />

förutsättning för en rationell användning <strong>av</strong> <strong>kommun</strong>ens stomnät men det var inte tillräckligt.<br />

RT R06<br />

2,5 gon V 0:-15<br />

RT 38 66<br />

2,5 gon V<br />

Figur 1. Ungefärliga tillämpningsområden för de två lokala koordinatsystemen<br />

De två lokala koordinatsystem (se figur 1) som använts utgjorde båda förtätningar <strong>av</strong> gamla<br />

regionala eller nationella system. <strong>Oskarshamns</strong> stad tillsammans med närmaste tätort söderut,<br />

Påskall<strong>av</strong>ik, täcktes in <strong>av</strong> ett koordinatsystem som grundade sig på RT 38 – alltså ett<br />

<strong>referenssystem</strong> som bygger på vinkelmätningar i trianglar i ett kransformigt nät, och som<br />

mättes och beräknades i samband med den andra rikstrianguleringen under första halvan <strong>av</strong><br />

1900-talet (Lantmäteriet, geodesi). Systemet förtätades i <strong>Oskarshamns</strong> <strong>kommun</strong> och räknades<br />

om 1966 och betecknades RT 38 66 2,5 gon V. Det andra lokala koordinatsystemet hade RT<br />

1


90 eller, mer preciserat, regionsystemet RT R06 (Småland) som grund, och utgjorde alltså en<br />

förtätning <strong>av</strong> mätningar gjorda under den tredje rikstrianguleringen mellan 1967 och 1982.<br />

Systemet RT R06 2,5 gon V 0:-15 täckte in resterande del <strong>av</strong> <strong>Oskarshamns</strong> <strong>kommun</strong> och<br />

inkluderade därmed de övriga omkring åtta tätorterna där primärkarta finns. De <strong>kommun</strong>ala<br />

stomnäten, liksom primärkartan, täckte endast in de tätbebyggda områdena (se<br />

stompunktskartan i figur 2), vilket innebär att den omkringliggande landsbygden fortfarande<br />

helt saknar <strong>kommun</strong>alt stomnät. Skillnaden mellan de lokala koordinatsystemen RT 38 66 och<br />

RT R06 var omkring 1,5 meter.<br />

Figur 2. Stompunktskarta för <strong>Oskarshamns</strong> <strong>kommun</strong>. Stomnäten täckte in <strong>Oskarshamns</strong> stad och de tätorter som<br />

omfattades <strong>av</strong> <strong>kommun</strong>ens primärkarta.<br />

2


Av praktiska skäl kommer här i rapporten fortsättningsvis benämningen ”<strong>Oskarshamns</strong> lokala<br />

koordinatsystem” att användas för koordinatsystemet RT 38 66 2,5 gon V, medan<br />

benämningen ”regionsystemet” kommer att syfta på koordinatsystemet som hade<br />

regionsystemet RT R06 som grund.<br />

1.2 Problem med två koordinatsystem<br />

Situationen med användning <strong>av</strong> två koordinatsystem parallellt var egentligen relativt<br />

hanterbar under lång tid, eftersom de två systemen var helt åtskilda och alltså inte överlappade<br />

varandra, och eftersom GNSS-tekniken ännu inte hade börjat användas. I fält ställde det<br />

faktum att det fanns två koordinatsystem till problem endast i de fall mätningsarbetet innebar<br />

att gränsen mellan systemen korsades, vilket inte hände särskilt ofta. Däremot var det<br />

problematiskt att hantera och skilja de olika koordinatsystemen åt i arbetet med <strong>kommun</strong>ens<br />

kartdatabas och GIS-programvaror. Det var också en begränsning att inte kunna lagra<br />

<strong>kommun</strong>täckande information i databasen.<br />

När GNSS-utrustning började användas var det <strong>av</strong> helt <strong>av</strong>görande betydelse att rätt<br />

koordinatsystem valdes för varje enskilt mätjobb. Här var det lätt att <strong>av</strong> misstag trycka fel i<br />

exempelvis fältdatorerna, vilket i sin tur innebar att risken fanns att fel koordinatsystem<br />

valdes och felaktiga projektionsparametrar därmed lästes in och användes. Tillsammans med<br />

vetskapen om den otillräckliga noggrannheten i Rix 95-parametrarna, som var de parametrar<br />

som användes för att direktprojicera de geodetiska koordinaterna i SWEREF 99 ner till de<br />

lokala koordinatsystemen, gjorde risken för sammanblandning <strong>av</strong> koordinatsystemen att en<br />

övergång till ett enhetligt <strong>referenssystem</strong> blev mer än välbehövlig. Rix 95-projektet och<br />

direktprojektionsmetoden beskrivs utförligare i <strong>av</strong>snittet 2.<br />

1.3 Varför gå över till SWEREF 99?<br />

Svårigheterna med att hantera två olika koordinatsystem i såväl fält som i kartdatabas- och<br />

GIS-sammanhang har redan nämnts, och de utgjorde naturligtvis en <strong>av</strong> orsakerna till att det<br />

var viktigt att byta till ett enhetligt <strong>referenssystem</strong>. Andra orsaker som direkt kunde ses, var<br />

att ny mätningsteknik krävde bättre kvalitet i stomnäten än vad de kunde ge, samt att<br />

dataöverföring mellan olika aktörer behövde underlättas. De möjligheter som GNSS-tekniken<br />

erbjuder kunde inte tillvaratas fullt ut med den bristande geometriska kvalitet som de gamla<br />

lokala stomnäten uppvisade. En kartläggning <strong>av</strong> de befintliga deformationerna i stomnäten var<br />

därför nödvändig, och i samband med denna geometriska renovering var det naturligt att gå<br />

över till ett <strong>referenssystem</strong> som är anpassat till global satellitn<strong>av</strong>igeringsteknik – svensk<br />

realisering: SWEREF 99. Att Lantmäteriet också <strong>plan</strong>erade för <strong>referenssystem</strong>sbyte sågs som<br />

ytterligare ett starkt skäl till att genomföra övergången.<br />

1.4 Varför byta <strong>referenssystem</strong> i <strong>plan</strong> men inte i höjd?<br />

Höjdsystemet som användes i <strong>Oskarshamns</strong> <strong>kommun</strong> var en dialekt <strong>av</strong> RH 70, vilken<br />

räknades om 1981 och därför fick benämningen RH 70/81. Systemet i höjd var liksom de<br />

<strong>plan</strong>a koordinatsystemen deformerat och behövde därför, även det, analyseras och därefter<br />

rätas upp. En övergång till det nya nationella höjdsystemet RH 2000 hade i samband med<br />

denna geometriska renovering varit lämplig. Det beslutades dock att skjuta höjdsystemsbytet<br />

3


på framtiden, för att på så sätt säkerställa att <strong>referenssystem</strong>sbytet i <strong>plan</strong> verkligen kunde<br />

genomföras. Det blev alltså till största delen en fråga om resurser men det fanns också en viss<br />

rädsla för den uppenbara risken att med differenser på ibland bara några centimeter blanda<br />

ihop eller förväxla höjdangivelser i det gamla höjdsystemet med höjdangivelser i det nya<br />

höjdsystemet. I <strong>plan</strong> var risken för förväxling <strong>av</strong> positionsangivelser inte så stor eftersom<br />

skillnaden mellan koordinaterna i de gamla koordinatsystemen och den nya projektionen i<br />

SWEREF 99 var så stor som omkring 140 mil.<br />

1.5 Mandat att genomföra övergången<br />

Redan på ett tidigt stadium ansåg <strong>kommun</strong>ens Kart- och <strong>plan</strong><strong>av</strong>delning att det var angeläget<br />

att söka central förankring för arbetet med att byta <strong>referenssystem</strong> till SWEREF 99, och även<br />

att få ett beslut på politisk nivå om att genomföra övergången. Detta var viktigt eftersom<br />

ingen i det läget hade full överblick över vilka de berörda aktörerna och datasamlingarna var.<br />

En förutsättning för en full satsning på genomförandet var arbetsro utan en massa störande<br />

och tidsödande resonemang och diskussioner med eventuella motarbetare. Visserligen<br />

bedömdes risken vara liten för en utveckling i den riktningen, men det kändes ändå skönt att<br />

låta ärendet få sin behandling först i Samhällsbyggnadsnämnden och därefter<br />

Kommunstyrelsen. Man ska vidare komma ihåg att de geografiska data som skulle<br />

transformeras har ett mycket betydande ekonomiskt värde – den <strong>kommun</strong>ala primärkartan<br />

utgör en databank som har byggts upp under flera decenniers tid – så ett väl förankrat beslut<br />

var en stor trygghet att ha i ryggen redan i projektets inledningsskede.<br />

4


2. Teori<br />

I detta <strong>av</strong>snitt beskrivs några viktiga geodetiska företeelser och begrepp som var centrala i<br />

arbetet och som snabbt kom i fokus.<br />

2.1 Rix 95<br />

Resultatet <strong>av</strong> Rix 95-projektets mätningar blev utgångspunkten för den detaljerade<br />

utvärdering <strong>av</strong> de lokala stomnätens geometrier som en övergång krävde i Oskarshamn. En<br />

allmän beskrivning <strong>av</strong> just det projektet är därför på sin plats och följer här.<br />

Rix 95-projektet startades 1995 och är ”ett nationellt projekt som syftar till att skapa goda<br />

transformationssamband mellan lokala (<strong>kommun</strong>ala) och nationella/globala <strong>referenssystem</strong>,<br />

för att underlätta utbyte <strong>av</strong> geografisk information och rationell användning <strong>av</strong> GNSS-teknik.”<br />

(Lantmäteriet, Rix 95, 2009). Det operativa ansvaret för projektets genomförande har<br />

Lantmäteriverket, och när det gällde själva mätningen var projektet slutfört i maj 2007.<br />

Återstod gjorde då slututjämningar och beräkning <strong>av</strong> transformationssamband (Andersson,<br />

2007).<br />

Genom Rix 95-mätningarna förtätas det <strong>plan</strong>a riksnätet med hjälp <strong>av</strong> GNSS-teknik, från<br />

omkring 3 800 punkter till drygt 9 000 punkter. De nya punkterna är valda så att de är<br />

lämpliga för GNSS-teknik; lättåtkomliga, utmed vägar och i närheten <strong>av</strong> tätorter. Eftersom<br />

syftet med projektet är att bland annat skapa anslutningar till befintliga <strong>kommun</strong>ala stomnät<br />

och det befintliga nationella höjdnätet, är många <strong>av</strong> de nya punkterna just befintliga<br />

<strong>kommun</strong>ala stompunkter eller höjdfixar i riks<strong>av</strong>vägningen.<br />

För att få anslutning till SWEREF 99 har vissa punkter (så kallade SWEREF-punkter) med 50<br />

kilometers punkt<strong>av</strong>stånd, bestämts relativt omkringliggande SWEPOS 1 -stationer. Genom att<br />

beräkna varje ny SWEREF-punkt utifrån dess sex till åtta närmaste SWEPOS-stationer, har<br />

en fritt utjämnad och jonosfärfri multistationslösning kunnat tas fram, som sedan har inpassats<br />

på SWEPOS-stationerna. SWEPOS-punkterna har därefter hållits fasta när övriga GNSSmätningar<br />

inom det aktuella beräkningsområdet utjämnats. För att få transformationssamband<br />

mellan olika lokala system och SWEREF 99 har slutligen inpassningar beräknats som så<br />

kallade direktprojektioner. Denna typ <strong>av</strong> projektioner har alltså spelat en <strong>av</strong>görande roll i Rix<br />

95-projektet och är därmed viktiga även i den här framställningen. I <strong>Oskarshamns</strong> <strong>kommun</strong>s<br />

<strong>referenssystem</strong>sbyte, liksom i andra <strong>kommun</strong>ers <strong>referenssystem</strong>sbyten, har nämligen Rix 95-<br />

sambanden varit helt centrala komponenter. Nästa stycke innehåller en mer utförlig<br />

beskrivning <strong>av</strong> direktprojektionsmetoden, men kort kan man säga att i direktprojiceringen<br />

beräknas parametrar för en Gauss-Krüger-projektion (transversell Mercator-projektion)<br />

mellan latitud och longitud i SWEREF 99 och <strong>plan</strong>a koordinater i det lokala systemet. I de<br />

fall det lokala systemet är felorienterat, kombineras direktprojektionen med en två- eller<br />

1 SWEPOS är ett nätverk <strong>av</strong> fasta referensstationer för GPS och GLONASS och möjliggör mätning med<br />

nätverks-RTK-teknik i Sverige. Fullt utbyggt kommer nätverket att bestå <strong>av</strong> omkring 120 referensstationer och<br />

täcka hela landet sånär som på Norrlands inland. Genom en tjänst baserad på nätverket distribuerar Lantmäteriet<br />

korrektionsdata som gör det möjligt att mäta i SWEREF 99 med centimeternoggrannhet. Nätverks-RTK och<br />

SWEREF 99 ger tillsammans de förutsättningar som krävs för att kunna redovisa positioner i lokala <strong>kommun</strong>ala<br />

koordinatsystem, utan att använda <strong>kommun</strong>ala stompunkter.<br />

5


tredimensionell Helmerttransformation som tar upp den aktuella vridningen (Engberg & Lilje,<br />

2006).<br />

2.2 Direktprojektion<br />

2.2.1 Bakgrund<br />

Transformationen från de gamla lokala koordinatsystemen över till det nya globala<br />

<strong>referenssystem</strong>et och den för <strong>Oskarshamns</strong> <strong>kommun</strong> aktuella kartprojektionen SWEREF 99<br />

16 30, innebar transformation från koordinatsystem med ursprung i rikets nationella<br />

<strong>referenssystem</strong> (RT 38 respektive RT 90), till ett nytt globalt anpassat <strong>referenssystem</strong> med en<br />

ny referensellipsoid – transformation från positioner på Bessels ellipsoid (Bessel 1841) som<br />

projicerats med Gauss-Krügers projektion, till positioner på ellipsoiden GRS 80.<br />

Transformationer <strong>av</strong> det slaget leder inte sällan till omständliga beräkningar med åtskilliga<br />

transformationssteg. Följande transformationskedja är typisk:<br />

( x, y) ( x,<br />

y) ↔ ( ϕ , λ) nationellt<br />

↔ ( X , Y,<br />

Z ) nationellt<br />

↔ ( X , Y,<br />

Z ) globalt<br />

↔ ( ϕ,<br />

λ) globalt<br />

lokalt<br />

↔ (1)<br />

nationellt<br />

Transformationer mellan två <strong>referenssystem</strong> utförs traditionellt med tredimensionell<br />

likformighetstransformation; vanligen benämnd sjuparametertransformation eller<br />

tredimensionell Helmerttransformation. Denna transformationsmetod definieras <strong>av</strong>:<br />

• tre translationer längs respektive axel: ∆ X , ∆ Y , ∆ Z<br />

• tre rotationer runt respektive koordinataxel: ω<br />

X<br />

, ω<br />

Y<br />

, ω<br />

Z<br />

• en skalförändring , som vanligen uttrycks som en skalkorrektion i ppm: δ<br />

En begränsning med sjuparametertransformationen är att även om transformationen endast<br />

omfattar <strong>plan</strong>a koordinater krävs ändå tillgång till höjdinformation. I transformationskedjan<br />

(1) kräver steget mellan ( ϕ, λ)<br />

och ( X , Y,<br />

Z ) tillgång till höjd över ellipsoiden. De<br />

transformerade geodetiska koordinaterna är alltså beroende <strong>av</strong> höjdinformation, vilken ofta<br />

inte finns när man i utgångsläget haft att göra med lägesangivelser i lokala <strong>kommun</strong>ala<br />

koordinatsystem. En möjlighet är då att sätta höjdvärdet lika med noll och använda<br />

sjuparametertransformationen i alla fall. Detta är en tänkbar lösning såvida det är möjligt att<br />

beräkna geodetiska koordinater (geodetisk latitud och longitud) från de lokala <strong>plan</strong>a<br />

koordinaterna (Engberg & Lilje, 2006).<br />

Liksom i många andra svenska <strong>kommun</strong>er saknade dock de lokala koordinatsystemen i<br />

<strong>Oskarshamns</strong> <strong>kommun</strong> rigorösa geodetiska definitioner. Visserligen fanns, som redan nämnts,<br />

släktskapet med rikets nationella <strong>referenssystem</strong> men lokala förtätningar, utvidgningar och<br />

sammanslagningar <strong>av</strong> olika stomnätsdelar hade med tiden gett upphov till sådana spänningar i<br />

stomnäten att beräkning <strong>av</strong> latitud och longitud från <strong>plan</strong>a koordinater i de lokala systemen<br />

inte kunde ge några tillförlitliga resultat. Därmed gick det heller inte att använda<br />

sjuparametertransformationen med det strikta matematiska tillvägagångssättet som<br />

exempelvis beräkningskedjan (1) fordrar.<br />

2.2.2 Direktprojektionsmetoden<br />

Med hjälp <strong>av</strong> en metod som utvecklats <strong>av</strong> B-G Reit på Lantmäteriets Geodesienhet, kan man i<br />

de flesta fall förkorta transformationskedjan i (1) till det mycket enkla uttrycket<br />

6


( x y) ( ϕ, λ) globalt<br />

, ↔ (2)<br />

lokalt<br />

Metoden kallas direktprojektionsmetoden och baseras på ett antagande att man ”givet ett<br />

geodetiskt datum A och ett <strong>plan</strong>t rektangulärt koordinatsystem <strong>av</strong> ett annat datum B, kan hitta<br />

en uppsättning projektionsparametrar (samma projektion som använts för de givna <strong>plan</strong>a<br />

koordinaterna med datum B används även här) för att definiera ett <strong>plan</strong>t koordinatsystem med<br />

datum A, som approximerar det <strong>plan</strong>a koordinatsystemet med datum B” (Reit, 1997, citatet är<br />

fritt översatt från engelska). Det återgivna antagandet beskriver den omvända riktningen i (2);<br />

det vill säga transformation från det tredimensionella globala <strong>referenssystem</strong>et till det <strong>plan</strong>a<br />

lokala koordinatsystemet. Transformationskedjan i uttrycket (3) förkortas alltså genom<br />

direktprojicering till uttrycket (4).<br />

( , λ) globalt<br />

↔ ( X , Y,<br />

Z ) globalt<br />

↔ ( X , Y,<br />

Z ) nationellt<br />

↔ ( ϕ,<br />

λ) nationellt<br />

↔ ( x,<br />

y) nationellt<br />

↔ ( x,<br />

y) lokalt<br />

ϕ (3)<br />

och ( , λ) globalt<br />

↔ ( x, y) lokalt<br />

ϕ (4)<br />

Direktprojektionsmetoden är inte någon matematiskt stringent metod och kan därför inte<br />

definiera några <strong>plan</strong>a system som matematiskt sett överensstämmer med de <strong>plan</strong>a system som<br />

de approximerar. Trots detta finns det tillämpningar där metoden är användbar eftersom de<br />

bristande överensstämmelserna i vissa fall kan accepteras (Engberg & Lilje, 2006).<br />

Grundtanken i den direktprojektionsmetod som Lantmäteriet (Reit, 2003) utformat är att det<br />

geodetiska globala systemet ska projiceras direkt till det lokala systemet. Detta åstadkoms<br />

genom att med hjälp <strong>av</strong> minsta-kvadrat-metoden skatta de fyra parametrarna λ<br />

0<br />

, k<br />

0<br />

, F<br />

N<br />

och<br />

FE<br />

i en inpassning baserad på en transversell Mercator-projektion med Gauss-Krügers<br />

formler. Här följer först en detaljerad beskrivning <strong>av</strong> denna Gauss konforma projektion (Reit,<br />

2009).<br />

Symboler och definitioner<br />

a<br />

f<br />

2<br />

e<br />

ϕ<br />

λ<br />

x<br />

y<br />

λ<br />

0<br />

k<br />

0<br />

ellipsoidens halva storaxel (semi-major axis)<br />

ellipsoidens <strong>av</strong>plattning (flattening)<br />

excentricitetskvadraten (first eccentricity squared)<br />

geodetisk latitud, positiv mot norr<br />

geodetisk longitud, positiv mot öster<br />

<strong>plan</strong> koordinat, positiv mot norr (grid coordinate)<br />

<strong>plan</strong> koordinat, positiv mot öster (grid coordinate)<br />

medelmeridianens longitud (longitude of the central meridian)<br />

skalfaktor på medelmeridianen (scale factor along the central meridian)<br />

∂ λ differens λ − λ0<br />

x-tillägg (false northing)<br />

x<br />

0<br />

y<br />

0<br />

y-tillägg (false easting)<br />

Alla vinklar (latitud, longitud o.s.v.) ska vara uttryckta i radianer. Observera att x-<br />

axeln pekar mot norr och y-axeln mot öster.<br />

7


Ur ellipsoidparametrarna a och f beräknas följande storheter<br />

e 2<br />

= f (2 − f )<br />

f<br />

n =<br />

(2 − f)<br />

a ⎛ 1<br />

â = ⎜1+<br />

n<br />

(1 + n) ⎝ 4<br />

2<br />

+<br />

1<br />

n<br />

64<br />

4<br />

+ ...<br />

Den konforma 2 latituden ϕ * beräknas genom<br />

⎞<br />

⎟<br />

⎠<br />

2<br />

4<br />

6<br />

( + Bsin<br />

ϕ + C sin ϕ + D sin + ...)<br />

ϕ*<br />

= ϕ − sinϕ<br />

cosϕ<br />

A ϕ<br />

(5)<br />

där koefficienterna A , B , C , och D erhålls ur formlerna:<br />

2<br />

A = e<br />

1<br />

B =<br />

6<br />

1<br />

C =<br />

120<br />

D =<br />

1260<br />

4 6<br />

( 5e − e )<br />

6 8<br />

( 104e − 45e + ...)<br />

1<br />

8<br />

( 1237e + ...)<br />

Storheterna ξ ′ och η′ definieras som<br />

ξ ′ = arctan(tanϕ<br />

* / cos( λ − λ0<br />

)) (6a)<br />

η′ = arctan h (cosϕ<br />

*sin( λ − λ0<br />

)) (6b)<br />

Då erhålls<br />

( ξ′ + β1<br />

sin 2ξ′<br />

cosh 2η′<br />

+ β2<br />

sin 4ξ′<br />

cosh 4η′<br />

+ β3<br />

sin 6ξ′<br />

cosh 6η′<br />

+ β 4 sin 8ξ′<br />

cosh 8η′<br />

+ ) 0<br />

( η′ + β1<br />

cos 2ξ′<br />

sinh 2η′<br />

+ β2<br />

cos 4ξ′<br />

sinh 4η′<br />

+ β3<br />

cos 6ξ′<br />

sinh 6η′<br />

+ β 4 cos8ξ′<br />

sinh 8η′<br />

+ ) 0<br />

x = k<br />

K +<br />

0â<br />

x<br />

y = k<br />

K +<br />

0â<br />

y<br />

(7a)<br />

(7b)<br />

β<br />

där koefficienterna 1, β2<br />

, β3<br />

och β 4<br />

beräknas ur<br />

1 2 2 5 3 41 4<br />

β1 = n − n + n + n + ...<br />

2 3 16 180<br />

13 2 3 3 557 4<br />

β2 = n − n + n + ...<br />

48 5 1440<br />

61 3 103 4<br />

β3 = n − n + ...<br />

240 140<br />

2 Äldre svensk litteratur benämner denna kvantitet isometrisk latitud. Idag används termen isometrisk latitud för<br />

storheten ψ = ln{ tan( π /4 + ϕ /2)[(1 − e sin ϕ)/(1<br />

+ e sin ϕ)]<br />

} . Den isometriska latituden beräknas ur den<br />

konforma latituden enligt formeln ψ = ln tan( π / 4 + ϕ * /2)<br />

. Jämför. John P. Snyder: Map Projections - A<br />

Working Manual, U.S. Geological Survey Professional Paper 1395.<br />

8<br />

e /2


49561<br />

β4 = n<br />

161280<br />

4<br />

+ ...<br />

Till inpassningen, som i det följande beskrivs med formler och redogörelse från (Reit, 2009),<br />

är sedan indata en uppsättning punkter (passpunkter) vars koordinater är kända i såväl det<br />

geodetiska som det lokala systemet. Koordinatparen latitud och longitud i geodetiska systemet<br />

λ<br />

x, y måste alltså finnas tillgängliga.<br />

( ϕ, ) G<br />

respektive <strong>plan</strong>a koordinater i lokala systemet ( ) L<br />

För att utföra en transversell Mercator-projektion behövs även värden för den använda<br />

ellipsoidens storaxel ( a ) och <strong>av</strong>plattning ( f ), vilka hämtas från det geodetiska systemet<br />

( ϕ, λ) G<br />

. Vidare behövs värden för longituden för medelmeridianen ( λ<br />

0<br />

), skalan på<br />

medelmeridianen ( k 0<br />

) och tilläggen i x- respektive y-led ( x<br />

0<br />

respektive y<br />

0<br />

).<br />

Nu betraktas x och y som funktioner <strong>av</strong> projektionsparametrarna på följande sätt:<br />

x = x λ , k , x , ) och y = y λ , k , x , ) . Taylorutveckling runt närmevärdena<br />

(<br />

0 0 0<br />

y0<br />

(<br />

0 0 0<br />

y0<br />

( λ<br />

0<br />

),( k 0<br />

),( x<br />

0<br />

) och ( y<br />

0<br />

) ger observationsekvationerna<br />

∂x<br />

∂x<br />

∂x<br />

∂x<br />

x + v x = x(( λ 0 ), (k 0 ), (x 0 ), (y 0 )) + ( ) 0 ∆λ 0 + ( ) 0 ∆k<br />

0 + ( ) 0 ∆x<br />

0 + ( ) 0 ∆y<br />

0 (8a)<br />

∂λ ∂k<br />

∂x<br />

∂y<br />

0<br />

∂y<br />

∂y<br />

∂y<br />

∂y<br />

+ v y = y(( λ 0 ), (k 0 ), (x 0 ), (y0<br />

)) + ( ) 0 ∆λ 0 + ( ) 0 ∆k<br />

0 + ( ) 0 ∆x<br />

0 + ( ) 0 ∆y<br />

(8b)<br />

∂λ ∂k<br />

∂x<br />

∂y<br />

y 0<br />

0<br />

0<br />

0<br />

0<br />

där ∆ λ 0 , ∆ k 0 , ∆x<br />

0 och ∆ y 0 är obekanta korrektioner till närmevärdena samt v<br />

x<br />

och<br />

v<br />

y<br />

är förbättringar till de observerade (kända) värdena x och y.<br />

0<br />

0<br />

0<br />

Genom att använda Gauss-Krügers formler som beskrivits ovan, kan uttryck för de<br />

partiella derivatorna härledas. Formlerna (7) ovan ger.<br />

x = k ˆ<br />

0<br />

a f ( ξ ´( λ0<br />

), η´(<br />

λ0<br />

)) + x0<br />

(9a)<br />

y = k a g( ξ ´( λ ), η´(<br />

λ )) +<br />

(9a)<br />

0<br />

ˆ<br />

0 0<br />

y0<br />

De partiella derivatorna blir<br />

∂x<br />

∂k<br />

0<br />

∂y<br />

∂k<br />

0<br />

= aˆ<br />

f<br />

= aˆ<br />

g<br />

∂x<br />

∂x0<br />

∂y<br />

∂x<br />

0<br />

= 1<br />

= 0<br />

∂x<br />

∂y<br />

0<br />

∂y<br />

∂y<br />

0<br />

= 0<br />

= 1<br />

(10a)<br />

(10a)<br />

∂x<br />

∂λ<br />

0<br />

= k<br />

0<br />

⎧ ∂f<br />

∂ξ´<br />

â⎨<br />

⎩ ∂ξ´<br />

∂λ 0<br />

∂f<br />

∂η´<br />

+<br />

∂η´<br />

∂λ<br />

0<br />

⎫<br />

⎬<br />

⎭<br />

(11a)<br />

∂y<br />

∂λ<br />

0<br />

= k<br />

0<br />

⎧ ∂g<br />

∂ξ´<br />

â⎨<br />

⎩ ∂ξ´<br />

∂λ 0<br />

∂g<br />

∂η´<br />

+<br />

∂η´<br />

∂λ<br />

0<br />

⎫<br />

⎬<br />

⎭<br />

(11a)<br />

9


Enligt ekv (7) och (9) fås<br />

4<br />

f ( ξ ′,<br />

η′<br />

) = ξ ′ + ∑ β sin 2 ξ ′ cosh 2iη′<br />

i<br />

i<br />

+ ... (12a)<br />

í = 1<br />

4<br />

g ( ξ ′,<br />

η ′)<br />

= η′<br />

+ ∑ β cos 2 ξ ′ sinh 2iη′<br />

i<br />

i<br />

+ ... (12b)<br />

í = 1<br />

För millimeterprecision räcker fyra termer i serieutvecklingen gott och väl (Reit, 2009).<br />

Formel (12) ger<br />

4<br />

∂f<br />

= 1+<br />

′ ′ +K<br />

∂ ′<br />

∑ 2i<br />

βi<br />

cos 2i<br />

ξ cosh 2iη<br />

ξ<br />

í = 1<br />

4<br />

∂f<br />

=<br />

′ ′ +K<br />

∂ ′<br />

∑ 2i<br />

βi<br />

sin 2i<br />

ξ sinh 2iη<br />

η<br />

í = 1<br />

∂g<br />

= −<br />

∂ξ<br />

′<br />

4<br />

∑<br />

í = 1<br />

∂g<br />

= 1+<br />

∂η′<br />

2i<br />

β sin 2i<br />

ξ ′ sinh 2iη′<br />

+K<br />

4<br />

∑<br />

í = 1<br />

i<br />

2i<br />

β cos 2i<br />

ξ ′ cosh 2iη′<br />

+K<br />

i<br />

(13a)<br />

(13b)<br />

(14a)<br />

(14b)<br />

En jämförelse mellan ekvation (13) och ekvation (14) visar att (13a) och (14b) är identiska.<br />

Samma sak gäller (13b) och (14a), d.v.s.<br />

∂f<br />

∂g<br />

=<br />

∂ξ′ ∂η′<br />

och<br />

∂f<br />

∂g<br />

= −<br />

∂η′ ∂ξ′<br />

Detta samband är generellt och gäller för alla konforma <strong>av</strong>bildningar. Det brukar benämnas<br />

Cauchy-Riemanns differentialekvationer (Reit, 2009).<br />

Till slut ger ekvationerna (6) efter viss formelhantering formlerna<br />

∂ξ<br />

′ sinϕ<br />

*cosϕ<br />

*sin( λ - λ0<br />

)<br />

= −<br />

2<br />

2 2<br />

∂λ<br />

sin ϕ * + cos ϕ *cos ( λ -<br />

0 λ0<br />

∂η′<br />

cosϕ<br />

*cos( λ - λ0<br />

)<br />

= −<br />

2<br />

2 2<br />

∂λ<br />

sin ϕ * + cos ϕ *cos ( λ -<br />

0 λ0<br />

)<br />

)<br />

(15a)<br />

(15b)<br />

För att kunna ställa upp observationsekvationerna behövs även närmevärden för de<br />

obekanta inför det första iterationsvarvet. Tilläggen x<br />

0<br />

och y<br />

0<br />

ingår linjärt i (9)<br />

och kan därmed sättas till 0. Tester som gjorts visar att även λ<br />

0<br />

kan sättas till 0. För<br />

k<br />

0<br />

väljs värdet 1 (Reit, 2009).<br />

10


Med minsta-kvadrat-metoden löses korrektionerna till parametrarna ut, och läggs<br />

till nästa iterationsvarvs närmevärden. Vanligtvis konvergerar iterationerna sedan<br />

snabbt mot ett slutvärde (Reit, 2009).<br />

Under förutsättning att de <strong>plan</strong>a koordinaterna har sitt ursprung i en transversell Mercatorprojektion<br />

(Reit, 2009) och att det lokala systemet täcker en rimligt stor yta, kommer en<br />

direktprojektion <strong>av</strong> det globala systemet att ge koordinater i god överensstämmelse med det<br />

lokala systemet (Engberg & Lilje, 2006).<br />

I <strong>Oskarshamns</strong> <strong>kommun</strong> fanns förutsättningar för tillämpning <strong>av</strong> direktprojektionsmetoden<br />

enligt det utförande som precis beskrivits. Bland andra <strong>kommun</strong>er finns det dock exempel där<br />

metoden ensam inte kan tillämpas rakt <strong>av</strong> vid byte <strong>av</strong> <strong>referenssystem</strong>. Om det lokala systemet<br />

är alltför felorienterat i förhållande till det globala, måste direktprojektionen nämligen<br />

användas i kombination med en transformation. Gränsens för hur stor rotationen <strong>av</strong> det lokala<br />

<strong>plan</strong>a systemet i förhållande till det global geodetiska systemet får vara ligger omkring 1-5<br />

mgon (Reit, 2009). För direktprojektion i kombination med transformation har Lantmäteriet i<br />

huvudsak använt sig <strong>av</strong> två olika metoder eftersom möjligheterna att implementera metoderna<br />

varierar för olika programvaror. Den första metoden kombinerar direktprojektionen med en<br />

efterföljande Helmerttransformation, medan direktprojektionen i den andra metoden föregås<br />

<strong>av</strong> en sjuparametertransformation. Resultaten för de två metoderna skiljer sig inte så mycket<br />

åt vad gäller noggrannhet, men de ger däremot olika koordinatvärden. Vilken metod som<br />

används varierar därför och <strong>av</strong>görs i varje enskilt fall (Engberg & Lilje, 2006).<br />

2.3 Restfel<br />

En central term i detta arbete var begreppet restfel. Begreppet restfel, eller passfel, som det<br />

också ofta benämns, betecknar det beräknade eller uppskattade värdet <strong>av</strong> ett <strong>referenssystem</strong>s<br />

deformation i en given punkt. Förutsättningen för att ett restfel ska kunna tas fram, är att det<br />

finns ett deformationsfritt <strong>referenssystem</strong> att jämföra med; ett homogent <strong>referenssystem</strong> som<br />

får utgöra facit. I det här sammanhanget är det SWEREF 99 som utgör detta facit, medan de<br />

lokala stomnäten, vars geometrier ofta är bristfälliga, är de <strong>referenssystem</strong> som restfel ska<br />

beräknas för.<br />

Den bristande geometrin i de lokala stomnäten beror på att de ofta, utgående från punkter <strong>av</strong><br />

högre ordning, successivt har förtätats, utvidgats för att omfatta större områden eller<br />

tillkommit genom hopslagning <strong>av</strong> flera mindre stomnät till ett större. Denna utveckling har<br />

lett till att spänningar och deformationer uppstått, och stomnätens fjärrnoggrannhet<br />

försämrats. Noggrannheten mellan närliggande punkter har dock i stor utsträckning kunnat<br />

bevaras tack vare att mätning traditionellt skett relativt de närmast omkringliggande punkterna<br />

(Lantmäteriet, informationsblad, 2009).<br />

Arbetet med att ta fram restfel för ett <strong>referenssystem</strong> börjar med inmätning <strong>av</strong> så kallade<br />

passpunkter. Som tidigare nämnts (se <strong>av</strong>snitt 2.2) är en passpunkt en mätpunkt som är<br />

koordinatsatt i två olika <strong>referenssystem</strong>. I SWEREF 99-sammanhang handlar det om att mäta<br />

in punkter i lokala <strong>kommun</strong>ala stomnät, punkter som alltså redan har koordinater bestämda i<br />

lokala system, och förse dem med koordinater i SWEREF 99. När passpunkterna väl mätts in<br />

kan en inpassning göras och parametrar för transformation mellan de två <strong>referenssystem</strong>en<br />

beräknas. I stora drag är det detta förfarande som ligger bakom framtagandet <strong>av</strong> Rix 95-<br />

sambanden mellan lokala och nationella <strong>referenssystem</strong> (se kapitel 1.2.1.1). Med hjälp <strong>av</strong><br />

transformationssambanden kan koordinaterna i ett <strong>av</strong> systemen transformeras till det andra.<br />

11


Differenserna mellan passpunkternas två olika koordinatangivelser, vilka efter<br />

transformationen är angivelser i samma <strong>referenssystem</strong>, är det som benämns restfel. Restfelet<br />

i en passpunkt utgörs alltså <strong>av</strong> differensen mellan de kända koordinaterna i ett <strong>referenssystem</strong><br />

och de, till samma <strong>referenssystem</strong>, transformerade koordinaterna för just den punkten.<br />

2.3.1 Restfelsmodell<br />

Utifrån restfelen i ett antal passpunkter kan en restfelsmodell genereras, som sedan kan<br />

användas för att korrigera koordinater för stompunkter och andra geografiska data så att<br />

deformationer minimeras. Den metod som Lantmäteriet har använt för detta ändamål innebär<br />

att en heltäckande och icke-överlappande restfelsmodell skapas inom det konvexa höljet 3 för<br />

alla passpunkterna med hjälp <strong>av</strong> en algoritm som knyter samman passpunkterna så att det<br />

bildas ett nätverk <strong>av</strong> trianglar. Det finns olika algoritmer som skulle kunna nyttjas för denna<br />

triangelbildning, men med tanke på att ansträngningar har gjorts i ett tidigare skede för att få<br />

en så jämn fördelning <strong>av</strong> passpunkterna som möjligt och att noggrannheten mellan<br />

närliggande passpunkter är hög, är Delaunay-triangulering 4 den algoritm som ansetts bäst<br />

lämpad för ändamålet. En Delaunay-triangulering resulterar i trianglar som är så ”likvinkliga”<br />

som möjligt; trianguleringen blir så regelbunden som möjligt, vilket är optimalt när den<br />

restfelsmodell som trianglarna tillsammans utgör ska användas för att räta upp en bristande<br />

<strong>referenssystem</strong>sgeometri. I <strong>av</strong>snitt 3.10 redovisas utseendet för restfelsmodellerna i<br />

<strong>Oskarshamns</strong> <strong>kommun</strong>.<br />

Upprätningen går till på så sätt att korrektioner, eller restfel, för de kartdetaljer som täcks in<br />

<strong>av</strong> restfelsmodellen interpoleras fram med hänsyn till i vilken triangel respektive kartdetalj<br />

återfinns. Varje interpolerat restfel bestäms till storlek och riktning utifrån restfelen för den<br />

aktuella triangelns tre hörnpasspunkter. Interpoleringen är en så kallad<br />

gummidukstransformation där en affin transformation används, baserad på restfelen hos de tre<br />

passpunkterna som bildar varje triangel i restfelsmodellen. Mellan triangelpunkterna ger<br />

interpoleringen kontinuitet i de interpolerade värdena och i passpunkterna erhålls de exakta<br />

till-koordinaterna. I närheten <strong>av</strong> det konvexa höljet bildas ofta trianglar med väldigt spetsiga<br />

och trubbiga vinklar, vilket gör att trianglarna får ett platt utseende och därmed blir olämpliga<br />

för interpolering. För att undvika detta problem och problem med extrapolering <strong>av</strong> punkter<br />

utanför konvexa höljet, bör försök göras att även hitta passpunkter som ligger utanför<br />

tillämpningsområdet; exempelvis en <strong>kommun</strong>s primärkarteområde, så att punkter som ska<br />

interpoleras i största möjligaste mån ligger innanför konvexa höljet (Gtrans, 2008).<br />

2.3.2 Redovisning och presentation <strong>av</strong> restfel<br />

Restfel kan exempelvis redovisas grafiskt som restfelsvektorer eller i restfelsvariationsbilder.<br />

I vektorkartor redovisas restfelen till position, storlek och riktning, vilket gör att man snabbt<br />

kan få en överblick <strong>av</strong> i vilka områden stark korrelation mellan närliggande punkters restfel<br />

föreligger samt var <strong>av</strong>vikande restfel finns. Stark korrelation uppnås i områden där<br />

närliggande restfelsvektorer har liknande storlek och riktning, medan korrelationen blir svag<br />

där restfelen varierar i storlek och mer eller mindre spretar åt olika håll. Områden med svagt<br />

3 Konvexa höljet är den minsta konvexa mängd som innesluter en given mängd. I <strong>plan</strong>et är konvexa höljet för en<br />

punktmängd lika med en konvex polygon.<br />

4 Ordet triangulering har här inte något att göra med mätning, beräkning och utjämning <strong>av</strong> geodetiska nät, utan<br />

<strong>av</strong>ser just nätverk <strong>av</strong> trianglar.<br />

12


korrelerade restfel ger upphov till en deformation <strong>av</strong> <strong>plan</strong>et. Beroende på hur restfelen för en<br />

triangels hörnpasspunkter pekar, kan en deformation uppfattas på olika sätt. Om alla tre<br />

restfelen är riktade in mot triangeln föreligger en kontraktion, om alla restfelen är riktade ut<br />

från triangeln är deformationen en töjning och om alla restfelen är riktade åt samma håll<br />

samtidigt som de är lika stora handlar det om en translation och ingen deformation. Mellan<br />

dessa tre huvudgrupper finns det naturligtvis en mängd olika variationer (Gtrans, 2008).<br />

I Lantmäteriets transformationsverktyg Gtrans används ett särskilt mått: ”RMAX”, på den<br />

maximala deformationen i en triangel. ”RMAX” antar värdet noll för ingen deformation,<br />

större än noll för töjning och mindre än noll för kontraktion. Detta mått gör det möjligt att i ett<br />

visst område <strong>av</strong>göra om interpoleringen är bra eller dålig och även vilken storleksordning<br />

felet har.<br />

Lantmäteriet har även tagit fram ett verktyg i programmet Matlab för att underlätta analysen<br />

<strong>av</strong> deformationernas variationer inom stomnäten. Matlab-verktyget tar restfelsvektorerna som<br />

indata och åskådliggör efter bearbetning restfelsvariationerna i form <strong>av</strong> en färgbildskarta. För<br />

att beskriva hur denna bearbetning går till betraktar vi en godtycklig triangelsida i en <strong>av</strong> de<br />

Delaunay-trianglar som bygger upp restfelsmodellen. I sidans två ändpunkter finns då en<br />

passpunkt med tillhörande restfel (se figur 3). Följande beräkningsprocedur utförs (Kempe et<br />

al., 2006).<br />

1. Restfelsvektorn för en <strong>av</strong> de två passpunkterna flyttas längs triangelsidan till den andra<br />

restfelsvektorn (se figur 4).<br />

2. En vektorsubtraktion utförs, med resultanten (resulterande vektorn) som resultat (se<br />

figur 5).<br />

3. Värdet på resultantens längd d hanteras nu som en skalär och flyttas till mitten <strong>av</strong><br />

triangelsidan (se figur 6). Normalt viktas d omvänt proportionellt mot kvadratroten ur<br />

passpunkts<strong>av</strong>ståndet, det vill säga triangelsidans längd, men de absoluta differenserna<br />

kan i många fall också vara till stor hjälp i resfelsanalysen.<br />

Figur 3. Delauny-triangel med hörnpunkternas restfel<br />

redovisade som blå vektorer<br />

Figur 4. En <strong>av</strong> restfelsvektorerna förflyttas längs<br />

triangelsidan till en annan passpunkt.<br />

13


d<br />

d<br />

Figur 5. En vektorsubtraktion resulterar i en vektor<br />

med längden d.<br />

Figur 6. Punkten mitt på triangelsidan ges värdet d,<br />

vilket anger hur restfelen varierar.<br />

Om proceduren sedan utförs på samma sätt för resten <strong>av</strong> triangelsidorna i modellen erhålls en<br />

uppsättning d-värden som kan presenteras exempelvis i en färgbildskarta, där olika stora<br />

värden resulterar i olika färger. Med hjälp <strong>av</strong> dessa tematiska kartor kan man därmed<br />

visualisera hur restfelen varierar i stomnäten och snabbt få en överblick <strong>av</strong> var spänningar<br />

finns och var dessa eventuellt behöver kartläggas ytterligare.<br />

För exempel på restfelsvektorkartor och restfelsvariationsbilder; se den följande<br />

beskrivningen <strong>av</strong> restfelsanalys och transformation.<br />

14


3. Restfelsanalys och transformation<br />

3.1 Resultatet <strong>av</strong> Rix 95-mätningarna i Småland och Oskarshamn<br />

Utvärderingen <strong>av</strong> de lokala stomnätens geometriska kvalitet tog, som redan nämnts, sitt<br />

<strong>av</strong>stamp i resultatet <strong>av</strong> Lantmäteriets Rix 95-mätningar i Småland och <strong>Oskarshamns</strong> <strong>kommun</strong>.<br />

Dessa mätningar gjordes i oktober-november 2001 och analysen för den delen i projektet<br />

slutfördes i början <strong>av</strong> året efter. När det gäller <strong>Oskarshamns</strong> <strong>kommun</strong> mättes sammanlagt 37<br />

punkter i regionsystemet och 24 punkter i <strong>Oskarshamns</strong> lokala koordinatsystem. Rix 95-<br />

inpassningens grundmedelfel ( σ<br />

0<br />

, 2D) var 51 millimeter för regionsystemet och 41<br />

millimeter för <strong>Oskarshamns</strong> lokala koordinatsystem. De kvadratiska medelvärdena (RMS,<br />

2D) för restfelen i de två systemen var 51 millimeter för regionsystemet respektive 41<br />

millimeter för <strong>Oskarshamns</strong> lokala koordinatsystem. De största restfelen (2D) beräknades till<br />

173 millimeter i regionsystemet respektive 124 millimeter i <strong>Oskarshamns</strong> lokala<br />

koordinatsystem.<br />

Tabellen nedan redovisar de Rix 95-parametrar som Lantmäteriet hade beräknat för<br />

”transformation” mellan latitud och longitud i SWEREF 99 och <strong>plan</strong>a koordinater i respektive<br />

koordinatsystem i <strong>Oskarshamns</strong> <strong>kommun</strong>.<br />

Tabell 1. Rix 95-parametrarna för de lokala koordinatsystemen i <strong>Oskarshamns</strong> <strong>kommun</strong><br />

Regionsystemet<br />

<strong>Oskarshamns</strong> lokala koordinatsystem<br />

Projektionstyp Gauss-Krüger Gauss-Krüger<br />

Medelmeridian 15° 48’ 22.572” Öst Greenwich 15° 48’ 23.580” Öst Greenwich<br />

Latitud för origo 0° 0°<br />

Skalreduktionsfaktor 1.00000506 (+5.06ppm) 0.99999579 (-4.21ppm)<br />

FN, x-tillägg -663.6418 m -606.0353 m<br />

FE, y-tillägg +1500063.7352 m +1500080.3408 m<br />

Grundmedelfelet är ett precisionsmått och visar spridningen <strong>av</strong> restfelen. Måttet <strong>av</strong>ser<br />

medelfelet för mätningar med vikten ett (Fan, 2003). De kvadratiska medelvärdena är<br />

statistiska mått på storleken <strong>av</strong> kvantiteter som varierar. Särskilt användbara är dessa mått när<br />

det finns såväl positiva som negativa värden i mätserien, eftersom aritmetiska medelvärden i<br />

dessa fall inte ger så mycket information. Information om största restfelet kan även det vara<br />

intressant, i och för sig, men det säger inte så mycket om restfelen som helhet. För att kunna<br />

upptäcka spänningar mellan olika stomnätsdelar och <strong>av</strong>vikande restfel, räcker inte enbart<br />

statistiska mått, utan restfelen måste även redovisas i deras respektive geografiska lägen. För<br />

denna redovisning används, som redan beskrivits, exempelvis restfelsvektorkartor och<br />

restfelsvariationsbilder. Sådana kartor och bilder utgör en del <strong>av</strong> Lantmäteriets redovisningar<br />

till <strong>kommun</strong>erna <strong>av</strong> Rix 95-mätningarnas resultat, så den ”restfelsgrafiken” fanns att tillgå<br />

även för <strong>Oskarshamns</strong> <strong>kommun</strong> redan från utvärderingens början.<br />

15


Figur 7. Restfelsvektorer för Rix 95-mätningarna i<br />

regionsystemet<br />

De röda cirklarna (dubbelcirklar) visar områden med<br />

restfel som pekar åt olika håll, medan de blå cirklarna<br />

markerar områden med eventuella trender hos restfelen.<br />

Figur 8. Restfelsvektorer för Rix 95-mätningarna i<br />

<strong>Oskarshamns</strong> lokala koordinatsystem<br />

Efter att ha studerat de grafiska presentationerna <strong>av</strong> restfelen och deras variationer inom<br />

respektive koordinatsystems tillämpningsområde i <strong>kommun</strong>en, framstod det än mer klart att<br />

en ytterligare och noggrannare kartläggning och med största säkerhet även upprätning <strong>av</strong><br />

stomnätsgeometrierna var nödvändig. Restfelsvektorerna i figur 7 och 8 tillsammans med de<br />

tematiska redovisningarna <strong>av</strong> restfelens variationer i figur 9 och 10 visar visserligen inte på<br />

några direkt alarmerande förhållanden, bortsett från i ett antal väldigt tydliga fall, men<br />

följande iakttagelser bör noteras:<br />

1. Stora restfel – Figurerna 7 och 8 visar flera restfel som är omkring en decimeter stora,<br />

vilket är mycket jämfört med de flesta andra vektorerna.<br />

2. Närliggande vektorer med olika riktningar – Vektorerna i de röda cirklarna<br />

(dubbelcirklar) i figur 7 och 8 pekar åt mycket olika håll, trots att de ligger nära<br />

varandra. I figur 9 och 10 ger detta upphov till mer eller mindre tydliga röda<br />

färgpartier.<br />

3. Eventuella trender – De blå cirklarna i figur 7 och 8 visar några områden där<br />

vektorernas riktningar är förhållandevis lika. Det skulle eventuellt kunna vara fråga<br />

om trender hos restfelen.<br />

4. Stomnätsdelar som inte är representerade – Detta är framförallt tydligt när det gäller<br />

<strong>Oskarshamns</strong> lokala koordinatsystem.<br />

5. Antalet inmätta punkter är väldigt litet i förhållande till stomnätens alla punkter<br />

(jämför med figur 2), vilket gör att passpunkts<strong>av</strong>stånden blir alltför stora för att de<br />

spänningar som finns i stomnäten ska kunna fångas upp. Lantmäteriets<br />

16


ekommendationer är därför att passpunkts<strong>av</strong>stånden i områden med primärkarta inte<br />

bör överskrida 500-700 meter (Kempe & Lohász., 2008).<br />

På de iakttagelser som ovanstående numrerade noteringar beskriver, följde bedömningen att<br />

en förtätning <strong>av</strong> mätpunkterna var ofrånkomlig. Kartläggningen <strong>av</strong> restfelens variationer var<br />

helt enkelt än så länge alltför bristfällig. På grund <strong>av</strong> de stora passpunkts<strong>av</strong>stånden <strong>av</strong>gjordes<br />

att behovet <strong>av</strong> förtätning var generellt och behövde göras i alla <strong>kommun</strong>ens områden med<br />

primärkarta.<br />

Figur 9. Restfelsvariationer för Rix 95-mätningarna i<br />

regionsystemet<br />

Restfelsvariationsbilderna är uppbyggda <strong>av</strong> färger i skalan vitt till<br />

rött, där den vita färgen används i områden med inga eller mycket<br />

små restfelsvariationer och röd färg återfinns i de områden där<br />

restfelen varierar som mest till storlek och riktning. De olika<br />

färgerna representerar alltså olika stora restfelsvariationer eller<br />

olika stora d-värden, vilka beskrivits mer utförligt i <strong>av</strong>snitt 2.3.2.<br />

Figur 10. Restfelsvariationer för Rix 95-<br />

mätningarna i <strong>Oskarshamns</strong> lokala<br />

koordinatsystem<br />

3.2 Planering <strong>av</strong> kompletteringsmätningar<br />

Rix 95-projektet visade alltså att kompletteringsmätningar behövde göras för att den bristande<br />

stomnätsgeometrin skulle kunna kartläggas och därefter rätas upp. Det optimala hade givetvis<br />

varit att samtliga stomnätspunkter i <strong>kommun</strong>en mättes in, eftersom det ju är<br />

17


stomnätspunkterna som utgjort utgångspunkter för i princip all mätning för ajourhållandet <strong>av</strong><br />

primärkarteverket. Av praktiska skäl var en så omfattande mätinsats naturligtvis inte<br />

genomförbar och det handlade därför om att göra ett så bra urval <strong>av</strong> stomnätspunkterna som<br />

möjligt. De punkter som mättes in skulle väljas på sådant sätt att deras deformationer eller<br />

restfel tillsammans kunde antas representera primärkartans deformationer. Med detta syfte i<br />

åtanke kunde några övergripande riktlinjer för utväljandet <strong>av</strong> passpunkter dras upp. Stabila<br />

punkter som ofta använts skulle tas med eftersom deras restfel överförts till många<br />

kartdetaljer. Vidare skulle passpunkterna väljas så att stomnätens delar i så hög grad som<br />

möjligt omringades. Detta var viktigt för att få en bra triangelbildning i restfelsmodellen och<br />

därmed undvika att deformationer i en stomnätsdel påverkade korrektionen i intilliggande<br />

delar (Kempe et al., 2006). Eftersom bedömningen gjordes att en generell förtätning <strong>av</strong><br />

passpunkterna var lämplig, skulle också en jämn fördelning <strong>av</strong> passpunkterna över hela<br />

stomnätsytan eftersträvas i mätningarnas första etapp.<br />

Samtidigt som passpunkterna valdes enligt ovanstående riktlinjer, var det helt nödvändigt att<br />

också tänka på vilka praktiska förutsättningar för mätning som fanns vid de utvalda<br />

punkterna. För att erhålla tillförlitliga mätresultat med GNSS-teknik krävs framförallt god<br />

satellittillgänglighet, vilket kan vara svårt att uppnå på många platser i stomnätens utkanter.<br />

Det gällde alltså att i möjligaste mån <strong>plan</strong>era mätning på punkter där risken för flervägsfel var<br />

liten och där satellitsignalerna inte hindrades alltför mycket <strong>av</strong> höga hus eller träd.<br />

Rent praktiskt gick valet <strong>av</strong> passpunkter till på följande sätt. Samtliga Rix 95-punkter och<br />

några få redan mätta kontroll-/kompletteringspunkter markerades på kartor över <strong>kommun</strong>ens<br />

stomnät. Utifrån dessa punkter valdes sedan flera stomnätspunkter ut till passpunkter så att<br />

<strong>av</strong>ståndet mellan dem var en till en och en halv kilometer. I de fall en utvald punkt fanns med<br />

på <strong>kommun</strong>ens egna ortofoton eller i andra fritt tillgängliga högupplösta ortofotosamlingar,<br />

som exempelvis Googles virtuella jordglob Google Earth, gjordes utifrån detta material en<br />

snabb och grov bedömning <strong>av</strong> förutsättningarna för GNSS-mätning. I de områden som inte<br />

täcktes in <strong>av</strong> ortofoton skickades en person med god lokalkännedom ut på en<br />

rekognoseringsrunda innan mätningarna där påbörjades. Dessa mätningsförberedande<br />

åtgärder gjorde att svårmätta punkter i många fall redan vid skrivbordet kunde väljas bort och<br />

därmed slapp besökas med mätutrustning, vilket annars hade krävt mycket arbete och tagit tid<br />

till ingen nytta.<br />

3.3 Val <strong>av</strong> mätmetod och mätteknik<br />

När det första urvalet <strong>av</strong> stomnätspunkter hade gjorts var det dags att välja metod för<br />

inmätning <strong>av</strong> passpunkterna. Eftersom kartbasens framtida noggrannhet skulle bli direkt<br />

beroende <strong>av</strong> passpunkternas uppmätta restfel, var det <strong>av</strong> största betydelse att valet föll på en<br />

mätmetod som gjorde det möjligt att utföra mätningarna på ett kontrollerat sätt. Det var restfel<br />

<strong>av</strong>seende skillnaderna mellan till- och frånsystemet i transformationen och inte mätfel som<br />

skulle modelleras. De tänkbara metoderna innefattade antingen statisk mätning eller RTKmätning.<br />

Därför fanns det anledning att göra en liten granskning <strong>av</strong> dessa två olika tekniker<br />

för GNSS-mätning samt ta reda på förutsättningar, för- och nackdelar för respektive teknik.<br />

3.3.1 Statisk mätning<br />

Statisk mätning är en variant <strong>av</strong> relativ mätning som ofta bygger på att satellitgeometrin<br />

förändras under observationstiden. För att geometrin ska förändras tillräckligt mycket krävs<br />

att observationstiden är minst 20 minuter lång. Ju längre baslinjerna, det vill säga <strong>av</strong>stånden<br />

18


mellan mottagarna, är desto längre observationstid krävs för att kunna lösa periodobekanta<br />

och få noggrant bestämda positioner. Med enfrekvensmätning kan baslinjer upp till 30<br />

kilometer bestämmas på 45-60 minuters observationstid. För längre <strong>av</strong>stånd blir<br />

tvåfrekvensmätning ett kr<strong>av</strong> (HMK-Ge:GPS, 1993).<br />

Ett alternativ till vanlig statisk mätning är snabb statisk mätning med kortare observationstid.<br />

På 5-20 minuter kan baslinjer upp till 20 kilometer bestämmas. De beräkningsalgoritmer som<br />

då används liknar dem som används vid RTK-mätning för att bestämma periodobekanta.<br />

Jämfört med vanlig statisk mätning blir noggrannheten med den snabbare mätmetoden inte<br />

lika hög.<br />

I ett projekt <strong>av</strong> den omfattning som inmätning <strong>av</strong> passpunkter för ett <strong>referenssystem</strong>sbyte har,<br />

ställer en metod med statisk mätning extra stora kr<strong>av</strong> på god och noggrann <strong>plan</strong>ering innan<br />

mättillfällena. Det finns inget utrymme för misstag, vilket i den meningen gör metoden sårbar<br />

och mindre flexibel än exempelvis RTK-metoderna. Resultatet <strong>av</strong> mätningar i bra <strong>plan</strong>erade<br />

sessioner är dock ett nät med utmärkt kontrollerbarhet och goda möjligheter till ingående<br />

geodetiska analyser. Som utgångspunkter vid statisk mätning används punkter som är kända i<br />

det aktuella <strong>referenssystem</strong>et; i detta fall SWEREF 99. Lantmäteriet rekommenderar att man<br />

håller en mottagare fast på en känd punkt under hela mätningen; exempelvis den ursprungliga<br />

utgångspunkten eller närmast tillgängliga Rix 95-punkt.<br />

Det normala är att beräkning <strong>av</strong> mätdata från statisk mätning – bärvågsmätningar – sker i<br />

efterhand i program som klarar baslinje- och nätutjämningsberäkningar. Jämfört med andra<br />

positionsbestämningsmetoder ger statisk mätning den högsta noggrannheten. Beroende på<br />

observationstid, baslinjelängd och beräkningsalgoritmer, är medelfelet mellan några<br />

millimeter och några centimeter (Lilje et al., 2007).<br />

Vid statisk mätning är mätinsatserna väldigt intensiva och resurskrävande dels<br />

personalmässigt och dels utrustningsmässigt. För att statisk mätning ska kunna genomföras på<br />

ett rationellt sätt rekommenderar Lantmäteriet att antalet GNSS-mottagare inte är mindre än<br />

fyra (Kempe et al., 2006). Med en eller två operatörer per mottagare inses då lätt att den<br />

dagliga mätverksamheten i exempelvis en <strong>kommun</strong> kan bli lidande; särskilt om flera dagar i<br />

rad måste tas i anspråk. Å andra sidan blir mätinsatserna mer koncentrerade och inte så<br />

utdragna tidsmässigt som vid nyttjande <strong>av</strong> RTK-metoder.<br />

3.3.2 RTK-mätning<br />

Vid RTK 5 -mätning tillämpas relativ bärvågsmätning i realtid. En referensstation med känd<br />

position skickar över bärvågsdata eller korrektioner som kombineras med bärvågsdata från<br />

den rörliga mottagaren. Fixlösning fås om de periodobekanta kan lösas; det vill säga om<br />

mottagaren kan initialiseras. Detta kan göras på lite olika sätt men idag används nästan alltid<br />

så kallad flygande bestämning <strong>av</strong> periodobekanta, vilket kräver tillgång till minst fem<br />

satelliter. Flygande bestämning <strong>av</strong> periodobekanta kan utföras även när mottagaren är i<br />

rörelse.<br />

5 RTK står för Real Time Kinematic<br />

19


Initialiseringstiden varierar med antal satelliter, satellitgeometri, <strong>av</strong>stånd till referensstationen<br />

och mottagarens kvalitet. Normal tidsåtgång kan vara allt mellan ett tiotal sekunder till några<br />

minuter. Om den fixlösning som erhållits vid initialiseringen tappas, måste initialiseringen<br />

göras om.<br />

RTK-mätning kan utföras med hjälp <strong>av</strong> enkelstations-RTK eller nätverks-RTK. Det som<br />

skiljer de olika teknikerna åt är framförallt det antal referensstationer som används. Till<br />

nätverks-RTK brukar en tjänst vara kopplad, vilken möjliggör flera samtidiga användare.<br />

Precis som benämningen antyder, används data vid enkelstations-RTK från enbart en<br />

referensstation i taget. I RTK-användningens inledningsfas las det ingen större vikt vid<br />

referensstationens och den rörliga mottagarens olika jonosfäriska förhållanden. Slutlösningen<br />

<strong>av</strong> de periodobekanta var en ren L1-lösning, även om både L1 och L2 användes för att snabba<br />

upp initialiseringen. Detta medförde att räckvidden var begränsad till omkring tio kilometer.<br />

Numera utnyttjas en jonosfärsfri linjärkombination <strong>av</strong> båda frekvenserna, för att reducera<br />

jonosfärens inverkan vid längre <strong>av</strong>stånd. Tack vare denna utveckling <strong>av</strong> tekniken kan<br />

enkelstations-RTK användas med <strong>av</strong>stånd upp till 40 kilometer mellan referensstation och den<br />

rörliga mottagaren.<br />

För att kunna använda enkelstations-RTK behövs alltså minst två RTK-utrustningar; en fast<br />

eller tillfälligt etablerad referensstation och en rörlig mottagare. Om tillgång finns till<br />

enkelstations-RTK via en tjänst, vilket är fullt möjligt om än inte så vanligt förekommande,<br />

med en eller flera referensstationer, räcker det för en användare <strong>av</strong> tjänsten att ha en RTKutrustning.<br />

Förutom en rörlig mottagare och tillgång till minst en referensstation, behövs även<br />

en radiomodem- eller GSM-länk för överföring <strong>av</strong> data mellan utrustningarna (Lilje et al.,<br />

2007).<br />

Om ett antal fasta referensstationer används tillsammans i ett nätverk och inte enbart en i<br />

taget, kan en yttäckande felmodell beräknas, som korrigerar för de fel som uppstår när GNSSsignalerna<br />

går genom atmosfären och för de fel som beror på felaktiga positioner för<br />

satelliternas banor. På detta sätt kan vanlig RTK-data kompletteras och noggrannheten<br />

förbättras, vilket i sin tur innebär att <strong>av</strong>ståndet mellan referensstationerna kan ökas från de 20<br />

till 30 kilometer som gäller i fallet med enkelstations-RTK, till upp mot 70 kilometer, utan att<br />

noggrannheten försämras eller initialiseringstiden förlängs oacceptabelt mycket (Lilje et al.,<br />

2007).<br />

Modellen med flera referensstationer i ett nätverk kallas nätverks-RTK och har bland annat de<br />

stora fördelarna att täckningsområdet saknar skarvar och att den ger kvalitetskontrollerade<br />

data.<br />

Det finns olika sätt att skicka ut data från en nätverks-RTK-tjänst på. I Lantmäteriets nätverk<br />

<strong>av</strong> referensstationer, SWEPOS (se <strong>av</strong>snitt 2.1), låter man data från flera referensstationer<br />

tillsammans skapa en felmodell för området som täcks <strong>av</strong> de berörda referensstationerna. Med<br />

kännedom om en rörlig mottagares approximativa position kan observationsdata från en<br />

referensstation i närheten korrigeras med hjälp <strong>av</strong> felmodellen. Den korrigerade informationen<br />

sägs komma från en virtuell referensstation (VRS 6 ) eftersom även en geometrisk korrektion<br />

6 VRS – Virtual Reference Station<br />

20


<strong>av</strong> informationen har gjorts så att den ser ut att vara insamlad vid den approximativa position<br />

som först skickades från den rörliga mottagaren. Oftast bygger detta datautbyte på<br />

tvåvägs<strong>kommun</strong>ikation – företrädesvis GSM 7 eller GPRS 8 – där data skickas mellan den<br />

rörliga mottagaren och en server med programvara för nätverks-RTK (von Malmborg, 2006).<br />

När den rörliga GNSS-mottagaren väl har upprättat förbindelse med servern på nätverks-<br />

RTK-tjänstens driftledningscentral, sker mätning relativt den virtuella referensstationen. I den<br />

meningen är nätverks-RTK alltså en polär mätteknik (Jansson, 2007).<br />

För inmätning <strong>av</strong> passpunkter inför ett <strong>referenssystem</strong>sbyte har Lantmäteriet tagit fram några<br />

riktlinjer för att tillförlitliga mätresultat ska kunna erhållas vid användning <strong>av</strong> RTK-teknik. På<br />

varje passpunkt ska enligt dessa rekommendationer två mätserier göras med tio mätningar i<br />

varje. Mellan varje mätning i en serie ska en ominitialisering göras så att de periodobekanta<br />

bestäms på nytt. Vidare ska mätserie ett och två göras vid olika tidpunkter på dygnet för att<br />

satellitkonstellationen ska hinna ändras. En annan anledning till att två mätserier bör göras är<br />

att risken för systematiska fel från exempelvis centrering eller satellitsystemet då minskar.<br />

Dålig centrering undviks för övrigt lättast genom att ha stabila uppställningar där GNSSantennen<br />

placeras på stativ (Kempe, 2006).<br />

RTK-mätning är jämfört med statisk mätning en väldigt lättillgänglig teknik tack vare att<br />

utrustningsbehoven är relativt små och att koordinater för inmätta punkter fås direkt i fält,<br />

utan komplicerade efterberäkningar <strong>av</strong> mätdata. I stora projekt blir utmaningen istället att i<br />

<strong>plan</strong>eringsskedet bedöma behovet <strong>av</strong> personal- och instrumentresurser i en väldigt lång och<br />

utdragen process. Utan en noggrann <strong>plan</strong>ering – i princip ner på dagsnivå – riskerar annars<br />

projektet att stanna <strong>av</strong>. Detta är sannolikt en större risk än om statiskt mätning tillämpas.<br />

3.3.3 Förutsättningarna i Oskarshamn<br />

Ur noggrannhetssynpunkt spelade det inte så stor roll om statisk mätning eller RTK-mätning<br />

valdes för inmätning <strong>av</strong> restfelspunkter. Under förutsättning att mätmetoden utformats så att<br />

den ger säkra mätresultat, kan nämligen kr<strong>av</strong>et på centimeternoggrannhet uppfyllas o<strong>av</strong>sett<br />

vilken <strong>av</strong> de båda mätteknikerna som används. Istället var det tillgång på resurser som<br />

<strong>av</strong>gjorde vilken mätmetod och teknik som kom att användas.<br />

<strong>Oskarshamns</strong> <strong>kommun</strong> var den enda <strong>kommun</strong>en i regionen som riktigt hade kommit igång<br />

med arbetet för en övergång till SWEREF 99, så någon naturlig samarbetspart i<br />

mätningsarbetet fanns inte. Eftersom <strong>kommun</strong>ens tillgång på GNSS-mottagare dessutom var<br />

begränsad, blev en RTK-metod därmed ett tämligen självklart val. Mätning med nätverks-<br />

RTK fanns det sedan tidigare vana med och eftersom <strong>kommun</strong>en saknade tillgång till egen<br />

referensstation fanns det ingen anledning att inte fortsätta med just nätverks-RTK-teknik.<br />

3.4 Personal, utrustning och upplägg<br />

Som redan nämnts blir mätinsatserna tidsmässigt väldigt krävande när RTK-teknik används.<br />

För att säkerställa att projektet med inmätning <strong>av</strong> restfelspunkter hela tiden gick framåt, valde<br />

7 GSM – Global System for Mobile Communications<br />

8 GPRS – General Packet Radio Service<br />

21


<strong>Oskarshamns</strong> <strong>kommun</strong> att hyra in en nätverks-RTK-utrustning för användning enbart i dessa<br />

inmätningar. På det sättet behövde aldrig brist på mätinstrument bromsa upp vare sig<br />

mätningarna för <strong>referenssystem</strong>sbytet eller <strong>kommun</strong>ens ordinarie mätverksamhet.<br />

Vidare gjordes det upp ett schema som dag för dag visade vem <strong>av</strong> den tillgängliga<br />

mätpersonalen som skulle mäta. Detta arbetsschema var till stor hjälp inte bara för att driva<br />

projektet framåt, utan även för att säkerställa att det inhyrda instrumentet inte behövde hyras<br />

längre än absolut nödvändigt. Målsättningen var givetvis att hyrutrustningen skulle ”gå varm”<br />

och inte stå outnyttjad i onödan.<br />

En förutsättning för att <strong>Oskarshamns</strong> <strong>kommun</strong> skulle kunna göra mätningarna i egen regi, var<br />

att <strong>kommun</strong>ens vardagliga verksamhet inte blev så lidande att den inte kunde fungera<br />

tillfredsställande. Naturligtvis kan inte en satsning <strong>av</strong> det format som ett <strong>referenssystem</strong>sbyte<br />

utgör, göras utan påfrestningar på andra håll i verksamheten, men genom en sund och rimlig<br />

prioritering bland mätuppdragen kunde <strong>kommun</strong>ens mätpersonal tillsammans med<br />

förstärkning från GIS-funktionen ändå driva SWEREF-projektet parallellt med den löpande<br />

mätverksamheten.<br />

3.5 Mätningarnas första etapp<br />

När det väl var dags att sätta igång med själva passpunktsinmätningen fanns det några saker<br />

som det behövde fästas extra uppmärksamhet vid redan på kontoret. För det första<br />

kontrollerades satellittillgängligheten så att varje dags mätning i möjligaste mån kunde styras<br />

till tidpunkter med god tillgänglighet. Det verktyg som användes för att ta reda på denna<br />

information var SWEPOS satellitprediktionstjänst som fanns och fortfarande finns tillgänglig<br />

på Internet. Tillsammans med mätpunkternas placering och omgivning i terrängen var<br />

satellittillgängligheten naturligtvis helt <strong>av</strong>görande för hur väl GNSS-mätning skulle låta sig<br />

göras. Därför blev mätrutterna i hög grad bestämda utifrån kunskap om detta. I <strong>plan</strong>eringen <strong>av</strong><br />

mätrutterna togs även kr<strong>av</strong>et på förändrad satellitgeometri mellan de två mätserierna på varje<br />

punkt, i beaktande.<br />

Mätningarna gick sedan till på det sätt som beskrivits under <strong>av</strong>snitt 3.3.2 om mätning med<br />

nätverks-RTK-teknik. De <strong>plan</strong>erade restfelspunkterna besöktes med en GNSS-mottagare, en<br />

fältdator samt ett stabilt stativ. Dubbla mätserier på varje punkt gjordes med minst tio<br />

mätningar i varje serie och ny initialisering mellan varje mätning. I övrigt följdes i så stor<br />

utsträckning som möjligt Lantmäteriets råd för parametervärden rörande GNSS-mätning med<br />

höga noggrannhetskr<strong>av</strong> – mätning på så kallad ”förhöjd nivå 2”. På den nivån rekommenderar<br />

Lantmäteriet bland annat att inte färre än sju GNSS-satelliter, inte högre PDOP-värden än två,<br />

inte mindre än 90 procents fri sikt ner till 15 graders elevationgräns samt inte längre<br />

initialiseringstider än en minut, ska accepteras (Norin et al., 2007). Till detta fördes också<br />

riktlinjen att inte acceptera alltför höga värden på GNSS-utrustningens interna kvalitetstal.<br />

Råden var en bra hjälp i försöken att få så tillförlitliga och bra mätningsresultat som möjligt<br />

men de var inte alltid helt lätta att följa. Detta gällde särskilt i områden där passpunkter låg i<br />

gränsen mellan samlad bebyggelse i ett förhållandevis öppet landskap och skogsområden.<br />

I mätetapp ett mättes sammanlagt 95 passpunkter, vilka utgjorde drygt två procent <strong>av</strong><br />

stomnätens alla punkter. Av dessa 95 punkter hörde 39 till regionsystemet och 56 till<br />

<strong>Oskarshamns</strong> lokala koordinatsystem. De mätta punkterna i regionsystemet var fördelade på<br />

åtta tätorter eller platser med samlad bebyggelse, medan punkterna som mättes in i<br />

22


<strong>Oskarshamns</strong> lokala system till omkring 66 procent återfanns i <strong>Oskarshamns</strong> stad och till<br />

resterande 44 procent i en tätort och ett par mindre platser med samlad bebyggelse.<br />

Som tidigare nämnts hade passpunkterna valts så att stomnätens delar i största möjligaste mån<br />

omringades. Detta innebar att det även i mindre orter mättes några punkter i stomnätets<br />

ytterkanter och några inom stomnätet, trots att passpunkterna därmed kom relativt nära<br />

varandra och bildade ett tätare nät än i större orter med större stomnätsdelar. Det viktiga var<br />

inte att nätet <strong>av</strong> inmätta stompunkter överallt blev lika tätt, utan att passpunkternas restfel på<br />

ett tillfredsställande sätt verkligen representerade de spänningar som fanns i stomnätet. Om<br />

detta skulle vara möjligt var det nödvändigt att restfelsmodellens trianglar fick ett så optimalt<br />

utseende som möjligt (se <strong>av</strong>snitt 2.3.1).<br />

3.6 Fortsatt analys<br />

Under mätetapp ett gjordes fortlöpande analyser <strong>av</strong> de utförda mätningarna för att bedöma om<br />

strategin i några områden direkt skulle revideras eller om mätningarna på de <strong>plan</strong>erade<br />

punkterna verkade kunna ge en bra kartläggning <strong>av</strong> stomnätens bristande geometrier.<br />

Analyskedjan började nu med att medeltalet för varje mätserie på varje punkt beräknades,<br />

varefter varje enskild mätnings <strong>av</strong>vikelse från medeltalet kontrollerades. Alla mätningar med<br />

<strong>av</strong>vikelser större än fyra centimeter från medeltalet plockades bort. I vissa fall behövde många<br />

mätningar plockas bort, vilket gjorde att medeltalsberäkningen till slut innehöll så få<br />

överbestämningar att det statistiska underlaget blev alltför bristfälligt. Antalet mätningar inom<br />

varje mätserie ökades därför därefter till närmare 15. De två seriernas medeltal för varje punkt<br />

jämfördes och accepterades endast om differensen hamnade under fyra centimeter. Att just<br />

fyra centimeter användes som ett riktvärde var ett resultat <strong>av</strong> en bedömning som Lantmäteriet<br />

gjort utifrån noggrannheten hos den mätteknik och den mätmetod som nyttjades (Odolinski &<br />

Sunna, 2009). Om de två medeltalen inte <strong>av</strong>vek mer än tillåtet från varandra, beräknades ett<br />

nytt medeltal <strong>av</strong> dem. Detta sista medeltal utgjorde resultatet <strong>av</strong> mätningarna på punkten,<br />

angivet i latitud och longitud i <strong>referenssystem</strong>et SWEREF 99.<br />

Nästa steg i analyskedjan omfattade transformation <strong>av</strong> de mätta punkternas positioner, från<br />

latitud och longitud i SWEREF 99 till x- och y-koordinater i <strong>Oskarshamns</strong> <strong>kommun</strong>s två<br />

lokala koordinatsystem. Transformationen gjordes med hjälp <strong>av</strong> Rix 95-parametrarna och var<br />

därmed en direktprojektion <strong>av</strong> positioner från ett geodetiskt globalt <strong>referenssystem</strong> till ett<br />

<strong>plan</strong>t lokalt system (se <strong>av</strong>snitt 2.2.2). Resultatet <strong>av</strong> direktprojiceringen var en uppsättning<br />

koordinatpar som kunde jämföras med de redan kända lokala koordinaterna för punkterna i<br />

stomnäten. De beräknade koordinatdifferenserna mellan kända och transformerade<br />

koordinater utgjorde de mätta stompunkternas restfel, vilka presenterades i såväl<br />

restfelsvektorkartor som restfelsvariationsbilder.<br />

För regionsystemet resulterade den <strong>plan</strong>a inpassningen som gjordes för <strong>av</strong> få fram restfelen i,<br />

ett grundmedelfel (<br />

0<br />

σ , 2D) på 48 millimeter, ett kvadratiskt medelfel (RMS, 2D) på 49<br />

millimeter och ett största restfel på 173 millimeter. I <strong>Oskarshamns</strong> lokala koordinatsystem<br />

blev restfelens grundmedelfel 47 millimeter, det kvadratiska medelfelet 47 millimeter och det<br />

största restfelet 124 millimeter. Passpunktsinmätningarna ledde alltså till att restfelens<br />

variation och spridning minskade något i regionsystemet, vilket var förväntat eftersom de nya<br />

inmätningarna inte hade gjorts i helt nya områden utan mest var förtätningar <strong>av</strong> tidigare<br />

mätningar (Rix 95-mätningarna). Att motsvarande mått inte minskade i <strong>Oskarshamns</strong> lokala<br />

koordinatsystem, berodde på att de nytillkomna restfelen här i många fall markant skilde sig<br />

23


från de tidigare relativt få restfelen från Rix 95-mätningarna. Det var nämligen ofta nya<br />

områden det hade gjorts mätningar i under mätetapp ett, så det blev inte på samma sätt som i<br />

regionsystemet fråga om en ren förtätning <strong>av</strong> Rix 95-mätningarna.<br />

Figur 11. Restfelsvektorer efter mätetapp ett i regionsystemet. Cirklarna används på samma sätt som i figur 7<br />

och 8.<br />

24


Figur 12. Restfelsvektorer efter mätetapp ett i <strong>Oskarshamns</strong> lokala koordinatsystem<br />

3.6.1 Regionsystemet<br />

I figur 11 syns restfelsvektorerna för de mätta punkterna fram till och med mätetapp ett i<br />

regionsystemet. Som figuren visar är vektorerna koncentrerade till ett antal något sånär<br />

friliggande områden. Detta är tydligare i figur 11 än det var i figur 7, som visade resultatet<br />

efter Rix 95-mätningarna. De friliggande områdena är de tätorter eller platser med samlad<br />

bebyggelse där stomnät finns utlagt. De enstaka restfelsvektorerna som syns mellan tätorterna<br />

har beräknats utifrån mätningar som gjordes i Rix 95-projektet, på punkter som antingen är<br />

rikstriangelpunkter eller punkter som tillkom i Rix 95-projektet; så kallade Rix 95-punkter.<br />

Eftersom <strong>kommun</strong>ens primärkarta endast täcker in den samlade bebyggelsen, var dessa<br />

enstaka restfelspunkter inte så intressanta i arbetet för <strong>referenssystem</strong>sbytet. Det som istället<br />

analyserades var restfelen på punkter i tätorterna. I de flesta fall hade stomnätet lagts ut, mätts<br />

25


in och utjämnats under en och samma mätinsats för hela tätorten, vilket innebar att den<br />

geometriska kvaliteten överlag var relativt god. Restfelen på punkter i samma tätort hade ofta<br />

liknande storlekar och riktningar, vilket g<strong>av</strong> upphov till tydliga trender i vektorkartorna (se<br />

vektorerna inom de blå cirklarna i figur 11). Detta var just ett <strong>av</strong> de förhållanden som<br />

analysen <strong>av</strong> Rix 95-mätningarna hade visat behövde undersökas närmare (se punkt 3 i <strong>av</strong>snitt<br />

3.1). Här och var fanns det dock restfel som stack ut och <strong>av</strong>vek från övriga restfel i den<br />

närmaste omgivningen. Dessa <strong>av</strong>vikelser fanns det också anledning att undersöka närmare.<br />

Restfelsvariationsbilder togs fram och i figur 13 syns restfelsvariationerna efter mätningarnas<br />

första etapp i regionsystemet. Utifrån i huvudsak figur 11 och 13 kunde nu följande<br />

sammanfattning <strong>av</strong> iakttagelser göras:<br />

1. De stora restfelen från analysen <strong>av</strong> Rix 95-mätningarna fanns kvar och <strong>av</strong>vek<br />

fortfarande från sin omgivning. Några <strong>av</strong>vikande restfel hade tillkommit. Se<br />

exempelvis de skarpt röda områdena i figur 13.<br />

2. Områdena med olika riktningar hos restfelsvektorerna hade blivit fler. Se de röda<br />

cirklarna (dubbelcirklar) i figur 11.<br />

3. Områdena med tydliga trender hos restfelsvektorerna hade blivit fler. Se de blå<br />

cirklarna i figur 11.<br />

4. De hittills gjorda mätningarna representerade och omringade väl stomnätets olika<br />

delar.<br />

5. Passpunkts<strong>av</strong>stånden var omkring 500 meter och borde inte försämra möjligheterna att<br />

räta upp den bristande stomnätsgeometrin.<br />

Av dessa iakttagelser var det framförallt de stora restfelen och områdena med spretande<br />

vektorer som behövde analyseras noggrannare. I restfelsvektorkartan var alltså de röda<br />

cirklarna (dubbelcirklar) mest intressanta och i restfelsvariationsbilderna de skarpt röda<br />

områdena. Med detta som utgångspunkt fortsatte analysen.<br />

26


Figur 13. Restfelsvariationer efter mätetapp ett i regionsystemet. Förstoringarna visar de två mest kritiska<br />

områdena; Misterhult i det övre utsnittet och Kristdala i det undre.<br />

I det övre utsnittet i figur 13 syns tydligt hur en punkt i det röda området har en stor<br />

restfelsvektor som pekar åt sydväst. De blå punkterna markerar punkternas inmätta lägen och<br />

vektorerna pekar mot punkternas lägen enligt de kända koordinaterna i regionsystemet. För de<br />

två närliggande punkterna åt nordväst är restfelen små och pekar åt väst och nord-nordväst. I<br />

ett begränsat område finns det alltså stora variationer hos restfelen. Visserligen ligger just<br />

dessa tre punkter utanför <strong>kommun</strong>ens primärkarteområde och behövdes inte för<br />

transformation <strong>av</strong> kartbilden, men punkten med det största restfelet var mätt i Rix 95-<br />

projektet, vilket gjorde det intressant att ändå titta lite närmare på punkterna. Den första<br />

analysen ledde därmed fram till ett beslut om en kontrollmätning <strong>av</strong> i första hand punkten<br />

med det stora restfelet. Med den punktens status säkert fastställd skulle det sedan gå att<br />

<strong>av</strong>göra om fler mätningar behövde göras på punkter runtomkring.<br />

27


För Kristdala, det andra området som krävde vidare analys och förmodligen fler mätningar,<br />

visade den första mätetappens resultat på flera oregelbundenheter i variationerna hos<br />

restfelen. I det undre utsnittet i figur 13 ses visserligen med lite god vilja dels en grupp<br />

restfelsvektorer med nord-sydlig riktning, och dels en grupp med öst-västlig riktning, men<br />

trenderna går in i varandra och dessutom finns några tydligt <strong>av</strong>vikande restfelsvektorer. Redan<br />

från början <strong>av</strong> analysen <strong>av</strong> området stod det klart att det var ett helt samhälles<br />

stomnätsdeformationer som behövde kartläggas ytterligare, men denna kompletterande<br />

mätinsats behövde dock någonting att ta <strong>av</strong>stamp ifrån. Istället för att helt oreflekterat<br />

genomföra förtätningsmätningar blev strategin att först försöka ta reda på om det som såg ut<br />

att vara <strong>av</strong>vikande restfel, verkligen var det, och därefter ta ställning till hur omfattande<br />

kompletteringsmätningarna behövde göras.<br />

Genom att mäta in ytterligare några stomnätspunkter i direkt anslutning till punkterna med<br />

stora och <strong>av</strong>vikande restfel, skulle det relativt snabbt och enkelt gå att <strong>av</strong>göra om dessa restfel<br />

trots <strong>av</strong>vikelserna skulle ingå i restfelsmodellen. Om de omgivande punkterna skulle uppvisa<br />

likartade restfel skulle det tyda på att det inte var fråga om någon isolerad <strong>av</strong>vikelse och att<br />

restfelet därmed skulle tas med i modellen. Omgivande punkter med restfel helt olika det<br />

restfel som skulle kontrolleras, skulle däremot visa att det verkligen handlade om en <strong>av</strong>vikelse<br />

som möjligtvis inte var representativ för den omgivande kartbilden, och därmed kanske inte<br />

skulle vara med och bygga upp restfelsmodellen. För att <strong>av</strong>göra om en rubbad stomnätspunkt<br />

skulle inkluderas i restfelsmodellen eller inte, skulle alltså en bedömning behöva göras <strong>av</strong> hur<br />

mycket punkten använts vid mätning efter att dess läge förändrats. Detta skulle naturligtvis<br />

vara nästan omöjligt att <strong>av</strong>göra utan någon kännedom om <strong>kommun</strong>ens mätningsverksamhet i<br />

ett historiskt perspektiv. Kunskapen hos den personal som jobbat länge med mätning i<br />

<strong>kommun</strong>en skulle därför komma väl till pass här.<br />

Kunskap om hur <strong>kommun</strong>ens stomnät byggts upp och använts var också väldigt viktig i<br />

bedömningen <strong>av</strong> var insatserna skulle koncentreras när det gällde den andra mätetappen och<br />

den vidare kartläggningen <strong>av</strong> det aktuella samhällets stomnätsdeformationer. Efter<br />

mätningarna i första mätetappen gjordes därför en efterforskning i <strong>kommun</strong>ens kart- och<br />

mätningsarkiv för att få information om hur tätortens stomnät lagts ut och allt eftersom byggts<br />

ut. För de andra mindre tätorterna hade stomnätet i de flesta fall etablerats vid ett enda<br />

tillfälle, men för det nu aktuella samhället hade stompunkterna lagts ut i tre olika etapper;<br />

visserligen under ett begränsat antal år men ändå i tre olika mätinsatser och <strong>av</strong> tre olika<br />

konsultfirmor. Möjligen kunde detta vara en förklaring till de lite större restfelsvariationerna i<br />

detta samhälle jämfört med övriga tätorter i <strong>kommun</strong>en. I vilket fall som helst kändes det<br />

angeläget att analysera och mäta extra grundligt i de delar <strong>av</strong> samhället där de gamla<br />

polygonnätskartorna visade att alla tre konsultfirmor hade varit och lagt ut stompunkter. Det<br />

var även i dessa områden de största oregelbundenheterna återfanns efter första etappens<br />

mätningar, även om det var svårt att få en fullständig bild <strong>av</strong> deformationerna utifrån det fåtal<br />

mätningar som gjordes då.<br />

Sammantaget resulterade alltså analysen efter första mätetappens mätningar i det andra <strong>av</strong> de<br />

två kritiska områdena i figur 13 i, att några punkter skulle mätas in i direkt anslutning till de<br />

två till tre punkter som hade mest <strong>av</strong>vikande restfel, samt att en generell förtätning <strong>av</strong><br />

mätningarna i samhället i övrigt skulle genomföras. I den förtätningen skulle särskild<br />

koncentration läggas på de områden i samhällets centrala delar där stompunkterna lagts ut <strong>av</strong><br />

olika mätfirmor vid olika tillfällen.<br />

28


I övriga tätorter som omfattades <strong>av</strong> regionsystemet var som sagt överensstämmelsen mellan<br />

restfelsvektorerna relativt god samtidigt som de inmätta passpunkterna låg så pass tätt att<br />

stomnätets geometrier därmed fick anses vara kartlagda. Den enda ytterligare mätinsats som<br />

ansågs behövas förutom de två redan nämnda, var därför en inmätning <strong>av</strong> ett par punkter i<br />

Fårbo och Figeholm, de två tätorterna söder om det första <strong>av</strong> de två uppförstorade röda<br />

områdena i figur 13. Restfelsvariationsbilden ger intrycket att det är ett område där<br />

restfelsvektorerna spretar ganska kraftigt åt två olika håll. I själva verket handlar det om två<br />

tätorter vars stomnät byggts upp vid olika tillfällen och därmed uppvisar lite olika geometrier.<br />

För att göra <strong>av</strong>gränsningen mellan de två stomnätsdelarna lite tydligare och kontrollera<br />

<strong>av</strong>vikande punkter, skulle därför ett par passpunktsinmätningar ytterligare vara lämpliga.<br />

3.6.2 <strong>Oskarshamns</strong> lokala koordinatsystem<br />

För det andra <strong>av</strong> de två lokala koordinatsystemen i <strong>kommun</strong>en blev angreppssättet för den<br />

fortsatta analysen lite annorlunda. Detta koordinatsystem täckte, som redan nämnts, in<br />

<strong>Oskarshamns</strong> stad och tätorten söderut, Påskall<strong>av</strong>ik. Mellan Oskarshamn och Påskall<strong>av</strong>ik<br />

finns ett område där <strong>kommun</strong>en inte har någon primärkarta och alltså inte heller något utbyggt<br />

stomnät. Stomnätet i Oskarshamn stad är det stomnät som <strong>av</strong> naturliga skäl genomgått de<br />

största förändringarna under årens lopp. Det har gjorts förtätningar och utbyggnader i flera<br />

etapper, små nät har slagits ihop till större och så vidare. För att kunna göra en vettig analys i<br />

denna del <strong>av</strong> <strong>Oskarshamns</strong> lokala koordinatsystem krävdes därför områdesvisa bedömningar<br />

samtidigt som en helhetsgranskning inte kunde utelämnas helt. För tätbebyggelsen kring och i<br />

Påskall<strong>av</strong>ik kunde däremot analysen göras på motsvarande sätt som för tätorterna i<br />

regionsystemet. Här fanns nämligen inga direkt kritiska områden men passpunkterna var<br />

alltför få, så förtätningsmätningar behövde göras.<br />

Figur 14. Restfelsvariationer efter mätetapp ett i<br />

<strong>Oskarshamns</strong> lokala koordinatsystem<br />

Figur 15. Restfelsvariationer efter mätetapp ett i<br />

<strong>Oskarshamns</strong> lokala koordinatsystem – Absoluta<br />

differenser<br />

29


Analysen tog sin utgångspunkt från restfelsvektorkartorna och restfelsvariationsbilderna<br />

tillsammans med <strong>kommun</strong>ens primärkarta, för att få en överblick över var de förväntade<br />

respektive icke förväntade större restfelsvariationerna påträffades. Trots att vektorerna i<br />

variationskartorna (se figur 14 och 15) överlag ser ut att spreta väldigt mycket, var det<br />

nämligen möjligt att med hjälp <strong>av</strong> de olika kartorna tillsammans, i analysen delvis <strong>av</strong>gränsa<br />

stomnätsdelar och områden där trender hos restfelen kunde ses eller borde kunna ses. Detta<br />

var möjligt eftersom primärkartan g<strong>av</strong> information om stadsbilden med dess uppdelning i<br />

exploaterade och icke exploaterade områden samt om stomnätets utsträckning, medan<br />

analysverktygen i form <strong>av</strong> restfelsredovisningar talade om var i stadsbilden<br />

passpunktsinmätningar hade gjorts och var spänningar och deformationer i stomnätet verkade<br />

finnas (se figur 16). Den bristande överensstämmelsen mellan restfelen i olika närliggande<br />

områden kunde därmed i många fall enkelt förstås utan närmare vetskap om stomnätets<br />

uppbyggnad historiskt sett, genom en titt i analysmaterialet där orsakerna till olikheterna<br />

kunde ses. Exempelvis kunde det mellan närliggande områden finnas naturliga barriärer;<br />

kanske i form <strong>av</strong> skogspartier eller bergshöjder, som gjorde det omöjligt att ansluta områdena<br />

till varandra på ett tillfredsställande sätt.<br />

Figur 16. Exempel på hur olika kartor och bilder kunde kombineras i restfelsanalysen. En transparent<br />

georefererad restfelsvariationsbild kombineras här med en karta som bland annat visar <strong>kommun</strong>ala<br />

stomnätspunkter (de röda punkterna).<br />

30


Den historiska bakgrunden till stomnätets uppbyggnad var dock viktig att ta i beaktande i<br />

analysarbetet. För det första behövde de slutsatser som dragits utifrån enbart mätresultaten<br />

stämmas <strong>av</strong> mot den historiska kunskap om stomnätet och dess brister, som fanns hos dem<br />

som jobbat länge med geodetisk mätning i <strong>kommun</strong>en. Den andra anledningen till att det<br />

historiska perspektivet behövde lyftas fram var att detta kanske kunde ge svar på förhållanden<br />

som inte upptäckts genom de dittills gjorda mätningarna. Med anledning <strong>av</strong> detta bjöds därför<br />

<strong>kommun</strong>ens förre kart- och mätchef in för att ge sin bild <strong>av</strong> stomnätens kvalitet i <strong>kommun</strong>en<br />

och även utifrån denna historiska kunskap föreslå var de vidare mätinsatserna borde<br />

koncentreras. I stora drag stämde analysresultaten och den historiska kunskapen väl överens. I<br />

de områden där sämre och motsägelsefulla mätresultat erhållits under årens lopp, visade<br />

analysen också ofta på spänningar i stomnätet. Det fanns dock exempel på brister som<br />

analysen fångat in men som inte tidigare upptäckts genom mätning med stompunkter som<br />

utgångspunkter. Likaså kom det på motsvarande sätt fram sådana värdefulla historiska fakta<br />

om stomnätskvaliteten och stomnätens uppbyggnad som analysen inte kunnat tala om.<br />

Analysen <strong>av</strong> det andra <strong>av</strong> de två lokala stomnäten hade därmed så långt resulterat i följande<br />

iakttagelser och bedömningar:<br />

1. Några tydligt <strong>av</strong>vikande restfel kunde hittas.<br />

2. Det var svårt att i vektorkartor och variationsbilder identifiera trender hos<br />

restfelsvektorerna och områden där trender borde finnas. Det var även svårt att hitta<br />

områden med uppenbart olika riktningar hos närliggande restfel. Restfelsgrafiken<br />

behövde kombineras med <strong>kommun</strong>ala primärkartan för att olika stomnätsdelar skulle<br />

kunna <strong>av</strong>gränsas och eventuella trender hittas. I figur 12 visar de blå cirklarna dock på<br />

sådant som skulle kunna vara trender eller delar <strong>av</strong> trender hos restfelens riktningar.<br />

Den röda cirkeln (dubbelcirkel) markerar ett delområde där restfelsvektorer i samma<br />

stomnätsdel uppvisar olika riktningar och storlekar.<br />

3. En del stomnätsdelar var sämre representerade än andra.<br />

4. Trender hos restfelen i tätorten Påskall<strong>av</strong>ik kunde ses men fler passpunkter med restfel<br />

behövdes.<br />

5. Historiska bakgrunden till stomnätsuppbyggnaden var till hjälp i arbetet med att<br />

<strong>plan</strong>era fortsatta mätningsåtgärder.<br />

Slutsatsen från analys och insamlande <strong>av</strong> historisk stomnätskunskap blev att en i det närmaste<br />

generell förtätning <strong>av</strong> första mätetappens mätningar i <strong>Oskarshamns</strong> stad var nödvändig.<br />

Eftersom stomnätet genomgått så många och olika förändringar under årens lopp, och därför<br />

uppvisade spänningar och deformationer <strong>av</strong> olika slag i olika delar, krävdes dock att de<br />

områdesvisa trender eller antydningar till trender hos restfelsvektorerna som kunnat påvisas,<br />

också i förtätningsmätningarna behandlades var för sig. Det var även viktigt att i den fortsatta<br />

analysen få utrett om det som <strong>av</strong> restfelsvariationsbilderna att döma verkade vara <strong>av</strong>vikande<br />

restfel verkligen var det. En del extra mätningar för att klargöra detta bedömdes alltså<br />

behövas. I Påskall<strong>av</strong>ik behövdes flera mätningar för att identifieringen <strong>av</strong> trender bland<br />

restfelen där skulle kunna göras på ett tillförlitligt sätt.<br />

31


3.7 Planering <strong>av</strong> andra mätetappens mätningar<br />

Valet <strong>av</strong> passpunkter som skulle mätas in i andra mätetappens mätningar styrdes givetvis i<br />

mångt och mycket <strong>av</strong> det analysen <strong>av</strong> första mätetappens mätningar gett vid handen. Bortsett<br />

från det var tankegången dock precis densamma som tidigare när det gäller kr<strong>av</strong> på goda<br />

mätförhållanden, tillförlitliga mätresultat och bra möjligheter till optimal triangelbildning vid<br />

uppbyggandet <strong>av</strong> restfelsmodellen.<br />

Mätningarna i andra mätetappen skulle göras <strong>av</strong> olika anledningar i olika områden. I vissa<br />

områden behövdes mätningar för att <strong>av</strong>vikande restfel skulle kontrolleras, medan en allmän<br />

förtätning <strong>av</strong> passpunktsinmätningarna var nödvändig i andra delar <strong>av</strong> stomnäten där<br />

spänningarna ännu inte var tillräckligt kartlagda. När det gäller <strong>plan</strong>ering <strong>av</strong> nya passpunkter<br />

var det mest i det senare fallet sådan var nödvändig eftersom de <strong>av</strong>vikande restfelen kunde<br />

kontrolleras ganska lätt med hjälp <strong>av</strong> kontrollmätningar eller kompletteringsmätningar på<br />

första bästa punkter i direkt anslutning till punkten med det <strong>av</strong>vikande restfelet. I områdena i<br />

Oskarshamn stad där passpunktsinmätningarna skulle förtätas generellt <strong>plan</strong>erades<br />

passpunkterna så att passpunkts<strong>av</strong>ståndet reducerades till mellan 500 och 800 meter. Däremot<br />

blev passpunkts<strong>av</strong>stånden inte lika långa i Kristdala; den tätort i regionsystemet där en<br />

förtätning <strong>av</strong> passpunktsinmätningarna var nödvändig. Det viktiga var dock återigen inte att<br />

<strong>av</strong>ståndet mellan passpunkterna var detsamma överallt, utan att den bristande<br />

stomnätsgeometrin kunde rätas upp.<br />

3.8 Mätningarnas andra etapp<br />

I mätetapp två genomfördes mätningarna nästan uteslutande på samma sätt som tidigare.<br />

Sånär som på vid något enstaka mättillfälle kunde GNSS-mätning göras, och när<br />

satellittillgängligheten inte var tillräcklig togs traditionell mätteknik i bruk för att mäta in<br />

passpunkterna excentriskt. Sammanlagt mättes 62 passpunkter in i andra mätetappen. Av<br />

dessa var 49 stompunkter i <strong>Oskarshamns</strong> stad och tätorten Påskall<strong>av</strong>ik, det vill säga punkter i<br />

<strong>Oskarshamns</strong> lokala koordinatsystem. Resterande inmätningar gjordes i några <strong>av</strong> tätorterna i<br />

regionsystemet, och då främst i Kristdala, där en generell förtätning behövdes (se <strong>av</strong>snitt<br />

3.6.1). I <strong>Oskarshamns</strong> lokala koordinatsystem hade därmed totalt 105 passpunkter mätts in<br />

och i regionsystemet 52 passpunkter.<br />

3.9 Analysens slutfas<br />

Analysen efter den andra mätetappen följde samma arbetsgång som vid analysen efter de<br />

första passpunktsinmätningarna. Medeltalsbildningen g<strong>av</strong> koordinater i SWEREF 99 och efter<br />

transformation med hjälp <strong>av</strong> Rix 95-parametrarna erhölls en uppsättning koordinater som<br />

kunde jämföras med de redan kända koordinaterna i de två gamla lokala koordinatsystemen.<br />

Den inpassning i <strong>plan</strong> som gjordes visade att restfelen i regionsystemet hade ett<br />

grundmedelfel (<br />

0<br />

σ , 2D) på 40 millimeter, ett kvadratiskt medelfel (RMS, 2D) på 43<br />

millimeter och ett största restfel på 95 millimeter. För <strong>Oskarshamns</strong> lokala koordinatsystem<br />

g<strong>av</strong> motsvarande inpassning grundmedelfelet 44 millimeter, kvadratiska medelfelet 44<br />

millimeter och största restfelet 124 millimeter. Jämfört med resultatet från mätningarna i<br />

mätetapp ett minskade alltså spridnings- och variationsmåtten i såväl regionsystemet som<br />

<strong>Oskarshamns</strong> lokala koordinatsystem (se tabell 2 och 3). Detta var ett troligt resultat eftersom<br />

mätningarna för båda koordinatsystemen under mätetapp två företrädesvis var<br />

förtätningsmätningar och inte mätningar i nya stomnätsdelar.<br />

32


Tabell 2. Passpunktsinmätningar gjorda i regionsystemet<br />

Rix95 Mätetapp 1 Mätetapp 2<br />

Antal punkter 37 39 13<br />

Gundmedelfel, σ<br />

0<br />

[mm] 51 48 40<br />

Kvadratiskt medelfel,<br />

RMS, 2D [mm]<br />

51 49 43<br />

Största restfel [mm] 173 173 95<br />

Tabell 3. Passpunktsinmätningar gjorda i <strong>Oskarshamns</strong> lokala koordinatsystem<br />

Rix95 Mätetapp 1 Mätetapp 2<br />

Antal punkter 24 56 49<br />

Gundmedelfel, σ<br />

0<br />

[mm] 41 47 44<br />

Kvadratiskt medelfel,<br />

RMS, 2D [mm]<br />

41 47 44<br />

Största restfel [mm] 124 124 124<br />

Nya restfelsvektorkartor (se figur 17 och 18) och restfelsvariationsbilder (se kommande<br />

<strong>av</strong>snitt) togs fram och det blev möjligt att analysera hur andra mätetappens mätningar hade<br />

förändrat bilden och förhoppningsvis resulterat i en bättre kartläggning <strong>av</strong> stomnätets<br />

geometri.<br />

33


Figur 17. Restfelsvektorer efter mätetapp två i regionsystemet<br />

34


Figur 18. Restfelsvektorer efter mätetapp två i <strong>Oskarshamns</strong> lokala koordinatsystem<br />

3.9.1 Regionsystemet<br />

I restfelsvariationsbilderna efter etapp två <strong>av</strong> mätningarna hade de kritiska områdena kunnat<br />

reduceras tack vare att en del <strong>av</strong>vikande punkter kunnat tas bort och förtätningsmätningarna<br />

åstadkommit en bättre kartläggning <strong>av</strong> stomnätets deformationer.<br />

35


Figur 19. Restfelsvariationer efter mätetapp två i regionsystemet. Förstoringen visar det mest kritiska området;<br />

Kristdala.<br />

Figur 20. Restfelsvariationer efter mätetapp två i <strong>Oskarshamns</strong> lokala koordinatsystem – Absoluta differenser<br />

36


Figurerna 19 och 20 visar restfelsvariationerna efter mätetapp två i regionsystemet. I figur 19<br />

har värdet på variationerna viktats omvänt proportionellt mot respektive passpunkts<strong>av</strong>stånd,<br />

vilket är det som oftast görs och som använts hittills i den här framställningen, medan figur 20<br />

presenterar variationerna utifrån de absoluta differenserna (se punkt 3 i <strong>av</strong>snitt 2.3.2). Att<br />

även de absoluta differenserna användes i analysen, berodde på att detta ofta är till god hjälp<br />

och på ett bra sätt kan komplettera analysen utifrån det normala presentationssättet. Ibland<br />

redovisas restfelsvariationerna nämligen alltför kraftigt eller alltför svagt när de viktas med<br />

hänsyn till passpunkts<strong>av</strong>ståndet. Exempelvis kan ett stort passpunkts<strong>av</strong>stånd göra att en stor<br />

variation redovisas alltför svagt, medan små variationer kan förstärkas obefogat mycket på<br />

grund <strong>av</strong> ett litet passpunkts<strong>av</strong>stånd. I sådana fall är en redovisning som använder de absoluta<br />

differenserna värdefull att stämma <strong>av</strong> mot.<br />

Den mest markanta skillnaden mellan de två olika presentationssättet i detta fall är att en<br />

punkt söder om de två tätorterna Fårbo och Figeholm i mitten österut, har ett kraftigare rött<br />

område kring sig i figur 20 än i figur 19. Denna punkt var en <strong>av</strong> de punkter som stack ut något<br />

efter mätningarna i mätetapp ett. I mätetapp två var tanken att punkten skulle kontrollmätas<br />

och att ytterligare någon punkt i nära anslutning skulle mätas in, för att möjliggöra jämförelse<br />

<strong>av</strong> restfelen. Kontrollmätningen <strong>av</strong> punkten resulterade i samma restfel som tidigare, och<br />

några punkter i direkt anslutning att mäta in med satellitteknik fanns inte. Eftersom punkten<br />

ligger helt i utkanten <strong>av</strong> stomnätet, har en stabil markering och inte kunde antas ha rubbats ur<br />

sitt ursprungliga läge, beslutades därför att punktens restfel kunde tas med i den kommande<br />

restfelsmodellen. För övrigt var stomnätets deformationer i dessa två tätorter kartlagda.<br />

Norr om Fårbo och Figeholm fanns efter första mätetappen ett restfel som också <strong>av</strong>vek från<br />

omgivningen. Efter en enklare kontrollmätning med samma resultat som tidigare, kunde<br />

slutsatsen dras att det rör i mark som markerar punkten, under årens lopp, mer och mer tryckts<br />

ut och bort från den väg som det är nerslaget intill. Detta var ett förhållande som under<br />

mätningarnas gång även kunde konstateras på liknande sätt för andra punkter vid äldre eller<br />

hårt trafikerade vägar. Den stora restfelsvektorn kunde alltså inte anses vara representativ för<br />

de omgivande restfelsvektorerna och utelämnades därför i den fortsatta processen för att ta<br />

fram en bra restfelsmodell för användning vid transformationen. Detta kunde göras utan<br />

vidare just eftersom inga kartdetaljer i primärkartan mätts in med utgångspunkt från den<br />

punkten. Samtidigt med att denna punkt utelämnades, togs även ett beslut att inte heller ta<br />

med de två andra punkterna som låg i direkt anslutning. Anledningen till det var att dessa<br />

punkter tillhörde en grupp punkter i stomnätet som dels låg utanför tätbebyggelse och därför<br />

inte hade använts som utgångspunkter vid inmätning <strong>av</strong> kartdetaljer, och dels inte hade<br />

tillkommit på ett sådant sätt att deras lägen var helt tillförlitliga.<br />

Den tredje och sista tätorten som kompletteringsmätningar i mätetapp två hade gjorts i var<br />

Kristdala. Analysen utifrån restfelsvariationsbilderna för denna tätort visade att mätningarna i<br />

mätetapp två hade gjort att de trender i nord-sydlig respektive öst-västlig riktning som kunde<br />

ses i analysmaterialet efter första mätetappens mätningar, blev förstärkta och tydligare.<br />

Fortfarande spretade restfelen i flera <strong>av</strong>seenden men det förhållandevis stora antalet<br />

passpunktsinmätningar gjorde att passpunkts<strong>av</strong>stånden var korta och kartläggningen <strong>av</strong><br />

stomnätsgeometrin mycket bättre än tidigare. Dessutom hade de punkter som efter första<br />

mätetappens mätningar hade restfel som <strong>av</strong>vek från de närliggande restfelen, kontrollerats och<br />

i vissa fall kunnat tas bort, medan de i andra fall fick vara kvar.<br />

37


En sammanfattning <strong>av</strong> resultaten från analysen <strong>av</strong> mätningarnas andra etapp i regionsystemet,<br />

ser ut så här:<br />

1. Några <strong>av</strong> de stora och <strong>av</strong>vikande restfelen hade kunnat tas bort eftersom de inte var<br />

representativa för övriga restfel i omgivningen, eller för att de inte skulle fylla någon<br />

funktion vid transformationen <strong>av</strong> primärkartan.<br />

2. Områden med restfel i olika riktningar fanns kvar, men de var bättre kartlagda nu, med<br />

kortare passpunkts<strong>av</strong>stånd och <strong>av</strong>vikande restfel kontrollerade.<br />

3. Tack vare fler mätningar hade trenderna som tidigare kunnat skönjas, nu blivit<br />

förstärkta och tydligare.<br />

3.9.2 <strong>Oskarshamns</strong> lokala koordinatsystem<br />

Figur 21. Restfelsvariationer efter mätetapp två i<br />

<strong>Oskarshamns</strong> lokala koordinatsystem<br />

Figur 22. Restfelsvariationer efter mätetapp två i<br />

<strong>Oskarshamns</strong> lokala koordinatsystem – Absoluta<br />

differenser<br />

Resultatet <strong>av</strong> passpunktsinmätningarna i mätetapp två i <strong>Oskarshamns</strong> lokala koordinatsystem<br />

kan ses i figurerna 21 och 22. Den stora skillnaden jämfört med restfelsvariationsbilderna<br />

efter första mätetappen, var utan tvekan antalet restfelsvektorer. Drygt 70 procent <strong>av</strong><br />

mätningarna i mätetapp två gjordes i <strong>Oskarshamns</strong> lokala koordinatsystem, vilket innebar en<br />

kraftig förtätning <strong>av</strong> inmätningarna. Denna generella förtätning var också huvudmålet med<br />

andra mätetappens mätningar i detta koordinatsystem. Stomnätet hade under lång tid<br />

genomgått så många förändringar <strong>av</strong> olika slag att en noggrann kartläggning genom en stor<br />

mängd mätningar var nödvändig. Efter andra mätetappens mätningar fanns de röda områdena<br />

i restfelsvariationsbilderna kvar men de hade blivit mindre och något fler. En del <strong>av</strong> dem dök<br />

upp först i mätningarnas andra etapp, vilket var naturligt med tanke på att de nya mätningarna<br />

38


fångade in nya spänningar och deformationer. En jämförelse mellan första och andra<br />

mätetappens restfelsvariationsbilder med absoluta differenser (jämför figur 22 med figur 15)<br />

visar dock att bilderna är ljusare för andra mätetappen än för första, i de områden där de flesta<br />

<strong>av</strong> andra mätetappens inmätningar gjordes. Restfelens variationer minskade alltså när nätet<br />

förtätades kraftigt, vilket tydde på att kartläggningen <strong>av</strong> spänningarna hade förbättrats.<br />

Fortfarande var det här och var restfel som stack ut, men de <strong>av</strong>vikande punkterna som låg<br />

inom <strong>kommun</strong>ens primärkarteområden hade kontrollerats och deras lägen var säkerställda.<br />

Anledningen till att de var kvar även efter andra mätetappens mätningar berodde därför på att<br />

de bedömdes vara viktiga punkter i kartläggningen <strong>av</strong> den bristande stomnätsgeometrin och<br />

därmed borde ingå i restfelsmodellen. Avvikande restfel utanför de områden som omfattades<br />

<strong>av</strong> primärkartan fick vara kvar i analysprocessen eftersom de inte störde som de nu såg ut. I<br />

vissa fall kunde de istället vara till hjälp i arbetet med att hitta trender eller grupper <strong>av</strong> restfel<br />

med liknande restfel.<br />

För tätorten Påskall<strong>av</strong>ik söder om Oskarshamn hade genom andra mätetappens mätningar en<br />

ren kombinerad kontroll- och förtätningsmätning gjorts. Det var ingen större mängd punkter<br />

som hade mätts in, men de mätningar som hade gjorts var till hjälp i bedömningen <strong>av</strong><br />

<strong>av</strong>vikande restfel samt eventuella mönster hos restfelens storlek och utseende.<br />

Analysen efter andra mätetappens mätningar i <strong>Oskarshamns</strong> lokala koordinatsystem, kunde<br />

sammanfattas på följande sätt.<br />

1. Avvikande restfel hade kontrollerats och en bedömning om de skulle vara kvar eller<br />

inte, hade gjorts.<br />

2. Såväl trender som olika riktningar hos restfelen hade blivit lättare att identifiera.<br />

3. Mätningar på ett stort antal nya passpunkter hade gjorts, vilket resulterat i att<br />

stomnätets olika delar var bättre representerade nu än tidigare, och att<br />

passpunkts<strong>av</strong>stånden förkortats till en mer tillfredsställande storlek. I<br />

restfelsvariationsbilderna kunde den förbättrade kartläggningen ses genom att bilderna<br />

var ljusare i områden där många passpunkter mätts in.<br />

3.10 Upprättande <strong>av</strong> restfelsmodeller<br />

Efter analysen <strong>av</strong> mätetapp två togs beslutet att inga fler passpunktsinmätningar behövde<br />

göras. Kartläggningen <strong>av</strong> stomnätens geometriska brister ansågs alltså vara färdig. Ett sådant<br />

beslut fattas någon gång under processen när tillräckligt många förtätningsmätningar har<br />

gjorts, de <strong>av</strong>vikande restfelen undersökts och passpunkternas lägen är sådana att en bra<br />

triangelbildning kan fås till stånd. Ett perfekt resultat är omöjligt att få och det är för de flesta<br />

även en fråga om hur mycket tid och resurser som kan läggas ner i projektet, men om<br />

passpunkts<strong>av</strong>stånden inte är för långa och tydliga trender och mönster hos restfelen har<br />

analyserats fram i stomnätens olika delar, så kan beslutet att <strong>av</strong>sluta inmätningarna tas.<br />

Tillsammans med passpunkterna från Rix 95-projektet hade sammanlagt 86 passpunkter mätts<br />

in i regionsystemet och 134 passpunkter i <strong>Oskarshamns</strong> lokala system. Projektet för byte <strong>av</strong><br />

<strong>referenssystem</strong> i <strong>Oskarshamns</strong> <strong>kommun</strong> kom alltså att sammanlagt innebära användning <strong>av</strong><br />

220 <strong>av</strong> <strong>kommun</strong>ens stompunkter, vilka utgjorde drygt 5,5 procent <strong>av</strong> totala antalet<br />

stompunkter.<br />

39


Nästa steg i arbetet med att ta fram transformationssamband var därmed upprättandet <strong>av</strong><br />

restfelsmodeller för de båda <strong>kommun</strong>ala stomnäten. För det ändamålet användes modulen<br />

TRIAD i Lantmäteriets transformationsverktyg Gtrans. TRIAD tar passpunkternas restfel som<br />

indata och bygger upp en restfelsmodell genom att generera en så kallad Delaunaytriangulering<br />

inom det konvexa höljet för alla passpunkterna (se <strong>av</strong>snitt 2.3.1).<br />

Figur 23. Samtliga restfelsvektorer i regionsystemet<br />

Figur 23 och 24 innehåller restfelsvektorkarta samt restfelsmodell för regionsystemet.<br />

Eftersom vektorkartan i princip bara innehåller vektorer samlade till ett antal ”öar” i<br />

<strong>kommun</strong>en, blir trianguleringen, som utgör restfelsmodell i figur 24, inte särskilt tydlig i varje<br />

enskild ”ö”. Det är dock möjligt att delvis bedöma trianglarnas lämplighet för interpolering. I<br />

restfelsmodellen ses bland annat att de flesta ”platta” trianglarna förekommer mellan<br />

vektorkoncentrationerna i regionsystemet, och eftersom <strong>kommun</strong>en inte har någon<br />

primärkarta där, har detta ingen betydelse. Inne i tätorterna blir trianguleringen överlag helt<br />

acceptabel, även om ett och annat mindre lämpligt utseende på trianglarna förekommer på<br />

sina håll. De olämpliga trianglarna påverkar resultatet men i områden där restfelen är små<br />

eller väl korrelerade, blir effekten begränsad och oftast obetydlig.<br />

40


Kartdetaljer som låg utanför de punkter som skulle omringa stomnätsdelarna, hamnade i<br />

trianglar där en eller ibland till och med två hörnpasspunkter låg långt ifrån den aktuella<br />

kartdetaljen. Detta var naturligtvis inte helt bra men gjorde i de flesta fall inte så mycket,<br />

eftersom interpoleringen utförs linjärt i triangeln och därmed låter de passpunkter som ligger<br />

närmast kartdetaljen i fråga, påverka mest.<br />

Figur 24. Restfelsmodell utifrån samtliga restfel i regionsystemet<br />

41


Figur 25. Samtliga restfelsvektorer i <strong>Oskarshamns</strong> lokala koordinatsystem<br />

42


Figur 26. Restfelsmodell utifrån samtliga restfel i <strong>Oskarshamns</strong> lokala koordinatsystem<br />

43


I figurerna 25 och 26 presenteras en restfelsvektorkarta och en restfelsmodell för<br />

<strong>Oskarshamns</strong> lokala koordinatsystem. Även i detta system g<strong>av</strong> restfelsvektorernas ibland<br />

oregelbundna geometri upphov till en och annan ”platt” triangel också i primärkarteområden.<br />

Tillsammans med motsättningar på grund <strong>av</strong> låg inbördes noggrannhet mellan<br />

restfelsvektorerna kunde det vara berättigat att i vissa områden ifrågasätta<br />

restfelsinterpolering utifrån en restfelsmodell skapad i TRIAD. Generellt sett fick<br />

trianguleringen dock betraktas som fullt godkänd i de delar där den skulle användas, vilket<br />

ledde till beslutet att, även för detta koordinatsystem inkludera restfelsmodellen vid<br />

framtagandet <strong>av</strong> transformationssambanden.<br />

3.11 Framtagande <strong>av</strong> transformationssamband<br />

Ett transformationssamband som involverar restfelsinterpolering består <strong>av</strong> två delar; dels<br />

transformationsparametrar och dels en restfelsmodell. För <strong>Oskarshamns</strong> <strong>kommun</strong> togs två<br />

sådana samband fram; ett för vart och ett <strong>av</strong> de två gamla koordinatsystemen. De<br />

transformationsparametrar som användes var de som Lantmäteriet tagit fram i Rix 95-<br />

projektet. Dessa parametrar var, som redan nämnts, behäftade med stora brister på grund <strong>av</strong><br />

att underlaget innehöll alltför få passpunkter. Otillräckligheterna som fanns skulle dock<br />

kompenseras med en bra restfelsmodell. Om transformationsparametrarna hade varit bättre<br />

hade däremot kr<strong>av</strong>en på restfelsmodellen inte behövt ställas lika högt. Så länge slutresultatet<br />

blir bra spelar det alltså inte så stor roll på vilken <strong>av</strong> de två ingående delarna i sambandet det<br />

läggs mest resurser för att få delen bra.<br />

Transformationssambanden togs fram med hjälp <strong>av</strong> Gtrans. Rix 95-parametrarna för<br />

respektive lokalt system kopplades ihop med motsvarande restfelsmodell, så att en<br />

transformationsfil som både rätade upp och transformerade till SWEREF 99 16 30 skapades.<br />

En transformation <strong>av</strong> geografiska data i exempelvis regionsystemet innebar med de nya<br />

sambanden därmed, att informationen först rätades upp, sedan transformerades till SWEREF<br />

99 och slutligen räknades över till kartprojektionen SWEREF 99 16 30. För vart och ett <strong>av</strong><br />

sambanden skapades även dess inverterade motsvarighet, som skulle kunna användas för att<br />

transformera från SWEREF 99 16 30 tillbaka till de gamla koordinatsystemen. Det var viktigt<br />

att denna möjlighet fanns eftersom den <strong>kommun</strong>ala kartverksamheten ofta innebär att<br />

information insamlad och redovisad i olika koordinatsystem ska kombineras och redovisas i<br />

samma koordinatsystem, eller att ny information i nya koordinatsystem ska tillföras och<br />

användas i äldre projekt där koordinatsystemet är gammalt.<br />

3.12 Test och kontroll <strong>av</strong> transformationssamband<br />

Efter att transformationssambanden var skapade var det dags att använda och kontrollera dem<br />

i en transformation <strong>av</strong> kartdata. En enkel metod för att snabbt skapa påhittade punkter som<br />

kan ersätta riktiga kartdetaljer under exempelvis ett testförfarande, genererar ett regelbundet<br />

rutnät <strong>av</strong> punkter över det aktuella transformationsområdet (Svanholm, 2000). Denna metod<br />

användes också i Oskarshamn, med en punkttäthet på 1000 meter i regionsystemet och 500<br />

meter i <strong>Oskarshamns</strong> lokala koordinatsystem. När rutnäten för de båda koordinatsystemen<br />

genererats täckte de sina respektive tillämpningsområden plus ytterligare en bit <strong>av</strong><br />

omgivningen. Punkterna i rutnäten transformerades med hjälp <strong>av</strong> de transformationssamband<br />

som tagits fram, och erhöll då var och en ett interpolerat restfel från restfelsmodellen. Nästa<br />

steg var att återigen transformera det genererade rutnätet, fast denna gång utelämna<br />

restfelsinterpoleringen i transformationen. På detta sätt togs två punktuppsättningar fram där<br />

skillnaden mellan uppsättningarna utgjordes <strong>av</strong> de interpolerade restfelen i varje punkt.<br />

44


Genom att beräkna koordinatdifferenserna och låta dessa till storlek och riktning representeras<br />

<strong>av</strong> vektorer i en vektorkarta, var det möjligt att åskådliggöra restfelsinterpoleringens påverkan<br />

på varje punkt. Det som syns i figur 27 och figur 30 kan alltså sägas vara skillnaden mellan att<br />

transformera med och utan restfelsmodell.<br />

Figur 27. Differensvektorer hos punkter i ett rutnät med punkttätheten 1000<br />

meter i regionsystemet. Vektorerna representerar skillnaden i resultat mellan<br />

transformation med och utan restfelsinterpoleringen.<br />

45


Figur 28. Restfelsmodell för regionsystemet med passpunkternas<br />

restfel redovisade<br />

Figur 29. Del <strong>av</strong> regionsystemets tillämpningsområde. De tunnare<br />

differensvektorerna har ett tydligt samband med de kraftigare<br />

restfelsvektorerna för passpunkterna i trianglarnas hörn.<br />

46


Genom att jämföra de interpolerade restfelen på de påhittade punkterna i rutnäten, med<br />

passpunkternas restfel i restfelsmodellerna (se figurerna 28 och 31), gick det att göra en grov<br />

förhandsbedömning <strong>av</strong> resultatet och den synbara kvaliteten i den kommande<br />

primärkartetransformationen i <strong>Oskarshamns</strong> <strong>kommun</strong>. Om en transformation med<br />

restfelsinterpolering går som <strong>plan</strong>erat tilldelas kartdetaljerna interpolerade restfel vars<br />

storlekar och riktningar har stora likheter med restfelen hos passpunkterna som bygger upp de<br />

trianglar i restfelsmodellen som kartdetaljerna finns i. Dessa likheter mellan de påhittade<br />

punkternas restfel och passpunkternas restfel kunde skönjas i figurerna 27 och 28 respektive<br />

30 och 31 men helt tydligt blev det först i en figur med större skala och där<br />

differensvektorerna kombinerades med restfelsmodellen. Ett exempel på hur detta kunde se ut<br />

visas i figur 29, där en liten del <strong>av</strong> regionsystemets tillämpningsområde åskådliggörs. De<br />

kraftigare vektorerna är restfelen för passpunkterna, medan de tunnare representerar restfelen<br />

hos de påhittade kartdetaljerna.<br />

Figur 30. Differensvektorer hos punkter i ett rutnät<br />

med punkttätheten 500 meter i <strong>Oskarshamns</strong> lokala<br />

koordinatsystem. Vektorerna representerar skillnaden i<br />

resultat mellan transformation med och utan<br />

restfelsinterpoleringen.<br />

47


Figur 31. Restfelsmodell för <strong>Oskarshamns</strong> lokala koordinatsystem med<br />

passpunkternas restfel redovisade<br />

Ett annat sätt att testa de framtagna sambanden på, som också användes i Oskarshamn,<br />

innebar att kartdetaljer i primärkartan mättes in i det nya koordinatsystemet. De koordinater<br />

som då erhölls jämfördes sedan med koordinaterna för samma detaljer i det gamla<br />

koordinatsystemet, efter att en transformation <strong>av</strong> dessa hade gjorts till den nya<br />

kartprojektionen i SWEREF 99. Istället för att hitta på nya kartdetaljer i form <strong>av</strong> exempelvis<br />

ett rutnät <strong>av</strong> punkter, gjordes här alltså en jämförelse mellan inmätta och transformerade<br />

kartdetaljer. I såväl x-led som y-led var de kvadratiska medelfelen (RMS) 16 millimeter och<br />

den största radiella <strong>av</strong>vikelsen 49 millimeter. Jämförelsen visade att de flesta differenser var i<br />

samma storleksordning som det förväntade felet i GNSS-mätningarna, vilket gjorde det svårt<br />

48


att dra några mer långtgående slutsatser utifrån dessa differenser. För att ytterligare jämföra<br />

koordinaterna från mätningarna med de transformerade koordinaterna kunde även en<br />

minstakvadratinpassning gjorts för att minimera differenserna. De nya minimerade<br />

skillnaderna hade på detta sätt kunnat påvisa eventuella systematiska brister i differenserna,<br />

genom att de i så fall hade <strong>av</strong>vikit mer eller mindre från de rena koordinatdifferenserna<br />

(Svanholm, 2000). Eftersom antalet inmätta punkter var få, kändes det dock inte vettigt att<br />

göra denna jämförelse i Oskarshamn.<br />

3.13 Transformation <strong>av</strong> <strong>Oskarshamns</strong> <strong>kommun</strong>s primärkarta<br />

Efter att transformationssambanden var testade och kontrollerade, var det dags att<br />

transformera <strong>kommun</strong>ens hela primärkarta från de gamla lokala koordinatsystemen till<br />

kartprojektionen SWEREF 99 16 30 i det nya och globalt anpassade <strong>referenssystem</strong>et<br />

SWEREF 99. I samarbete med leverantören <strong>av</strong> <strong>kommun</strong>ens kartsystem hade<br />

transformationsverktyget Tekis-Trans tagits fram, som använde Lantmäteriets program Gtrans<br />

som transformationsmotor och transformerade informationen i den Oracle-databas som<br />

geografisk data lagrades i.<br />

Förutsättningarna var att transformationsmetoden skulle klara <strong>av</strong> att transformera en databas<br />

där information från två olika koordinatsystem lagrats i samma objektsklasser. Detta var<br />

möjligt tack vare att en områdesindelning kunde skapas, som skilde de två gamla<br />

koordinatsystemen från varandra. Valet <strong>av</strong> transformationssamband vid transformationen <strong>av</strong><br />

en kartdetalj kunde alltså göras beroende på var kartdetaljen återfanns geografiskt. En<br />

kartdetalj inom regionsystemets tillämpningsområde transformerades alltså med<br />

transformationssambandet som gällde just regionsystemet, medan en kartdetalj i <strong>Oskarshamns</strong><br />

lokala koordinatsystem transformerades med det samband som gällde inom det<br />

koordinatsystemets tillämpningsområde. Att denna metod kunde användas vid<br />

transformationen berodde på att de två koordinatsystemen hade tillämpningsområden som inte<br />

överlappade varandra.<br />

Den stora delen under själva transformationsarbetet var alltså bytet <strong>av</strong> referens- och<br />

koordinatsystem för primärkartan; det vill säga informationen i kartdatabasen. Det fanns dock<br />

även geografisk information i en mängd filer <strong>av</strong> olika slag och format. Även vid<br />

transformation <strong>av</strong> dessa kunde Tekis-Trans användas, tack vare att Safe Softwares<br />

transformationsprogram FME, med alla dess möjligheter att hantera en mängd olika filformat,<br />

integrerades i transformationsprocessen och kunde anropas <strong>av</strong> Tekis-Trans.<br />

49


4. Diskussion och slutsatser<br />

Det här examensarbetet har bestått i att beskriva <strong>Oskarshamns</strong> <strong>kommun</strong>s arbete med att byta<br />

<strong>referenssystem</strong> till SWEREF 99 och därmed också gå över från två gamla koordinatsystem till<br />

en kartprojektion i det nya <strong>referenssystem</strong>et. Ur geodetisk synvinkel var det ingen ny mark<br />

som bröts; studier <strong>av</strong> olika sätt att transformera kartdetaljer och svårigheter i samband med<br />

det, har gjorts tidigare. Det är även så att Lantmäteriet, på det här området, har tagit fram<br />

omfattande och bra riktlinjer för <strong>kommun</strong>er att arbeta utifrån. Dokumentationen som gjorts<br />

här har därför mer redovisat hur en <strong>kommun</strong>, helt konkret med utgångspunkt från det material<br />

som redan finns, kan genomföra ett <strong>referenssystem</strong>sbyte. Siktet har även varit inställt på att<br />

beskriva den analysprocess <strong>av</strong> mätresultat som Lantmäteriet normalt utför för <strong>kommun</strong>ernas<br />

räkning. Detta arbete har lett till kontakt med viktiga geodetiska begrepp och företeelser, och<br />

också inneburit att många ställningstaganden med geodetisk anknytning har behövt göras.<br />

De frågeställningar och reflektioner som kommit fram under arbetets gång har redan, tidigare<br />

i rapporten, diskuterats i sitt sammanhang, men här följer en sammanfattning och<br />

kompletterad diskussion kring en del <strong>av</strong> det som berörts.<br />

4.1 Valet <strong>av</strong> mätteknik<br />

I <strong>av</strong>snittet om valet <strong>av</strong> mätteknik och mätmetod (3.3) behandlas de två teknikerna statisk<br />

mätning och RTK-mätning. Det blir klart att det finns såväl för- som nackdelar med båda<br />

teknikerna men att båda kan användas vid inmätning <strong>av</strong> passpunkter, såvida den mätmetod<br />

man väljer klarar <strong>av</strong> att ge tillförlitliga resultat. I Oskarshamn föll valet på en mätmetod som<br />

använder nätverks-RTK-teknik. Anledningen till detta var att det hela blev en resursfråga i<br />

flera <strong>av</strong>seenden. Det handlade om tillgång på instrument, beräkningsprogram och<br />

mätpersonal. Dessutom passade den flexibilitet och lättillgänglighet som nätverks-RTKtekniken<br />

erbjuder bra med tanke på kr<strong>av</strong>et att den löpande mätverksamheten i <strong>kommun</strong>en<br />

också skulle fungera och kunna drivas parallellt med projektet för <strong>referenssystem</strong>sbytet.<br />

En bättre framförhållning och <strong>plan</strong>ering för övergång till SWEREF 99 på länsnivå, hade dock<br />

inneburit ett helt annat läge med andra resurser och möjligheter. Om flera <strong>av</strong> <strong>Oskarshamns</strong><br />

<strong>kommun</strong>s grann<strong>kommun</strong>er också hade varit i färd med att byta <strong>referenssystem</strong>, är det mycket<br />

möjligt att en gemensam satsning på en mätmetod med statisk mätteknik hade varit att<br />

föredra. Statisk mätning ger trots allt en helt annan kontrollerbarhet än RTK-teknik och<br />

resulterar i ett nätutjämnat GNSS-mätt stomnät. Dessutom är det möjligt att tidsmässigt<br />

koncentrera mätinsatserna mer med statisk mätning än med RTK-mätning. Finns det<br />

önskemål om att använda statisk mätning vid inmätning <strong>av</strong> passpunkter vid exempelvis ett<br />

<strong>referenssystem</strong>sbyte i någon annan <strong>kommun</strong>, bör därför alla tänkbara intressenter på ett tidigt<br />

stadium kontaktas så att en gemensam strategi kan tas fram för hur samordning och samarbete<br />

kan etableras.<br />

4.2 Analysarbete och historisk stomnätskunskap<br />

Som redan nämnts gjorde <strong>Oskarshamns</strong> <strong>kommun</strong> själv analysarbetet <strong>av</strong> resultatet från<br />

passpunktsinmätningarna. De möjligheterna har självklart inte alla <strong>kommun</strong>er och normalt tar<br />

Lantmäteriet hand om denna del <strong>av</strong> <strong>referenssystem</strong>sbytet för <strong>kommun</strong>ernas räkning. Det kan<br />

dock ändå vara på sin plats att beröra några <strong>av</strong> de fördelar som finns med att på egen hand<br />

50


sköta hela kedjan från förberedelserna inklusive <strong>plan</strong>eringen <strong>av</strong> passpunkterna, via<br />

mätningarna med analysprocessen och fram till upprättandet <strong>av</strong> de slutgiltiga<br />

transformationssambanden. I denna utvidgade redogörelse <strong>av</strong> och diskussion kring<br />

analysarbetet kommer även nyttan <strong>av</strong> kunskapen om <strong>kommun</strong>ens stomnät och geodetiska<br />

mätningsarbete ur ett historiskt perspektiv att lyftas fram.<br />

Beskrivningen <strong>av</strong> analysprocessen visar att det, tack vare att analyserna <strong>av</strong> mätdata kunde<br />

göras fortlöpande under mätningarnas gång, blev möjligt att på ett tidigt stadium se hur<br />

kartläggningen <strong>av</strong> den bristande geometrin fortskred. Ganska snart efter att mätningar inom<br />

ett område hade påbörjats, kunde det bedömas hur omfattande de fortsatta mätinsatserna där<br />

skulle behöva göras. Att analyserna gjordes <strong>av</strong> <strong>kommun</strong>en själv möjliggjorde också att den<br />

stora kunskap om stomnäten och den erfarenhet <strong>av</strong> mätningsarbete i <strong>kommun</strong>en som fanns, på<br />

ett utmärkt sätt kunde tas tillvara. Det blev möjligt att diskutera och resonera om<br />

analysresultat direkt när de tagits fram och ofta kompletterade den historiska kunskapen<br />

analyserna och gjorde dem bättre. Den historiska kunskapen om hur etableringen <strong>av</strong><br />

stomnätens olika delar gått till var exempelvis något som många gånger kunde förklara annars<br />

svårförståeliga analysresultat. I vissa fall ledde kunskapen om stomnäten och deras svagheter<br />

till att fler mätningar gjordes än vad som från början var <strong>plan</strong>erat. Genom att ha<br />

analysprocessen i egen regi kunde upplägget för mätningarna förbättras i olika <strong>av</strong>seenden<br />

redan innan en hel mätetapp var genomförd. Vid flera tillfällen blev det även möjligt att redan<br />

under första mätetappen fånga in trender och mönster hos restfel, som det annars hade krävts<br />

en andra mätetapp och nya tidskrävande besök i fält för att kartlägga. På detta sätt kunde<br />

mätnings- och analysförfarandet göras mer flexibelt och dynamiskt än om mätdata skulle ha<br />

skickats iväg till Lantmäteriet för analys.<br />

Framställningen redogör återkommande för nyttan <strong>av</strong> att ha det historiska perspektivet och<br />

kunskapen med under genomförandet <strong>av</strong> <strong>referenssystem</strong>sbytet. Detta gällde även den<br />

historiska kunskapen om <strong>kommun</strong>ens mätningstekniska arbete, som kom väl till pass vid valet<br />

<strong>av</strong> passpunkter och generellt när det gällde arbetet i fält. Exempelvis underlättades valet <strong>av</strong><br />

bra passpunkter betydligt, tack vare att det fanns god lokalkännedom och dessutom vetskap<br />

om och mångårig erfarenhet <strong>av</strong> mätning på de utlagda stompunkterna.<br />

4.3 Hjälpmedel i analysprocessen<br />

De analysverktyg som användes i Oskarshamn var de som även Lantmäteriet använder när det<br />

analyserar restfel på passpunkter i samband med <strong>kommun</strong>ala <strong>referenssystem</strong>sbyten. Gtrans,<br />

Lantmäteriets egen produkt, användes tillsammans med ett Matlab-skript, som även det<br />

utvecklats <strong>av</strong> Lantmäteriet, för att ta fram restfelsvektorkartor och restfelsvariationsbilder.<br />

Dessa grafiska presentationer utgjorde de viktigaste underlagen i analysprocessen, och med<br />

dem som grund togs beslut om var mätinsatser skulle sättas in och koncentreras.<br />

Beskrivningen <strong>av</strong> Oskarhamns <strong>kommun</strong>s <strong>referenssystem</strong>sbyte visar att det alltså är fullt<br />

möjligt och genomförbart att i princip enbart utifrån tolkning <strong>av</strong> dessa grafiska presentationer<br />

<strong>av</strong>göra på vilka stompunkter mätningar ska göras för att få till stånd en restfelsmodell som<br />

med tillräcklig kvalitet ska kunna användas som upprätningsmodell i samband med<br />

transformationen till SWEREF 99.<br />

Den fulla nyttan <strong>av</strong> de grafiska presentationerna <strong>av</strong> restfelen erhölls dock när de även<br />

kombinerades tillsammans med <strong>kommun</strong>alt kartmaterial på olika sätt. Exempelvis kunde<br />

primärkartan eller ett ortofoto tillsammans med en överlagrad transparent georefererad<br />

restfelsvariationsbild bli ett mycket kraftfullt verktyg i arbetet med att försöka förklara<br />

51


estfelens utseende i olika delar <strong>av</strong> stomnäten. De kritiska områdena i restfelsvariationsbilden<br />

kunde i många fall direkt förstås genom en titt på ortofotot, eftersom detta <strong>av</strong>slöjade<br />

företeelser som exempelvis naturliga hinder och barriärer. Kombinationen <strong>av</strong> kartmaterial<br />

eller flygbildsmaterial och grafiskt presenterat analysmaterial g<strong>av</strong> en överblick över trender<br />

och mönster hos restfelen och erbjöd även en ökad förståelse för varför stomnäten såg ut som<br />

de gjorde och hur de under årens lopp hade förändrats och utvecklats på olika sätt i olika<br />

områden. Listan över vettiga och användbara sätt att kombinera presentationer <strong>av</strong><br />

analysresultaten med kartinformation på, kan säkerligen göras lång. Egentligen är det väl,<br />

liksom i så många andra fall, endast fantasin som sätter gränserna. En bra idé kan vara att<br />

förena kunskap i geodesi med kunskap i GIS och kartframställning för att försöka skapa<br />

meningsfulla och ändamålsenliga kartprodukter som kan användas för att ytterligare förbättra<br />

och förenkla analysarbetet.<br />

4.4 Transformationens kvalitet<br />

Det övergripande projektmålet för byte <strong>av</strong> <strong>referenssystem</strong> i <strong>Oskarshamns</strong> <strong>kommun</strong> var att ett<br />

byte skulle genomföras på ett sådant sätt att kartbasens information efter bytet hade en så hög<br />

kvalitet som möjligt. Förutsättningen för att uppnå detta var självklart att de<br />

transformationssamband som togs fram verkligen kunde användas för att dels transformera<br />

från lokala koordinatsystem till en kartprojektion i ett nytt <strong>referenssystem</strong>, dels för att räta upp<br />

en primärkarta som under drygt ett halvt sekel byggts upp genom mätning med stomnätens<br />

punkter som utgångspunkter.<br />

I <strong>av</strong>snittet om test och kontroll <strong>av</strong> transformationssambanden (<strong>av</strong>snitt 3.12) beskrivs ett par<br />

olika sätt för hur kartinformationens kvalitet kunde bedömas efter transformationen till<br />

SWEREF 99. Figurerna 27 och 30 visar skillnaden mellan att transformera med och utan<br />

restfelsmodell i respektive lokalt koordinatsystem. Detta sätt att åskådliggöra<br />

restfelsmodellens inverkan vid transformationen ger visserligen svar på hur den bristande<br />

geometrin hos kartbasens information rätas upp i olika delar <strong>av</strong> primärkartan, men några<br />

direkta mått på hur väl den lägesbundna informationen stämmer med verkligheten efter<br />

transformationen ges inte på detta sätt. Det som dock kan ses är som sagt att interpolationen<br />

med restfelsmodellen verkligen ger ett synligt samband mellan restfelsvektorerna hos<br />

passpunkterna och restfelen hos kartdetaljerna.<br />

Det andra sättet som transformationens resultat kontrollerades på beskrivs också i <strong>av</strong>snitt<br />

3.12, och innebar att ett försök gjordes för att ta fram ett slags facit genom inmätning <strong>av</strong><br />

kartdetaljer i det nya koordinatsystemet. De inmätta kartdetaljernas koordinatvärden<br />

jämfördes med de transformerade koordinaterna, vilket g<strong>av</strong> en uppfattning om<br />

transformationsresultatet. Jämförelsen visade att transformation med restfelsinterpolering g<strong>av</strong><br />

koordinatvärden som överlag stämde väl överens med de inmätta värdena. Tyvärr fanns inte<br />

resurserna för att genomföra det antalet testinmätningar som hade behövts för att jämförelsen<br />

och kontrollen skulle kunna betraktas som riktigt rigorös. Inmätningarna fick mer ses som ett<br />

litet antal stickprov i de områden där stomnätens bristande geometrier varit svårast att<br />

kartlägga. Bland de kartdetaljer som mättes in blev den största radiella <strong>av</strong>vikelsen från de<br />

transformerade värdena 49 millimeter, vilket, trots det bristfälliga statistiska materialet, får ses<br />

som en antydan om hur stora <strong>av</strong>vikelserna i primärkarteområdet kan vara. Detta resultat i<br />

<strong>Oskarshamns</strong> <strong>kommun</strong> liknar resultaten i andra <strong>kommun</strong>er där <strong>referenssystem</strong>sbytet också har<br />

genomförts.<br />

52


Utöver bristerna i kontrollen <strong>av</strong> transformationsresultatet tillkommer ytterligare ett par<br />

osäkerhetsfaktorer som bör lyftas fram i bedömningen <strong>av</strong> transformationens kvalitet. För det<br />

första var de kartdetaljer som kontrollmättes väl markerade punkter i stomnäten och därmed<br />

förmodligen betydligt stabilare än de flesta andra kartdetaljer i primärkartan. Det är<br />

exempelvis inte säkert att <strong>av</strong>vikelserna för kontrollmätta gränspunkter skulle ha blivit<br />

desamma (Svanholm, 2002). Den andra osäkerheten i transformationsresultatet ligger i det<br />

faktum att resultat från GNSS-mätningar under flera år blandats med resultat från traditionella<br />

mätningar i kartdatabasen, utan att detta markerats på något sätt. Vid transformationen kunde<br />

därför ingen separat hantering <strong>av</strong> de GNSS-mätta kartdetaljerna göras, vilket självklart<br />

medförde försämrad kvalitet. Resultaten <strong>av</strong> de utförda GNSS-mätningarna var ju inte<br />

behäftade med de geometriska brister som stomnäten hade, och behövde därför inte<br />

transformeras till det nya <strong>referenssystem</strong>et med inblandning <strong>av</strong> någon restfelsmodell. Så här i<br />

efterhand inses alltså lätt att behandlingen <strong>av</strong> mätdata från GNSS-mätningar i kartdatabasen<br />

borde ha gjorts med mer eftertanke, så att dessa mätningar hade kunnat lyftas ut och<br />

transformerats separat, utan användning <strong>av</strong> restfelsmodell.<br />

53


Litteraturförteckning<br />

Böcker, rapporter och föredrag<br />

Andersson, Bengt (2007) Projektsammanfattning <strong>av</strong> RIX 95 samt Lantmäteriets officiella<br />

transformationer. Föredraget, som hölls under konferensen MätKart 07 i Uppsala 9-11 maj<br />

2007 (arrangör: SKMF - Sveriges Kart- och Mätningstekniska Förening), återfinns i<br />

konferensens föredragsdokumentation, Session 4a, föredrag 2, sidorna 1-4.<br />

Fan, Huaan (2003) Theory of Errors and Least Squares Adjustment. Stockholm. Institutionen<br />

för geodesi och fotogrammetri vid <strong>Kungliga</strong> Tekniska Högskolan (KTH).<br />

Jansson, Patric (2007) Behöver <strong>kommun</strong>erna underhålla sina stomnät?. Artikel i SINUS, nr.<br />

2:2007.<br />

Lantmäteriet (2008), Hjälpfilen till transformationsprogrammet Gtrans ver. 3.62. (Gtrans,<br />

2008)<br />

Lantmäteriet (1993), Handbok till mätningskungörelsen Geodesi, GPS. (HMK Ge:GPS, 1993)<br />

Gävle. Lantmäteriet.<br />

Lilje, Christina, Engfeldt, Andreas och Jivall, Lotti (2007). Introduktion till GNSS. Gävle.<br />

LMV-rapport 2007:11.<br />

Norin, Dan; Engfeldt, Andreas; Johansson, Daniel; Lilje Christina (2007) Kortmanual för<br />

mätning med SWEPOS Nätverks-RTK-tjänst. Utgåva 2. Gävle. LMV-rapport 2006:2.<br />

Odolinski, Robert och Sunna, Johan (2009) Detaljmätning med nätverks-RTK – en<br />

noggrannhetsundersökning. Gävle. Examensarbete i geodesi vid KTH i Stockholm och LMVrapport<br />

2009:2.<br />

Reit, Bo-Gunnar (1997) A simple way of introducing a global reference frame in surveying<br />

and mapping. Svensk lantmäteritidskrift, 34(264) 1997 sid. 87-90.<br />

Reit, Bo-Gunnar (2003) Connecting a local system to a geocentric reference frame. Intern<br />

artikel, Lantmäteriet, 2003-04-04.<br />

Reit, Bo-Gunnar (2009) Om geodetiska transformationer. Rapport som ska publiceras.Gävle.<br />

Lantmäteriet.<br />

Svanholm, Niklas (2000) Studier <strong>av</strong> deformationer vid byte <strong>av</strong> koordinatsystem. Gävle.<br />

Examensarbete i geodesi vid KTH i Stockholm och LMV-rapport 2000:4.<br />

von Malmborg, Helena (2006). Jämförelse <strong>av</strong> Epos och nätverks-DGPS. Gävle.<br />

Examensarbete i geodesi vid KTH i Stockholm och LMV-rapport 2006:5.<br />

54


Webplatser och elektroniska dokument<br />

Lantmäteriet, geodesi (geodesiinformation från Lantmäteriets hemsida) [www]. Hämtat från<br />

http://www.lantmateriet.se/geodesi. Hämtat 16 mars 2007.<br />

Lantmäteriet, Rix 95 (information om Rix 95-projektet från Lantmäteriets hemsida) [www].<br />

Hämtat från http://www.lantmateriet.se/rix95. Hämtat 16 mars 2009.<br />

RIX 95-projektet, Lantmäteriets hemsida [www]. Hämtat från<br />

http://www.lantmateriet.se/upload/filer/kartor/geodesi_gps_och_detaljmatning/aktuella_projekt/RIX_95/rix95inf<br />

ormation.pdf. Hämtat 16 mars 2009.<br />

Samband mellan nationella <strong>referenssystem</strong> och lokala/<strong>kommun</strong>ala system, Lantmäteriets<br />

hemsida [www]. Hämtat från<br />

http://www.lantmateriet.se/upload/filer/kartor/geodesi_gps_och_detaljmatning/aktuella_projekt/RIX_95/mer_om<br />

_samband.pdf. Hämtat 16 mars 2009.<br />

Engberg, Lars E. och Lilje, Mikael (2006). Direct Projection – An Efficient Approach for<br />

Datum Transformation of Plane Co-ordinates, [www]. Hämtat från<br />

http://www.lantmateriet.se/upload/filer/kartor/geodesi_gps_och_detaljmatning/Rapporter-<br />

Publikationer/Publikationer/PS_5_6_Engberg.pdf. Publicerat 13 december 2006. Hämtat 16 mars<br />

2009.<br />

Kempe, Christina; Alfredsson, Anders; Engberg, Lars E. och Lilje, Mikael (2006). Correction<br />

Model to Rectify Distorted Co-ordinate Systems, [www]. Hämtat från<br />

http://www.lantmateriet.se/upload/filer/kartor/geodesi_gps_och_detaljmatning/Rapporter-<br />

Publikationer/Publikationer/CORRECTION_MODEL_060707_korr.pdf. Publicerat 13 december 2006.<br />

Hämtat 16 mars 2009.<br />

Kempe, Christina och Géza Lohász (2008). Att ta fram en restfelsmodell, [www]. Hämtat från<br />

http://www.lantmateriet.se/upload/filer/kartor/geodesi_gps_och_detaljmatning/Nytt_referenssytem/Dokument/A<br />

tt_ta_fram_restfelsmodell_080924.pdf. Publicerat 2 oktober 2008. Hämtat 16 mars 2009.<br />

Korrektion med restfelsinterpolering. Nummer sex i Lantmäteriets serie informationsblad om<br />

införandet <strong>av</strong> SWEREF 99 och RH 2000 [www]. Hämtat från<br />

http://www.lantmateriet.se/upload/filer/kartor/geodesi_gps_och_detaljmatning/Nytt_referenssytem/Infoblad/info<br />

_blad-6.pdf. Hämtat 16 mars 2009.<br />

Gauss Conformal Projection (Transverse Mercator) , Lantmäteriets hemsida [www]. Hämtat<br />

från<br />

http://www.lantmateriet.se/upload/filer/kartor/geodesi_gps_och_detaljmatning/geodesi/Formelsamling/Gauss_C<br />

onformal_Projection.pdf. Publicerat 2005-08-31. Hämtat 16 mars 2009.<br />

55


Reports in Geographic Information Technology 2006-2009<br />

The TRITA-GIT Series - ISSN 1653-5227<br />

2006<br />

06-001 Uliana Danila. Corrective surface for GPS-levelling in Moldova. Master of Science thesis in geodesy<br />

No. 3089. Supervisor: Lars Sjöberg. TRITA-GIT EX 06-001. January 2006.<br />

06-002 Ingemar Lewén. Use of gyrotheodolite in underground control network. Master of Science thesis in<br />

geodesy No. 3090. Supervisor: Erick Asenjo. TRITA-GIT EX 06-002. January 2006.<br />

06-003 Johan Tornberg. Felfort<strong>plan</strong>tningsanalys i GIS-projekt genom Monte Carlo-simulering. Master of<br />

Science thesis in geoinformatics. Supervisor: Mats Dunkars. TRITA-GIT EX 06-003. February 2006.<br />

06-004 Constantin-Oct<strong>av</strong>ian Andrei. 3D affine coordinate transformations. Master of Science thesis in<br />

geodesy No. 3091. Supervisor: Huaan Fan. TRITA-GIT EX 06-004. March 2006.<br />

06-005 Helena von Malmborg. Jämförelse <strong>av</strong> Epos och nätverks-DGPS. Master of Science thesis in geodesy<br />

No. 3092. Supervisor: Milan Horemuz. TRITA-GIT EX 06-005. March 2006.<br />

06-006 Lina Ståhl. Uppskattning <strong>av</strong> kloridhalt i brunnar - modellering och visualisering med hjälp <strong>av</strong> SAS-<br />

Bridge. Master of Science thesis in geoinformatics. Supervisor: Hans Hauska. TRITA-GIT EX 06-006.<br />

May 2006.<br />

06-007 Dimitrios Chrysafinos. VRS network design considerations applicable to the topology of the Hellenic<br />

Positioning System (HEPOS) stations. Master of Science thesis in geodesy No.3093. Supervisor: Lars<br />

Sjöberg. TRITA-GIT EX 06-007. May 2006.<br />

06-008 Tao Zhang. Application of GIS and CARE-W systems on water distribution networks. Master of<br />

Science thesis in geoinformatics. Supervisor: Mats Dunkars. TRITA-GIT EX 06-008. May 2006.<br />

06-009 Krishnasamy Satish Kumar. Usability engineering for Utö tourism information system. Master of<br />

Science thesis in geoinformatics. Supervisor: Mats Dunkars. TRITA-GIT EX 06-009. May 2006.<br />

06-010 Irene Rangle. High resolution satellite data for mapping Landuse/land-cover in the rural-urban fringe<br />

of the Greater Toronto area. Supervisor: Yifang Ban. TRITA-GIT EX 06-010. May 2006.<br />

06-011 Kazi Ishtiak Ahmed. ENVISAT ASAR for land-cover mapping and change detection. Supervisor:<br />

Yifang Ban. TRITA-GIT EX 06-011. May 2006.<br />

06-012 Jian Liang. Synergy of ENVISAT ASAR and MERIS data for landuse/land-cover classification.<br />

Supervisor: Yifang Ban. TRITA-GIT EX 06-012. May 2006.<br />

06-013 Assad Shah. Systematiska effecter inom Riks<strong>av</strong>vägningen. Master of Science thesis in geodesy<br />

No.3094. Supervisor: Tomas Egeltoft. TRITA-GIT EX 06-013. August 2006.<br />

06-014 Erik Trehn. GPS Precise Point Positioning – An Investigation In Reachable Accuracy. Master of<br />

Science thesis in geodesy No.3095. Supervisor: Milan Horemuz. TRITA-GIT EX 06-014. August 2006.


2007<br />

07-001 Carl Schedlich. Turn at the roundabout: A practical assessment of spatial representations in two<br />

different GPS interfaces from a pedestrian’s perspective. Bachelor of Science thesis in geoinformatics.<br />

Supervisor: Michael Le Duc. January 2007.<br />

07-002 Staffan Bengtsson. Förändringsanalys i ortofoton. Master of Science thesis in geoinformatics.<br />

Supervisor: Jonas Nelson and Patric Jansson. TRITA-GIT EX 07-002. March 2007.<br />

07-003 Joseph Addai. Qantification of temporal changes in metal loads – Moss data over 20 years. Master of<br />

Science thesis in geoinformatics. Supervisor: Katrin Grunfeld. March 2007.<br />

07-004 Stephen Rosewarne. Deformation study of the Vasa Ship. Bachelor of Science thesis in geodesy<br />

No.3097. Suppervisor: Milan Horemuz. March 2007.<br />

07-005 Naeim Dastgir. Processing SAR Data Using Range Doppler and Chirp Scaling Algorithms. Master of<br />

Science thesis in geodesy No.3096. Supervisor: Lars Sjöberg. April 2007.<br />

07-006 Torgny Israelsson and Youssef Shoumar. Motion Detection with GPS. Master of Science thesis in<br />

geodesy No.3098. Supervisor: Milan Horemuz. April 2007.<br />

07-007 Akjol Djenaliev. Multicriteria decision making and GIS for railroad <strong>plan</strong>ning in Kyrgyzstan. Master of<br />

Science thesis in geoinformatics. Supervisor: Hans Hauska. May 2007.<br />

07-008 Anna Hammar. Quality comparison of automatic 3D city house modelling methods from laser data.<br />

Master of Science thesis in geodesy No.3099. Supervisor: Milan Horemuz. May 2007.<br />

07-009 Md Ubydul Haque. Mapping malaria vector habitats in the dry season in Bangladesh using Spot<br />

imagery. Master of Science thesis in geoinformatics. Supervisor: Hans Hauska. May 2007.<br />

07-010 Jing Jiang. Analysis of the Suitable and Low-Cost Sites for Industrial Land Using Multi Criteria<br />

Evaluation: A Case of Panzhihua, China. Master of Science thesis in geoinformatics. Supervisor:<br />

Yifang Ban. June 2007.<br />

07-011 Ragh<strong>av</strong>endra Jayamangal. Quantification of coastal erosion along Spey Bay and the Spey River using<br />

photogrammetry and LiDAR imagery-derived DTMs. Master of Science thesis in geoinformatics.<br />

Supervisor: Yifang Ban and Jim Hansom. June 2007.<br />

07-012 Alicia E. Porcar Lahoz. An analysis of how geographical factors affect real estate prices. Master of<br />

Science thesis in geoinformatics. Supervisor: Yifang Ban. October 2007.<br />

07-013 Ebenezer Kwakye Bentum. Detection of land use change in the Accra Metropolitan Area from 1990<br />

to 2000. Master of Science thesis in geoinformatics. Supervisor: Hans Hauska. November 2007.<br />

07-014 Jesper Gradin. Aided inertial n<strong>av</strong>igation field tests using an RLG IMU. Master of Science thesis in<br />

geodesy No.3100. Supervisor: Milan Horemuz. December 2007.


2008<br />

08-001 Wan Wen. Road Roughness Detection by analyzing IMU Data. Master of Science thesis in<br />

geodesy No.3101. Supervisor: Milan Horemuz. Janaury 2008.<br />

08-002 Ilias Daras. Determination of a gr<strong>av</strong>imetric geoid model of Greece using the method of KTH. Master<br />

of Science thesis in geodesy No.3102. Supervisor: Huaan Fan and Kalliopi Papazissi. Janaury 2008.<br />

08-003 Anwar Negash Surur. Surveying, modelling and visualisation of geological structures in the<br />

Tunberget tunnel. Master of Science thesis in geoinformatics. Supervisor: Hans Hauska. February 2008.<br />

08-004 Helena Swanh. Noggrannhetskontroll <strong>av</strong> data 3D-modell från pulsskanner och fasmätningsskanner.<br />

Master of Science thesis in geodesy No.3103. Supervisor: Milan Horemuz. March 2008.<br />

08-005 Henrik Löfqvist. Inpassning <strong>av</strong> mätdata i inhomogena nät. Master of Science thesis in geodesy<br />

No.3104. Supervisor: Milan Horemuz and Anders Viden. April 2008.<br />

08-006 Cement Takyi. Lineament study in the Björkö area using satellite radar and airborne VLF data -<br />

possible tools for groundwater prospecting. Master of Science thesis in geoinformatics. Supervisor:<br />

Herbert Henkel. June 2008.<br />

08-007 Qarin Bånkestad. 3D data: Faktisk hantering för att kunna lagra, återfå och visualisera tredimensionell<br />

information. Master of Science thesis in geoinformatics. Supervisor: Hans Hauska. June 2008.<br />

08-008 Tadesse Tafesse Kebede. Development and implementation of filter algorithms and controllers to a<br />

construction machine simulator. Master of Science thesis in geodesy No.3105. Supervisor: Lars<br />

Sjöberg. September 2008.<br />

08-009 Gulilat Tegane Alemu. Assessments of effects on mixing different types of GPS antennas and<br />

receivers. Master of Science thesis in geodesy No.3106. Supervisor: Lars Sjöberg. September 2008.<br />

08-010 Jure Kop. Tests of New Solutions to the Direct and Indirect Geodetic Problems on the Ellipsoid.<br />

Master of Science thesis in geodesy No.3107. Supervisor: Lars Sjöberg. September 2008.<br />

08-011 Karoliina Kolehmainen. Monitoring and Analysis of Urban Land Cover Changes over Stockholm<br />

Region between 1986 and 2004 using Remote Sensing and Spatial Metrics. Master of Science thesis in<br />

geoinformatics. Supervisor: Yifang Ban. October 2008.<br />

08-012 Johan Wallin. Land-Cover Mapping in Stockholm Using Fusion of ALOS PALSAR and SPOT Data.<br />

Master of Science thesis in geoinformatics. Supervisor: Yifang Ban. October 2008.<br />

08-013 Mårten Korall. A comparison between GNSS, Laser scanning and Photogrammetry for volume<br />

calculation of aggregates and gr<strong>av</strong>el materials. Master of Science thesis in geodesy No.3108.<br />

Supervisor: Milan Horemuz. December 2008.


2009<br />

09-001 Ahmed Abdallah. Determination of a gr<strong>av</strong>imetric geoid if Sudan using the KTH method. Master of<br />

Science thesis in geodesy No.3109. Supervisor: Huaan Fan. Janaury 2009.<br />

09-002 Hussein Mohammed Ahmed Elhadi. GIS, a tool for p<strong>av</strong>ement management. Master of Science thesis<br />

in geoinformatics. Supervisor: Hans Hauska. February 2009.<br />

09-003 Robert Odolinski and Johan Sunna. Detaljmätning med nätverks-RTK – en<br />

noggrannhetsundersökning (Detail surveying with network RTK – an accuracy research). Master of<br />

Science thesis in geodesy No.3110. Supervisor: Clas-Göran Persson and Milan Horemuz. March 2009.<br />

09-004 Jenny Illerstam och Susanna Bosrup. Restfelshantering med Natural Neighbour och TRIAD vid byte<br />

<strong>av</strong> koordinatsystem i <strong>plan</strong> och höjd. Master of science thesis in geodesy No. 3111. Supervisor: Milan<br />

Horemuz and Lars Engberg. March 2009.<br />

09-005 Erik Olsson. Exporting 3D Geoinformation from Baggis Database to CityGML. Supervisors: Peter<br />

Axelsson and Yifang Ban. April 2009.<br />

09-006 Henrik Nilsson. Referenssystemsbyte i <strong>Oskarshamns</strong> <strong>kommun</strong> – en fallstudie (Change of reference<br />

systems in Oskarshamn – a case study). Master’s of Science thesis in geodesy No.3112. Supervisor:<br />

Huaan Fan. May 2009.


TRITA-GIT EX 09-06<br />

ISSN 1653-5227<br />

ISRN KTH/GIT/EX--09/006-SE

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!