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