23.07.2013 Views

Ved Askeengen Ved Askeengen og Harrestrup Ådal - VBN - Aalborg ...

Ved Askeengen Ved Askeengen og Harrestrup Ådal - VBN - Aalborg ...

Ved Askeengen Ved Askeengen og Harrestrup Ådal - VBN - Aalborg ...

SHOW MORE
SHOW LESS

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

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

Kortlægning <strong>og</strong><br />

Landinspektørstudiet 5.semester<br />

<strong>Aalborg</strong> Universitet – København<br />

afsætning<br />

~ <strong>Ved</strong> <strong>Askeengen</strong> <strong>og</strong> <strong>Harrestrup</strong> <strong>Ådal</strong> ~<br />

1<br />

Udarbejdet af:<br />

Andreas K. Jensen &<br />

Johan V. Eckhoff


TITELBLAD<br />

Titel: Kortlægning <strong>og</strong> afsætning – <strong>Ved</strong> <strong>Askeengen</strong> <strong>og</strong> <strong>Harrestrup</strong> <strong>Ådal</strong><br />

Tema: Kortlægning <strong>og</strong> afsætning<br />

Projektperiode:<br />

02-09-2011 til 06-12-2011<br />

Projektgruppe: L5-ball2011_01 Synopsis:<br />

Deltagere:<br />

Andreas Kudsk Jensen<br />

Johan Vibe Eckhoff<br />

Vejledere:<br />

Karsten Jensen<br />

Jens Juhl<br />

Peter Cederholm<br />

Oplagstal: 3<br />

Sideantal: 100 heraf 32 siders bilag<br />

Bilagsantal <strong>og</strong> -art: 44 heraf 32 på CD<br />

Afsluttet den 05-12-2011<br />

Rapportens indhold er frit tilgængeligt, men offentliggørelse (med kildeangivelse) må kun ske efter aftale med<br />

forfatterne.<br />

3<br />

Dette projekt omhandler temaet kortlægning <strong>og</strong> afsætning. Det<br />

består af fire faser, som omhandler forskellige kortlægnings- <strong>og</strong><br />

afsætningsmetoder. Endvidere vurderes kvaliteten af forskellige<br />

kortprodukter.<br />

Fase 1 resulterer i udarbejdelsen af et teknisk kort med<br />

højdeinformation samt en digital terrænmodel. Kortet er blevet til<br />

ved RTK-måling. Derudover er der indmålt 4 GI-planfikspunkter for<br />

at vurdere spændinger i det landsdækkende net. Der er endvidere<br />

indmålt 15 kontrolpunkter til brug ved kontrol af det tekniske kort<br />

<strong>og</strong> terrænmodellen.<br />

I 2. fase foretages afsætning af skel med tilhørende vej samt en<br />

bygning. Det afsatte kontrolleres, <strong>og</strong> nøjagtigheden af afsætningen<br />

vurderes. Endvidere indeholder fasen to udjævninger efter mindste<br />

kvadraters princip af et geometrisk nivellement samt et net, hvor<br />

koordinaterne til hjælpepunkter bestemmes. Fase 2 er baseret på<br />

både RTK-måling samt måling med totalstation.<br />

Fase 3 består af flere kortlægninger ved fot<strong>og</strong>rammetrisk opmåling<br />

herunder et teknisk kort med højdeinformation, en digital<br />

terrænmodel <strong>og</strong> to ortofotos. Der er foretaget en relativ <strong>og</strong> absolut<br />

orientering ud fra to billeder udleveret ved semesterstart. Det<br />

tekniske kort <strong>og</strong> terrænmodellen er kontrolleret med de 15<br />

kontrolpunkter fra fase 1.<br />

Fase 4 indeholder 11 kortsammenligninger med både kort<br />

produceret i fase 1 <strong>og</strong> 3 samt andre landsdækkende kort, ortofotos<br />

<strong>og</strong> højdemodeller. Kortene vurderes i forhold til en beregning af<br />

forventede spredninger.<br />

Det kan konkluderes, at n<strong>og</strong>le kortlægninger er foretaget med større<br />

nøjagtighed end andre, hvilket hovedsageligt skyldes<br />

opmålingsmetoden. Endvidere er nøjagtigheden <strong>og</strong>så et udtryk for<br />

erfaringen med kortlægningsarbejdet, herunder kan nævnes<br />

opmåling fot<strong>og</strong>rammetrisk. Der blev under projektarbejdet opdaget<br />

en nævneværdig forskel fra landsdækkende ortofotos til egne<br />

producerede ortofotos.


FORORD<br />

Projektet er udarbejdet på landinspektøruddannelsens 5. semester ved <strong>Aalborg</strong> Universitet København. Projektets<br />

overordnede tema er ”Kortlægning <strong>og</strong> afsætning” med fokus på n<strong>og</strong>le af de tekniske opgaver, som en privat<br />

praktiserende landsinspektør varetager. Projektet er udarbejdet i perioden fra 02-09-2011 til 06-12-2011.<br />

Formålet med semesteret er, at de studerende skal tilegne sig viden om teori <strong>og</strong> anvendelse af forskellige opmålings-<br />

<strong>og</strong> afsætningsmetoder herunder RTK-måling, afsætning med totalstation <strong>og</strong> fot<strong>og</strong>rammetrisk kortlægning. Derudover<br />

skal den studerende tilegne sig viden om lineær algebra samt udjævning efter mindste kvadraters princip. Den<br />

studerende skal endvidere tilegne sig kompetencer ved kombination af forskellige opmålingsmetoder samt vurdering<br />

af kvaliteten af de sammensatte opmålingsmetoder.<br />

Projektet er opdelt i fire faser, som hver især beskæftiger med forskellige opmålings- <strong>og</strong> kortlægningsmetoder.<br />

Projektet beskriver de valg <strong>og</strong> fravalg, som projektgruppen har foretaget undervejs i udarbejdelsen af de fire faser.<br />

Til udarbejdelse af projektet skal der rettes en særlig tak til de grundejere, som har stillet deres grund til rådighed ved<br />

opmåling med RTK under Fase 1 samt nivellement under Fase 2.<br />

Projektet er understøttet af tre kurser: Databearbejdning <strong>og</strong> vurdering (DBV), Dataindsamling <strong>og</strong> modellering (DIM) <strong>og</strong><br />

Lineær Algebra (LiA).<br />

Til kortlægninger <strong>og</strong> beregninger er anvendt følgende pr<strong>og</strong>rammer: GeoCAD, TMK, ImageStation, Matlab <strong>og</strong> diverse<br />

scripts udviklet af gruppen eller vejledere. I visse kort er COWI’s ortofoto DDOland 2010 anvendt, som<br />

baggrundsbillede.<br />

Projektrapporten indeholder 14.433 ord.<br />

4


INDHOLDSFORTEGNELSE<br />

TITELBLAD .................................................................................................................................................................. 3<br />

FORORD ...................................................................................................................................................................... 4<br />

INDLEDNING ............................................................................................................................................................... 7<br />

VURDERING AF RTK-MÅLING .....................................................................................................................................10<br />

FASE 1 – KORTLÆGNING VED RTK-MÅLING ...............................................................................................................12<br />

GI-PUNKTER OG INDBYRDES NØJAGTIGHED.............................................................................................................................. 12<br />

KRAVSPECIFIKATION FOR DETAILMÅLING ................................................................................................................................ 14<br />

NØJAGTIGHED AF MÅLINGERNE ............................................................................................................................................ 15<br />

FREMSTILLING AF TEKNISK KORT ........................................................................................................................................... 16<br />

FREMSTILLING AF DIGITAL TERRÆNMODEL .............................................................................................................................. 17<br />

VURDERING AF NØJAGTIGHED I TEKNISK KORT ......................................................................................................................... 21<br />

OPSAMLING ...................................................................................................................................................................... 22<br />

FASE 2 – AFSÆTNING .................................................................................................................................................23<br />

KRAVSPECIFIKATION ........................................................................................................................................................... 23<br />

AFSÆTNING AF SKEL OG VEJE ............................................................................................................................................... 23<br />

BYGNINGSAFSÆTNING ........................................................................................................................................................ 26<br />

OPSAMLING ...................................................................................................................................................................... 34<br />

FASE 3 – KORTLÆGNING VED FOTOGRAMMETRI .......................................................................................................35<br />

KRAVSPECIFIKATION ........................................................................................................................................................... 35<br />

BILLEDEKVALITET ............................................................................................................................................................... 36<br />

RELATIV ORIENTERING ........................................................................................................................................................ 37<br />

ABSOLUT ORIENTERING ....................................................................................................................................................... 39<br />

KONTROL AF ORIENTERING .................................................................................................................................................. 42<br />

FOTOGRAMMETRISK KORT ................................................................................................................................................... 44<br />

DTM-MODEL ................................................................................................................................................................... 46<br />

ORTOFOTO ....................................................................................................................................................................... 49<br />

OPSAMLING ...................................................................................................................................................................... 51<br />

FASE 4 – VURDERING AF FORSKELLIGE KORTLÆGNINGSMETODER ............................................................................52<br />

RTK-TK* KONTRA FOTO-TK* .............................................................................................................................................. 54<br />

RTK-TK* KONTRA ORTO-EGET* ........................................................................................................................................... 55<br />

RTK-TK* KONTRA ORTO-BLOM .......................................................................................................................................... 56<br />

RTK-TK* KONTRA COWI-DDO .......................................................................................................................................... 57<br />

COWI-DDO KONTRA ORTO-BLOM* ................................................................................................................................... 58<br />

COWI-DDH KONTRA RTK-DTM* ....................................................................................................................................... 59<br />

FOTO-DTM* KONTRA COWI DDH ...................................................................................................................................... 60<br />

FOTO-TK* KONTRA FOT ..................................................................................................................................................... 61<br />

ORTO-EGEN* KONTRA FOT ................................................................................................................................................. 62<br />

ORTO-BLOM* KONTRA FOT .............................................................................................................................................. 63<br />

RTK-TK* KONTRA FOT ...................................................................................................................................................... 64<br />

OPSAMLING ...................................................................................................................................................................... 65<br />

5


KONKLUSION .............................................................................................................................................................66<br />

KILDELISTE .................................................................................................................................................................68<br />

BILAGSLISTE ...............................................................................................................................................................69<br />

OVERSIGT OVER BILAG PÅ CD ............................................................................................................................................... 69<br />

BILAG – VURDERING AF RTK-MÅLING .................................................................................................................................... 71<br />

BILAG-VERIFICERING AF INSTRUMENTER ................................................................................................................................ 74<br />

BILAG – VURDERING AF TRANSFORMATION AF GI-PUNKTER ....................................................................................................... 77<br />

BILAG – VURDERING AF KONTROLMÅLING AF BYGNINGSDIMENSIONER ......................................................................................... 85<br />

BILAG – AFSÆTNING AF SKEL OG VEJE .................................................................................................................................... 87<br />

BILAG – DOKUMENTATIONSFIL MINIMUM FASTHOLDT UDJÆVNING ............................................................................................. 89<br />

BILAG – DOKUMENTATIONSFIL FASTHOLDT UDJÆVNING ............................................................................................................. 91<br />

BILAG – BINGOFIL .............................................................................................................................................................. 93<br />

BILAG - BEREGNING AF SOLVINKEL ........................................................................................................................................ 96<br />

BILAG – DOKUMENTATION FOR PASPUNKTER .......................................................................................................................... 97<br />

6


INDLEDNING<br />

Dette projekt består af fire faser, hvor der fokuseres på forskellige kortlægnings- <strong>og</strong> opmålingsmetoder samt en<br />

vurdering af kvaliteten af disse.<br />

Fase 1 omhandler undersøgelse af kortlægning ved RTK-måling, hvor nøjagtigheden af RTK-målinger undersøges ved<br />

indmåling af et mindre projektområde. Endvidere indmåles 4 GI-planfikspunkter for at undersøge spændinger i det<br />

landsdækkende fikspunktsnet. Fase 1 resulterer i to produkter: et teknisk kort i 2½D <strong>og</strong> en digital terrænmodel.<br />

Fase 2 omhandler afsætning, hvor en mindre udstykning med vej afsættes med GPS, en bygning afsættes med<br />

totalstation <strong>og</strong> sidst vurderes afsætningerne. Dertil hører en udjævning af et geometrisk nivellement foretaget fra fire<br />

kendte højdefikspunkter <strong>og</strong> udjævning af et net bestående af hjælpepunkter efter mindste kvadraters princip.<br />

Udjævningen af netmålingen skal ske ved brug af Anblok <strong>og</strong> udjævning af det geometriske nivellement sker ved brug<br />

af eget script udviklet til formålet. Fase 2 resulterer i: to udjævninger, to afsætninger <strong>og</strong> en kontrol af det afsatte.<br />

Fase 3 beskæftiger sig med kortlægning ved fot<strong>og</strong>rammetri, hvor forskellige kortlægningsmetoder ved fot<strong>og</strong>rammetri<br />

undersøges <strong>og</strong> vurderes. Fase 3 benytter to billeder henholdsvis 84B <strong>og</strong> 85B, som er blevet udleveret ved<br />

semesterstart. Disse billeder bliver orienteret relativt <strong>og</strong> absolut i forhold til UTM32/EUREF89 <strong>og</strong> anvendes derefter til<br />

at kortlægge samme område, som blev indmålt i fase 1. Fase 3 resulterer i: en modtagerkontrol af billederne, en<br />

relativ <strong>og</strong> absolut orientering, et teknisk kort i 2½D, en digital terrænmodel samt fremstilling af to ortofotos, hvor det<br />

ene udarbejdes med en digital terrænmodel produceret i fasen, <strong>og</strong> det andet tager udgangspunkt i BLOM’s<br />

højdemodel.<br />

Fase 4 undersøger forskellige kortlægningers nøjagtighed. Fase 4 anvender kortlægninger fra projektets andre faser<br />

samt n<strong>og</strong>le landsdækkende kortlægninger, ortofotos <strong>og</strong> højdemodeller til sammenligning af de forskellige<br />

kortlægnings- <strong>og</strong> opmålingsmetoder. Fase 4 resulterer i 11 kortsammenligninger, hvor forskellen vurderes i forhold til<br />

n<strong>og</strong>le beregnede forventede spredninger.<br />

Billede 1: Oversigt over projektområdet<br />

Der er anvendt en punktnummerstrategi for indmålinger jf. studievejledningen:<br />

”Fikspunkter (alle typer): Originale (lange) numre<br />

• Frie opstillinger: 1 - 99<br />

7


• Hjælpepunkter til afsætning: 300 - 399<br />

• Hjælpepunkter til terrestrisk måling: 400 - 499<br />

• Paspunkter (plan): 500 - 599<br />

• Paspunkter (højde): 600 - 699<br />

• Paspunkter (kombineret): 700 - 799<br />

• Detailpunkter terrestrisk måling: 1000 - 4999<br />

• Detailpunkter fladenivellement: 5000 – 8999”<br />

Til udarbejdelse af projektet er der anvendt en punktspredning, som er defineret efter (Jensen, 2005):<br />

= + <br />

2<br />

Nedenstående diagram viser opbygningen af de forskellige faser samt resultaterne fra disse.<br />

8<br />

(Juhl, 2011a, s. 6)


Figur 1: Projektopbygning<br />

9


Inden fase 1 skal det endvidere vurderes, hvor godt der måles med RTK i de to servicesystemer Leica Smartnet <strong>og</strong><br />

Gpsnet.dk. Denne vurdering danner grundlag for vurdering af spredninger i de efterfølgende faser, hvor kvaliteten af<br />

forskellige kortlægninger skal vurderes.<br />

VURDERING AF RTK-MÅLING<br />

Vurderingen af RTK-målingen er foretaget på baggrund af 20 definerede punkter, som er målt i alt 4 gange; 2 gange i<br />

GPSnet.dk <strong>og</strong> 2 gange i Leica Smartnet (Bilag – vurdering af RTK-måling).<br />

I Leica Smartnet er der foretaget flere beregninger, hvorefter punkter, som afviger mere end 3 ∙ √2 ∙ (,,), er<br />

sorteret fra. Derefter er beregningerne fortaget igen, indtil alle målingerne ligger indenfor fejlgrænsen. I alt ender det<br />

med, at der er 17 kontrolpunkter, som ikke er behæftet med grove fejl.<br />

For at lave en generel vurdering af, hvor godt der måles i de to referencesystemer samlet set, anvendes definitionen<br />

for punktspredning (Jensen, 2005), som beskrevet i indledningen. I Leica Smartnet ligger spredningen på _ =<br />

0.009 <strong>og</strong> = 0.012 .<br />

I GPSnet.dk er der ligeledes foretaget flere beregninger, hvor punkter, som afviger mere end 3 ∙ √2 ∙ (,,), er<br />

sorteret fra. I GPSnet.dk ligger spredningen på _ = 0.012 <strong>og</strong> = 0.013 . I alt er der 19 kontrolpunkter,<br />

som ligger indenfor fejlgrænsen.<br />

De målinger, som bliver sorteret fra, da de ligger uden for fejlgrænsen, er afhængige af flere faktorer. Dels var der<br />

forholdsvis meget blæst på dagen, hvorfor centreringen over punktet i visse tilfælde ikke har været optimal. Endvidere<br />

foregik kontrolopmålingen tæt på en bygning, hvilket kan have ført til multipath, hvor signalet muligvis kan have været<br />

reflekteret på bygninger eller andet i nærheden af opmålingen. Nedenstående billede viser et eksempel på multipath:<br />

Billede 2: Illustration af multipath<br />

Kilde: Kowoma.de, 2011<br />

I udregningen viser tendensen, at Leica Smartnet er mere præcis end GPSnet.dk, d<strong>og</strong> er der flere kontrolpunkter, som<br />

ligger til grund for spredningen i GPSnet.dk. Det vurderes, at der måles med en punktspredning på 11 mm i planen i<br />

10


egge net ud fra beregningen (,) = .. ≈ 0.011 . Derudover vurderes det, at der måles<br />

med en punktspredning i højden på 13 mm i begge net ud fra beregningen () = .. ≈<br />

0.013.<br />

<br />

11


FASE 1 – KORTLÆGNING VED RTK-MÅLING<br />

Formålet med fase 1 er at fremstille et teknisk kort <strong>og</strong> en digital terrænmodel, som kan anvendes til projektering af<br />

fremtidige tekniske anlæg i projektområdet. Opmålingen til det tekniske kort skal ske ved RTK-måling. Derudover skal<br />

der foretages en vurdering af det tekniske kort i forhold til de dobbeltmålte nærtliggende GI-planfikspunkter.<br />

Fase 1 består dermed af en vurdering af GI-planfikspunkternes nøjagtighed, en kravspecifikation til fremstilling af det<br />

tekniske kort, en beskrivelse af fremstillingen af det tekniske kort samt den digitale terrænmodel. Endvidere<br />

indeholder fase 1 en vurdering af nøjagtigheden af det tekniske kort <strong>og</strong> terrænmodellen i forhold til de nøjagtigheder,<br />

som blev opstillet under kravspecifikationen.<br />

GI-PUNKTER OG INDBYRDES NØJAGTIGHED<br />

Formålet med at indmåle GI-planfikspunkterne er at vurdere om, der er n<strong>og</strong>le spændinger i nettet. GI-fikspunkterne er<br />

blevet målt ind med enten GPS eller en kombination af GPS <strong>og</strong> totalstation i UTM32/EUREF89.<br />

Kort 1: Oversigt over planfikspunkter i nærheden af <strong>Harrestrup</strong>.<br />

Det viste sig, at der var få punkter, som omkranser området ved <strong>Harrestrup</strong>. Punkt 2-06-00005 blev fravalgt, da det var<br />

uegnet for GPS-måling, da det var omkranset af flere træer. Derudover var punkt 2-01-00002 målt ind tidligere,<br />

hvorfor det blev vurderet, at punkt 2-06-00005 ikke var nødvendigt. Punkt 1-11-00804 <strong>og</strong> 1-11-00807 viste sig at være<br />

tabtgået efter en søgning efter begge punkter. Til søgning af punkterne blev der anvendt afsætningsfunktionen i<br />

GPS’en.<br />

Punkt 1-04-00830 blev vurderet til at være utilgængeligt, da det lå på taget af Leica’s bygning. Der blev konfereret om<br />

punktet med en person fra Leica, som mente, at det kun er en referencestation. Punkt 1-04-00006 blev vurderet til at<br />

være udefinerbart, da det er på gavlen af Herlev Kirke, <strong>og</strong> det ikke var klart, hvor på gavlen punktet er. Desuden blev<br />

punkterne 2-01-00007 <strong>og</strong> 2-01-00880 fravalgt, da et andet punkt blev indmålt, som lå tættere på området ved<br />

<strong>Harrestrup</strong>.<br />

12


Følgende punkter blev indmålt: 2-01-00002, 2-01-00883, 2-04-00870 <strong>og</strong> 2-06-00009.<br />

2-01-<br />

00002<br />

2-01-<br />

00883<br />

2-06-<br />

00009<br />

2-04-<br />

00870<br />

Valdemar E Valdemar N Middel E Middel N Afvigelse<br />

E<br />

13<br />

Afvigelse<br />

N<br />

709153.647 6179837.402 709153.678 6179837.405 0.031 0.003<br />

711989.401 6180644.459 711989.400 6180644.456 -0.001 -0.003<br />

712879.280 6178459.564 712879.271 6178459.573 -0.011 0.011<br />

714902.506 6178706.280 714902.407 6178706.183 -0.099 -0.091<br />

Alle punkter er dobbeltmålte for at hindre grove fejl. N<strong>og</strong>le punkter er endda målt 3 gange, da der var<br />

uoverensstemmelse mellem de to første målinger. Der viste sig ikke n<strong>og</strong>et tydeligt mønster i afvigelserne fra<br />

Valdemar. D<strong>og</strong> var afvigelserne på E-koordinaten i n<strong>og</strong>le tilfælde større end på N-koordinaten.<br />

Punkt 2-04-00870 afviger meget. Punktet er defineret som midten af toppen på en skorsten, <strong>og</strong> det blev målt ind med<br />

totalstation ved hjælp af en fremskæring. Da det er svært at definere midten af skorstenen blev der målt en<br />

horisontalretningsvinkel på hver side af skorstenen, hvorefter vinklen blev midlet. Der blev i alt foretaget tre<br />

observationer fra tre etablerede hjælpepunkter, <strong>og</strong> skorstenen er målt med 1 sats. Koordinater er beregnet parvis i<br />

TMK, da TMK kun kan beregne fremskæring mellem to målinger. Derefter er de tre koordinatsæt til 2-04-00870 blevet<br />

midlet. Der er overensstemmelse mellem de tre beregninger, <strong>og</strong> der er målt indenfor 3 mgon mellem første <strong>og</strong> anden<br />

hovedstilling.<br />

Kort 2: Skitse over måling af GI-punkt 2-04-00870<br />

Punkt 2-01-00883 er vedligeholdt af LE34, da punktet sidder på deres bygning. Dette er sandsynligvis medvirkende til,<br />

at afvigelse for dette punkt bliver mindre end de resterende punkter.<br />

VURDERING AF TRANSFORMATION<br />

Det er forsøgt at transformere de indmålte GI-punkter ind over GI-punkterne fra Valdemar ved hjælp af 2Dtranslation,<br />

2D-transformation (Helmert) samt 2D-translation med skalering (Bilag – Vurdering af transformation af GIpunkter).


2D-translationen viser et stort residual for punkt 2-04-00870, hvilket stemmer overens med den tidligere vurdering af<br />

punktet. Punktet er enten dårligt indmålt eller <strong>og</strong>så stemmer koordinaterne fra Valdemar ikke overens med punktets<br />

nuværende placering. De tre beregninger af punktet stemmer som nævnt tidligere n<strong>og</strong>enlunde overens med højeste<br />

residual fra middelkoordinatet på 24 mm i E <strong>og</strong> 15 mm i N. Det tyder derfor på, at punktets koordinater ikke stemmer<br />

overens med Valdemar-koordinaterne.<br />

Derfor er punktet taget ud af beregningen, da forskellen er meget større end i de andre GI-fikspunkter. Der er herefter<br />

foretaget en 2D-translation, en 2D-translation med skalering samt en 2D-transformation (Helmert) af de resterende<br />

tre punkter over Valdermars koordinater. Det blev her konstateret, at residualerne efter en 2D-translation er relativt<br />

små (længste residualvektor på 25 mm). Samtidig var spredningen på vægtenheden 16 mm. Der skete en lille rykning,<br />

hvor tx=-0.006m <strong>og</strong> ty=-0.004m.<br />

Der blev på trods af den lille spredning på vægtenheden alligevel forsøgt med en 2D-transformation(Helmert) samt en<br />

2D-translation med skalering. Disse to gav n<strong>og</strong>enlunde samme resultat. Eneste afvigelse var længste residualvektor,<br />

som ændrede sig 1 mm. Dette er et udtryk for, at der ikke er n<strong>og</strong>en nævneværdig drejning. Skaleringen blev 1.000010,<br />

hvilket svarer til en ændring på 10 ppm. Der var derfor en forskel fra 2D-translation uden skalering på både<br />

residualvektorerne, spredning på vægtenheden samt konfidensellipserne, hvor residualerne nu er nede på et par<br />

millimeter.<br />

Dette betyder, at der skal forekomme en lille skalering for at få de indmålte punkter til at stemme bedre overens med<br />

Valdemar-koordinaterne <strong>og</strong> en lille ændring <strong>og</strong> . En drejning er d<strong>og</strong> unødvendig. D<strong>og</strong> er ændringerne meget<br />

små, <strong>og</strong> da RTK-målingen blev vurderet til ca. 1 cm i Leica Smartnet, kan det være en del af forklaringen på<br />

uoverensstemmelser mellem koordinaterne. Endvidere skal det tages i betragtning, at der <strong>og</strong>så er en spredning på<br />

KMS’s opmåling af punkterne. Med hensyn til skaleringen kan der gisnes om, hvorvidt der ved tidligere udjævning af<br />

GI-punkterne ikke er blevet reduceret til ellipsoiden, hvilket kunne føre til fejl i den størrelsesorden.<br />

KRAVSPECIFIKATION FOR DETAILMÅLING<br />

Detailopmålingen skal kunne benyttes til fremtidigt teknisk anlægsarbejde indenfor projektområdet.<br />

Detailopmålingen tager udgangspunkt i FOT-standarden (FOTdanmark, 2010) <strong>og</strong> skal tilpasses de aktuelle forhold i<br />

projektområdet <strong>og</strong> opmålingsmetoden. Nedenstående tabel indeholder de elementer, som er givet ved FOTstandarden.<br />

Indmålingen af de forskellige elementer er tilpasset de aktuelle forhold.<br />

I FOT-standarden nævnes geometrisk, l<strong>og</strong>isk <strong>og</strong> tematisk nøjagtighed samt fuldstændighed. Den geometriske<br />

nøjagtighed vil blive gennemgået senere. Krav til registrering <strong>og</strong> den l<strong>og</strong>iske nøjagtighed fremgår af skemaet nedenfor.<br />

Med hensyn til fuldstændighed forventes det som udgangspunkt, at der ikke bliver lavet n<strong>og</strong>en fejl grundet områdets<br />

lille udstrækning.<br />

14


Objektgruppe Objekttype Geometri Beskrivelse Bemærkning<br />

Bygninger Bygning Flade Fri mur over sokkel<br />

(Målt ved<br />

jordoverflade)<br />

15<br />

Bygninger har en<br />

mindste størrelse på<br />

10 m 2 . Garager<br />

registreres <strong>og</strong>så, hvis<br />

de er<br />

sammenhængende<br />

med bygningen.<br />

Carport Flade <strong>Ved</strong> tagudhæng Carporte registreres,<br />

hvis de er<br />

sammenhængende<br />

med bygningen.<br />

Trafik Vejkant Skillelinje<br />

Fortov Skillelinje<br />

Teknik Mast Punkt Højde mindst 3 m<br />

Nedløbsrist Punkt Rensebrønde <strong>og</strong><br />

nedløbsriste på<br />

parceller indmåles<br />

ikke.<br />

Brønddæksel Punkt<br />

Installationsskab Punkt Mindst 40 cm<br />

Diverse Punkt Eks. Brandhaner<br />

Natur Levende hegn Linje Mindst 2 meter langt<br />

Træ Punkt Træ defineres, som<br />

én stamme, der går i<br />

jorden.<br />

Trægruppe Polygon med<br />

punkter<br />

Hvis der er mere end<br />

3 træer, som ikke er<br />

på linje, registreres<br />

dette som<br />

trægruppe<br />

Plankeværk Linje Mindst 2 meter langt<br />

Hegn/Stendige Linje Mindst 2 meter langt Diger, som er under<br />

2 meter bredt <strong>og</strong> 50<br />

meter langt,<br />

registreres som hegn<br />

Det tekniske kort skal registreres i 2½D. Dette betyder, at alle målinger skal have en kote. Koten er registreret ved<br />

jordoverfladen, hvilket gælder for alle oplyste koter i kortet.<br />

NØJAGTIGHED AF MÅLINGERNE<br />

Det blev vurderet til at starte med, at der kunne måles med en præcision i både Leica Smartnet <strong>og</strong> GPSnet.dk på<br />

= 0.011 <strong>og</strong> = 0.013 . FOT-standarden er baseret på fot<strong>og</strong>rammetrisk måling. Heri defineres en<br />

nøjagtighed i forhold en GSD på 10 cm på = 0.10 <strong>og</strong> = 0.15 . Dermed vil punktspredningen for målingerne<br />

være = .<br />

√ =0.07m i forhold til den definerede punktspredning. Nøjagtigheden i plan <strong>og</strong> højde vurderes d<strong>og</strong> ud fra<br />

nøjagtigheden af RTK-måling, da der måles med GPS. Det vil sige, at hvert veldefineret registreret punkt skal leve op til<br />

disse nøjagtigheder.


Til terrænmodellen anvendes halvdelen af den oplyste nøjagtighed fra FOT-standarden til at finde den rette pilhøjde.<br />

Dette gøres, da detailpunkter er indmålt med GPS, hvilket indebærer, at nøjagtigheden er bedre end 15 cm. Samtidig<br />

ønskes det at reducere datamængden, hvorfor nøjagtigheden er højere end den mulige ved RTK-måling.<br />

NØJAGTIGHED AF MÅLEMETODER<br />

Da der ikke kun er anvendt enkeltpunkts RTK-måling, er det <strong>og</strong>så nødvendigt at opstille en nøjagtighed for de andre<br />

målemetoder, som er taget i brug under indmålingen af parcellerne. Der er brugt flugtlinjer, bueskæring samt<br />

ort<strong>og</strong>nalmåling.<br />

Flugtlinjerne er foretaget ud fra en antagelse om lige vægge. A priori spredningen på en flugtlinje vurderes ud fra et<br />

punkts geometri langs en linje. Der er derfor taget udgangspunkt i et standardafstandsforhold mellem de indmålte<br />

punkter <strong>og</strong> detailpunktet P, = <br />

= −0.9, som er beregnet ud fra en stikprøve på et af de indmålte huse i<br />

detailområdet.<br />

<br />

Ud fra dette beregnes a priori spredningen, som = √2 ∙ + 2 ∙ + 1 ∙ . Derfor bliver a priori spredningen for<br />

en flugtlinje = √2 ∙ −0.9 + 2 ∙ −0.9 + 1 ∙ 0.011 = 0.010. D<strong>og</strong> er der ikke taget højde for, hvor nøjagtigt<br />

flugtningen kan defineres. Derfor bliver der tillagt et bidrag på 2 cm. Derved bliver den nye a priori spredning<br />

= 2 ∙ + 2 ∙ + 1 ∙ <br />

+ = 0.022. Da de fleste hushjørner antages at være baseret på skæringen<br />

mellem to flugtlinjer, skal a priori spredningen for skæringen <strong>og</strong>så beregnes. Det antages, at alle flugtlinjer er<br />

foretaget med lige stor nøjagtighed, så = . Dermed bliver spredningen på to flugtlinjer = <br />

<br />

=<br />

<br />

∗.<br />

<br />

= 0.022 (Cederholm & Jensen, 2006).<br />

Bueskæring er anvendt til punkter, hvor det ikke er muligt at lave flugtlinjer. Bueskæringer foretages med<br />

tilnærmelsesvis rette vinkler, <strong>og</strong> spredningen for bueskæring beregnes ved hjælp af = ∗( <br />

å<br />

16<br />

<br />

<br />

ø<br />

Variansen for de forskellige parametre optræder 2 gange, da bueskæringer konstrueres ud fra en måling til et<br />

hushjørne fra to punkter målt med RTK. Det defineres, at der måles med en nøjagtighed på 1 cm med stålmålebånd,<br />

da midtpunktet på GPS-stokken skal defineres. Endvidere kan et hushjørne defineres indenfor 0.5 cm, da det er<br />

forholdsvis nye huse. Derfor giver spredningen på bueskæringen = 0.016 (Cederholm & Jensen, 2006).<br />

Sidste målemetode, som kan anvendes, er ort<strong>og</strong>onalmåling. Her anvendes samme udgangspunkt, som ved flugtlinjer,<br />

hvilket gav en spredning på 0.022 m. D<strong>og</strong> skal der tillægges en nøjagtighed på måling med stålmålebånd, som sættes<br />

til 1 cm. Dette betyder, at = 2 ∙ + 2 ∙ + 1 ∙ <br />

+ + å = 0.024 (Cederholm & Jensen,<br />

2006).<br />

FREMSTILLING AF TEKNISK KORT<br />

Det tekniske kort er fremstillet ved RTK-måling. <strong>Ved</strong> RTK-målingen er flere målemetoder anvendt herunder<br />

bueskæring, flugtlinjer <strong>og</strong> ort<strong>og</strong>onalmåling. <strong>Ved</strong> ort<strong>og</strong>onalmåling er der fortaget opmåling med målebånd langs<br />

flugtlinje til at definere eksempelvis et hjørne på en carport.<br />

<br />

)<br />

.


De indmålte punkter konverteres i TMK, <strong>og</strong> punkter med en høj PDOP-værdi bliver sorteret fra. Det blev vurderet, at<br />

træer <strong>og</strong> levende hegn ikke skal have samme nøjagtighed, som eksempelvis hushjørner, <strong>og</strong> derfor godt kan have en<br />

højere DOP-værdi. Derfor blev n<strong>og</strong>le punkter bibeholdt på trods af en lidt for høj PDOP-værdi. Der var n<strong>og</strong>le målinger<br />

af hushjørner (4 punkter), som havde en for høj PDOP-værdi. Disse målinger blev foretaget igen for sikre, at husene fik<br />

en højere nøjagtighed.<br />

Kortet rummer visse undtagelser fra FOT-standarden, da indmåling er tilpasset området <strong>og</strong> RTK-måling som tidligere<br />

nævnt. Bygninger er målt som fri mur over sokkel i stedet for ved tagudhæng. Derudover er det fravalgt at registrere<br />

garager <strong>og</strong> carporte, som ikke er sammenhængende med en permanent bygning. Fortov er registeret ved skillelinjer<br />

med anden brug. Det vil sige, at der både forekommer en skillelinje ved vej/fortov <strong>og</strong> ved fortov/parcel. Desuden er<br />

størrelsen på indmålte installationsskabe tilpasset de størrelser, som befinder sig i området. Der er derfor valgt en<br />

mindre størrelse på skabe end FOT-standardens. Desuden registreres brandhaner, hvilket ikke fremgår af FOTstandarden.<br />

Trægrupper registreres på parceller, såfremt at der er mere end tre træer, som ikke står på linje. I FOTstandarden<br />

registreres trægrupper ikke på parceller (FOTDanmark, 2010, 3.1.1-3.5.12.2).<br />

Kort 3: Teknisk kort (2½D)<br />

FREMSTILLING AF DIGITAL TERRÆNMODEL<br />

Højderne til generering af en digital terrænmodel er givet af detailmålingen, hvor de sorterede indmålte 265 punkter<br />

anvendes. Fra indmålingen er der endvidere lavet brudlinjer, som bruges til beregningen, ligesom der bruges<br />

udeladelsespolygon (Kort 4).<br />

17


Kort 4: Oversigt over indmålte punkter, omfangs- <strong>og</strong> udeladelsespolygon samt brudlinjer<br />

D<strong>og</strong> blev der sorteret et punkt fra, som viste sig at have en fejlagtig højde. Da højden i FOT har en nøjagtighed på<br />

15cm, vurderes det, at højdemodellen skal have den halve nøjagtighed, da den er målt med RTK. Dette kunne gøres<br />

mere nøjagtigt grundet GPS’ens nøjagtighed. Men for at gøre datamængden mindre vælges en højere spredning.<br />

Derfor skal spredningen være 7.5 cm i forhold til de generede trekanter. Der er blevet forsøgt med forskelle pilhøjder,<br />

men sidst blev en pilhøjde på 15 cm valgt, da spredning (7.6 cm) tilnærmelsesvist er lig med den valgte spredning.<br />

Nedenstående boks viser trekantsrapporten for genereringen af TIN-modellen.<br />

Trekantsrapport:<br />

145 trekanter er genereret. Apriori beregnet antal: 145<br />

Areal af alle trekanter 2D ( ): 5300.496639<br />

3D ( ): 5310.004924<br />

Volumen under alle trekanter ( ): 129962.485513 (til kote: 0.000)<br />

Min / max kote til trekanter: 23.784 / 25.253<br />

Spredning på kotepunkter. der ikke indgår i trekantsnettet: 0.076<br />

Antal kotepunkter til denne beregning: 171<br />

93 punkter ud af 264 benyttes til trekanter ( 35%)<br />

Under genereringen udtyndes antallet af trekanter således, at der kun er det nødvendige antal. Der er generet 145<br />

trekanter. Der er benyttet 93 punkter ud af de 264 oprindelige punkter til at generere trekanterne.<br />

18


Kort 5: Genereret trekantsmodel<br />

Dernæst skal der generes n<strong>og</strong>le højdekurver ud fra trekantsmodellen. Der er valgt at danne fuldoptrukne linjer for<br />

hele <strong>og</strong> halve meter <strong>og</strong> stiplede linjer for de resterende højder delelige med 0.25 m. Derefter er der påsat n<strong>og</strong>le<br />

kotetekster samt reduceret i datamængden således, at der kun er 371 punkter i forhold til de oprindelige 1166<br />

punkter, der blev genereret.<br />

Kort 6: Højdekurver ud fra trekantsmodellen<br />

19


VUDERING AF NØJAGTIGHED AF DIGITAL TERRÆNMODEL<br />

Terrænmodellen kontrolleres ved hjælp af de veldefinerede kontrolpunkter, som <strong>og</strong>så bliver anvendt til at kontrollere<br />

nøjagtigheden af det tekniske kort. Da det ene kontrolpunkt ligger udenfor de genererede trekanter, er der kun 14<br />

kontrolpunkter. Dernæst beregnes spredningen for de interpolerede koter fra kontrolpunkterne i trekantsnettet ud<br />

fra differensen mellem terrænkoten <strong>og</strong> de interpolerede koter fra kontrolpunkterne.<br />

Dette giver følgende spredninger:<br />

Beregning af spredninger mm. på grundlag af MULTI-FIND-punkter (kort tekst3):<br />

Antal værdier i beregningen: . . . 14<br />

Maximum værdi (abs) før udvægtning: 0.225<br />

Antal værdier udvægtet: . . . . . 0<br />

Minimum værdi: . . . . . . . . . . -0.028<br />

Maximum værdi: . . . . . . . . . . 0.113<br />

Middel af værdier: . . . . . . . . 0.038<br />

Translation (middel tal) af værdier: 0.000<br />

Spredning på værdier: . . . . . . 0.057<br />

Beregning af spredninger mm. på grundlag af MULTI-FIND-punkter (kort tekst3)<br />

reduceret med middel tal:<br />

Antal værdier i beregningen: . . . 14<br />

Maximum værdi (abs) før udvægtning: 0.225<br />

Antal værdier udvægtet: . . . . . 0<br />

Minimum værdi: ( -0.028) . . . -0.066<br />

Maximum værdi: ( 0.113) . . . 0.075<br />

Translation (middel tal) af værdier: 0.038<br />

Spredning på værdier: . . . . . . 0.043<br />

Den maksimale værdi er valgt ud fra en betragtning af, at spredning på trekantsmodellen var 7.6 cm. Derfor er det<br />

valgt, at alle værdier, som overstiger = 3 ∙ = 3 ∙ 7.6 = 0.225 er behæftet med grove fejl, <strong>og</strong> derfor skal<br />

sorteres fra. Ingen af de interpolerede koter ligger over denne værdi. Det ses, at nøjagtigheden på de interpolerede<br />

koter er 5.7 cm før transformation, mens nøjagtigheden efter transformation er 4.3 cm. Dermed overholdes<br />

nøjagtigheden på 7.5 cm.<br />

20


Kort 7: Kort over trekantsmodellen <strong>og</strong> de 15 kontrolpunkter<br />

VURDERING AF NØJAGTIGHED I TEKNISK KORT<br />

Vurdering af nøjagtigheden af det tekniske kort sker dels ud fra en kontrolopmåling med GPS af n<strong>og</strong>le veldefinerede<br />

punkter dels ud fra kontrol af bygningers dimensioner (Bilag – Vurdering af kontrolmåling af bygningsdimensioner).<br />

Der er foretaget målinger til 15 riste <strong>og</strong> brønddæksler, som omkranser området. Disse er både centralsymmetriske<br />

samt veldefinerede i både terræn <strong>og</strong> fra ortofotos. De er desuden veldefinerede i højden <strong>og</strong> i planen. De 15 punkter er<br />

indmålt fire gange – to gange med GPSnet <strong>og</strong> to gange med Leica Smartnet. De er målt ind ved at placere GPS-stokken<br />

i midten af dækslet eller risten. PDOP ligger i alle tilfælde under 5. De fire målinger er midlet, <strong>og</strong> koordinaterne er<br />

sammenlignet med koordinaterne fra detailmålingen. Ud fra differenserne er der beregnet en spredning. Dette giver<br />

en punktspredning på = <br />

<br />

.<br />

=<br />

√ =0.009m <strong>og</strong> en spredning på højden på <br />

= = 0.022. Dette<br />

betyder, at der er målt mere nøjagtigt end først antaget i planen. Formodningen var, at kortets nøjagtighed var<br />

= 0.011 . D<strong>og</strong> er højden ikke så nøjagtig som først antaget, da den forventede spredning på højden var<br />

= 0.013 , hvilket blev konstateret under afsnittet Vurdering af RTK-måling.<br />

Det tyder på, at der er sket en systematisk fejl i målingen af højden, da alle afvigelserne ligger mellem 14 <strong>og</strong> 32 mm i<br />

den samme retning. Dette kan tyde på, at der er sket en fejl i beregning af korrektion mellem referencestationen <strong>og</strong><br />

roveren, da fejlen er entydig.<br />

Endvidere er målingerne mere nøjagtige end FOT-standardens nøjagtighed både i planen <strong>og</strong> i højden.<br />

Punktspredningen i FOT-standarden lå på = 0.07m.<br />

I en anden kontrol vurderes nøjagtigheden for de objekter, som er fremkommet ved flugtmåling. Her er der foretaget<br />

en kontrol af samtlige bygningsdimensioner. Kontrollen indebærer, at der er foretaget målinger med stålmålebånd i<br />

marken, hvorefter disse er sammenlignet med mål af samme dimensioner foretaget i GeoCAD (Bilag – Vurdering af<br />

21


kontrolmåling af bygningsdimensioner). Afvigelserne er ikke entydige i samme retning, hvilket kan betyde, at<br />

stålmålebåndet har været strakt for meget eller for lidt.<br />

Der er anvendt forskellige metoder til at konstruere husene i det tekniske kort. D<strong>og</strong> er den hyppigste metode<br />

flugtmåling. Som nævnt tidligere er spredningen for et skæringspunkt mellem to flugtlinjer = 0.022 .<br />

Afvigelsen mellem målingerne i marken <strong>og</strong> GeoCAD betegnes = − − + − . <strong>Ved</strong> hjælp af<br />

fejlforplantningsloven udledes udtrykket for spredningen på afvigelsen .<br />

hvor . Dette betyder. at den forventede spredning for afvigelsen mellem målingen i marken <strong>og</strong><br />

målingen i GeoCAD er = (√2 ∙ 0.022) + 0.01 = 0.033. D<strong>og</strong> er den beregnede spredning på forskellen<br />

mellem målingen i GeoCAD <strong>og</strong> målingen i marken 0.049 m.<br />

Dette vil sige, at opmålingen af husene ikke er foretaget med den forventede nøjagtighed, hvilket kan skyldes, at<br />

vurderingen af nøjagtighed af, hvor godt flugtmålingerne er foretaget, er lavere, end opmålingen reelt er foretaget<br />

med. Samtidig kan unøjagtigheden <strong>og</strong>så skyldes dårlig geometri, da GPS’ens placering kan have gjort, at en del af<br />

satellitterne på himlen ikke kunne ses. Nøjagtigheden kunne være blevet bedre ved at foretage indmålingen af diverse<br />

hushjørner med totalstation.<br />

OPSAMLING<br />

Første fase indeholdt en fremstilling af et teknisk kort <strong>og</strong> en digital terrænmodel, som skulle kunne bruges i<br />

forbindelse med projektering eller lignende. Opmålingen skete ved hjælp af RTK. Nøjagtigheden af det tekniske kort<br />

blev vurderet til 0.9cm for veldefinerede punkter, hvilket var lidt bedre end forventet. I sammenligning af koter viste<br />

der sig en ensidig fejl, der bevirkede, at spredningen var højere end forventet.<br />

Kontrolopmåling af bygningsdimensioner viste, at der ved flugtmåling ikke blev målt med den forventede nøjagtighed.<br />

Nøjagtigheden af terrænmodellen blev vurderet til enten 5.7 cm eller 4.3 cm henholdsvis før <strong>og</strong> efter transformation.<br />

Disse spredninger skal anvendes videre i fase 4 under kortsammenligningen. I fase 4 vil det tekniske kort blive<br />

vurderet i forhold til n<strong>og</strong>le andre kortprodukter.<br />

22


FASE 2 – AFSÆTNING<br />

Formålet med fase 2 er at foretage en afsætning af en mindre udstykning med tilhørende veje samt en<br />

bygningsafsætning, hvor der afsættes et bestemt antal modullinjer. Til bygningsafsætningen hører udjævning af et net<br />

af hjælpefikspunkter til afsætning i <strong>Harrestrup</strong> <strong>Ådal</strong> samt udjævning af et geometrisk nivellement.<br />

Fase 2 består af en kravspecifikation, en afsætning af skel <strong>og</strong> veje samt en bygningsafsætning. Til de forskellige<br />

afsætninger hører en beskrivelse af den geometriske konstruktion, afsætningen i marken samt en vurdering af de<br />

forskellige afsætningers nøjagtigheder i forhold til den opstillede kravspecifikation.<br />

KRAVSPECIFIKATION<br />

AFSÆTNING AF SKEL OG VEJE<br />

Afsætning skal foretages i GPSnet.dk. Tidligere blev spredningen på en RTK-måling i planet i GPSnet.dk beregnet til at<br />

være = 0.012 .<br />

Der skal foretages afsætning af 50 punkter, som skal repræsentere både skelpunkter samt stationeringslinjer til veje.<br />

Forskellen mellem kontrolopmålingen af de afsatte punkter <strong>og</strong> designkoordinaterne bør ikke overskride =<br />

= ±3 ∙ + = ±3 ∙ √2 ∙ = ±0.051 , <strong>og</strong> spredningen på vægtenheden i EN bør ligge på omkring<br />

√2 ∙ = 0.017 (Jensen, 2005).<br />

Afsætningen af skel <strong>og</strong> veje foretages ved RTK-måling, hvor punkterne afsættes i planen.<br />

BYGNINGSAFSÆTNING<br />

Der skal foretages en udjævning af et net, som består af fire hjælpepunkter henholdsvis 100, 190, 300 <strong>og</strong> 390. Disse er<br />

i første omgang målt ind med RTK, <strong>og</strong> udjævningen sker på baggrund af en supplerende terrestrisk måling.<br />

Højderne til de fire hjælpepunkter bestemmes ved hjælp af geometrisk nivellement, hvor der nivelleres ud fra 4<br />

højdefikspunkter. Dernæst udjævnes disse ved hjælp af mindste kvadraters metode i et script, skrevet til formålet.<br />

Bygningshjørnerne afsættes ved hjælp af modullinjer en meter fra hushjørnerne. Afsætningen markeres med søm i<br />

træpæle. Det forventes, at bygningsafsætning har en nøjagtighed på = 0.002 .<br />

Bygningsafsætningen foretages ved måling med totalstation.<br />

AFSÆTNING AF SKEL OG VEJE<br />

GEOMETRISK KONSTRUKTION<br />

Afsætningen tager udgangspunkt i en geometrisk konstruktion, hvor koordinaterne til de punkter, der ønskes afsat, er<br />

konverteret fra lokalt system til UTM32/EUREF89. Afsætningsdata er konstrueret i GeoCAD, hvor koordinaterne er<br />

importeret, konverteret til UTM <strong>og</strong> roteret på plads i forhold til, hvor det er muligt at afsætte i <strong>Harrestrup</strong> <strong>Ådal</strong>. Det er<br />

23


desuden vigtigt, at afsætningsdata er afstandskorrigeret i forhold til UTM32/EUREF89. Følgende formel kan bruges til<br />

at beregne afstandskorrektionen (Jensen 2005):<br />

=<br />

<br />

(1 − 10 )1 + 10 <br />

PPMsys blev ud fra GeoCAD vurderet til 153 ppm. Da Ballerup ligger uden for cylinderens zonegrænser skal<br />

afstandskorrektionen være positiv. Samtidig skal der korrigeres for ellipsoidehøjden, hvilket gør afstandene mindre.<br />

PPMn blev beregnet ud fra = ()<br />

∙ 10<br />

<br />

= ()<br />

∙ 10 = 10. Dermed blev den samlede<br />

afstandskorrektion:<br />

= 143<br />

Nedenstående kort viser modellen skaleret <strong>og</strong> konverteret til UTM32/EUREF89.<br />

Kort 8: Kort over designkoordinater<br />

AFSÆTNING<br />

I marken blev alle punkter markeret ved hjælp af træpæle, hvor midten af træpælen angav det afsatte punkt. Derefter<br />

blev alle 50 punkter kontrolleret ved at foretage en kontrolmåling, som blev foretaget en time efter afsætningen fandt<br />

sted. Der er afsat flere skel, d<strong>og</strong> er der <strong>og</strong>så afsat centerlinjer til n<strong>og</strong>le af de tilhørende veje.<br />

Nedenstående kort viser de 50 afsatte punkter i <strong>Harrestrup</strong> <strong>Ådal</strong>.<br />

24


Kort 9: Kort over afsatte punkter<br />

KONTROL OG VURDERING<br />

Kontrollen indebærer, at de kontrolmålte punkter sammenlignes med designkoordinaterne (Bilag – Afsætning af skel<br />

<strong>og</strong> veje). Ingen af de afsatte punkter overstiger fejlgrænsen på ±0.051 i E eller N, hvilket <strong>og</strong>så ses af nedenstående<br />

diagram.<br />

Afvigelse i mm<br />

35<br />

30<br />

25<br />

20<br />

15<br />

10<br />

5<br />

0<br />

±d<br />

Af nedenstående figur ses fejlvektorer for de afsatte punkter. Spredningen på vægtenheden er 0.013 <strong>og</strong> ligger derfor<br />

lidt under den forventede spredning.<br />

25<br />

Δe i mm<br />

Δn i mm


BYGNINGSAFSÆTNING<br />

Figur 2: Beregning af koordinatdifferencer<br />

<strong>Ved</strong> bygningsafsætningen er der taget udgangspunkt i n<strong>og</strong>le netpunkter, som er etableret til formålet. Til<br />

netpunkterne foretages der et geometrisk nivellement samt en netmåling for at få n<strong>og</strong>le nøjagtige koordinater.<br />

Netpunkterne er forud for netmålingen blevet indmålt med GPS, hvilket ligger til grund for den endelige beregning af<br />

de plane koordinater.<br />

GEOMETRISK NIVELLEMENT<br />

Der er foretaget dobbeltnivellement fra fire kendte højdefikspunkter henholdsvis 2-01-09011, 2-01-09130, 2-01-09136<br />

<strong>og</strong> 2-01-09216 til fire hjælpepunkter henholdsvis 300, 190, 390 <strong>og</strong> 100.<br />

26


Der blev noteret følgende højdeforskelle:<br />

Kort 10: Skitse over nivellement<br />

Fra Til Højdeforskel Afstand<br />

(km)<br />

9011 300 -0.247 0.955<br />

300 9011 0.245 0.969<br />

9130 300 -2.909 0.761<br />

300 9130 2.898 0.761<br />

9136 9216 -4.474 0.506<br />

9216 9136 4.469 0.499<br />

9216 190 1.475 0.666<br />

190 9216 -1.479 0.564<br />

190 390 2.916 0.127<br />

390 190 -2.915 0.109<br />

390 100 -2.159 0.102<br />

100 390 2.16 0.101<br />

100 300 0.974 0.079<br />

300 100 -0.975 0.083<br />

27


Der er derefter foretaget to udjævninger af nivellementer i to scripts henholdsvis<br />

minimum_fastholdt_udjævning_nivellement.m <strong>og</strong> fastholdt_udjævning_nivellement.m (Bilag 16, Oversigt over bilag<br />

på CD), som er udviklet til formålet. I det første script beregnes en minimum fastholdt udjævning, hvor kun et<br />

fikspunkt bliver fastholdt med en høj vægt. Her undersøges det, hvor godt dobbeltnivellementerne passer overens, <strong>og</strong><br />

det kan dermed undersøges, hvor godt der nivelleres. Efterfølgende er de andre GI-højdefikspunkter taget med i<br />

beregningen for at undersøge, hvor godt de passer overens<br />

Udjævning er foretaget ved hjælp af mindste kvadraters metode, hvor løsningen er bestemt ved = ( ∙ ∙ ) ∙<br />

∙ ∙ , hvor A er designmatricen, C er vægtmatricen, <strong>og</strong> b indeholder observationerne fra nivellementer <strong>og</strong> koter til<br />

højdefikspunkterne. Vægtene i C-matricen til nivellementerne er beregnet på baggrund af spredningen på et<br />

enkeltnivellement ∆ = √ ∙ , <strong>og</strong> vægtene til højdefikspunkter er beregnet på baggrund af en estimeret<br />

nøjagtighed.<br />

<strong>Ved</strong> minimum fastholdt beregning blev nivellementet i sig selv undersøgt for grove fejl. I dette tilfælde blev 2-01-<br />

09011 fastholdt. <strong>Ved</strong> en a priori kilometerspredning på 2.8 mm gav det en spredning på vægtenheden på 1.6. D<strong>og</strong><br />

viste de normaliserede residualer, at der ikke var n<strong>og</strong>en grove fejl. A priori spredningen blev justeret så spredningen<br />

på vægtenheden gav ca. 1, hvilket betød, at a priori spredning blev justeret til ca. 4.4 mm. Spredningerne, som<br />

fremgår af kovariansmatricen, er på op til 4.8mm (Bilag – Dokumentationsfil Minimum fastholdt udjævning).<br />

Derefter blev der foretaget en udjævning, hvor alle højdefikspunkter blev inddraget med vægt, <strong>og</strong><br />

kilometerspredningen fra den minimum fastholdte udjævning blev benyttet. Først blev højdefikspunkternes<br />

nøjagtighed sat til 2 mm. Her fremgik det, at 2-01-09130 lå udenfor 3 gange spredningen, <strong>og</strong> fejlen forplantede sig<br />

endvidere til et af hjælpepunkterne. Spredningen på vægtenheden gav 1.5. Derfor blev nøjagtigheden for det<br />

pågældende højdefikspunkt sat til 1 m, hvilket betyder, at punktet stort set udelades af beregningerne. Dette giver en<br />

spredning på vægtenheden i nærheden af 1, hvilket bliver accepteret. Endvidere er de normaliserede residualer alle<br />

under 3, <strong>og</strong> kovariansmatricen har spredninger på lige under 3 mm for alle punkterne (bilag – Dokumentationsfil<br />

fastholdt udjævning).<br />

Følgende højder fremkommer af udjævningen:<br />

Punkt Kote (DVR90)<br />

2-01-09011 25.674<br />

2-01-09136 26.696<br />

2-01-09216 22.225<br />

2-01-09130 28.334<br />

300 25.430<br />

190 23.701<br />

390 26.616<br />

100 24.456<br />

Disse bliver anvendt til koordinatfilen, som skal bruges til bygningsafsætningen.<br />

NETUDJÆVNING<br />

Der er foretaget en netudjævning for at få indbyrdes nøjagtige koordinater til hjælpepunkter i forbindelse med<br />

bygningsafsætning. Fra seks frie opstillinger er der målt til de fire hjælpepunkter 100, 190, 300 <strong>og</strong> 390. Til disse<br />

punkter er der opgivet koordinater i UTM32/EUREF89, som er målt med GPS.<br />

28


Kort 11: Skitse over netmåling<br />

Netudjævningen er foretaget i TMK ved hjælp af anblok, hvor der er foretaget transformation <strong>og</strong> udjævning af<br />

observationerne. I beregningen foretages to udjævninger; én med <strong>og</strong> én uden GPS-koordinaterne. I resultaterne af<br />

den ene udjævning kan målingernes indbyrdes overensstemmelse undersøges, da GPS-målingerne udelades. Denne<br />

beregning er dermed et udtryk for den relative nøjagtighed, <strong>og</strong> dermed hvor godt der er blevet målt. I denne<br />

udjævning indgår en af de frie opstillinger som referencekoordinater. I den anden indgår GPS-observationerne som<br />

referencekoordinater, <strong>og</strong> dermed kan det udledes, hvor godt målingerne kommer til at passe overens med GPSobservationerne<br />

<strong>og</strong> den absolutte orientering.<br />

Udjævning relativt<br />

Netudjævning blev indledningsvist foretaget med udgangspunkt i a priori spredninger på 0.001 m for<br />

modelkoordinaterne <strong>og</strong> 0.01 m for referencekoordinater (GPS-koordinater) i planet. I udjævningen, hvor den relative<br />

nøjagtighed undersøges, viser der sig en spredning på vægtenheden på 1.6, selvom de normaliserede residualer ikke<br />

viser n<strong>og</strong>en grove fejl. Derfor blev a priori spredningen på modelkoordinaterne sat til 0.0015 m i planet. Dette gav en<br />

spredning på vægtenheden på 1.07, hvilket blev vurderet til acceptabelt.<br />

Spredning i E <strong>og</strong> N på de beregnede koordinater ligger i dette tilfælde på ca. 1 mm i alle tilfælde. For<br />

højdeudjævningen viser der sig en spredning på vægtenheden på 0.35. A priori spredningen skal derfor korrigeres til<br />

0.5mm før spredningen på vægtenheden giver 1. Spredningen på højderne er på 0.3 mm i alle tilfælde.<br />

29


Punkt rE rE’ rN rN’ s_E s_N rH rH’<br />

m<br />

m<br />

300 0.000 0.0 0.0009 0.8 0.0011 0.0010 -0.0005 -1.2<br />

100 -0.0001 -0.1 -0.0011 -1.0 0.0010 0.0012 0.0002 0.4<br />

190 -0.0001 -0.1 0.0010 0.8 0.0013 0.0010 0.0001 0.2<br />

390 0.0002 0.1 -0.0008 -0.8 0.0010 0.0013 0.0002 0.6<br />

Udjævning absolut<br />

I udjævningen, hvor GPS-koordinaterne indgår, giver beregningen spredning på vægtenheden på 1.02 i planen med de<br />

korrigerede a priori spredninger. Ingen af de normaliserede residualer ligger over 3. Spredningen i E <strong>og</strong> N ligger nu på<br />

mellem 5 <strong>og</strong> 8 mm, som er et udtryk for den absolutte nøjagtighed. Dette skyldes den højere unøjagtighed i opmåling<br />

med GPS. Udjævningen er derfor forløbet stort set efter planen. I højden er spredningen på vægtenheden omkring 1,<br />

<strong>og</strong> spredningen på punkterne er på små 6 mm.<br />

Punkt rE rE’ rN rN’ s_E s_N rH rH’<br />

m<br />

m<br />

300 -0.0045 -0.6 -0.0081 -1.0 0.0061 0.0059 0.0039 0.4<br />

100 0.0068 0.8 0.0017 0.2 0.0051 0.0068 -0.0001 0.0<br />

190 0.0000 0.0 0.0095 1.1 0.0074 0.0055 -0.0055 -0.6<br />

390 -0.0023 -0.3 -0.0031 -0.4 0.0054 0.0074 0.0017 0.2<br />

Netmålingen var i planen ikke så nøjagtig som først antaget. Grunden til at a priori spredningen skulle sættes højere<br />

kan skyldes, at et af punkterne lå i højt græs. Samtidig blev der brugt en prismehøjde på 40 cm for at kunne se<br />

prismet, hvilket kan have forrykket midtpunktet. D<strong>og</strong> viste det sig, at højderne havde langt højere relativ nøjagtighed<br />

end først antaget. I det geometriske nivellement blev der nivelleret over alle hjælpepunkter. Disse koter bliver brugt i<br />

forbindelse med bygningsafsætning, men det trigonometriske nivellement har d<strong>og</strong> vist sig at være så nøjagtigt, at det<br />

havde været rigeligt kun at foretage geometrisk nivellement til et at hjælpepunkterne, <strong>og</strong> derefter bruge højder fra<br />

det trigonometriske nivellement.<br />

Koordinaterne fra udjævningen bliver dernæst anvendt til at konstruere bygningen.<br />

GEOMETRISK KONSTRUKTION<br />

Det blev vurderet, at den geometriske konstruktion kunne foretages på to måder, da udjævningen resulterede i to sæt<br />

koordinater henholdvis lokalt <strong>og</strong> i UTM32/EUREF89.<br />

De lokale koordinater blev derefter anvendt, da disse ikke skulle afstandskorrigeres for ppmn <strong>og</strong> ppmsys ligesom vist<br />

tidligere under afsætning af skel <strong>og</strong> vej. Dermed svarer afstande mellem de lokale koordinater til afstande i<br />

virkeligheden.<br />

De lokale koordinater blev indlæst i GeoCAD, <strong>og</strong> fik derefter højderne fra det geometriske nivellement. Derefter blev<br />

et lokalt koordinatsystem defineret. Bygning blev placeret efter systemets northing- <strong>og</strong> easting akse, således at<br />

bygningen ligger parallelt med koordinatsystemet. Dette letter markarbejdet, da modullinjer derefter kan afsættes<br />

med en lille afvigelse i den ønskede retning. Bygningens dimensioner er 20x80 m.<br />

30


Kort 12: Kort over designkoordinater til bygningsafsætning<br />

Efter den geometriske konstruktion var fuldendt blev punkterne eksporteret til en koordinatfil, som blev indlæst i<br />

totalstationen.<br />

AFSÆTNING<br />

Før afsætningen blev der foretaget en testberegning for at se, hvor mange opstillinger, der skal til for at afsætte de<br />

otte punkter. Testberegningen viste, at det ved opstilling midt i området er muligt at holde spredningen i E <strong>og</strong> N under<br />

de 2 mm, som er angivet i kravspecifikationen. Dette er under forudsætning af, at centrering af prismet anses som<br />

værende ubetydelig, da der måles til miniprisme.<br />

31


Figur 3: Testberegning af opstilling før afsætning<br />

Derefter blev der afsat otte punkter i marken, som markerede modullinjerne for hvert bygningshjørne. Punkter blev<br />

sat 1 m i northing- <strong>og</strong> easting-retningen fra bygningshjørnet, således at der kunne spændes snor mellem punkter, som<br />

vist på nedenstående afsætningsplan. Der blev foretaget to opstillinger; én hvor punkterne blev afsat <strong>og</strong> én, hvor<br />

punkterne blev kontrolleret. Endvidere er der lavet en testberegning, hvor hjælpepunkternes spredninger er sat til 0<br />

for at beregne punkternes relative nøjagtighed i forhold til opstillingen. I testberegningen bliver midlet 1.7 mm på<br />

detailpunkterne relativt i forhold til opstillingen, hvilket betyder, at = 0.0017.<br />

Der blev foretaget følgende korrektioner for atmosfæriske forhold i marken <strong>og</strong> følgende konfigurationsparametre blev<br />

anvendt under afsætningen:<br />

Parameter Korrektion<br />

HzAccOri 0.030 g<br />

PosAccTarget 0.010 m<br />

HtAccTarget 0.002 m<br />

Temperatur 8 grader<br />

Luftfugtighed 80 %<br />

Atmosfærisk ppm -7.2<br />

Der blev derudover taget højde for refraktionskoefficienten på 0.13.<br />

32


KONTROL OG VURDERING<br />

Kort 13: Afsætningsplan af bygningsafsætning<br />

Opstillingen vurderes på baggrund af de konfigurationsparametre, som blev indstillet inden afsætningen.<br />

Opstillingerne må afvige følgende:<br />

• Spredningen på orienteringen< to gange HzAccOri<br />

• Spredningen på E <strong>og</strong> N< to gange PosAccTarget<br />

• Spredningen på højden< to gange HtAccTarget<br />

Spredningerne på begge opstillinger er begge under ovenstående fejlgrænser for grove fejl, hvilket betyder, at<br />

opstillingerne er godkendt (Bilag 13, Oversigt over bilag på CD).<br />

33<br />

(Jensen, 2010a)<br />

Punkt ∆E (mm) ∆N (mm) rE(mm) rN (mm)<br />

1.301 -1 0 -1 0<br />

2.301 -2 0 -2 0<br />

1.302 1 0 1 0<br />

2.302 -1 0 -1 0<br />

1.303 2 -1 2 -1<br />

2.303 1 0 1 0<br />

1.304 1 -2 1 -2


2.304 -1 1 -1 1<br />

Middel tE=0 tN≈0<br />

De korrigerede residualer er beregnet ud fra følgende formler:<br />

= ∆ − <strong>og</strong> = ∆ − <br />

Som det ses i tabellen er der ingen translationsparametre, hvilket betyder, at der ikke sker n<strong>og</strong>en ændring i<br />

residualerne. Punktspredningen for afsætningen vurderes til at være 2 mm på baggrund af testberegningen. Den<br />

relative spredning vurderes til at være ca. 1.7 mm.<br />

Grovfejlsgrænser er på den baggrund beregnet således (Jensen, 2010b):<br />

= = ±3 ∙ √2 √<br />

34<br />

= ±3 ∙ √2 0.002<br />

√8<br />

= 0.003<br />

= = ±3 ∙ √2 ∙ = ±3 ∙ √2 ∙ 0.0017 = 0.007<br />

Da der ikke er n<strong>og</strong>en translationsparametre, overholdes denne fejlgrænse. For residualerne gælder det desuden, at de<br />

alle ligger under fejlgrænsen.<br />

OPSAMLING<br />

Der er foretaget to afsætninger. Afsætning af skel <strong>og</strong> veje skete ved RTK-måling. <strong>Ved</strong> afsætningen var forskellen<br />

mellem designkoordinater <strong>og</strong> de kontrolmålte punkter alle indenfor fejlgrænsen på 0.051 m. Den beregnede<br />

spredning på forskellen var desuden under den forventede.<br />

I forbindelse med bygningsafsætningen blev et net etableret. Nettet, som på forhånd var indmålt med GPS, blev<br />

indmålt med totalstation, <strong>og</strong> observationerne blev udjævnet ved hjælp af Anblok for at få relativt gode koordinater. Til<br />

fastsættelse af koter til netpunkterne blev der foretaget et geometrisk nivellement, <strong>og</strong> en udjævning i et script blev<br />

foretaget. Den konstruerede bygning blev afsat i et lokalt system, hvor en efterfølgende kontrolopmåling viste, at alle<br />

residualerne lå indenfor fejlgrænsen.


FASE 3 – KORTLÆGNING VED FOTOGRAMMETRI<br />

Målet med fase 3 er en fremstilling af et teknisk kort med højdeinformation, en digital terrænmodel <strong>og</strong> to ortofotos<br />

ved fot<strong>og</strong>rammetrisk kortlægning. De fire kortlægningstyper er alle udarbejdet ud fra de to billeder 84B <strong>og</strong> 85B, som<br />

er blevet udleveret ved semesterstart.<br />

Fase 3 består af en kravspecifikation, en beskrivelse af billedkvaliteten, en relativ <strong>og</strong> absolut orientering, en<br />

beskrivelse af det fot<strong>og</strong>rammetriske kort, en beskrivelse af den digitale terrænmodel <strong>og</strong> en beskrivelse af<br />

fremstillingen af de to ortofotos.<br />

KRAVSPECIFIKATION<br />

Kravspecifikationen nedenfor beskriver krav til det tekniske kort, højdemodellen samt til de ortofotos, der genereres.<br />

Krav, der knytter sig til relativ <strong>og</strong> absolut orientering, vil blive gennemgået i de afsnit, der omhandler disse.<br />

TEKNISK KORT<br />

Med udgangspunkt i FOT-standarden vil der blive gennemgået n<strong>og</strong>le overordnede krav til registrering af det tekniske<br />

kort. Ligesom i Fase 1 stilles der krav til registreringen af objekter, men de vil d<strong>og</strong> variere lidt i forhold til den tidligere<br />

beskrivelse som følge af en anderledes registrering (se tabel). Alle objekter registreres i 2½D i UTM32/EUREF89.<br />

DVR90.<br />

Som udgangspunkt vælges at tage kote fra terræn ved siden af objekt eller på selve objektet, hvis det er<br />

sammenfaldende med terrænet i højden. Eventuelle afvigelser fremgår af beskrivelse.<br />

Objektgruppe Objekttype Geometri Beskrivelse Bemærkning<br />

Bygninger Bygning Flade Tagudhæng (Samt<br />

højder herfra)<br />

Carport Flade <strong>Ved</strong> tagudhæng<br />

(samt højder herfra)<br />

35<br />

Bygninger har en<br />

mindste størrelse på<br />

10 m 2 . Garager<br />

registreres <strong>og</strong>så, hvis<br />

de er<br />

sammenhængende<br />

med bygningen.<br />

Carporte registreres,<br />

hvis de er<br />

sammenhængende<br />

med bygningen.<br />

Trafik Vejkant Skillelinje<br />

Fortov Skillelinje<br />

Teknik Mast Punkt<br />

Nedløbsrist Punkt Rensebrønde <strong>og</strong><br />

nedløbsriste på<br />

parceller indmåles<br />

ikke.<br />

Brønddæksel Punkt<br />

Diverse Punkt Eks. Brandhaner<br />

Natur Levende hegn Linje Mindst 2 meter langt<br />

Træ Punkt Træ defineres, som<br />

én stamme, der går i


Trægruppe Polygon med<br />

punkter<br />

36<br />

jorden.<br />

Hvis der er mere end<br />

3 træer, som ikke er<br />

på linje, registreres<br />

dette som<br />

trægruppe<br />

Plankeværk Linje Mindst 2 meter langt<br />

Hegn/Stendige Linje Mindst 2 meter langt Diger, som er under<br />

2 meter bredt <strong>og</strong> 50<br />

meter langt,<br />

registreres som hegn<br />

På billederne kan det være sværere at se, hvad der befinder sig på jorden end ved RTK-opmåling. Eksempelvis kan det<br />

være svært at kende forskel på træer <strong>og</strong> buske. D<strong>og</strong> er kendskabet til området så stort, at der ikke må være n<strong>og</strong>et, der<br />

afviger.<br />

Med hensyn til fuldstændighed afviger denne kravspecifikation fra den i fase 1 ved, at der kun er krav til de objekter,<br />

som er synlige på billederne. Igen stilles der høje krav til fuldstændighed grundet dets ringe udbredelse. Der vil derfor<br />

ikke være n<strong>og</strong>en tolerance i forhold til manglende objekter, som er synlige på de to fotos.<br />

DTM-MODEL<br />

Den genererede DTM-model skal indeholde hele det overordnede projektområde <strong>og</strong> skal dermed kunne bruges til<br />

generering af et ortofoto. Derudover indgår DTM-modelen i fase 4 ved en kontrol af højdemodeller. Den forventede<br />

spredning vil blive beregnet længere nede. Højdemodellen vil endvidere blive kontrolleret i forhold til de 15<br />

kontrolpunkter, som blev indmålt under fase 1.<br />

ORTOFOTO<br />

Der skal genereres tre små ortofotos, som viser opløsningen ved forskellige pixelstørrelser. Dette skal indgå i en<br />

vurdering af, hvilken pixelstørrelse, der er den optimale pixelstørrelse for de overordnede ortofotos. Det skal resultere<br />

i to overordnede ortofotos. Ét ortofoto genereres ved hjælp af den beregnede højdemodel, <strong>og</strong> det andet beregnes ud<br />

fra BLOM’s højdemodel. Der er som udgangspunkt ikke n<strong>og</strong>le krav til de genererede ortofotos, men de vil d<strong>og</strong> indgå i<br />

flere kortsammenligninger i fase 4.<br />

BILLEDEKVALITET<br />

Kalibreringsrapporten tilhørende billede 84B <strong>og</strong> 85B er udarbejdet den 8. maj 2008. Billeder er taget med kameraet<br />

UltraCam X., som har en kamerakonstant på 100.5 mm.<br />

Billederne er taget tirsdag den 7. april 2009, hvilket er udregnet ud fra bingofilen (Bilag – Bingofil). De nøjagtige<br />

tidspunkter for billedernes optagetidspunkt er for 84B 07/04-2009 kl. 11:35:32.737 <strong>og</strong> 85B 07/04-2009 kl.<br />

11:35:35.142. Der skal lægges to timer til sekundtallet i Bingofilen, da GPS er i forhold til Greenwich Mean Time, hvor<br />

Danmark er +1 i tidszone. Desuden er billedet taget, mens det var sommertid, hvilket der <strong>og</strong>så skal tages højde for.<br />

På dette tidspunkt var solvinklen i området 37.4°, hvilket fremgår af bingo-filen. Tidspunkter i tabellen svarer d<strong>og</strong> ikke<br />

til sommertid. Derfor trækkes en time fra ovenstående tidspunkter for at se solens vinkel på billedernes<br />

optagetidspunkt.


Overlappet skal være i nærheden af 60 %, hvilket er det planlagte overlap. Det reelle overlap for hvert billede ses af<br />

nedenstående tabel.<br />

84B 85B<br />

Venstre 56.3% 55.0%<br />

Højre 55.6% 55.9%<br />

Pixelstørrelsen beregnes ud fra en afstand målt diagonalt i det ene billede ved hjælp af ERviewer <strong>og</strong> GeoCAD.<br />

Afstanden er taget mellem to riste målt i GeoCAD <strong>og</strong> den skrå pixelafstand målt i ERviewer. Dette giver følgende:<br />

= 527.614<br />

= 0.0329<br />

16030.59<br />

Dette betyder, at billeder har en Ground Sampling Distance (GSD) på 3.29 cm, hvilket repræsenterer afstanden for en<br />

pixel i virkeligheden.<br />

Billederne er endvidere kontrolleret for skyer, dis, skygger, pixelfejl <strong>og</strong> hotspots. Der blev ikke fundet n<strong>og</strong>et<br />

utilfredsstillende i søgningen i billeder for disse fejl. Der eksisterer selvfølgelig skygger, men ikke n<strong>og</strong>le nævneværdige.<br />

Billederne er utroligt skarpe, hvilket <strong>og</strong>så underbygges af den lave pixelstørrelse. Billedernes hist<strong>og</strong>rammer viser<br />

endvidere en symmetrisk kurve, hvilket underbygger, at der ikke er n<strong>og</strong>en hotspots. Begge billeder har farver indenfor<br />

hele farvespektret i RGB. Farverne fordeler sig n<strong>og</strong>enlunde ligeligt.<br />

Billede 3: Billede hist<strong>og</strong>rammer (84b tv. <strong>og</strong> 85B th.)<br />

RELATIV ORIENTERING<br />

Den relative orientering er blevet foretaget i pr<strong>og</strong>rammet ImageStation Digital Mensuration. Billederne skal flyttes<br />

relativt i forhold til hinanden for at fjerne x <strong>og</strong> y-parallakser i billederne. En ændring i x vil have en direkte indvirkning<br />

på z. Dette sker for at danne en fot<strong>og</strong>rammetrisk model, hvor man kan se stereo.<br />

37


Billede 4: Illustration af en flytning i x’<br />

Pr<strong>og</strong>rammet danner en billedpyramide, hvor billederne bliver komprimeret til færre pixels for lettere at kunne<br />

håndtere datamængden i pr<strong>og</strong>rammet.<br />

Billede 5: Illustration af billedpyramide<br />

I den relative orientering fastholdes det venstre billede, mens det højre billede flyttes <strong>og</strong> indpasses. <strong>Ved</strong><br />

overbestemmelse sker indpasningen efter mindste kvadraters princip.<br />

KRAV TIL RELATIV ORIENTERING<br />

Den forventede nøjagtighed på den relative orientering beregnes ud fra flere fejlbidrag. Dels er der en restfejl fra<br />

kalibreringen på kameraet. Denne vurderes til 1 µm per billede. Derudover er der fejl i måling af y-parallakser, som<br />

sættes til 1/3 pixel i den manuelle måling <strong>og</strong> 1/5 pixel i den automatiske ved måling på veldefinerede punkter. Da alle<br />

punkter ikke forventes at være lige veldefinerede ganges fejlbidraget på måling af y-parallakse med 1.5. Derudover er<br />

der et refraktionsbidrag, som er afhængigt af flyvehøjden. Refraktionsbidraget beregnes som en differens mellem to<br />

stråler svarende til 1 µm ganget med flyvehøjden på ca. 500 meter. Dette giver et bidrag på 0.5 µm. Det vil sige, at den<br />

forventede spredning på den relative orientering manuelt <strong>og</strong> automatisk er:<br />

38


, = 1 ∙ 2 + 1<br />

<br />

∙ 7.2 ∙ 1.5 + 0.5<br />

3 = 3.9<br />

, = 1 ∙ 2 + 1<br />

<br />

∙ 7.2 ∙ 1.5 + 0.5<br />

5 = 2.6<br />

Der er foretaget en relativ orientering manuelt <strong>og</strong> automatisk.<br />

Manuel orientering<br />

39<br />

(Juhl, 2011b)<br />

I pr<strong>og</strong>rammet er der manuelt udvalgt fem sammenknytningspunkter både i højden <strong>og</strong> i planen, som måles i hvert<br />

billede for at danne den bedste indpasning. Resultatet af den manuelle relative orientering gav efter måling af<br />

sammenknytningspunkter i begge billeder en RMS på 0.6 µm <strong>og</strong> en y-parallakse på 1.4 µm, hvilket vurderes som<br />

acceptabelt.<br />

Automatisk orientering<br />

I den automatisk relative orientering vælger pr<strong>og</strong>rammet selv n<strong>og</strong>le sammenknytningspunkter i højden <strong>og</strong> planet, som<br />

måles i hvert billede. Pr<strong>og</strong>rammet finder automatisk de steder, hvor punkter stemmer bedst overens <strong>og</strong> laver derefter<br />

en indpasning af billederne.<br />

Pr<strong>og</strong>rammet laver en matching af punkter. som derved giver en spredning ved RMS. Resultatet af den automatiske<br />

relative orientering gav en spredning ved RMS på 0.7 µm. Derudover resulterede orienteringen i en y-parallakse på 2.2<br />

µm.<br />

Resultat<br />

Det kan konkluderes, at begge orienteringer overholder den forventede spredning både manuelt <strong>og</strong> automatisk.<br />

Resultatet kan variere mellem 1.5-4 µm (Juhl, 2011b).<br />

Den automatiske orientering er valgt til videre bearbejdning på trods af resultatet er dårligere end den manuelle<br />

orientering. Den er valgt på baggrund af en vurdering af, at der i den automatiske orientering er brugt flere punkter i<br />

indpasningen <strong>og</strong> derfor mindre sandsynlighed for at lave fejl.<br />

ABSOLUT ORIENTERING<br />

I den absolutte orientering indpasses billederne efter verdenskoordinater. Dette gøres ved at transformere den<br />

relative model over et antal RTK-indmålte paspunkter. Den absolutte orientering er <strong>og</strong>så lavet i pr<strong>og</strong>rammet<br />

ImageStation Digital Mensuration. Orienteringen er foretaget ved at registrere paspunkternes placering i begge<br />

billeder samtidig.


Billede 6: Illustration af proces ved absolut orientering<br />

KRAV TIL ABSOLUT ORIENTERING<br />

Der er både en forventet spredning i planen <strong>og</strong> i højden. Den forventede nøjagtighed på den absolutte orientering i<br />

planen beregnes ud fra flere fejlbidrag. Dels er der en restfejl fra kalibreringen på kameraet. Denne vurderes til 1 µm<br />

per billede. Derudover er der fejl i måling af veldefinerede punkter med en pixelstørrelse på 3.5 cm, som sættes til 1/6<br />

pixel. Så er der refraktionsbidrag, som i dette tilfælde er 1.5 µm, som er afhængig af flyvehøjden (ca.500m), hvilket<br />

giver en refraktionsbidrag på 0.8µm. Paspunkterne giver endvidere et fejlbidrag, som svarer til GPS’ens nøjagtighed på<br />

ca. 1 cm. Sidst er der en definitionsafvigelse mellem fot<strong>og</strong>rammetri <strong>og</strong> GPS, som sættes til 1.5cm. Dette kan d<strong>og</strong><br />

variere i forhold til, hvor veldefineret punktet er. Spredningen skal ses i forhold til målforholdet, som beregnes til:<br />

åℎ:<br />

, =<br />

∙ <br />

ø()<br />

<br />

∙ <br />

<br />

<br />

= =4564. Det vil sige, at målforholdet er 1:4564.<br />

.<br />

+ <br />

∙ 3.5 + 0.8 ∙<br />

<br />

+ 1 + 1.5 ≈ 2 (Juhl, 2011b)<br />

I højden er der samme fejlbidrag som i planen bortset fra, at der skal tages hensyn til højde/basisforholdet.<br />

Basisforholdet beregnes ud fra koordinatafstanden mellem de to optagecentre i billederne:<br />

= (712649.67 − 712575.965) + (6179398.676 − 6179283.302) =136.9<br />

Flyvehøjden bliver beregnet ud fra målforholdet <strong>og</strong> kamerakonstanten på 100.5 mm:<br />

ℎø =<br />

Dermed bliver højdebasisforholdet ℎø = <br />

4542<br />

1 ∙<br />

,ø =<br />

10000<br />

2<br />

<br />

∙ 3.35<br />

.<br />

4564 ∙ 100.5<br />

1000<br />

= 3.35.<br />

40<br />

= 459 <br />

+ 1<br />

<br />

∙ 3.5 ∙ 3.35 + 0.8 ∙<br />

6 4542<br />

<br />

∙ 3.35 + 1<br />

10000 + 1.5 = 3.1


PASPUNKTER<br />

41<br />

(Juhl, 2011b)<br />

Til den absolutte orientering blev der udvalgt otte højde- <strong>og</strong> planpaspunkter, som er dobbeltmålte <strong>og</strong> midlet for at<br />

finde de bedste koordinater. D<strong>og</strong> var der flere punkter, hvor dobbeltmålingen gav en stor afvigelse, <strong>og</strong> dermed var<br />

kun en af målingerne gældende. Derfor var det kun den ene måling, som indgik med mindre vægt end de<br />

dobbeltmålte. Derudover blev der målt to højdepaspunkter, hvoraf det ene var dobbeltmålt. Punkterne blev udvalgt<br />

på forhånd for at sikre, at de kunne ses i modellen. Paspunkter er lagt således, at det giver en god bestemmelse af<br />

Billede 7: Oversigt over paspunkters<br />

placering<br />

planen <strong>og</strong> højden indenfor bebyggede områder <strong>og</strong> en god bestemmelse af<br />

højden udenfor bebyggede områder, hvilket fremgår af nedenstående billede.<br />

Paspunkterne er nummeret således, at højdepaspunkter har numre fra 600-<br />

699, <strong>og</strong> plan- <strong>og</strong> højdepaspunkter har numre fra 700-799. Paspunkternes<br />

specifikke placering <strong>og</strong> beskrivelse fremgår af Bilag – Dokumentation for<br />

paspunkter. Paspunkterne er udvalgt, så de ligger i udkanten af modellen d<strong>og</strong><br />

med undtagelse af paspunkterne nederst i billedet, hvor det har været svært at<br />

finde n<strong>og</strong>et veldefineret. De er valgt, så det meste af modellen bliver korrigeret<br />

for flyets drejninger (ϕ,κ <strong>og</strong> ω).<br />

Resultat<br />

Før den absolutte orientering kan laves, skal en tekst-fil indlæses med<br />

paspunkternes koordinater samt spredninger, så orienteringen kan foretages.<br />

Punkterne er i ImageStation målt flere gange for at få en bedre nøjagtighed.<br />

Som nævnt under kravspecifikationen er den forventede nøjagtighed for<br />

planen <strong>og</strong> højden henholdsvis 2 <strong>og</strong> 3.2 cm. Nøjagtigheden af den absolutte<br />

orientering er både et resultat af den relative orientering <strong>og</strong> orienteringen efter<br />

paspunkter. Den relative orientering resulterede i en y-parallakse på 2.2 µm.<br />

Endvidere blev der givet følgende værdier:<br />

RMS<br />

X=0.013 m<br />

Y=0.011 m<br />

Z=0.021 m<br />

Resultatet i ImageStation angives ikke i spredninger, men i Root Mean Square (RMS).<br />

= ∑ <br />

<br />

<br />

<br />

= ∑ <br />

<br />

For at omregne RMS til spredning skal der derfor omregnes ved hjælp af en faktor k<br />

= ∙ → ∑ <br />

<br />

− = ∙ ∑ <br />

<br />

<br />

1 1 <br />

→ = ∙ → =<br />

− −


Derved bliver spredningen i den absolutte orientering:<br />

= √ ∙ <br />

I dette tilfælde er m=8 i planen <strong>og</strong> m=10 i højden. N er i planen 2 <strong>og</strong> 3 i højden.<br />

Dermed bliver spredninger på X, Y <strong>og</strong> Z:<br />

= <br />

∙ .<br />

=0.015m<br />

= <br />

∙ .<br />

=0.013m<br />

= <br />

∙ .<br />

=0.025m<br />

Og spredningen i planen:<br />

_ = 0.015 + 0.013 = 0.014 <br />

2<br />

Det vil sige, at de forventede spredninger fra kravspecifikationen på 2 <strong>og</strong> 3.1 cm overholdes. Derved er den absolutte<br />

orientering gået bedre end forventet, hvilket ses, da de beregnede spredninger er lavere end de forventede<br />

spredninger.<br />

De to billeder har følgende drejninger i Omega, Phi <strong>og</strong> Kappa:<br />

ω ϕ κ<br />

Billede 84 -0.033gon -0.407gon -135.718gon<br />

Billede 85 -0.008gon -0.427gon -135.673gon<br />

KONTROL AF ORIENTERING<br />

Den absolutte orientering kan sammenlignes med COWI’s orientering, som det fremgår af bingo-filen. Endvidere<br />

kontrolleres den i forhold til en beregning af den ydre orientering, som er foretaget i et Matlab-script for billede 84B<br />

<strong>og</strong> 85B.<br />

84B Absolut orientering COWI’s orientering Ydre Orientering<br />

X0 712649.670 712649.258 712649.576<br />

Y0 6179398.676 6179399.615 6179399.117<br />

Z0 485.188 484.999 485.309<br />

ω -0.033 gon -0.051 gon 0.031 gon<br />

ϕ -0.407 gon -0.419 gon -0.392 gon<br />

κ -135.718 gon -135.655 gon -135.711 gon<br />

I den ydre orientering er der beregnet en spredning på de ubekendte. Dette gøres ved hjælp af scriptet opretning.m<br />

(Bilag 22, Oversigt over bilag på CD), hvor spredninger på den ydre orientering beregnes ud fra kovariansmatricen. Da<br />

der samtidig er givet en spredning på COWI’s orientering fra bingofilen, kan der beregnes en fejlgrænse for, hvor<br />

meget afvigelsen må være mellem COWI’s <strong>og</strong> den ydre orientering. Den forventede spredning på forskellen beregnes<br />

42


som: ∆(,,,,,) = (,,,,,) + (,,,,,) . Fejlgrænsen kan dermed opstilles som:<br />

±3 ∙ ∆(,,,,,).<br />

Kontrol Cowi kontra ydre Spredning (ydre) Spredning(COWI) Fejlgrænse<br />

X0 -0.318 0.42 0.03 1.263<br />

Y0 0.498 0.42 0.03 1.263<br />

Z0 -0.310 0.16 0.03 0.488<br />

ω -0.082 0.007 0.004 0.0242<br />

ϕ -0.027 0.007 0.004 0.0242<br />

κ 0.056 0.002 0.01 0.0306<br />

Spredningen på vægtenheden giver i scriptet 0.016 mm, hvilket er et udtryk for, hvor godt der måles i billedet.<br />

Den ydre orientering af billede 84B passer i X0, Y0 <strong>og</strong> Z0 sammen med COWI’s orientering, da afvigelserne alle er under<br />

fejlgrænsen. D<strong>og</strong> er afvigelserne i ω, φ <strong>og</strong> κ alle over fejlgrænsen, hvilket betyder, at der er en uoverensstemmelse i<br />

optagecentrenes position <strong>og</strong> flyets drejninger <strong>og</strong> vipninger.<br />

Og for billede 85B:<br />

85B Absolut orientering COWI’s orientering Ydre orientering<br />

X0 712575.965 712575.654 712575.400<br />

Y0 6179283.302 6179284.39 6179284.197<br />

Z0 483.672 483.47 483.599<br />

ω -0.008 gon -0.0452 gon -0.005 gon<br />

ϕ -0.427 gon -0.4287 gon -0.469 gon<br />

κ -135.673 gon -135.6093 gon -135.667 gon<br />

Igen benyttes scriptet til at beregne spredningerne for den ydre orientering af billede 85B, <strong>og</strong> der opstilles fejlgrænser.<br />

Kontrol Cowi kontra ydre Spredning (Ydre) Spredning(COWI) Fejlgrænse<br />

X0 0.254 0.09 0.03 0.285<br />

Y0 0.193 0.11 0.03 0.342<br />

Z0 -0.129 0.036 0.03 0.141<br />

ω -0.0402 0.002 0.004 0.0134<br />

ϕ 0.0403 0.002 0.004 0.0134<br />

κ 0.0577 0.001 0.01 0.0301<br />

Spredning på vægtenheden bliver 0.004 mm.<br />

Den ydre orientering af billede 85B passer i X0, Y0 <strong>og</strong> Z0 sammen med COWI’s orientering, da afvigelserne alle er under<br />

fejlgrænserne. D<strong>og</strong> er afvigelserne i ω, φ <strong>og</strong> κ alle over fejlgrænsen, hvilket betyder, at der er en uoverensstemmelse i<br />

flyets drejninger <strong>og</strong> vipninger.<br />

Grunden til, at ω, φ <strong>og</strong> κ i begge billeder ikke passer sammen med COWI’s orientering, er, at der ikke skal en stor<br />

drejning eller vipning af optagecenteret til at rykke meget på disse parametre.<br />

Det kan ud fra koordinaterne konkluderes, at den ydre orientering er acceptabel med hensyn til optagecentrets<br />

position, da forskellen på den ydre orientering <strong>og</strong> COWI’s orientering med hensyn til , <strong>og</strong> Z er indenfor<br />

fejlgrænserne, som er beregnet på baggrund af scriptet opretning.m <strong>og</strong> COWI’s orientering, som givet i bingofilen. D<strong>og</strong><br />

er forskellen i flyets vipninger , <strong>og</strong> over fejlgrænsen.<br />

43


Der blev undervejs konstateret en stor afvigelse mellem den absolutte orientering <strong>og</strong> COWI’s orientering. Denne fejl<br />

beror på en forkert indtastning af værdier for forskydning af kameraets midtpunkt. Dette betyder, at kameraets<br />

midtpunkt bliver flyttet i den forkerte retning, hvilket giver en stor afvigelse i forhold til COWI’s orientering. D<strong>og</strong> vil<br />

fejlen ikke have indflydelse på de efterfølgende produkter, da drejninger af flyet ikke er store, <strong>og</strong> fejlen samtidig er<br />

den samme for begge billeder.<br />

FOTOGRAMMETRISK KORT<br />

Det tekniske kort udarbejdet ved hjælp af fot<strong>og</strong>rammetri skal tage udgangspunkt i FOT-standarden, ligesom tilfældet<br />

var under fremstilling af det tekniske kort i fase 1. Det tekniske kort skal fremstilles i 2½D, hvilket betyder, at det<br />

tekniske kort skal være i 2D med en kote til hvert indmålt punkt. Det fot<strong>og</strong>rammetriske kort skal være fremstillet over<br />

samme område, som det tekniske kort fremstillet under fase 1. Når man ser nærmere på kravspecifikation vil det<br />

fremgå, at det stort set er det samme tillæg til FOT-standarden, som danner udgangspunkt for begge tekniske kort.<br />

D<strong>og</strong> er der små ændringer, som beror på den anden type af registrering. FOT-standarden henvender sig hovedsaligt til<br />

tekniske kort dannet ved fot<strong>og</strong>rammetrisk opmåling.<br />

GEOMETRISK KONSTRUKTION<br />

Inden indmålingen af terrængenstande begynder, skal der i første omgang dannes en kodetabel, således at de<br />

indmålte data defineres med en kode for hver objekttype. Dette gøres, så data fra Microsystem kan konverteres til<br />

GeoCAD-format. Data indmåles i første omgang i ImageStation Stereo Display. Derefter eksporteres data som en dxffil,<br />

hvorefter data bliver konverteret <strong>og</strong> indlæst i GeoCAD for at færdiggøre det tekniske kort med de rette signaturer.<br />

Data indmåles i stereo i forhold til den absolutte orientering, så punkterne får de rigtige UTM32/EUREF89-koordinater<br />

<strong>og</strong> koter i DVR90. Dette indebærer, at det er vigtigt at være opmærksom på, at krydset i ImageStation er placeret<br />

enten på det objekt, som skal måles i højden eller i terræn for genstande, som skal indmåles ved terræn.<br />

Nedenstående kort er resultatet af den fot<strong>og</strong>rammetriske opmåling:<br />

44


Kort 14: Fot<strong>og</strong>rammetrisk teknisk kort<br />

KONTROL OG VURDERING<br />

Det tekniske kort skal kontrolleres i forhold til de 15 kontrolpunkter, som blev indmålt med RTK under fase 1. Den<br />

forventede spredning i planen er afhængig af flere faktorer. For det første er der et bidrag fra spredningen på den<br />

absolutte orientering på 0.014 m. Derudover er der et bidrag for, hvor godt der måles i stereo, hvilket afhænger af,<br />

hvor god man er til at se stereo, hvor detaljeret man er omkring indmålingen af punkter, <strong>og</strong> hvor erfaren man er til<br />

måling i stereo. Dette bidrag sættes til 1, som skal omsættes til virkeligheden ved at multiplicere med GSD’en. Den<br />

forventede nøjagtighed af det fot<strong>og</strong>rammetriske tekniske kort i planen bliver dermed:<br />

= (1 ∙ 0.0327) + 0.014 = 0.036 <br />

45<br />

(Juhl, 2011b)<br />

Den forventede spredning i højden er afhængig af den absolutte orienterings bidrag i højden på 0.025 m. Endvidere er<br />

spredningen i højden afhængig af, hvor god man er til at se stereo, hvor detaljeret man er omkring indmålingen af<br />

punkter, <strong>og</strong> hvor erfaren man er til måling i stereo. Dette bidrag sættes igen til 1 <strong>og</strong> skal multipliceres med<br />

højde/basisforholdet på 3.5 cm.<br />

_ø = (1 ∙ 0.035) + 0.025 = 0.043 <br />

(Juhl, 2011b)<br />

FOT specificerer, at deres kort har en nøjagtighed på 10 cm ved en GSD på 10 cm. Tidligere blev de 10 cm delt med √2<br />

for at få en spredning tilsvarende de andre. Derfor divideres GSD’en på 3.27 cm med √2 for at få en estimeret<br />

spredning på FOT-kort med tilsvarende GSD.<br />

_ = 0.0327<br />

= 0.023<br />

√2


I højden er nøjagtigheden ganget med faktor 1.5, hvilket betyder, at GSD’en skal ganges med denne værdi<br />

(FOTDanmark, 2010).<br />

0.0327 ∙ 1.5<br />

_ø = = 0.035<br />

√2<br />

Spredninger på det tekniske kort er beregnet ved hjælp af afstande til indmålte kontrolpunkter foretaget i GeoCAD<br />

(Bilag 21, Oversigt over bilag på CD). Der er både foretaget direkte sammenligninger, men det er <strong>og</strong>så foretaget en 1D<br />

<strong>og</strong> 2D translation. Resultatet af translationerne er:<br />

Spredninger Højde Plan<br />

Før translation σh=0.251 m σEN=0.068 m<br />

Efter translation σh=0.113 m σEN=0.033 m<br />

I højden har translationen den indflydelse, at højden bliver flyttet 22.5 cm. Der har været en ensidig fejl i placeringen<br />

af målemærket, hvilket har resulteret i store afvigelser mellem de kontrollerede brønddæksler <strong>og</strong> riste.<br />

Dette er både højere end den forventede spredning <strong>og</strong> nøjagtighedskravet fra FOT-standarden. Unøjagtigheden<br />

skyldes stor usikkerhed ved placering af målemærket på terræn. Endvidere har den relative <strong>og</strong> absolutte orientering<br />

en indvirkning, da målingen afhænger af, hvor godt modellen er indpasset i højden.<br />

I planen undersøges spredninger først uden translation. Denne er 6.8 cm, hvilket er over den forventede værdi. Derfor<br />

forsøges en 2D-translation, hvilket bevirker, at spredningen i planen bliver 3.3 cm. Spredningen er under den<br />

forventede spredning, d<strong>og</strong> er den højere end den estimerede FOT-nøjagtighed.<br />

DTM-MODEL<br />

Højdemodellen genereres ud fra den absolutte orientering i pr<strong>og</strong>rammet Image Station Automatic Elevation.<br />

Højdemodellen er en digital terrænmodel, hvor terrænkoter automatisk registreres i pr<strong>og</strong>rammet.<br />

KRAV TIL DTM-MODEL<br />

DTM-modellen skal genereres for hele overlappet mellem de to billeder med en gridstørrelse på 1.5 m, da dette svarer<br />

til gridstørrelsen i de højdemodeller, der skal kontrolleres med i fase 4.<br />

I nøjagtigheden af den endelige højdemodel tages der afsæt i FOT-specifikationen. Her ganges GSD med en faktor 1.5<br />

for at få den ønskede nøjagtighed i cm. For en GSD på 3.27 cm giver det en nøjagtighed på 4.91 cm. D<strong>og</strong> tages der ikke<br />

hensyn til denne beregning. Den forventede spredning på de genererede punkter i forbindelse med DTM-modellen<br />

kan udregnes som:<br />

= 0.00015 ∙ ℎø = 0.00015 ∙ 459 = 0.069 m<br />

46<br />

(Juhl, 2011b)<br />

Denne nøjagtighed svarer til den nøjagtighed, som vil blive beregnet ud fra de forskellige bidrag herunder<br />

kalibreringen af kameraet, bidraget fra den absolutte orientering, hvor godt der måles i stereo <strong>og</strong> refraktionen.


Kort 15: Fot<strong>og</strong>rammetrisk digital terrænmodel<br />

GEOMETRISK KONSTRUKTION<br />

Terræntypen sættes til ’Flat’, da det er et meget fladt terræn. Det har betydning for længden af de epipolære linier<br />

(Epipolar Line Distance) <strong>og</strong> dermed, hvad der tillades at blive kaldt terræn. Samtidig blev Smoothing Weight sat til 1.5,<br />

hvilket har betydning for udjævning af terrænet <strong>og</strong> dermed, hvor meget af det ”lokale” der udelades.<br />

Der opstår en punktsky, som er terrænmålte punkter. Punkterne har forskellige farver alt efter, hvor godt de er<br />

bestemt. Orange betyder, at det er godt bestemt. Lilla betyder, at det er kritisk <strong>og</strong> blå, at det er dårligt bestemt. De<br />

genererede højder er orange på nær hustagene, som i n<strong>og</strong>le tilfælde er lilla.<br />

47


I GeoCAD importeres data som en que-fil, <strong>og</strong> derefter skabes en TIN-model ud fra de målte punkter for at opnå en<br />

spredning, som fremgår af kravspecifikationen. Der er dannet en udtyndet TIN-model med pilehøjde på 15 cm, hvilket<br />

giver en spredning på 0.056 på de punkter, som ikke indgår i nettet.<br />

KONTROL OG VURDERING AF HØJDEMODEL<br />

Højderne kontrolleres med de 15 kontrolmålte punkter fra fase 1. De er blevet interpoleret til TIN-modellen for at få<br />

en kotedifference. Det har ført til en beregning af en spredning på TIN-modellen både ved direkte afvigelser <strong>og</strong><br />

translation. Afvigelserne giver en spredning, som begge fremgår af nedenstående skema.<br />

Som det fremgår, overholder TIN-modellen den forventede spredning, som blev udregnet tidligere til 6.8 cm. D<strong>og</strong><br />

overholder modellen ikke den estimerede FOT-nøjagtighed hverken ved direkte afvigelser eller translation. En metode<br />

48


til at overholde FOT-standarden kunne være at sætte pilhøjden ned. Endvidere kunne den højere spredning skyldes, at<br />

der er flere usikre punkter omkring husene, som giver en større unøjagtighed i højdemodellen.<br />

ORTOFOTO<br />

Ortofotoerne genereres i pr<strong>og</strong>rammet ImageStation Base Rectifier. Her benyttes de absolut orienterede billeder <strong>og</strong><br />

højdemodellen til at genere ortofotos. Pixelstørrelsen defineres, hvorefter ortofotoet genereres. Ortofotos generes<br />

kun ud fra et af billederne.<br />

Der generes tre ortofotos over et mindre område med pixelstørrelser på 6 cm, 4 cm <strong>og</strong> 3 cm for at vurdere, hvilken<br />

pixelstørrelse, der giver den bedste opløsning. Billederne er taget med en pixelstørrelse på ca. 3.3 cm.<br />

Billede 8: 6 cm pixels over område ved teknisk kort<br />

49


Billede 9: 4 cm pixels over område ved teknisk kort<br />

Billede 10: 3 cm pixels over område ved teknisk kort<br />

Der skal herefter generes to endelige ortofotos. Det ene ortofoto bliver genereret ud fra DTM-modellen dannet<br />

fot<strong>og</strong>rammetrisk, hvilket indebærer en opretning af ortofotoet ud fra de højder, som er generet i DTM-modellen. Det<br />

andet ortofoto genereres ud fra BLOM’s højdemodel.<br />

50


Når man sammenligner de forskellige ortofotos ses der en nævneværdig ændring i opløsningen mellem 6 <strong>og</strong> 4 pixels.<br />

Desuden ses der en lille forandring i opløsning fra 4 til 3 pixels. Det mest optimale ville derfor være at vælge i<br />

nærheden af 3 pixels i pixelstørrelse.<br />

D<strong>og</strong> skal det endelige ortofoto konstrueres således, at størrelsen på filen ikke overstiger 500 Mb, da GeoCAD ellers<br />

ikke kan konvertere ortofotoet til ECW-fil, hvilket er væsentligt nemmere at håndtere. Derfor vælges den endelige<br />

pixelstørrelse til 4.2, da dette giver den laveste pixelstørrelse i forhold til størrelsen af filen.<br />

Højdemodellen <strong>og</strong> dermed ortofotoet bliver dårligere jo længere væk fra nadir-punktet, man er. Dette sker, da<br />

orienteringen kun er sket ud fra to billeder <strong>og</strong> ikke en aerotriangulation af et større område. <strong>Ved</strong> en aerotriangulation<br />

ville fejlene i kanten af modellen blive oprettet af det næste billede i modellen. D<strong>og</strong> ville der stadig være en fejl i<br />

yderkanten af blocken. Desuden opstår der problemer i nærheden af hustage, hvor højdemodellen er dårligere<br />

korreleret. Denne fejl kunne afhjælpes ved hjælp at brudlinjer. D<strong>og</strong> er dette tidskrævende.<br />

Der ses øjensynligt en forskel mellem de to ortofotos lavet på de to højdemodeller. Objekter på det ortofoto, der er<br />

genereret ud fra egne højder, er meget forvrængede, hvorimod det, der er genereret med BLOM’s højdedata, er mere<br />

lige. En endelig sammenligning af de to ortofotos sker i fase 4.<br />

OPSAMLING<br />

I fase 3 har det ud fra to digitale billeder været hensigten at lave n<strong>og</strong>le forskellige produkter, der kan bruges ved<br />

projektering. Der er digitaliseret et teknisk kort, lavet en højdemodel <strong>og</strong> produceret to ortofotos.<br />

Det tekniske kort er digitaliseret efter, at der er foretaget en relativ <strong>og</strong> en absolut orientering af billederne. Det har<br />

ført til digitalisering af objekterne, som er sammenlignet med tidligere RTK-indmålte kontrolpunkter. Det viser en<br />

nøjagtighed, som er dårligere end den forventede i højden. Dette kan skyldes manglende erfaring i digitaliseringen.<br />

Endvidere viste spredningen i planen en dårligere nøjagtighed end forventet. D<strong>og</strong> var fejlen efter en translation under<br />

fejlgrænsen. Dette kan begrundes med unøjagtighed i forbindelse med den absolutte orientering.<br />

DTM-modellen er genereret automatisk ud fra de to billeder. Dens nøjagtighed vil blive vurderet i fase 4. Sidst er der<br />

genereret to ortofotos. N<strong>og</strong>le prøveberegninger viste, at en pixelstørrelse på ca. 3 cm giver det bedste resultat i<br />

forhold til billederne, men at der ikke er den store forskel mellem 3 <strong>og</strong> 4 cm pixels. D<strong>og</strong> er der af praktiske årsager<br />

valgt at beregne de endelige ortofotos med en pixelstørrelse på 4.2 cm. Det har resulteret i to ortofotos, hvor der i det<br />

ene ortofoto er brugt højder fra den genererede DTM-model, <strong>og</strong> i det andet er brugt BLOM’s højdedata. Det ses<br />

umiddelbart, at højdemodellerne har stor indvirkning på kvaliteten af ortofotoet.<br />

51


FASE 4 – VURDERING AF FORSKELLIGE KORTLÆGNINGSMETODER<br />

Formålet med denne fase er at vurdere forskellige kortprodukter indbyrdes <strong>og</strong> spredningen på deres<br />

koordinatdifferencer. Kravet til sammenligning fremgår af specifikationen for nøjagtigheder. Stjerne bag et navn<br />

angiver, at kortet er produceret i en tidligere fase.<br />

Fase 4 består af en specifikation af nøjagtigheder <strong>og</strong> 11 kortsammenligninger mellem forskellige kortlægninger, som<br />

vist i nedenstående tabel. Til hver sammenligning fremgår en beskrivelse af sammenligningen, et søjlediagram over<br />

spredninger samt et residualplot.<br />

RTK-TK* RTK- Foto-TK* Foto- Orto- Orto- FOT/TK3 COWI - COWI-<br />

DTM*<br />

DTM* Egen* BLOM*<br />

DDH DDO<br />

RTK-TK* x x x X x<br />

RTK-<br />

DTM*<br />

x<br />

Foto-TK* X<br />

Foto-<br />

DTM*<br />

x<br />

Orto-<br />

Egen*<br />

X<br />

Orto-<br />

BLOM*<br />

X x<br />

SPECIFIKATION AF NØJAGTIGHEDER<br />

Differencerne må vurderes i forhold til en spredning beregnet på baggrund af den forventede nøjagtighed af de to<br />

kort. Vurdering i planen bliver foretaget ud fra veldefinerede punkter for eksempel riste eller brønddæksler. Dette<br />

gøres, da definitionen af for eksempel huse varierer for de forskellige kortlægninger. En vurdering af de beregnede<br />

spredninger gøres på baggrund af nedenstående formel, der beskriver den forventede spredning:<br />

<br />

, = + Nøjagtighed for kortprodukter fremgår af nedenstående tabel. Spredning for RTK-TK, RTK-DTM, Foto-TK <strong>og</strong> Foto-DTM<br />

er beregnet tidligere på baggrund af differencer til kontrolpunkter.<br />

Plan Kote<br />

COWI’s højdemodel DDH - 25cm<br />

FOT-TK3 7cm 7.5cm<br />

COWI’s DDO<br />

25<br />

= 18<br />

√2<br />

RTK-TK* 1.1cm 1.3cm<br />

RTK-DTM* - 5.7cm<br />

Foto-TK* 6.8cm 25.1cm<br />

Foto-DTM* - 6.5cm<br />

Orto-Egen* 4.2cm -<br />

Orto-BLOM* 4.0cm -<br />

I forhold til spredningerne i egne ortofoto’s er den maksimale middelpunktsfejl valgt som spredning. Dette gøres med<br />

bevidsthed om, at spredningerne på punkterne i ortofotoet varierer afhængigt af, hvor man befinder sig i billedet.<br />

52


Dette indebærer, at når man er i nadir-punkter indgår kun spredningen fra absolut orientering. D<strong>og</strong> skal der vælges en<br />

spredning på punkterne, således at der kan beregnes en forventet spredning.<br />

Spredningen på ortofoto’s beregnes ud fra:<br />

= + (()) <br />

√<br />

() =<br />

+ ∙ <br />

2 ∙ <br />

√2<br />

53<br />

(Geoforum Danmark, 2005)<br />

Spredning σdhm er et udtryk for spredningen på højdemodellen. Desuden skal bidraget fra højdemodellen divideres<br />

med kvadratrod to for at lave spredningen om til den valgte definition fra indledningen. Da spredningen er forskellig<br />

på ortofotoet kan det defineres i et interval:<br />

For det ortofoto, som er produceret efter egen højdemodel, er intervallet:<br />

0.025 ≤ ≤ 0.042<br />

For det ortofoto, som er produceret efter BLOM’s højdemodel, er intervallet:<br />

0.025 ≤ ≤ 0.040<br />

Det vælges som sagt at vurdere ud fra den højeste middelpunkstfejl, som er henholdsvis 42 mm <strong>og</strong> 40 mm. Dette er<br />

de spredninger, der skal vurderes ud fra i forhold til kontrollen i Fase 4.


RTK-TK* KONTRA FOTO-TK*<br />

<strong>Ved</strong> sammenligning af de to tekniske kort, der er<br />

fremstillet ved hjælp af henholdsvis RTK-måling <strong>og</strong><br />

fot<strong>og</strong>rammetri, ses specielt en afvigelse i koten. Dette<br />

stemmer overens med tidligere beregninger af<br />

nøjagtigheden af det tekniske kort fremstillet på basis af<br />

fot<strong>og</strong>rammetri. <strong>Ved</strong> translation ses en stor forbedring<br />

af nøjagtigheden for koten, hvilket tyder på en<br />

systematisk fejl. Der kan spekuleres i, om det er den<br />

absolutte orientering, det er gået galt med, eller om det<br />

er digitaliseringen af det tekniske kort. D<strong>og</strong> kunne det<br />

tyde på, at det er manglende erfaring i digitaliseringen i<br />

stereo, da den absolutte orientering ser ud til at være<br />

acceptabel.<br />

I planen har der ved den direkte sammenligning vist sig<br />

overensstemmelse mellem den forventede <strong>og</strong><br />

beregnede spredning. D<strong>og</strong> er der ved translation blevet endnu lavere spredning. Residualplottet viser da <strong>og</strong>så n<strong>og</strong>et,<br />

der kunne tyde på en lille systematisk fejl i planen. Der er indgået 34 punkter i beregning, <strong>og</strong> der har ikke vist sig<br />

n<strong>og</strong>en grove fejl i n<strong>og</strong>en af sammenligningerne.<br />

Residalplot<br />

Direkte koordinatsammenligning Translation<br />

54<br />

Plan Kote<br />

Antal punkter 34 34<br />

Antal grove fejl 0 0


RTK-TK* KONTRA ORTO-EGET*<br />

I det tekniske kort fra RTK-indmålingen <strong>og</strong> ortofotoet, som<br />

er baseret på den fot<strong>og</strong>rammetriske højdemodel, er de<br />

plane koordinater blevet sammenlignet. Der ses en god<br />

sammenhæng ved den direkte sammenligning.<br />

Det ses endvidere som forventet, at spredningen er<br />

mindre end ved sammenligningen mellem det tekniske<br />

kort fra GPS-målingen <strong>og</strong> det tekniske kort, som er<br />

digitaliseret på baggrund af fot<strong>og</strong>rammetri.<br />

<strong>Ved</strong> translation ses en lille forbedring. På residualplottet<br />

med direkte koordinatsammenligning ses da <strong>og</strong>så n<strong>og</strong>et,<br />

der minder om en mindre systematisk fejl. Denne fejl<br />

kunne bero på fejl i den absolutte positionering.<br />

Residalplot<br />

Direkte koordinatsammenligning Translation<br />

55<br />

Plan<br />

Antal sammenligninger 33<br />

Antal grove fejl 0


RTK-TK* KONTRA ORTO-BLOM<br />

Forventningerne til ortofotoet baseret på BLOM’s<br />

højdedata var en smule større end til det ortofoto, der<br />

blev genereret ud fra egen højdemodel.<br />

Resultatet var da <strong>og</strong>så lidt bedre efter translation, men<br />

hvor den direkte sammenligning gav det samme resultat.<br />

I sammenligningen blev der sorteret to målinger fra, som lå<br />

over fejlgrænsen, hvilket betyder, at spredningen for<br />

ortofotoet, som er dannet på baggrund af egen<br />

højdemodel, er beregnet på baggrund af flere punkter.<br />

Dette kan godt have en minimal indflydelse på<br />

spredningerne, som ses i de to kortsammenligninger.<br />

Residalplot<br />

Direkte koordinatsammenligning Translation<br />

56<br />

Plan<br />

Antal sammenligninger 31<br />

Antal grove fejl 2


RTK-TK* KONTRA COWI-DDO<br />

<strong>Ved</strong> sammenligning af det tekniske kort med COWI’s<br />

ortofoto ses større systematisk fejl. COWI har da <strong>og</strong>så sat<br />

en stor fejlgrænse på deres ortofoto, men fejlen ligger<br />

stadigvæk udenfor denne grænse.<br />

Det tekniske kort er tidligere blevet sammenlignet med<br />

kontrolopmålte brønddæksler <strong>og</strong> riste under fase 1,<br />

hvorfra man kan deducere, at den store fejl må bero på<br />

COWI’s ortofoto. Der ses ved translation en spredning<br />

bedre end det forventede, hvilket kan tyde på, at<br />

ortofotoet er forskubbet. Dette kan udledes, da<br />

koordinaterne passer godt sammen, hvis koordinater fra<br />

det tekniske kort flyttes hen over ortofotoets.<br />

Der kan spekuleres i om aerotriangulationen måske ikke<br />

er foregået helt hensigtsmæssigt i forhold til den<br />

absolutte orientering af COWI’s ortofoto. Dette udledes af, at aerotriangulation er foregået over større områder,<br />

således at den passer godt til hele landet, mens der måske forekommer lokale udfald, hvorimod det tekniske kort er<br />

meget præcist bestemt i forhold til et større antal paspunkter i området.<br />

Residalplot<br />

Direkte koordinatsammenligning Translation<br />

57<br />

Plan<br />

Antal sammenligninger 27<br />

Antal grove fejl 0


COWI-DDO KONTRA ORTO-BLOM*<br />

<strong>Ved</strong> sammenligning med COWI’s ortofoto ses der igen en<br />

stor spredning. Denne gang sammenlignet med ortofoto<br />

baseret på BLOM-højder. Det underbygger tesen om, at<br />

der er en systematisk fejl. Det ses, at spredningen er større<br />

end, da COWI’s ortofoto blev sammenlignet med det<br />

tekniske kort fra RTK-opmåling.<br />

Den større unøjagtighed kan skyldes, at det er sværere at<br />

definere midtpunktet i orto-BLOM, end det var tidligere<br />

var med det tekniske kort, da der her kunne snappes til<br />

punkterne i GeoCAD.<br />

D<strong>og</strong> lader det til, at fejlen beror på COWI’s ortofoto, da<br />

det tekniske kort er kontrolleret op mod n<strong>og</strong>le<br />

kontrolmålte riste <strong>og</strong> brønddæksel under fase 1, hvor<br />

spredningen på disse i planen var under 1 cm. Endvidere er orto-BLOM kontrolleret med det tekniske kort, hvilket igen<br />

gav en lavere spredning end forventet. Derfor kan det konstateres, at da både COWI’s ortofoto <strong>og</strong> orto-BLOM begge<br />

er kontrolleret med det tekniske kort, <strong>og</strong> kun orto-BLOM passer overens med det tekniske kort, som er konstateret at<br />

være korrekt, hvorfor COWI’s kort må være mindre nøjagtigt.<br />

Derudover ses det, at translation har en bedre spredning end forventet.<br />

Der var en enkelt grov fejl, som blev frasorteret.<br />

Residalplot<br />

Direkte koordinatsammenligning Translation<br />

58<br />

Plan<br />

Antal sammenligninger 21<br />

Antal grove fejl 1


COWI-DDH KONTRA RTK-DTM*<br />

<strong>Ved</strong> sammenligning af de to højdemodeller; én fra COWI<br />

<strong>og</strong> én, der er genereret ud fra RTK-opmålingen, ses det, at<br />

spredningen er lidt bedre end den forventede. <strong>Ved</strong><br />

translation bliver den endnu bedre.<br />

På residualplottet er forskellene vist med farver. Turkis<br />

<strong>og</strong> gul er derfor under 1xspredning. Lys-rød <strong>og</strong> lys-blå er<br />

under 2xspredning, <strong>og</strong> blå <strong>og</strong> rød er under 3xspredning.<br />

De mørkeblå <strong>og</strong> mørkerøde felter ligger over<br />

3xspredning. Antallet af grove fejl passer godt med de<br />

0.3 %, som ligger over 3xspredning, da 3xspredningen<br />

gerne skulle svare til 99.7 % af alle målinger, som ikke er<br />

behæftet med grove fejl<br />

Højdesammenligning<br />

Direkte kotesammenligning Translation<br />

59<br />

Kote<br />

Antal sammenligninger 5312<br />

Antal grove fejl 19


FOTO-DTM* KONTRA COWI DDH<br />

<strong>Ved</strong> sammenligning af højdemodel genereret ved<br />

fot<strong>og</strong>rammetri <strong>og</strong> COWI’s højdemodel fås som<br />

udgangspunkt en lavere spredning end forventet. Dette<br />

ses både direkte <strong>og</strong> ved translation.<br />

D<strong>og</strong> ses det <strong>og</strong>så, at ca. 10 % af sammenligningerne er<br />

sorteret fra efter at være gået over fejlgrænsen. Dette kan<br />

dels begrundes med, at højdemodellen, som er genereret<br />

ved hjælp af Fot<strong>og</strong>rammetri, havde n<strong>og</strong>le usikkerheder -<br />

især hvor der ligger bygninger- <strong>og</strong> dels, at data, der er<br />

modtaget fra COWI, er uredigerede, <strong>og</strong> derfor <strong>og</strong>så kan<br />

indeholde fejl.<br />

På højdesammenligningerne ses de punkter, som er<br />

sorteret fra på grund af grovfejlsgrænsen som røde<br />

plamager.<br />

Højdesammenligning<br />

Direkte kotesammenligning Translation<br />

60<br />

Kote<br />

Antal sammenligninger 82363<br />

Antal grove fejl 8910


FOTO-TK* KONTRA FOT<br />

I sammenligning mellem teknisk kort fra fot<strong>og</strong>rammetri<br />

<strong>og</strong> FOT’s tekniske kort er der en klart bedre spredning i<br />

planen end forventet.<br />

Samtidig ses det, at der ikke sker n<strong>og</strong>en forbedring ved<br />

translation, tværtimod. Dette beror på, at der indgår to<br />

ubekendte i beregningen for spredningen. Dette<br />

indebærer en højere spredning, da nævneren i brøken<br />

får en højere værdi. Der er sorteret en grov måling fra.<br />

For koten er den bedre end forventet. Den forventede<br />

spredning for koten er da <strong>og</strong>så sat højt, hvilket beror på<br />

sammenligning mellem det tekniske kort fra fase 3 <strong>og</strong> de<br />

15 kontrolpunkter, hvor det er spredningen fra den<br />

direkte sammenligning, som danner udgangspunkt for<br />

den forventede spredning.<br />

Residalplot<br />

Direkte koordinatsammenligning Translation<br />

61<br />

Plan Kote<br />

Antal sammenligninger 30 30<br />

Antal grove fejl 1 0


ORTO-EGEN* KONTRA FOT<br />

Ortofotoet baseret på egen højdemodel sat op mod FOT<br />

viser en bedre spredning end forventet, hvilket er en<br />

indikation på, at den absolutte orientering er sket<br />

acceptabelt, da orienteringen af billederne passer godt<br />

sammen med det landsdækkende kortværk. <strong>Ved</strong><br />

translation sker der en lille forbedring.<br />

Dette passer overens med tidligere sammenligninger, der<br />

viser, at ortofotoet har en god nøjagtighed. Samtidig<br />

passer FOT <strong>og</strong>så ind med nærværende kort.<br />

Residalplot – Direkte sammenligning<br />

Translation<br />

62<br />

Plan<br />

Antal sammenligninger 57<br />

Antal grove fejl 0


ORTO-BLOM* KONTRA FOT<br />

Her sættes ortofoto baseret på BLOM-højder op mod FOT.<br />

Resultaterne for den direkte sammenligning er stort set de<br />

samme, som da FOT blev sat op mod ortofoto baseret på<br />

egne højder, hvilket igen er en indikation på, at den<br />

absolutte orientering af billederne er sket med en<br />

acceptabel nøjagtighed.<br />

D<strong>og</strong> har translationen givet en mærkbar forskel, hvilket<br />

kunne tyde på en systematisk fejl, hvilket residualplottet<br />

<strong>og</strong>så viser.<br />

Residalplot – Direkte sammenligning<br />

Translation<br />

63<br />

Plan<br />

Antal sammenligninger 71<br />

Antal grove fejl 2


RTK-TK* KONTRA FOT<br />

Sammenligning af RTK teknisk kort <strong>og</strong> FOT viser ved<br />

direkte sammenligning en spredning, der lægger sig tæt<br />

op ad den forventede. <strong>Ved</strong> translation bliver denne<br />

spredning endnu bedre. Det fremgår af residualplottet,<br />

at der er forekommet n<strong>og</strong>le mindre systematiske fejl,<br />

hvilket fremgår af, at pilene peger i samme retning.<br />

Residalplot<br />

Direkte koordinatsammenligning Translation<br />

64<br />

Plan Kote<br />

Antal sammenligninger 29 29<br />

Antal grove fejl 1 0


OPSAMLING<br />

I fase 4 blev forskellige kortprodukter vurderet. Her indgik det tekniske kort <strong>og</strong> højdemodel fra fase 1 samt teknisk<br />

kort, højdemodel <strong>og</strong> to ortofotos fra fase 3. Derudover indgik COWI DDO, COWI DDH <strong>og</strong> FOT’s tekniske kort.<br />

Overordnet set passede de forventede spredninger på forskellen meget godt med de beregnede. D<strong>og</strong> var der to<br />

kortsammenligninger, som gav en større spredning end forventet. Den ene kortsammenligning var det tekniske kort<br />

fra fase 3, hvor højden viste en større afvigelse end først antaget. Den anden var COWI’s ortofoto, som viste en<br />

systematisk fejl i planen, som tydede op, at ortofotoet generelt var skubbet i en nordlig retning.<br />

Førstnævnte fejl kan bero på dårlig placering af målemærket. Derudover kan der være et bidrag fra den absolutte<br />

orientering, hvis den ikke har været så nøjagtig. Med hensyn til COWI’s ortofoto går denne fejl igen i alle<br />

sammenligninger. Det kunne tyde på fejl i aerotriangulationen eller generel unøjagtighed. Dette underbygges af<br />

sammenligningerne med egne ortofotos, som er kontrolleret op mod det tekniske kort fra fase 1, der endvidere er<br />

kontrolleret ud fra 15 kontrolpunkter. Disse sammenligninger gav alle spredninger under det forventede.<br />

65


KONKLUSION<br />

Projektet omhandler forskellige tekniske opgaver, som vedrører landinspektørers arbejde. Dette gennemgås i fire<br />

faser, som omhandler forskellige tilgange til opmåling <strong>og</strong> afsætning. Fase 1 omhandler kortlægning ved RTK-måling <strong>og</strong><br />

resulterer i et teknisk kort <strong>og</strong> en digital højdemodel. Fase 2 omhandler afsætning ved hjælp af RTK <strong>og</strong> totalstation. I<br />

fase 3 er kortlægningen foretaget ved hjælp af fot<strong>og</strong>rammetri, <strong>og</strong> der bliver i denne fase <strong>og</strong>så genereret en<br />

højdemodel <strong>og</strong> to ortofotos. Den sidste fase er en opsamlingsfase, hvor forskellige kortlægningsmetoder<br />

sammenlignes i forhold til nøjagtigheden, herunder indgår andre eksterne kortværker såsom FOT, DDOland 2010 <strong>og</strong><br />

COWI’s DDH i analysen.<br />

I den første fase er hovedformålet at udarbejde et teknisk kort i 2 1 /2D <strong>og</strong> en digital højdemodel. Kravet til kortet er, at<br />

det skal kunne bruges til projektering af fremtidige tekniske anlæg. Som en indledning til opmålingen er de<br />

nærliggende GI-planfikspunkter kontrolleret for at undersøge, hvor godt de er bestemt i forhold til deres RTK-målte<br />

koordinater. Denne test viste en lille forskel på n<strong>og</strong>le cm i x <strong>og</strong> y <strong>og</strong> en lille skalering på 10ppm. Der kunne være krav<br />

fra rekvirenten om, at indmålingen skulle knyttes til GI-fikspunkterne. Derfor skulle det vurderes, om en<br />

transformation af detailmålte koordinater skulle foretages for at knytte detailopmålingen til det eksisterende<br />

fikspunktsnet, men i dette tilfælde blev det vurderet til ikke at komme på tale. Transformationen skulle være en 2D<br />

translation med skalering. Yderligere ændringer havde været unødvendige.<br />

Det tekniske kort blev til på baggrund af en kravspecifikation udarbejdet med udgangspunkt i FOT’s kravspecifikation.<br />

Grundet de ændrede vilkår for indmåling <strong>og</strong> den højere nøjagtighed ved RTK-måling blev der udarbejdet et tillæg til<br />

FOT-standarden. Heri blev beskrevet de øgede krav til den geometriske nøjagtighed <strong>og</strong> ændring i registrering. For<br />

eksempel skal bygninger registreres ved fri mur over sokkel <strong>og</strong> ikke ved tagudhæng som i FOT-standarden. Det<br />

tekniske kort blev efterfølgende undersøgt i forhold til n<strong>og</strong>le indmålte kontrolpunkter. Kortet havde en bedre<br />

nøjagtighed i planen end den forventede med en spredning på 9 mm, <strong>og</strong> i højden fremkom en ensidig fejl, hvilket<br />

betød, at spredningen blev større end den forventede med en beregnet spredning på 22 mm. Kontrolpunkterne var<br />

veldefinerede <strong>og</strong> er ikke udtryk for hele kortets spredning. <strong>Ved</strong> indmåling af bygninger ved bueskæring eller flugtmål<br />

må der forventes en højere spredning, hvilket bunder i en større unøjagtighed blandt andet på grund af<br />

definitionsvanskeligheder. Endvidere er der skøn forbundet med indmåling af buske, træer m.m. Den digitale<br />

terrænmodel blev <strong>og</strong>så kontrolleret i forhold til de indmålte kontrolpunkter. <strong>Ved</strong> interpolering af højderne kunne der<br />

beregnes en spredning på 5.7 cm før <strong>og</strong> 4.3 cm efter translation. Denne spredning er under den valgte i<br />

kravspecifikationen. Der kan skrues på parametrene i genereringen af tin-modellen, hvilket indebærer, at<br />

nøjagtigheden enten kan være dårligere eller bedre alt efter, hvilken kvalitet rekvirenten ønsker.<br />

Anden fase vedrører afsætning, <strong>og</strong> indeholder to hovedopgaver. Dels skal n<strong>og</strong>le udstykninger <strong>og</strong> centerlinjer på n<strong>og</strong>le<br />

veje afsættes, <strong>og</strong> dels skal en bygning afsættes med større præcision. Kravet til afsætning af skel <strong>og</strong> veje er et par cm,<br />

hvilket betyder, at det kan afsættes ved RTK-måling. Punkterne afmærkes med pæle. Forud for afsætningen er der<br />

sket en geometrisk konstruktion, så punkterne er lagt på plads i UTM32/EUREF89 med den korrekte<br />

afstandskorrektion, så punkterne kan afsættes ved RTK. Afsætningen er forløbet således, at der ikke er n<strong>og</strong>en af<br />

punkternes differencer til designkoordinaterne, som ligger udenfor tre gange den forventede spredning. Spredningen<br />

er stort set som forventet.<br />

<strong>Ved</strong> afsætning af bygningen blev der stillet større krav til nøjagtigheden. Kravspecifikationen lød på en nøjagtighed på<br />

et par millimeter relativt set. Afsætningen blev derfor foretaget med totalstation. Forud for afsætningen blev der<br />

foretaget en netmåling <strong>og</strong> netudjævning, som var med til at give et meget præcist net. Endvidere blev højder givet ved<br />

et geometrisk nivellement, hvor højderne blev udjævnet i et særligt pr<strong>og</strong>rammeret Matlab-script. Sidst skete en<br />

geometrisk konstruktion, hvor bygningen blev placeret indenfor fikspunkterne. <strong>Ved</strong> afsætningen blev nøjagtigheden<br />

overholdt med residualer, der lå indenfor fejlgrænsen. Samtidig gav kontrolmålingen ingen translation.<br />

66


I tredje fase skal der laves et teknisk kort baseret på fot<strong>og</strong>rammetri. Endvidere udarbejdes en digital højdemodel <strong>og</strong> to<br />

ortofotos. Det tekniske kort bygger igen på FOT-standarden, men med et tillæg der beskriver de ændrede vilkår. Det<br />

tekniske kort bygger på billeder med en højere opløsning end dem, der bruges til udarbejdelse af FOT’s tekniske kort.<br />

For at være i stand til at digitalisere det tekniske kort skal de udleverede fotos først orienteres relativt <strong>og</strong> absolut. Den<br />

relative orientering blev forsøgt både manuelt <strong>og</strong> automatisk. Det var muligt at skabe en orientering, som var under<br />

fejlgrænsen i begge tilfælde. D<strong>og</strong> blev den automatiske orientering valgt til den videre behandling. Til den absolutte<br />

orientering blev nøje udvalgte paspunkter RTK-indmålt. Paspunkterne var valgt til at give den bedste geometri i<br />

forhold til at bestemme optagecenterets position <strong>og</strong> flyets krængninger ω, ϕ <strong>og</strong> κ.<br />

Orienteringen gjorde det muligt at digitalisere det tekniske kort i stereo. Det tekniske kort blev kontrolleret i forhold til<br />

kontrolpunkter i projektområdet. For koten var der en ensidig fejl, som bevirkede en stor spredning. Efter translation<br />

var der stadigvæk en for stor fejl. Denne fejl skyldes sandsynligvis manglende erfaring med digitaliseringen. I planen<br />

var nøjagtigheden efter transformation under den forventede. Med de orienterede billeder var det desuden muligt at<br />

generere en højdemodel. Denne højdemodel bestående af en punktsky blev eksporteret til GeoCAD <strong>og</strong> en TIN-model<br />

blev beregnet med en tilstrækkelig pilhøjde, så det passede med kravspecifikationen. Sidst skulle der konstrueres to<br />

ortofotos. Det ene blev genereret med højder beregnet ud fra orienteringen. Den anden var fra BLOM’s højdemodel.<br />

Der sås en tydelig forbedring i ortofotoet med BLOM’s mere nøjagtige højdemodel.<br />

I den fjerde fase blev kortene fra de forskellige faser undersøgt. Her indgik det tekniske kort <strong>og</strong> højdemodellen fra fase<br />

1 samt teknisk kort, højdemodel <strong>og</strong> to ortofotos fra fase 3. Derudover indgik COWI DDO, COWI DDH <strong>og</strong> FOT’s tekniske<br />

kort. Overordnet set passede de forventede spredninger på forskellen godt overens med de beregnede. D<strong>og</strong> viste det<br />

tekniske kort fra fase 3 en større afvigelse end først antaget, hvor højden ikke passede så godt med de andre kort.<br />

Desuden viste COWI’s ortofoto en systematisk fejl i planen, hvor det lignede, at ortofotoet generelt var skubbet i en<br />

nordlig retning sammenlignet med de andre kort.<br />

67


KILDELISTE<br />

Jensen. K., 2005, Landmåling i Teori <strong>og</strong> Praksis, 2. udgave, <strong>Aalborg</strong> Universitet, <strong>Aalborg</strong><br />

Juhl. J., 2011a, Studievejledning 5.sem, <strong>Aalborg</strong> Universitet, <strong>Aalborg</strong><br />

Cederholm. P. & Jensen. K., 2006, GPS måling af utilgængelige detailpunkter, Landinspektøren 3, København<br />

FOTdanmark, 2010, Specifikation FOT 4.1, Version: FOT 4.1 20100328, FOTdanmark. Lokaliseret på<br />

http://www.fotdanmark.dk/Materiale/Files/FOTspecifikation/FOT4.1+alt_dk.pdf (09-09-2011)<br />

Juhl. J., 2011b, Forventet fot<strong>og</strong>rammetrisk nøjagtighed (spredning), <strong>Aalborg</strong> Universitet, <strong>Aalborg</strong><br />

Geoforum Danmark, 2005, Specifikation for ortofotos, 2.udgave, Geoforum Danmark, København. Lokaliseret på<br />

http://www.geoforum.dk/Files/Filer/Rapporter/Orto2004_2.udg_version_1.pdf (01-11-2011)<br />

COWI, 2008, DDOland 2008, COWI, Lyngby. Lokaliseret på ftp://ftp.gishotel.dk/pub/Konferencer_udstillinger/MIKon2008-PP/MI-Konf%20-%20DDOland2008_MapInfokonference_2008.pdf<br />

(10-11-2011)<br />

COWI, 2011, DDH - Danmarks Digitale Højdemodel, COWI, Lyngby. Lokaliseret på<br />

http://www.cowi.dk/menu/project/Ge<strong>og</strong>rafiskinformation<strong>og</strong>it/Kortlaegning<strong>og</strong>dataindsamling/Laserscanning/Pages/d<br />

dhdanmarksdigitalehoejdemodel.aspx (10-11-2011)<br />

68


BILAGSLISTE<br />

OVERSIGT OVER BILAG PÅ CD<br />

CD Nr. Filnavn på CD Beskrivelse<br />

Fase 1<br />

1 110920-Detail-samlet.txt Detailopmåling<br />

2 110920-Detail-samlet.l<strong>og</strong> L<strong>og</strong>fil fra konvertering af<br />

detailpunkter<br />

3 110920-Detail-Samlet.asc Original GeoCAD-asciifil<br />

med detailpunkter<br />

4 Teknisk kort Mappe med filer til teknisk<br />

kort<br />

5 DTM Mappe med filer til<br />

højdemodel<br />

6 Indmåling af GI-punkter Mappe med filer fra<br />

indmålingen af GI-punkter<br />

7 Kontrolpunkter Mappe med filer fra<br />

indmåling af<br />

kontrolpunkter<br />

Fase 2<br />

8 Afs.skel+vej.sot Original fil Afsætning<br />

skel+vej<br />

9 Bygningsafsætning.sot Original fil<br />

bygningsafsætning<br />

10 Skel_<strong>og</strong>_veje.bnr GeoCAD-binær<br />

designkoordinater<br />

11 Bygningsafsætning_design.asc GeoCAD-asciifil<br />

designkoordinater<br />

12 Afsætning af skel <strong>og</strong> veje Mappe med filer fra<br />

afsætning af skel <strong>og</strong> veje<br />

13 Bygningsafsætning Mappe med filer fra<br />

bygningsafsætning<br />

14 Testberegning Mappe med filer fra<br />

testberegning af<br />

bygningsafsætning<br />

15 Nedudjævning Mappe med filer fra<br />

netudjævning<br />

16 Nivellement Mappe med filer fra<br />

nivellement<br />

Fase 3<br />

17 Kalibreringsrapport.pdf Kalibreringsrapport for<br />

billede 84B <strong>og</strong> 85B<br />

18 Absolut_Relativ orientering Mappe med<br />

modelrapporter fra relativ<br />

<strong>og</strong> absolut orientering<br />

19 DTM Mappe med filer fra<br />

højdemodel<br />

20 Ortofotos Mappe med filer fra<br />

ortofotos<br />

21 Teknisk Kort Mappe med filer fra<br />

69


Fase 4<br />

teknisk kort<br />

22 Ydre Orientering Mappe med filer fra ydre<br />

orientering<br />

23 Spredninger.xlsx Excel-fil over forventede<br />

spredninger til fase 4<br />

24 COWIDDH-FotoDTM Mappe med filer fra<br />

sammenligning<br />

25 COWIDDH-RTKDTM Mappe med filer fra<br />

sammenligning<br />

26 COWIDDO-OrtoBLOM Mappe med filer fra<br />

sammenligning<br />

27 COWIDDO-RTKTK Mappe med filer fra<br />

sammenligning<br />

28 FotoTK-FOT Mappe med filer fra<br />

sammenligning<br />

29 OrtoBLOM-FOT Mappe med filer fra<br />

sammenligning<br />

30 OrtoBLOM-RTKTK Mappe med filer fra<br />

sammenligning<br />

31 OrtoEgen-FOT Mappe med filer fra<br />

sammenligning<br />

32 OrtoEgen-RTKTK Mappe med filer fra<br />

sammenligning<br />

33 RTKTK-FOT Mappe med filer fra<br />

sammenligning<br />

34 RTKTK-FotoTK Mappe med filer fra<br />

sammenligning<br />

70


BILAG – VURDERING AF RTK-MÅLING<br />

Spredningerne til RTK-målingerne er beregnet ud fra formler fra GPS-undervisningen:<br />

Spredningen på afvigelsen mellem to uafhængige målinger:<br />

,,∆ = ∑ , , ℎ ∆<br />

20<br />

Alle målinger antages at være uafhængige <strong>og</strong> målt med samme præcision. Derfor kan forplantningsloven anvendes, så<br />

det generelle udtryk ser ud som følger:<br />

Fejlgrænsen er dermed:<br />

Spredningen på en måling kan bestemmes ud fra:<br />

71<br />

<br />

,,∆ = √2 ∙ (,,)<br />

,,∆ = 3 ∙ √2 ∙ (,,)<br />

(,,) = ,,∆<br />

√2<br />

Leica Smartnet<br />

Kontrolpunkt Δe i m Δn i m Δh i m<br />

KTR. 1 0.010 -0.021 0.086<br />

KTR. 2 0.001 0.004 0.029<br />

KTR. 3 -0.035 0.045 0.129<br />

KTR. 4 -0.015 0.005 0.010<br />

KTR. 5 -0.001 0.002 0.031<br />

KTR. 6 -0.005 0.022 0.018<br />

KTR. 7 -0.008 0.015 0.021<br />

KTR. 8 0.003 0.016 -0.008<br />

KTR. 9 0.010 0.096 -0.215<br />

KTR. 10 -0.002 -0.008 0.003<br />

KTR. 11 -0.008 -0.005 -0.014<br />

KTR. 12 -0.008 -0.002 -0.025<br />

KTR. 13 -0.003 -0.004 -0.013<br />

KTR. 14 -0.010 -0.002 0.007<br />

KTR. 15 0.001 0.021 0.025<br />

KTR. 16 -0.007 0.015 -0.017<br />

KTR. 17 -0.007 0.021 0.021<br />

KTR. 18 -0.007 0.040 0.028<br />

KTR. 19 0.000 0.041 0.013<br />

KTR. 20 -0.002 0.017 -0.006


Spredningen for<br />

målinger Leica<br />

Smartnet<br />

(1. beregning)<br />

,,∆ (,,) ,,∆ E 0.010 0.007 0.031<br />

N 0.030 0.021 0.089<br />

H 0.062 0.044 0.185<br />

Derefter blev målinger sorteret fra <strong>og</strong> spredningerne blev beregnet igen. Dette blev gjort to gange. Følgende punkter<br />

blev sorteret fra KTR.1, KTR.3 <strong>og</strong> KTR. 9 Derefter blev følgende spredninger beregnet:<br />

Spredningen for<br />

målinger Leica<br />

Smartnet<br />

(3. beregning)<br />

,,∆ (,,) ,,∆ E 0.006 0.0017 0.017<br />

N 0.004 0.012 0.012<br />

H 0.018 0.051 0.052<br />

Derefter blev lignende beregninger foretaget for GPSnet.dk<br />

GPSnet.dk<br />

Kontrolpunkt Δe i m Δn i m Δh i m<br />

KTR. 1 0.025 0.032 -0.016<br />

KTR. 2 0.005 0.013 -0.017<br />

KTR. 3 0.014 0.025 -0.005<br />

KTR. 4 0.017 0.019 -0.027<br />

KTR. 5 -0.008 -0.009 -0.002<br />

KTR. 6 -0.004 -0.024 0.007<br />

KTR. 7 0.016 -0.020 -0.001<br />

KTR. 8 0.019 -0.034 -0.005<br />

KTR. 9 0.016 0.010 -0.015<br />

KTR. 10 0.006 0.002 -0.037<br />

KTR. 11 0.012 0.000 -0.008<br />

KTR. 12 0.012 -0.018 -0.017<br />

KTR. 13 0.014 0.015 0.023<br />

KTR. 14 0.021 0.002 0.009<br />

KTR. 15 0.019 -0.005 -0.011<br />

KTR. 16 0.018 -0.013 -0.031<br />

KTR. 17 0.029 0.007 0.010<br />

KTR. 18 0.013 -0.018 0.026<br />

KTR. 19 0.013 -0.025 -0.007<br />

KTR. 20 0.025 0.032 -0.016<br />

72


Spredningen for<br />

målinger GPSnet.dk<br />

,,∆ (,,) ,,∆<br />

E 0.016 0.012 0.049<br />

N 0.021 0.015 0.063<br />

H 0.019 0.013 0.056<br />

Til slut midles koordinaterne, som i beregningerne ikke er faldet udenfor grænsen 3 ∙ √2 ∙ (,,):<br />

Midlede<br />

koordinat<br />

er<br />

Leica<br />

Smartnet<br />

E N H GPSnet.dk E N H<br />

KTR. 1 712846.608 6181482.283 34.666<br />

KTR. 2 712849.654 6181494.631 34.800 712849.647 6181494.649 34.780<br />

KTR. 3 712856.917 6181511.981 34.799<br />

KTR. 4 712857.871 6181518.833 34.809 712857.866 6181518.850 34.801<br />

KTR. 5 712862.615 6181529.480 34.848 712862.608 6181529.484 34.839<br />

KTR. 6 712867.249 6181545.683 34.828 712867.250 6181545.681 34.817<br />

KTR. 7 712874.149 6181570.898 34.771 712874.148 6181570.881 34.758<br />

KTR. 8 712875.572 6181571.838 34.781 712875.558 6181571.849 34.785<br />

KTR. 9 712898.021 6181599.328 35.052<br />

KTR. 10 712887.275 6181623.290 34.653 712887.272 6181623.290 34.667<br />

KTR. 11 712882.696 6181631.179 34.774 712882.690 6181631.171 34.775<br />

KTR. 12 712879.873 6181616.257 34.764 712879.879 6181616.258 34.775<br />

KTR. 13 712869.217 6181582.974 34.684 712869.216 6181582.974 34.679<br />

KTR. 14 712864.373 6181567.883 34.802 712864.376 6181567.891 34.811<br />

KTR. 15 712859.511 6181552.359 34.796 712859.508 6181552.352 34.784<br />

KTR. 16 712854.542 6181536.980 34.796 712854.552 6181536.983 34.797<br />

KTR. 17 712849.743 6181521.883 34.778 712849.746 6181521.881 34.766<br />

KTR. 18 712848.706 6181519.017 34.786 712848.703 6181519.021 34.775<br />

KTR. 19 712844.715 6181506.197 34.818 712844.719 6181506.193 34.802<br />

KTR. 20 712840.050 6181491.488 34.836 712840.044 6181491.474 34.795<br />

73


BILAG-VERIFICERING AF INSTRUMENTER<br />

Totalstation<br />

Instrument: AAU 76887<br />

Dato: 15/09-11<br />

Job: 20110915-kontrol<br />

1) Udstyr til signalering af punkter<br />

OK<br />

2) Indstilling af parametre<br />

OK<br />

3) Elektronisk libelle<br />

OK<br />

4) Dåselibelle<br />

OK<br />

5) Laserlod<br />

Rykker sig 1-2 mm ved omdrejning. Vurderes til ikke at have større betydning for opgaven.<br />

6) Kollimationsfejl<br />

OK<br />

7) Horisontalakseskævhed<br />

: 288.5291<br />

: 88.5290<br />

288.5291 − (88.5290 + 200)<br />

= = 0.00005<br />

2<br />

= 61.5323<br />

= 20<br />

= 115<br />

= 115.15<br />

= arctan 0.15<br />

= 0.0075<br />

20<br />

sin (0.0075)<br />

= arctan <br />

= 0.0078<br />

tan (100 − 61.5323)<br />

Dette vurderes til at være for højt. Derfor foretages en anden kontrol.<br />

= 61.5336<br />

74


Dette er OK.<br />

8) Vertikalkredsens indeksfejl<br />

Dette er OK.<br />

= 20<br />

= 115<br />

= 115.10<br />

= arctan 0.10<br />

= 0.005<br />

20<br />

sin (0.005)<br />

= arctan <br />

= 0.001<br />

tan (100 − 61.3623)<br />

= 95.9316<br />

= 304.0703<br />

400 − 95.9316 − 304.0703<br />

= = −0.0001<br />

2<br />

9) Udstyr til bestemmelse af instrument <strong>og</strong> sigteskivehøjde:<br />

OK<br />

10) Kompensatorens funktionsområde<br />

= 0.07 (med uret)<br />

= −0.07 (mod uret)<br />

= 0.07 (med uret)<br />

= −0.07 (mod uret)<br />

11) Udstyr til bestemmelse af lufttryk <strong>og</strong> temperatur<br />

OK<br />

12) Distancemålingens enhed<br />

1. hoved<br />

2. hoved<br />

Middel<br />

stilling<br />

stilling<br />

V Sd V Sd V Sd<br />

Pkt. 2 102.744 22.437 297.256 22.437 102.744 22.437<br />

Pkt. 3 100.395 51.101 299.604 51.101 100.395 51.101<br />

Pkt. 4 99.587 77.102 300.414 77.102 99.586 77.102<br />

Pkt. 5 99.747 118.526 300.255 118.525 99.746 118.526<br />

Pkt. 6 99.818 146.095 300.181 146.095 99.818 146.095<br />

75


Nivelleringsinstrument<br />

Nivelleringsinstrumenter Leica Sprinter 150+ er blevet verificeret digitalt.<br />

Der blev rettet følgende fejl:<br />

Fejl -7.1”<br />

Højdedifference 0.9mm<br />

Endvidere blev stadie <strong>og</strong> libelle verificeret til at være ok.<br />

76


BILAG – VURDERING AF TRANSFORMATION AF GI-PUNKTER<br />

2D-TRANSLATION AF FIRE GI-PUNKTER<br />

*************************************************************************************<br />

2D TRANSFORMATION trans2d version 2005.10.27<br />

mattrans version 2007.08.14<br />

2 translationer 2011-09-23 14:17:13<br />

*************************************************************************************<br />

FÆLLESPUNKTER<br />

-------------------------------------------------------------------------------------<br />

Fællespunkter i FRA system: (c:\users\andreas\documents\matlab\gps\fra2.txt)<br />

x [m] y [m]<br />

1 709153.678 6179837.405<br />

2 711989.400 6180644.456<br />

3 712879.269 6178459.575<br />

4 714902.407 6178706.189<br />

Fællespunkter i TIL system: (c:\users\andreas\documents\matlab\gps\til2.txt)<br />

x [m] y [m]<br />

1 709153.647 6179837.402<br />

2 711989.401 6180644.459<br />

3 712879.280 6178459.564<br />

4 714902.506 6178706.280<br />

TRANSFORMATIONSPARAMETRE<br />

-------------------------------------------------------------------------------------<br />

tx : 0.020 m<br />

ty : 0.020 m<br />

sigma 0 : 0.052 m<br />

TRANSFORMATIONSFORMLER<br />

-------------------------------------------------------------------------------------<br />

|xTil| |xFra| |tx|<br />

| | = | | + | |<br />

|yTil| |yFra| |ty|<br />

TRANSFORMEREDE FÆLLESPUNKTER (TIL system)<br />

-------------------------------------------------------------------------------------<br />

x [m] Res [m] Spred [m] y [m] Res [m] Spred [m]<br />

1 709153.698 0.051 0.037 6179837.425 0.023 0.037<br />

2 711989.420 0.019 0.037 6180644.476 0.017 0.037<br />

3 712879.289 0.009 0.037 6178459.595 0.031 0.037<br />

4 714902.427 -0.079 0.037 6178706.209 -0.071 0.037<br />

77


2D-TRANSFORMATION (HELMERT) AF FIRE GI-PUNKTER<br />

*************************************************************************************<br />

2D TRANSFORMATION trans2d version 2005.10.27<br />

mattrans version 2007.08.14<br />

2 translationer. 1 rotation. 1 skala (Helmert) 2011-09-23 13:29:31<br />

*************************************************************************************<br />

FÆLLESPUNKTER<br />

-------------------------------------------------------------------------------------<br />

Fællespunkter i FRA system: (c:\users\andreas\documents\matlab\gps\fra2.txt)<br />

x [m] y [m]<br />

1 709153.678 6179837.405<br />

2 711989.400 6180644.456<br />

3 712879.269 6178459.575<br />

4 714902.407 6178706.189<br />

78


Fællespunkter i TIL system: (c:\users\andreas\documents\matlab\gps\til2.txt)<br />

x [m] y [m]<br />

1 709153.647 6179837.402<br />

2 711989.401 6180644.459<br />

3 712879.280 6178459.564<br />

4 714902.506 6178706.280<br />

TRANSFORMATIONSPARAMETRE<br />

-------------------------------------------------------------------------------------<br />

a : 1.0000156223<br />

b : 0.0000166794<br />

tx : 91.962 m<br />

ty : -108.396 m<br />

rot : 0.001 gon (afledt af a.b)<br />

skala : 1.000016 (afledt af a.b)<br />

sigma 0 : 0.037 m<br />

TRANSFORMATIONSFORMLER<br />

-------------------------------------------------------------------------------------<br />

|xTil| |a -b| |xFra| |tx|<br />

| | = | | * | | + | |<br />

|yTil| |b a| |yFra| |ty|<br />

TRANSFORMEREDE FÆLLESPUNKTER (TIL system)<br />

-------------------------------------------------------------------------------------<br />

x [m] Res [m] Spred [m] y [m] Res [m] Spred [m]<br />

1 709153.643 -0.004 0.032 6179837.380 -0.022 0.032<br />

2 711989.396 -0.005 0.021 6180644.491 0.032 0.021<br />

3 712879.315 0.035 0.021 6178459.591 0.027 0.021<br />

4 714902.481 -0.025 0.030 6178706.243 -0.037 0.030<br />

79


2D-TRANSLATION AF TRE GI-PUNKTER<br />

*************************************************************************************<br />

2D TRANSFORMATION trans2d version 2005.10.27<br />

mattrans version 2007.08.14<br />

2 translationer 2011-09-23 14:20:55<br />

*************************************************************************************<br />

FÆLLESPUNKTER<br />

-------------------------------------------------------------------------------------<br />

Fællespunkter i FRA system: (c:\users\andreas\documents\matlab\gps\fra2-870.txt)<br />

x [m] y [m]<br />

1 709153.678 6179837.405<br />

2 711989.400 6180644.456<br />

3 712879.269 6178459.575<br />

Fællespunkter i TIL system: (c:\users\andreas\documents\matlab\gps\til2-870.txt)<br />

x [m] y [m]<br />

1 709153.647 6179837.402<br />

2 711989.401 6180644.459<br />

3 712879.280 6178459.564<br />

TRANSFORMATIONSPARAMETRE<br />

-------------------------------------------------------------------------------------<br />

tx : -0.006 m<br />

ty : -0.004 m<br />

sigma 0 : 0.016 m<br />

TRANSFORMATIONSFORMLER<br />

-------------------------------------------------------------------------------------<br />

|xTil| |xFra| |tx|<br />

| | = | | + | |<br />

|yTil| |yFra| |ty|<br />

80


TRANSFORMEREDE FÆLLESPUNKTER (TIL system)<br />

-------------------------------------------------------------------------------------<br />

x [m] Res [m] Spred [m] y [m] Res [m] Spred [m]<br />

1 709153.672 0.025 0.013 6179837.401 -0.001 0.013<br />

2 711989.394 -0.007 0.013 6180644.452 -0.007 0.013<br />

3 712879.263 -0.017 0.013 6178459.571 0.007 0.013<br />

2D-TRANSFORMATION (HELMERT) AF TRE GI-PUNKTER<br />

*************************************************************************************<br />

2D TRANSFORMATION trans2d version 2005.10.27<br />

mattrans version 2007.08.14<br />

2 translationer. 1 rotation. 1 skala (Helmert) 2011-09-23 14:22:53<br />

*************************************************************************************<br />

FÆLLESPUNKTER<br />

-------------------------------------------------------------------------------------<br />

81


Fællespunkter i FRA system: (c:\users\andreas\documents\matlab\gps\fra2-870.txt)<br />

x [m] y [m]<br />

1 709153.678 6179837.405<br />

2 711989.400 6180644.456<br />

3 712879.269 6178459.575<br />

Fællespunkter i TIL system: (c:\users\andreas\documents\matlab\gps\til2-870.txt)<br />

x [m] y [m]<br />

1 709153.647 6179837.402<br />

2 711989.401 6180644.459<br />

3 712879.280 6178459.564<br />

TRANSFORMATIONSPARAMETRE<br />

-------------------------------------------------------------------------------------<br />

a : 1.0000100734<br />

b : 0.0000009536<br />

tx : -1.279 m<br />

ty : -62.932 m<br />

rot : 0.000 gon (afledt af a.b)<br />

skala : 1.000010 (afledt af a.b)<br />

sigma 0 : 0.004 m<br />

TRANSFORMATIONSFORMLER<br />

-------------------------------------------------------------------------------------<br />

|xTil| |a -b| |xFra| |tx|<br />

| | = | | * | | + | |<br />

|yTil| |b a| |yFra| |ty|<br />

TRANSFORMEREDE FÆLLESPUNKTER (TIL system)<br />

-------------------------------------------------------------------------------------<br />

x [m] Res [m] Spred [m] y [m] Res [m] Spred [m]<br />

1 709153.649 0.002 0.004 6179837.401 -0.001 0.004<br />

2 711989.399 -0.002 0.003 6180644.463 0.004 0.003<br />

3 712879.279 -0.001 0.004 6178459.561 -0.003 0.004<br />

82


2D-TRANSLATION MED SKALERING AF TRE GI-PUNKTER<br />

*************************************************************************************<br />

2D TRANSFORMATION trans2d version 2005.10.27<br />

mattrans version 2007.08.14<br />

2 translationer. 1 skala 2011-09-23 14:25:30<br />

*************************************************************************************<br />

FÆLLESPUNKTER<br />

-------------------------------------------------------------------------------------<br />

Fællespunkter i FRA system: (c:\users\andreas\documents\matlab\gps\fra2-870.txt)<br />

x [m] y [m]<br />

1 709153.678 6179837.405<br />

2 711989.400 6180644.456<br />

3 712879.269 6178459.575<br />

Fællespunkter i TIL system: (c:\users\andreas\documents\matlab\gps\til2-870.txt)<br />

x [m] y [m]<br />

1 709153.647 6179837.402<br />

2 711989.401 6180644.459<br />

3 712879.280 6178459.564<br />

TRANSFORMATIONSPARAMETRE<br />

-------------------------------------------------------------------------------------<br />

tx : -7.172 m<br />

ty : -62.253 m<br />

skala : 1.000010<br />

sigma 0 : 0.004 m<br />

TRANSFORMATIONSFORMLER<br />

-------------------------------------------------------------------------------------<br />

|xTil| |xFra| |tx|<br />

| | = skala * | | + | |<br />

|yTil| |yFra| |ty|<br />

83


TRANSFORMEREDE FÆLLESPUNKTER (TIL system)<br />

-------------------------------------------------------------------------------------<br />

x [m] Res [m] Spred [m] y [m] Res [m] Spred [m]<br />

1 709153.650 0.003 0.004 6179837.403 0.001 0.002<br />

2 711989.400 -0.001 0.002 6180644.462 0.003 0.003<br />

3 712879.278 -0.002 0.003 6178459.559 -0.005 0.003<br />

84


BILAG – VURDERING AF KONTROLMÅLING AF BYGNINGSDIMENSIONER<br />

Adresse Opmålt i marken Opmålt i GeoCAD Afvigelse<br />

<strong>Askeengen</strong> 48 12.99 12.805 0.185<br />

4.075 4.074 0.001<br />

1.395 1.544 -0.149<br />

0.571 0.571 0<br />

4.415 4.405 0.01<br />

9.41 9.351 0.059<br />

4.43 4.403 0.027<br />

5.282 5.352 -0.07<br />

1.395 1.406 -0.011<br />

0.705 0.705 0<br />

12.985 12.915 0.07<br />

8.345 8.316 0.029<br />

Udestue 4.22 4.212 0.008<br />

4.11 4.066 0.044<br />

4.21 4.189 0.021<br />

Askengen 50 12.958 12.959 -0.001<br />

8.35 8.397 -0.047<br />

12.969 13.064 -0.095<br />

8.325 8.449 -0.124<br />

Udestue 3.8 3.74 0.06<br />

7.49 7.434 0.056<br />

3.78 3.677 0.103<br />

<strong>Askeengen</strong> 52 11.91 11.91 0<br />

14.465 14.468 -0.003<br />

14.36 14.203 0.157<br />

8.94 8.947 -0.007<br />

Carport 8.28 8.243 0.037<br />

Udestue 5.88 5.86 0.02<br />

2.935 2.917 0.018<br />

<strong>Askeengen</strong> 54 11.885 11.735 0.15<br />

14.52 14.387 0.133<br />

8.96 8.984 -0.024<br />

8.355 8.353 0.002<br />

2.9 2.886 0.014<br />

6.115 6.034 0.081<br />

Carport 6 6.179 -0.179<br />

3.2 3.21 -0.01<br />

<strong>Askeengen</strong> 56 14.48 14.504 -0.024<br />

8.93 8.918 0.012<br />

9.22 9.208 0.012<br />

4.87 4.872 -0.002<br />

5.28 5.236 0.044<br />

13.8 13.781 0.019<br />

85


Alle de målinger, som er markeret med rød, overstiger den maksimale afvigelse vurderet ud fra apriorispredningen på<br />

afvigelsen mellem målingen i marken <strong>og</strong> målingen GeoCAD = (√2 ∙ 0.022) + 0.01 = 0.033. Derfor sorteres<br />

alle målinger, som ligger over = ±3 ∙ √2 ∙ = 0.140 fra.<br />

Derefter beregnes spredningen på afvigelserne mellem målingerne i GeoCAD <strong>og</strong> målingerne i marken til at være<br />

= ∑ <br />

<br />

<br />

= 0.049.<br />

Tilfældige fejl betragtes som betydningsløse, da alle mål med stålmålebånd er foretaget med omhu med længder<br />

under 50 m (Jensen, 2005).<br />

86


BILAG – AFSÆTNING AF SKEL OG VEJE<br />

Punkt Δe i mm Δn i mm<br />

122 -2 -13<br />

140 -2 -9<br />

179 19 -11<br />

129 -18 -23<br />

128 13 -13<br />

178 7 -29<br />

138 5 -4<br />

127 4 -33<br />

137 0 -16<br />

177 7 -6<br />

126 7 10<br />

136 -14 -6<br />

125 -7 -17<br />

135 -9 1<br />

176 9 -8<br />

100 -3 -2<br />

175 -21 -9<br />

102 -8 19<br />

174 4 -10<br />

103 1 11<br />

173 4 16<br />

104 -11 14<br />

153 3 -7<br />

165 7 13<br />

172 -22 2<br />

101 -3 26<br />

171 8 7<br />

106 -7 4<br />

105 -20 27<br />

120 8 10<br />

166 -3 -13<br />

167 -16 25<br />

152 11 2<br />

158 -13 -15<br />

116 -11 23<br />

147 5 0<br />

168 -23 12<br />

169 -9 -7<br />

146 -4 -14<br />

170 13 -27<br />

119 -6 -4<br />

118 -19 0<br />

111 6 -1<br />

110 0 16<br />

87


109 6 -2<br />

108 6 -4<br />

107 17 19<br />

124 -8 3<br />

123 -10 -3<br />

Derefter beregnes spredning på E <strong>og</strong> spredningen på N ud fra formlen = = ∑ <br />

.<br />

Spredning E N<br />

0.011 0.015<br />

Dette giver en spredning på = <br />

<br />

= . . <br />

<br />

88<br />

= 0.013.


BILAG – DOKUMENTATIONSFIL MINIMUM FASTHOLDT UDJÆVNING<br />

Matlab-script:minimum_fastholdt_udjaevningsscript.m<br />

L5-projekt - Kortlægning <strong>og</strong> afsætning.<br />

Gr. 1 - Andreas Kudsk Jensen & Johan V. Eckhoff.<br />

Udjævning af nivellement, hvor 1 højdefikspunkt; 2-01-09011 er fastholdt.<br />

Der er foretaget dobbeltnivellementer til følgende hjælpepunkter; 100, 190, 300 <strong>og</strong> 390,<br />

Parametre vedr. nivellement:<br />

Kilometerspredning a priori (sigma_k)= 0.00400 m/km.<br />

Spredninger på enkeltnivellementer er beregnet jf. Jensen 2005.<br />

Punkt Udjævnede højder<br />

Pnr. m.<br />

9011 25.672<br />

9130 26.690<br />

9136 22.218<br />

9216 28.329<br />

300 25.426<br />

190 23.696<br />

390 26.611<br />

100 24.452<br />

Residualvektor:<br />

Fra Til<br />

9011 300 0.001 m<br />

300 9011 0.001 m<br />

9130 300 0.005 m<br />

300 9130 0.006 m<br />

9136 9216 0.003 m<br />

9216 9136 0.002 m<br />

9216 190 0.002 m<br />

190 9216 0.002 m<br />

190 390 -0.001 m<br />

390 190 -0.000 m<br />

390 300 -0.001 m<br />

300 390 -0.000 m<br />

300 100 0.000 m<br />

100 300 0.001 m<br />

9011 NaN 0.000 m<br />

Spredningen på vægtenheden:<br />

1.0971<br />

Antal af ligninger <strong>og</strong> ubekendte:<br />

15<br />

8<br />

Antal overbestemmelser:<br />

7<br />

Spredninger i kovarians matricen:<br />

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000<br />

0.0000 0.0048 0.0042 0.0030 0.0030 0.0035 0.0033 0.0032<br />

0.0000 0.0042 0.0042 0.0030 0.0030 0.0035 0.0033 0.0032<br />

0.0000 0.0030 0.0030 0.0041 0.0030 0.0030 0.0030 0.0030<br />

0.0000 0.0030 0.0030 0.0030 0.0030 0.0030 0.0030 0.0030<br />

0.0000 0.0035 0.0035 0.0030 0.0030 0.0035 0.0033 0.0032<br />

0.0000 0.0033 0.0033 0.0030 0.0030 0.0033 0.0033 0.0032<br />

0.0000 0.0032 0.0032 0.0030 0.0030 0.0032 0.0032 0.0032<br />

89


Normaliserede residualer:<br />

Fra Til<br />

9011 300 0.344<br />

300 9011 0.344<br />

9130 300 2.128<br />

300 9130 2.128<br />

9136 9216 1.190<br />

9216 9136 1.190<br />

9216 190 0.861<br />

190 9216 0.861<br />

190 390 -0.491<br />

390 190 -0.491<br />

390 300 -0.530<br />

300 390 -0.530<br />

300 100 0.593<br />

100 300 0.593<br />

9011 NaN 0.000<br />

90


BILAG – DOKUMENTATIONSFIL FASTHOLDT UDJÆVNING<br />

Matlab-script:fastholdt_udjaevningsscript.m<br />

L5-projekt - Kortlægning <strong>og</strong> afsætning.<br />

Gr. 1 - Andreas Kudsk Jensen & Johan V. Eckhoff.<br />

Udjævning af nivellement fra 4 højde fikspunkter; 2-01-09011, 2-01-09130, 2-01-09136 <strong>og</strong> 2-01-09216.<br />

Der er foretaget dobbeltnivellementer til følgende hjælpepunkter; 100, 190, 300 <strong>og</strong> 390.<br />

Parametre vedr. nivellement:<br />

Kilometerspredning (sigma_k)= 0.00400 m/km.<br />

Spredninger på enkeltnivellementer er beregnet jf. Jensen 2005.<br />

Spredninger på højdefikspunkter er bestemt til 0.003 m. Derudover er spredningen på højdefikspunkt<br />

2-01-09130 bestemt til 1 m, da det under beregninger blev fundet til at være dårligt<br />

Punkt Udjævnede højder<br />

Pnr. m.<br />

9011 25.674<br />

9130 26.696<br />

9136 22.225<br />

9216 28.334<br />

300 25.430<br />

190 23.701<br />

390 26.616<br />

100 24.456<br />

Residualvektor:<br />

Fra Til<br />

9011 300 0.003 m<br />

300 9011 -0.001 m<br />

9130 300 0.006 m<br />

300 9130 0.005 m<br />

9136 9216 0.003 m<br />

9216 9136 0.002 m<br />

9216 190 0.001 m<br />

190 9216 0.003 m<br />

190 390 -0.001 m<br />

390 190 -0.000 m<br />

390 300 -0.001 m<br />

300 390 -0.000 m<br />

300 100 0.000 m<br />

100 300 0.001 m<br />

9011 NaN 0.002 m<br />

9130 NaN 0.000 m<br />

9136 NaN -0.002 m<br />

9216 NaN 0.017 m<br />

Spredningen på vægtenheden:<br />

1.0334<br />

Antal af ligninger <strong>og</strong> ubekendte:<br />

18<br />

8<br />

Antal overbestemmelser:<br />

10<br />

Spredninger i kovarians matricen:<br />

0.0026 0.0011 0.0013 0.0020 0.0020 0.0018 0.0019 0.0020<br />

0.0011 0.0023 0.0018 0.0015 0.0015 0.0016 0.0016 0.0015<br />

0.0013 0.0018 0.0022 0.0018 0.0018 0.0019 0.0019 0.0018<br />

0.0020 0.0015 0.0018 0.0038 0.0028 0.0025 0.0026 0.0027<br />

0.0020 0.0015 0.0018 0.0028 0.0028 0.0025 0.0026 0.0027<br />

91


0.0018 0.0016 0.0019 0.0025 0.0025 0.0027 0.0026 0.0025<br />

0.0019 0.0016 0.0019 0.0026 0.0026 0.0026 0.0027 0.0027<br />

0.0020 0.0015 0.0018 0.0027 0.0027 0.0025 0.0027 0.0027<br />

Normaliserede residualer:<br />

Fra Til<br />

9011 300 0.949 m<br />

300 9011 -0.311 m<br />

9130 300 2.193 m<br />

300 9130 2.193 m<br />

9136 9216 1.130 m<br />

9216 9136 1.104 m<br />

9216 190 0.344 m<br />

190 9216 1.371 m<br />

190 390 -0.725 m<br />

390 190 -0.235 m<br />

390 300 -0.764 m<br />

300 390 -0.310 m<br />

300 100 0.396 m<br />

100 300 0.803 m<br />

9011 NaN 1.382 m<br />

9130 NaN 0.029 m<br />

9136 NaN -1.094 m<br />

9216 NaN 0.017 m<br />

92


BILAG – BINGOFIL<br />

# **********************************************************************<br />

# AEROoffice V5.1e 2008-07-24-Alpha<br />

# Copyright by IGI mbH, 1996-2008<br />

#<br />

# Dongle-ID: 07D5-00CF-0209-A551<br />

# Owner: COWI A/S<br />

#<br />

# Datamanager Outputfile<br />

# 23-04-2009 09:58:21<br />

#<br />

# **********************************************************************<br />

#<br />

# Project: 090407_1_CKP<br />

# Projectfile: L:\Fly-data\2009\AeroControl\DK\090407_1_CKP\090407_1_CKP.aop<br />

#<br />

# Event Marks: L:\Fly-data\2009\AeroControl\DK\090407_1_CKP\work\090407_1_CKP_XS1.aom<br />

#<br />

# Format Type: GIP Bingo<br />

#<br />

# Sensor-Leverarm: 0.161m 0.000m 0.426m ()<br />

#<br />

# Boresight Alignment:<br />

# Boresight from Project: 090318_2_CKP<br />

# Roll : 0.3298°<br />

# Pitch : -0.1842°<br />

# Yaw : -0.1100°<br />

#<br />

# Meridian Convergence corrected<br />

#<br />

# Local Coordinate System:<br />

# UTM- WGS84 - DVR90<br />

# Selected Zone: 32N<br />

#<br />

# **********************************************************************<br />

#<br />

# Infos from the postprocessing l<strong>og</strong>file:<br />

# AEROoffice V5.1e 2008-07-24-Alpha .<br />

# Dongle-ID: 07D5-00CF-0209-A551 .<br />

# Owner: COWI A/S .<br />

# 08-04-2009 13:09:29 .<br />

# Header of imported GPS File .<br />

# Project: 090407_1_CKP .<br />

# Pr<strong>og</strong>ram: GrafNav Version 8.10.2110 .<br />

# Profile: IGI AEROCTRL .<br />

# Source: GPS Epochs(Combined) .<br />

# ProcessInfo: Run (2) by MBE on 04/08/2009 at 13:07:15 .<br />

# Datum: WGS84, (processing datum) .<br />

# Master 1: Name VirA, Status ENABLED .<br />

# GPS-Leverarm: 1.135m -0.037m -1.184m .<br />

# GPS-Leverarm: variable .<br />

# Mount-Center-Leverarm: 0.115m 0.000m 0.302m .<br />

#<br />

# **********************************************************************<br />

#<br />

#<br />

# Selected Units:<br />

# Angular Units: Grad/Gon (0..400)<br />

# Length Units: Meter<br />

#<br />

# Format:<br />

# ID Time Easting Northing Height Phi Omega Kappa dPos dPhi/Omega dKappa<br />

#<br />

93


# Output of event data<br />

# File will contain 59 online Events<br />

#<br />

# LINE 004 111000<br />

35 206428.121 711970.332 6178371.203 483.716 0.4472 0.0440 64.5088 0.03 0.004 0.006<br />

36 206430.025 712041.166 6178486.866 484.570 0.4445 0.0414 64.3534 0.03 0.004 0.006<br />

37 206431.820 712108.056 6178595.643 485.346 0.4412 0.0962 64.1895 0.03 0.004 0.006<br />

38 206433.726 712178.771 6178710.936 485.372 0.4996 0.1071 64.2453 0.03 0.004 0.006<br />

39 206435.616 712249.199 6178824.643 484.405 0.4023 -0.0094 64.4188 0.03 0.004 0.006<br />

40 206437.520 712320.448 6178938.697 482.156 0.4516 0.1031 64.3967 0.03 0.004 0.006<br />

41 206439.426 712391.873 6179052.584 480.193 0.4354 0.0836 64.3401 0.03 0.004 0.006<br />

42 206441.378 712465.227 6179168.758 479.179 0.4539 0.0995 64.2679 0.03 0.004 0.006<br />

43 206443.221 712534.802 6179277.669 478.641 0.4107 0.0885 64.4261 0.03 0.004 0.006<br />

44 206445.126 712606.938 6179389.346 479.511 0.4413 0.0383 64.2004 0.03 0.004 0.006<br />

45 206447.015 712678.470 6179499.282 481.276 0.4534 0.0780 64.3209 0.03 0.004 0.006<br />

46 206448.919 712750.872 6179609.322 482.951 0.4303 0.0216 64.2379 0.03 0.004 0.006<br />

47 206450.825 712823.520 6179718.997 484.606 0.4057 0.0115 64.3404 0.03 0.004 0.006<br />

48 206452.825 712899.831 6179834.129 485.788 0.4415 0.0673 64.4424 0.03 0.004 0.006<br />

49 206454.716 712971.864 6179943.377 485.793 0.4036 0.0370 64.3591 0.03 0.004 0.006<br />

50 206456.619 713044.184 6180053.853 486.431 0.4068 0.0406 64.2308 0.03 0.004 0.006<br />

51 206458.523 713116.519 6180164.115 487.598 0.4226 0.0120 64.0717 0.03 0.004 0.006<br />

52 206460.524 713192.580 6180279.271 488.985 0.4454 0.0031 64.3457 0.03 0.004 0.006<br />

53 206462.418 713264.653 6180387.762 489.761 0.4198 -0.0083 64.6628 0.03 0.004 0.006<br />

54 206464.378 713338.814 6180500.285 488.719 0.4572 -0.0185 64.4943 0.03 0.004 0.006<br />

# LINE 003 111000<br />

55 206826.141 711113.576 6181734.951 484.812 0.0743 0.5183 -20.9386 0.03 0.004 0.006<br />

56 206828.544 711236.276 6181692.261 484.838 0.0089 0.5451 -21.1146 0.03 0.004 0.006<br />

57 206831.043 711363.382 6181646.705 484.153 0.0615 0.3634 -21.1595 0.03 0.004 0.006<br />

58 206833.541 711490.628 6181600.848 484.600 0.0518 0.4889 -20.8163 0.03 0.004 0.006<br />

59 206835.946 711614.129 6181557.188 481.742 0.0015 0.5297 -21.2623 0.03 0.004 0.006<br />

60 206838.444 711741.927 6181511.579 479.877 0.0557 0.4155 -21.1105 0.03 0.004 0.006<br />

61 206840.834 711862.945 6181467.527 481.545 0.0714 0.5027 -21.0267 0.03 0.004 0.006<br />

62 206843.364 711990.607 6181421.354 482.496 0.0427 0.5036 -20.8528 0.03 0.004 0.006<br />

63 206845.831 712114.658 6181376.246 483.082 0.0836 0.3919 -20.9486 0.03 0.004 0.006<br />

64 206848.361 712241.583 6181329.724 485.350 0.0318 0.5522 -21.1503 0.03 0.004 0.006<br />

65 206850.845 712366.371 6181284.793 486.608 0.0651 0.4231 -21.1834 0.03 0.004 0.006<br />

66 206853.358 712492.535 6181238.874 486.445 0.1026 0.3422 -20.9039 0.03 0.004 0.006<br />

67 206855.841 712617.431 6181193.771 484.202 0.0220 0.5094 -20.7357 0.03 0.004 0.006<br />

68 206858.230 712737.967 6181151.185 482.276 0.0292 0.5117 -21.0231 0.03 0.004 0.006<br />

69 206860.729 712864.342 6181107.824 483.334 0.0617 0.3951 -21.0636 0.03 0.004 0.006<br />

70 206863.227 712990.356 6181064.451 483.437 0.0344 0.3943 -21.0820 0.03 0.004 0.006<br />

71 206865.727 713116.259 6181021.238 482.293 -0.0124 0.4471 -20.9671 0.03 0.004 0.006<br />

72 206868.242 713242.174 6180977.845 481.998 0.0467 0.4708 -20.9865 0.03 0.004 0.006<br />

73 206870.831 713371.476 6180934.304 482.222 0.0108 0.4576 -21.0297 0.03 0.004 0.006<br />

# LINE 004 111000<br />

74 207309.438 713351.076 6180525.861 486.070 -0.4819 -0.0525 -135.6591 0.03 0.004 0.010<br />

75 207311.750 713280.485 6180414.023 484.119 -0.4445 -0.0786 -135.6471 0.03 0.004 0.010<br />

76 207314.045 713211.073 6180302.723 483.475 -0.4559 -0.0896 -135.7161 0.03 0.004 0.010<br />

77 207316.433 713139.537 6180186.717 483.705 -0.4287 -0.0692 -135.8042 0.03 0.004 0.010<br />

78 207318.746 713070.269 6180074.579 483.589 -0.3657 -0.0218 -135.6524 0.03 0.004 0.010<br />

79 207321.041 713002.047 6179963.446 481.551 -0.4937 -0.1234 -135.5693 0.03 0.004 0.010<br />

80 207323.447 712930.334 6179847.055 481.151 -0.4391 -0.0074 -135.4533 0.03 0.004 0.010<br />

81 207325.742 712861.769 6179736.144 481.645 -0.5378 -0.0825 -135.6186 0.03 0.004 0.010<br />

82 207328.130 712789.803 6179620.909 484.244 -0.4305 -0.0363 -135.9052 0.03 0.004 0.010<br />

83 207330.441 712719.560 6179509.712 485.999 -0.4047 -0.0391 -135.7710 0.03 0.004 0.010<br />

84 207332.737 712649.258 6179399.615 484.999 -0.4193 -0.0514 -135.6551 0.03 0.004 0.010<br />

85 207335.142 712575.654 6179284.390 483.470 -0.4287 -0.0452 -135.6093 0.03 0.004 0.010<br />

86 207337.437 712505.923 6179174.552 482.996 -0.4719 -0.0617 -135.6427 0.03 0.004 0.010<br />

87 207339.843 712433.324 6179059.249 482.532 -0.4952 -0.0712 -135.6373 0.03 0.004 0.010<br />

88 207342.152 712363.537 6178948.223 483.014 -0.3662 -0.0308 -135.9181 0.03 0.004 0.010<br />

89 207344.433 712294.162 6178838.651 483.583 -0.4786 -0.0549 -135.7456 0.03 0.004 0.010<br />

94


90 207346.838 712220.913 6178722.973 484.359 -0.4004 -0.0431 -135.9113 0.03 0.004 0.010<br />

91 207349.135 712151.008 6178612.684 484.455 -0.4313 0.0001 -135.3428 0.03 0.004 0.010<br />

92 207351.539 712077.821 6178497.042 485.418 -0.4979 -0.0718 -135.6601 0.03 0.004 0.010<br />

93 207353.834 712008.006 6178386.569 487.377 -0.4282 -0.0013 -135.8367 0.03 0.004 0.010<br />

95


BILAG - BEREGNING AF SOLVINKEL<br />

Solvinkler er beregnet vha. af en udregning på http://aa.usno.navy.mil/data/docs/AltAz.php.<br />

Det gav følgende tabel:<br />

Astronomical Applications Dept.<br />

U.S. Naval Observatory<br />

Washington. DC 20392-5420<br />

(NO NAME GIVEN)<br />

o . o .<br />

E 12 23. N55 24<br />

Altitude and Azimuth of the Sun<br />

Apr 7. 2009<br />

Zone: 1h East of Greenwich<br />

Altitude Azimuth<br />

(E of N)<br />

h m o o<br />

04:30 -8.4 64.8<br />

05:00 -4.4 71.1<br />

05:30 0.2 77.3<br />

06:00 4.1 83.5<br />

06:30 8.2 89.6<br />

07:00 12.5 95.8<br />

07:30 16.6 102.2<br />

08:00 20.7 108.7<br />

08:30 24.7 115.6<br />

09:00 28.4 122.8<br />

09:30 31.8 130.4<br />

10:00 34.9 138.6<br />

10:30 37.4 147.2<br />

11:00 39.5 156.4<br />

11:30 40.8 166.0<br />

12:00 41.5 175.8<br />

12:30 41.5 185.8<br />

13:00 40.7 195.6<br />

13:30 39.2 205.1<br />

14:00 37.1 214.2<br />

14:30 34.5 222.8<br />

15:00 31.4 230.9<br />

15:30 27.9 238.4<br />

16:00 24.2 245.6<br />

16:30 20.2 252.4<br />

17:00 16.1 258.9<br />

17:30 11.9 265.3<br />

18:00 7.7 271.5<br />

18:30 3.6 277.6<br />

19:00 -0.8 283.8<br />

19:30 -4.9 290.0<br />

20:00 -8.8 296.4<br />

96


BILAG – DOKUMENTATION FOR PASPUNKTER<br />

L5-projekt – Kortlægning <strong>og</strong> afsætning<br />

Paspunkter model 84B – 85B<br />

Skitse Beskrivelse Koordinater (UTM. DVR90)<br />

Midte af rist<br />

Sydlig midte mellem to vejstriber<br />

Midte af rist<br />

97<br />

E<br />

712595.993<br />

N<br />

6179451.489<br />

Z<br />

28.367<br />

E<br />

712595.147<br />

N<br />

6179445.424<br />

Z<br />

28.330<br />

E<br />

712510.865<br />

N<br />

6179310.533<br />

Z<br />

25.405<br />

700<br />

701<br />

702


Skitse Beskrivelse Koordinater (UTM. DVR90)<br />

Sydlig midte mellem to vejstriber<br />

Samling mellem sten<br />

Samling mellem sten<br />

98<br />

E<br />

712514.684<br />

N<br />

6179310.645<br />

Z<br />

25.413<br />

E<br />

712750.553<br />

N<br />

6179156.499<br />

Z<br />

22.518<br />

E<br />

712750.050<br />

N<br />

6179155.630<br />

Z<br />

22.530<br />

703<br />

704<br />

705


Skitse Beskrivelse Koordinater (UTM. DVR90)<br />

Samling mellem sten<br />

Midte af vejstribe ud for skilt<br />

Midte af sti ved knæk<br />

99<br />

E<br />

712856.747<br />

N<br />

6179290.685<br />

Z<br />

24.018<br />

E<br />

712856.303<br />

N<br />

6179290.907<br />

Z<br />

24.028<br />

E<br />

712389.168<br />

N<br />

6179386.624<br />

Z<br />

25.277<br />

706<br />

707<br />

600


Skitse Beskrivelse Koordinater (UTM. DVR90)<br />

Midte af sti ved knæk<br />

100<br />

E<br />

712494.978<br />

N<br />

6179468.933<br />

Z<br />

28.732<br />

602

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!