<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EpostStil17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="NO-BOK" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Hi,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span lang="EN-US">The Norwegian Mapping Authority is working on getting our old reference frame, NGO1948, implemented In proj.
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">The transformation we want to implement is the transformation from/to EUREF89 – NGO1948<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">We already have a solution that works, but we have been advised to discuss the use of other libraries to solve our problem.
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Some background information: <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">When the Norwegian Mapping Authority introduced the new geodetic datum EUREF89, it was necessary to find a proper transformation technique between the deprecated geodetic datum NGO1948 to the EUREF89. NGO1948 was significant
deformed, hereby it was computed and realized in different areas of the country. A consequence of this is huge planar gaps in some counties and municipalities. The gaps are upon 2-3 meters on borders. Description and evaluation of the method are further documented
in the articles see [OMathisen2002] and [OMathisen2003]. The processing is done in three steps:
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">1. Selection of point pairs <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">2. Deterministic step: 2D Helmert transformation
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">3. Statistic step: Smoothing Least Squared Collocation<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">In the solution we have now we have used an open source state of the art matrix library (EIGEN3) to compute the Helmert parameters with LSC. This process will demand inverting and other matrix operations on large matrices
(40x40).<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">The size of the matrices are also dynamic.
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">At this point there is no matrix library in proj that work with large matrices. Does anyone have a suggestion on how we can solve this? Can we install external packages or do we have to write our own matrix library, so
that the library is a part of the source code?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:NO-BOK">Best Regards<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:NO-BOK"> <o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:NO-BOK">Karoline A. Skaar<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:NO-BOK">Norwegian Mapping Authority<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:NO-BOK">Tlf./mob.: 32 11 82 16/907 40 302<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:NO-BOK">E-post: <a href="mailto:karoline.skar@kartverket.no">
<span style="color:#0563C1">karoline.skar@kartverket.no</span></a><o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:NO-BOK"> <o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:NO-BOK"><a href="http://www.kartverket.no/"><span style="color:blue">www.kartverket.no</span></a><o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:NO-BOK"> <o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#1F497D;mso-fareast-language:NO-BOK"><img border="0" width="89" height="56" style="width:.925in;height:.5833in" id="Picture_x0020_1" src="cid:image001.jpg@01D691A5.764EDEC0" alt="logo_sign"></span><span style="mso-fareast-language:NO-BOK"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:NO-BOK"> <o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>