<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 12 (filtered medium)">
<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;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
pre
{mso-style-priority:99;
mso-style-link:"HTML - vooraf opgemaakt Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";}
span.HTML-voorafopgemaaktChar
{mso-style-name:"HTML - vooraf opgemaakt Char";
mso-style-priority:99;
mso-style-link:"HTML - vooraf opgemaakt";
font-family:Consolas;}
p.style4, li.style4, div.style4
{mso-style-name:style4;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
p.style5, li.style5, div.style5
{mso-style-name:style5;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
span.E-mailStijl21
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page Section1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.Section1
{page:Section1;}
-->
</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=NL link=blue vlink=purple style='word-wrap: break-word;-webkit-nbsp-mode: space;
-webkit-line-break: after-white-space'>
<div class=Section1>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Mijn coordinaten stelsel is<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Natuurlijk dat van OSM<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Graden minuten en seconden, of de decimale versie daarvan<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Op de aangegeven link van martijn vindt je<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Nog veel meer op <o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><a href="http://mathforum.org/library/drmath/view/51879.html">http://mathforum.org/library/drmath/view/51879.html</a><o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Gert<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<div>
<div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm'>
<p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>Van:</span></b><span
style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>
talk-nl-bounces@openstreetmap.org [mailto:talk-nl-bounces@openstreetmap.org] <b>Namens
</b>Martijn van Exel<br>
<b>Verzonden:</b> maandag 22 oktober 2007 22:26<br>
<b>Aan:</b> mlinden@zeelandnet.nl; OpenStreetMap NL discussion list<br>
<b>Onderwerp:</b> Re: [OSM-talk-nl] [algorithme gezocht]<o:p></o:p></span></p>
</div>
</div>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Ik ben blij dat er iemand is die het wel echt weet; na het
posten van dit bericht begon ik te twijfelen aan mijn formulering en ben
sindsdien aan het graven in geheugen en bronnen, maar mijn intellectuele
vermogens lijken een bermudadriehoek te hebben waar het gaat om projecties en
coördinaatsystemen...<o:p></o:p></p>
<div>
<p class=MsoNormal>Sorry voor de verwarring! <o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<div>
<p class=MsoNormal>Martijn<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><o:p> </o:p></p>
<div>
<div>
<p class=MsoNormal>Op 22-okt-2007, om 22:03 heeft Milo van der Linden het
volgende geschreven:<o:p></o:p></p>
</div>
<p class=MsoNormal><br>
<br>
<o:p></o:p></p>
<p class=MsoNormal>Martijn van Exel schreef: <o:p></o:p></p>
<pre>Voor afstandberekening in het RD-stelsel kun je volgens mij inderdaad <o:p></o:p></pre><pre>veilig pythagoras gebruiken, omdat het een km2-grid is met vaste <o:p></o:p></pre><pre>afstand tussen de gridlijnen en rechte hoeken. Maar als je <o:p></o:p></pre><pre>cartesische coördinaten hebt (zoals Gert heeft) zou je die eerst <o:p></o:p></pre><pre>moeten omrekenen naar RD en dan krijg je ook met dergelijke <o:p></o:p></pre><pre>trigonometrie te maken.<o:p></o:p></pre><pre> <o:p></o:p></pre>
<p class=MsoNormal>Dit is niet helemaal juist, cartesisch betekent namelijk
"plat". Gert zijn coordinaten kunnen dus ook gebruik maken van
dezelfde stelling van pythagoras. <br>
<br>
Feitelijk zijn er 2 soorten coordinaatstelsels op GIS gebied: <br>
<br>
1) Geografische coordinaat systemen <br>
en <br>
2) "geprojecteerde" coordinaat systemen. <br>
<br>
In een coordinaat systeem uit 1) geldt dat trigonometrie of "great
circle" berekeningen dienen te worden toegepast.<br>
Bij 2) kun je altijd volstaan met de stelling van pythagoras, maar afhankelijk
van de projectie zullen er afwijkingen optreden in afstand of in oppervlakte.
Er zijn namelijk zogenaamde "Equal Area" (Oppervlakte getrouw) en
"Equidistant" (Afstands getrouw) projecties en mengvormen van deze
twee. Vliegroute kaarten zijn b.v. equidistant omdat afstand hier belangrijker
is dan oppervlakte. Thematische vlakken kaarten zijn vaak oppervlakte getrouw. <br>
<br>
Tevens is het gebied waarvoor de projectie specifiek is ingericht essentieel.
In Nederland snijdt de projectie bijvoorbeeld door de aarde in een specifieke
straal rondom Amersfoort. Dit betekent dat je vlak bij Amersfoort een negatieve
afwijking hebt, (in Zeeland en Groningen een positieve) van enkele centimeters.
Op de snijring is de afwijking het kleinste.<br>
<br>
Indien de afstandsberekening bedoeld is voor b.v. routeplanning zou ik daarom
adviseren om de extra moeite te nemen om de "great circle" berekening
door te voeren. Deze is dan automatische bij benadering gelijk aan de GPS
afstand.<br>
<br>
(Pfff, 't zat diep weg in de uithoeken van mijn RAM, maar als HTS
Geo-Informaticus/Geodeet moest ik hier even op reageren)<br>
<br>
<br>
<o:p></o:p></p>
<pre>Martijn<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Op 22-okt-2007, om 20:40 heeft Just van den Broecke het volgende <o:p></o:p></pre><pre>geschreven:<o:p></o:p></pre><pre><o:p> </o:p></pre><pre> <o:p></o:p></pre>
<blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><pre>Volgens mij heet de formule "great circle distance" (zoals c#<o:p></o:p></pre><pre>voorbeeld). Zie <a
href="http://en.wikipedia.org/wiki/Great-circle_distance">http://en.wikipedia.org/wiki/Great-circle_distance</a><o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Pythagoras kan m.i. alleen met zgn "equidistante projectie" (bijv.<o:p></o:p></pre><pre>pixels evenredig aan meters in x en y richting) zoals RD op bijv.<o:p></o:p></pre><pre>topokaarten in NL.<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Ik heb wel een aantal implementaties bijv. in Java of JavaScript.<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Just<o:p></o:p></pre><pre><o:p> </o:p></pre><pre><o:p> </o:p></pre><pre><o:p> </o:p></pre><pre>Martijn van Exel wrote:<o:p></o:p></pre><pre> <o:p></o:p></pre>
<blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><pre>Dat maakt het een stuk complexer. Je moet je afvragen of je dat <o:p></o:p></pre><pre>ervoor<o:p></o:p></pre><pre>over hebt, maar enfin, ik vind de volgende uitwerking:<o:p></o:p></pre><pre><o:p> </o:p></pre><pre><a
href="http://www.codeproject.com/csharp/distancebetweenlocations.asp">http://www.codeproject.com/csharp/distancebetweenlocations.asp</a><o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Succes ermee;<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Martijn<o:p></o:p></pre><pre><o:p> </o:p></pre><pre><o:p> </o:p></pre><pre>Op 22-okt-2007, om 13:29 heeft Gert Gremmen het volgende geschreven:<o:p></o:p></pre><pre><o:p> </o:p></pre><pre> <o:p></o:p></pre>
<blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><pre>Dank je Martijn,<o:p></o:p></pre><pre>maar nu in lengte en breedte graden....<o:p></o:p></pre><pre>en over de aardbol, hoewel dat<o:p></o:p></pre><pre>binnen Nederland wel ongeveer hetzelfde<o:p></o:p></pre><pre>is als op een plat vlak.<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Regards,<o:p></o:p></pre><pre><o:p> </o:p></pre><pre><o:p> </o:p></pre><pre><o:p> </o:p></pre><pre>Je bedoelt de hemelsbrede afstand?<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Dat is de implementatie van Pythagoras:<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>(in c#)<o:p></o:p></pre><pre>public double Distance(Point pt1, Point pt2)<o:p></o:p></pre><pre>{<o:p></o:p></pre><pre>return Math.Sqrt(Math.Pow(pt2.X.ToDouble() - pt1.X.ToDouble(),<o:p></o:p></pre><pre>2) +<o:p></o:p></pre><pre>Math.Pow(pt2.Y.ToDouble() - pt1.Y.ToDouble(), 2));<o:p></o:p></pre><pre>}<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Voor afstand over de weg kun je een webservice aanroepen, er zijn<o:p></o:p></pre><pre>verschillende partijen die dat aanbieden (Locatienet, Geodan, ...)<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Martijn<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Op 22-okt-2007, om 13:04 heeft Gert Gremmen het volgende geschreven:<o:p></o:p></pre><pre><o:p> </o:p></pre><pre> <o:p></o:p></pre>
<blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><pre>Voor het schrijven van een eenvoudige applicatie<o:p></o:p></pre><pre>zoek ik een algorithme om de afstand tussen<o:p></o:p></pre><pre>twee coordinaten ( uit OSM natuurlijk)<o:p></o:p></pre><pre>te berekenen. (in Nederland)<o:p></o:p></pre><pre>Het hoeft niet sub-meter nauwkeurig....<o:p></o:p></pre><pre>Wie heeft een voorbeeld in universele code?<o:p></o:p></pre><pre><o:p> </o:p></pre><pre><o:p> </o:p></pre><pre>Regards,<o:p></o:p></pre><pre><o:p> </o:p></pre><pre> Gert Gremmen<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>_______________________________________________<o:p></o:p></pre><pre>Talk-nl mailing list<o:p></o:p></pre><pre><a
href="mailto:Talk-nl@openstreetmap.org">Talk-nl@openstreetmap.org</a> <a
href="mailto:Talk-nl@openstreetmap.org"><mailto:Talk-nl@openstreetmap.org></a><o:p></o:p></pre><pre><a
href="http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl">http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl</a><o:p></o:p></pre><pre> <o:p></o:p></pre></blockquote>
<pre>_______________________________________________<o:p></o:p></pre><pre>Talk-nl mailing list<o:p></o:p></pre><pre><a
href="mailto:Talk-nl@openstreetmap.org">Talk-nl@openstreetmap.org</a> <a
href="mailto:Talk-nl@openstreetmap.org"><mailto:Talk-nl@openstreetmap.org></a><o:p></o:p></pre><pre><a
href="http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl">http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl</a><o:p></o:p></pre><pre><o:p> </o:p></pre><pre>_______________________________________________<o:p></o:p></pre><pre>Talk-nl mailing list<o:p></o:p></pre><pre><a
href="mailto:Talk-nl@openstreetmap.org">Talk-nl@openstreetmap.org</a> <a
href="mailto:Talk-nl@openstreetmap.org"><mailto:Talk-nl@openstreetmap.org></a><o:p></o:p></pre><pre><a
href="http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl">http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl</a><o:p></o:p></pre><pre> <o:p></o:p></pre></blockquote>
<pre>--------------------------------------------------------------------- <o:p></o:p></pre><pre>---<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>_______________________________________________<o:p></o:p></pre><pre>Talk-nl mailing list<o:p></o:p></pre><pre><a
href="mailto:Talk-nl@openstreetmap.org">Talk-nl@openstreetmap.org</a><o:p></o:p></pre><pre><a
href="http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl">http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl</a><o:p></o:p></pre><pre> <o:p></o:p></pre></blockquote>
<pre>-- <o:p></o:p></pre><pre><o:p> </o:p></pre><pre>--Just<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Just van den Broecke <a
href="mailto:just@justobjects.nl">just@justobjects.nl</a><o:p></o:p></pre><pre>Just Objects B.V. tel +31 65 4268627 Skype: justb4<o:p></o:p></pre><pre>The Netherlands <a
href="http://www.justobjects.nl">http://www.justobjects.nl</a><o:p></o:p></pre><pre><o:p> </o:p></pre><pre><o:p> </o:p></pre><pre><o:p> </o:p></pre><pre><o:p> </o:p></pre><pre><o:p> </o:p></pre><pre><o:p> </o:p></pre><pre>_______________________________________________<o:p></o:p></pre><pre>Talk-nl mailing list<o:p></o:p></pre><pre><a
href="mailto:Talk-nl@openstreetmap.org">Talk-nl@openstreetmap.org</a><o:p></o:p></pre><pre><a
href="http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl">http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl</a><o:p></o:p></pre><pre> <o:p></o:p></pre></blockquote>
<pre>_______________________________________________<o:p></o:p></pre><pre>Talk-nl mailing list<o:p></o:p></pre><pre><a
href="mailto:Talk-nl@openstreetmap.org">Talk-nl@openstreetmap.org</a><o:p></o:p></pre><pre><a
href="http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl">http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl</a><o:p></o:p></pre><pre><o:p> </o:p></pre><pre> <o:p></o:p></pre>
<p class=MsoNormal style='margin-bottom:12.0pt'><o:p> </o:p></p>
<div>
<p class=MsoNormal style='margin-bottom:12.0pt'>-- <o:p></o:p></p>
<table class=MsoNormalTable border=0 cellspacing=0 cellpadding=0 width=400
style='width:300.0pt'>
<tr style='height:58.5pt'>
<td width=15 rowspan=2 valign=bottom style='width:11.25pt;padding:0cm 0cm 0cm 0cm;
height:58.5pt'></td>
<td width=0 colspan=2 valign=top style='width:.3pt;padding:0cm 0cm 0cm 0cm;
height:58.5pt'>
<p class=style4>Milo van der Linden <br>
<a href="skype:milovanderlinden?add">skype: milovanderlinden</a><br>
<a href="mailto:mlinden@zeelandnet.nl">mlinden@zeelandnet.nl</a><br>
<a href="mailto:milovanderlinden@gmail.com">milovanderlinden@gmail.com</a><br>
<a href="mailto:milo@3dsite.nl">milo@3dsite.nl</a><br>
<a href="http://www.3dsite.nl">http://www.3dsite.nl</a><o:p></o:p></p>
</td>
<td width=15 rowspan=2 valign=top style='width:11.25pt;padding:0cm 0cm 0cm 0cm;
height:58.5pt'>
<p class=MsoNormal> <o:p></o:p></p>
</td>
<td width=300 valign=top style='width:225.0pt;padding:0cm 0cm 0cm 0cm;
height:58.5pt'>
<p class=style5>De informatie in dit bericht reflecteert mijn persoonlijke
mening en niet die van een bedrijf of instantie. Aan de informatie kunnen
geen rechten worden ontleend. Indien dit bericht onderdeel is van een forum,
mailing-list of community dan gelden automatisch de bij het betreffende
medium behorende voorwaarden. <span lang=EN>The information in this message
reflects my personal opinion and not that of a company or public body. All
rights reserved.If this message is contained in a mailing-list or community,
the rights on the medium are automatically adapted.</span><o:p></o:p></p>
</td>
</tr>
<tr>
<td style='padding:0cm 0cm 0cm 0cm'></td>
<td style='padding:0cm 0cm 0cm 0cm'></td>
<td style='padding:0cm 0cm 0cm 0cm'></td>
</tr>
<tr height=0>
<td width=7 style='border:none'></td>
<td width=97 style='border:none'></td>
<td width=97 style='border:none'></td>
<td width=9 style='border:none'></td>
<td width=190 style='border:none'></td>
</tr>
</table>
</div>
<div>
<p class=MsoNormal>_______________________________________________<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal>Talk-nl mailing list<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><a href="mailto:Talk-nl@openstreetmap.org">Talk-nl@openstreetmap.org</a><o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><a
href="http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl">http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl</a><o:p></o:p></p>
</div>
</div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
</div>
</body>
</html>