konvertering sweref - rt90

[Leif]

någon som har ett lätt sätt att konvertera koordinater inmätta i sweref till rt90?

[Nils Ristare]

Lantmäteriet har en "kalkylator" på webben:

http://www.lantmateriet.se/templates/LMV_Enkelkoordinattransformation.aspx?id=11500

[Kepler]

Det finns även ett gratisprogram som heter GeoPos tillgängligt från Kartforum http://www.djuvfeldt.se/fusion/downloads.php?cat_id=3.

Själv har jag tagit fram ett Excelblad för dessa konverteringar: SETrans. Det ger exakt samma resultat som Lantmäteriets GTrans. Eftersom man matar in koordinaterna i ett Excelblad fungerar det utmärkt om man har många punkter. Om det är någon som är intresserad skickar jag det gärna. 

[Micke]

Är konverteringen komplicerad? Om man kan göra det enkelt i ett excellblad antar jag att konverteringen är rätt banal. Därmed borde det vara hyfsat enkelt att skriva ett program som t.ex. tar en tabavgränsad textfil med koordinater och som spottar ut en fil med koordinaterna i det nya systemet. Om någon kan berätta för mig exakt hur konvereteringen går till skulle jag nog kunna slänga i hop nåt i c++...

/Micke 

[Micke]

Jag hittade den här pdf:en och jag måste erkänna att matten inte verkar helt enkel...

/Micke

[Kepler]

Det är inte helt trivialt även om matematiken inte är så farlig (om man bortser från en matrisräkning som inte tillhör gymnasiematten). Det gäller mer att hålla tungan rätt i mun när man skriver in formlerna. Min Excelfil innerhåller ca 900 rader skrivna i Visual Basic for Applications. Om du vill får du gärna koden.

Det som tar tid är egentligen att förstå de bakomliggande principerna för beräkningproceduren. Det finns enklare ämnen än geodesi.

[Jan Owe]

Jag har gjort en rutin som konverterar från RT90 till Latitud och Longitud.
Rutinen bygger på det som står i dokumentet
http://www.lantmateriet.se/upload/filer/kartor/geodesi_gps_och_detaljmatning/geodesi/Formelsamling/Gauss_Conformal_Projection.pdf
Och på de konstanter som står i dokumentet
http://www.lantmateriet.se/upload/filer/kartor/geodesi_gps_och_detaljmatning/Nytt_referenssytem/Infoblad/info_blad-12_transf_swe99_rt90.pdf
Se vidare min bifogade pascalkod. Genom att följa formlerna på sidorna 2 och 3 i det förstnämnda dokumentet kan man istället konvertera åt andra hållet - och koden blir inte så lång :-).
/Jan

[tobbedj]

Hej

Jag är personen bakom programmet GeoPos som länkas till tidigare i tråden.

Om det programmet saknar funktionalitet så tar jag gärna emot synpunkter och önskemål om hur programmet kan göras bättre.

Kepler's excel-fil har sitt användningsområde, GeoPos ett annat. Men varför uppfinna hjulet på nytt när gratisalternativen redan finns?

/torbjörn

[BeeWee]

Man kan naturligtvis också köpa programmet Swetrans från Swegis. http://www.swegis.se/swedish/software/swetrans/swetrans.asp

[Kepler]

Jan Owe,

Det är inte riktigt så enkelt. Vad din algoritm gör är att konvertera plana koordinater (x, y) i referenssystemet RT90 till lat/lon i samma referenssystem. Dock används aldrig dessa lat/lon värden; istället anges de i referensystemen SWEREF 99 / WGS 84 - och då blir det något andra värden. Man måste därför transformera från RT90 till SWEREF 99.

GPS mottagaren räknar primärt ut lat/lon i referenssystem WGS 84 som i praktiken är det samma som SWEREF 99 (skillnaden är max ett par decimeter). Sedan kan man ställa in GPS mottagaren så att den räknar om och visar t.ex. plana koordinater i RT 90.




 

[Micke]

Men varför uppfinna hjulet på nytt när gratisalternativen redan finns?

Det är visserligen bättre att ge bort binära filer än att ta betalt för dem, men det bästa är om man gör programmet helt fritt och inte bara gratis. Om du släpper ditt program under en fri licens så behöver ingen uppfinna hjulet flera gånger.  Personligen har jag ingen nytta av ditt gratisprogram eftersom du bara gjort de körbara filerna tillgängliga för Windows som för mig är ett oanvändbart operativsystem. Jag kan dessutom inte lära mig särskilt mycket av ditt program, eftersom jag inte får tillgång till källkoden.

Min Excelfil innerhåller ca 900 rader skrivna i Visual Basic for Applications. Om du vill får du gärna koden.

Jag skulle gärna titta på koden.

/Micke

[tobbedj]

Om du släpper ditt program under en fri licens så behöver ingen uppfinna hjulet flera gånger.

Den synpunkten har inte dykt upp förrän nu när det gäller GeoPos. Det skulle dessutom vara svårt att realisera då flera funktioner i programmet bygger på inköpta tredjepartskomponenter som inte kan lämnas ut för fri användning.

Surt, sa räven...

[merrimack]

Jag har publicerat ett fritt och öppet (Creative Commons) bibliotek i .NET som konverterar kooridnater mellan SWEREF99, RT90 och WGS84. Det är fritt att använda för en egen applikation eller att porta koden till andra språk. Koden finns på min blog: http://blog.sallarp.com/translate-coordinates-between-rt90-wgs84-and-sweref99-using-net

[Micke]

Jag har publicerat ett fritt och öppet (Creative Commons) bibliotek i .NET som konverterar kooridnater mellan SWEREF99, RT90 och WGS84. Det är fritt att använda för en egen applikation eller att porta koden till andra språk. Koden finns på min blog: http://blog.sallarp.com/translate-coordinates-between-rt90-wgs84-and-sweref99-using-net

Schysst! Är det den här js-koden du utgick i från alltså?

http://mellifica.se/geodesi/gausskruger.js

/Micke

[Micke]

Jag har nu fått behov av att konvertera rt90 2.5 gon v till latitude och longitude och skrev därför ett tunt php omslag till det fina javascript biblioteket. Så här blev det:

Kodstycke: [Välj]

<?php
if(isset($_GET[type])) {
    $type = addslashes(strip_tags($_GET[type]));
}
else {
    $type = "rt90_2.5_gon_v";
}
$xcoord = addslashes(strip_tags($_GET[xcoord]));
$ycoord = addslashes(strip_tags($_GET[ycoord]));
$s = `rhino coord.js $type $xcoord $ycoord`;
$lat_lon = explode(',', trim($s, "\n"));
header('Content-Type: text/xml version=1.0 charset=UTF-8');
header('Content-Disposition: inline; filename=coord.xml');
echo "<?xml version=\"1.0\" encoding=\"UTF-8\"?>                                                                                                            
<!DOCTYPE coord [
<!ELEMENT coord (point) >
<!ELEMENT point (xcoord, ycoord, lat, long) >
<!ELEMENT xcoord (#PCDATA) >
<!ELEMENT ycoord (#PCDATA) >
<!ELEMENT lat (#PCDATA) >
<!ELEMENT long (#PCDATA) >
]>                                                                                                                                      
<coord>
    <point>
        <xcoord>$xcoord</xcoord>
        <ycoord>$ycoord</ycoord>
        <lat>$lat_lon[0]</lat>
        <long>$lat_lon[1]</long>
    </point>
</coord>\n";
?>
 

javascriptet (coord.js) som anropas med rhino ser ut så här:

Kodstycke: [Välj]

load("http://latlong.mellifica.se/latlong.js", "http://mellifica.se/geodesi/lat_lon_conv.js", "http://mellifica.se/geodesi/gausskruger.js");
swedish_params(arguments[0]);
var lat_lon = grid_to_geodetic(arguments[1], arguments[2]);
print( lat_lon[0] + ','+ lat_lon[1]);

Här kan man kolla på en demo:
http://mickenordin.se/coord/coord.php?type=rt90_2.5_gon_v&xcoord=6575900.000&ycoord=1575700.000
som type kan man ange:

• rt90_7.5_gon_v
• rt90_5.0_gon_v
• rt90_2.5_gon_v
• rt90_0.0_gon_v
• rt90_2.5_gon_o
• rt90_5.0_gon_o

Det går också köra utan type:
http://mickenordin.se/coord/coord.php?xcoord=6575900.000&ycoord=1575700.000

I så fall defaultar den till rt 90 2.5 gon v
Man får gärna använda denna lilla kod under samma villkor som js.biblioteket:
Creative Commons: Attribution-Noncommercial-Share Alike

/Micke

[er-bjoer]

Det finns även ett gratisprogram som heter GeoPos tillgängligt från Kartforum http://www.djuvfeldt.se/fusion/downloads.php?cat_id=3.

Själv har jag tagit fram ett Excelblad för dessa konverteringar: SETrans. Det ger exakt samma resultat som Lantmäteriets GTrans. Eftersom man matar in koordinaterna i ett Excelblad fungerar det utmärkt om man har många punkter. Om det är någon som är intresserad skickar jag det gärna. 

Hei Kepler!

Er du fortsatt villig til å dele dette excel arket ?
Jeg er veldig interessert dersom du har mulighet

//Erik

[Kepler]

Visst, inga problem. Du bör ha filen i din mail nu.

Mvh Kepler.

[perssoh1]

Hej Kepler

Skulle vara mycket tacksam ifall du fortfarande bidrar med ditt Excel-alster! Vill kunna göra om Eniros WGS84-koordinater till RT90. Med RT90 koordinater kan man sedan räkna ut en orts medeltamperatur med data från SMHI etc ...

Mycket tacksam! 

mvh Håkan

[john_g]

Visst, inga problem. Du bör ha filen i din mail nu.

Mvh Kepler.

Kan du skicka den till min också?  Skulle bli överlycklig för en sån grej.

/John

[Kepler]

Javisst, men jag kan inte se din e-post adress. Men du kan ju skicka den i ett meddelande.
/Ulf