[OSM-talk-nl] [algorithme gezocht]

Mon Oct 22 20:33:45 UTC 2007

Mijn coordinaten stelsel is

Natuurlijk dat van OSM

Graden minuten en seconden, of de decimale versie daarvan

Op de aangegeven link van martijn vindt je

Nog veel meer op 

http://mathforum.org/library/drmath/view/51879.html

Gert

Van: talk-nl-bounces at openstreetmap.org [mailto:talk-nl-bounces at openstreetmap.org] Namens Martijn van Exel
Verzonden: maandag 22 oktober 2007 22:26
Aan: mlinden at zeelandnet.nl; OpenStreetMap NL discussion list
Onderwerp: Re: [OSM-talk-nl] [algorithme gezocht]

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

Sorry voor de verwarring! 

Martijn

Op 22-okt-2007, om 22:03 heeft Milo van der Linden het volgende geschreven:

Martijn van Exel schreef: 

Voor afstandberekening in het RD-stelsel kun je volgens mij inderdaad  
veilig pythagoras gebruiken, omdat het een km2-grid is met vaste  
afstand tussen de gridlijnen en rechte hoeken. Maar als je  
cartesische coördinaten hebt (zoals Gert heeft)  zou je die eerst  
moeten omrekenen naar RD en dan krijg je ook met dergelijke  
trigonometrie te maken.

Dit is niet helemaal juist, cartesisch betekent namelijk "plat". Gert zijn coordinaten kunnen dus ook gebruik maken van dezelfde stelling van pythagoras. 

Feitelijk zijn er 2 soorten coordinaatstelsels op GIS gebied: 

1) Geografische coordinaat systemen 
en 
2) "geprojecteerde" coordinaat systemen. 

In een coordinaat systeem uit 1) geldt dat trigonometrie of "great circle" berekeningen dienen te worden toegepast.
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. 

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.

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.

(Pfff, 't zat diep weg in de uithoeken van mijn RAM, maar als HTS Geo-Informaticus/Geodeet moest ik hier even op reageren)

Martijn

Op 22-okt-2007, om 20:40 heeft Just van den Broecke het volgende  
geschreven:

	Volgens mij heet de formule "great circle distance" (zoals c#
	voorbeeld). Zie http://en.wikipedia.org/wiki/Great-circle_distance

	Pythagoras kan m.i. alleen met zgn "equidistante projectie" (bijv.
	pixels evenredig aan meters in x en y richting) zoals RD op bijv.
	topokaarten in NL.

	Ik heb wel een aantal implementaties bijv. in Java of JavaScript.

	Just

	Martijn van Exel wrote:

		Dat maakt het een stuk complexer. Je moet je afvragen of je dat  
		ervoor
		over hebt, maar enfin, ik vind de volgende uitwerking:

		http://www.codeproject.com/csharp/distancebetweenlocations.asp

		Succes ermee;

		Martijn

		Op 22-okt-2007, om 13:29 heeft Gert Gremmen het volgende geschreven:

			Dank je Martijn,
			maar nu in lengte en breedte graden....
			en over de aardbol, hoewel dat
			binnen Nederland wel  ongeveer hetzelfde
			is als op een plat vlak.

			Regards,

			Je bedoelt de hemelsbrede afstand?

			Dat is de implementatie van Pythagoras:

			(in c#)
			public double Distance(Point pt1, Point pt2)
			{
			return Math.Sqrt(Math.Pow(pt2.X.ToDouble() - pt1.X.ToDouble(),
			2) +
			Math.Pow(pt2.Y.ToDouble() - pt1.Y.ToDouble(), 2));
			}

			Voor afstand over de weg kun je een webservice aanroepen, er zijn
			verschillende partijen die dat aanbieden (Locatienet, Geodan, ...)

			Martijn

			Op 22-okt-2007, om 13:04 heeft Gert Gremmen het volgende geschreven:

				Voor het schrijven van een eenvoudige applicatie
				zoek ik een algorithme om de afstand tussen
				twee coordinaten ( uit OSM natuurlijk)
				te berekenen. (in Nederland)
				Het hoeft niet sub-meter nauwkeurig....
				Wie heeft een voorbeeld in universele code?

				Regards,

				 Gert Gremmen

				_______________________________________________
				Talk-nl mailing list
				Talk-nl at openstreetmap.org <mailto:Talk-nl at openstreetmap.org> <mailto:Talk-nl at openstreetmap.org> 
				http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl

			_______________________________________________
			Talk-nl mailing list
			Talk-nl at openstreetmap.org <mailto:Talk-nl at openstreetmap.org> <mailto:Talk-nl at openstreetmap.org> 
			http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl

			_______________________________________________
			Talk-nl mailing list
			Talk-nl at openstreetmap.org <mailto:Talk-nl at openstreetmap.org> <mailto:Talk-nl at openstreetmap.org> 
			http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl

		--------------------------------------------------------------------- 
		---

		_______________________________________________
		Talk-nl mailing list
		Talk-nl at openstreetmap.org
		http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl

	-- 

	--Just

	Just van den Broecke  just at justobjects.nl
	Just Objects B.V.     tel +31 65 4268627 Skype: justb4
	The Netherlands       http://www.justobjects.nl

	_______________________________________________
	Talk-nl mailing list
	Talk-nl at openstreetmap.org
	http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl

_______________________________________________
Talk-nl mailing list
Talk-nl at openstreetmap.org
http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl

-- 

	Milo van der Linden 
skype: milovanderlinden <skype:milovanderlinden?add> 
mlinden at zeelandnet.nl
milovanderlinden at gmail.com
milo at 3dsite.nl
http://www.3dsite.nl

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

_______________________________________________

Talk-nl mailing list

Talk-nl at openstreetmap.org

http://lists.openstreetmap.org/cgi-bin/mailman/listinfo/talk-nl

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.openstreetmap.org/pipermail/talk-nl/attachments/20071022/93f46b78/attachment.htm>