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 ...
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