henca

[henca]

Snyggt jobbat!

m v h Henrik

[henca]

Tack! De är nu inlagda som K301W01 och K301E02.

m v h Henrik

[henca]

Tack för positionerna! Inlagda som:
K509E01
K509S02
K509S03
K509N04
K509S05
K509N06
K509S07

m v h Henrik

[henca]

Sen nyhetsbrevets integration mot bbPress avaktiverades har det inte kommit några fler ”bekräfta-mail” för varje skrivet inlägg.

m v h Henrik

[henca]

Nåväl, det är lättare att ta bort än att lägga till 🙂

m v h Henrik

[henca]

Snabbt och snyggt jobbat!

Kanske hade vi redan Torsång sen tidigare, då jag kikade på
http://poi.fartkameraforum.se/www/poi_list.php?selected_database=bulor
fanns den i dubbel upplaga med exakt samma koordinater. Jag tog bort den ena nu, antagligen den som var inmatad sen tidigare då det fortfarande står att Torsång är vår nyaste vägbula.

m v h Henrik

[henca]

Trots att rutan inte syns och varken går att kryssa i eller ur så kommer det fortfarande ett bekräftelse-mail för varje inlägg man skriver på forumet.

m v h Henrik

[henca]

[quote quote=232]Diskussionen pågår redan i tråden om IGO men jag tycker att vi försöker hålla ett ämne per tråd så jag kastar in denna för att fortsätta.[/quote]

Bra initiativ! Då får resten av igo8-tråden handla om just igo8.

m v h Henrik

[henca]

När du skrev att du kunde konvertera en koordinat utan problem trodde jag att du menade att du hade programkod och formel för konvertering av en koordinat. Då menade jag att det skulle vara trivialt att modifiera programkoden till att loopa över flera koordinater. Att någon tredje part implementerat detta till en tjänst där vem som helst kan konvertera en koordinat hjälper oss inte så mycket, utan tillgång till källkod och formel kan vi inte modifiera eller skriva om koden till att köra batch-job.

m v h Henrik

[henca]

Ja, men även att jämföra koordinaterna mellan trafikverkets alla kameror och våra kameror kommer kräva en hel del manuellt arbete. Beakta t ex positione för den här fartkameran som vi fick in häromdagen:
https://www.google.com/maps/place/62%C2%B027’18.0%22N+17%C2%B009’45.4%22E/@62.455,17.1604213,17z/data=!3m1!4b1!4m5!3m4!1s0x0:0x0!8m2!3d62.455!4d17.16261

Genom att bara studera kartan skulle jag gissa att precisionen på den positionsinmätningen har ett fel som i meter överstiger hur många meter bred vägen är. Men även om positionen skulle vara 10, 20 eller rentav 50 meter fel duger den mer än väl till att låta våra GPSer varna för dess position.

Låt oss anta att trafikverkets positioner är mer exakta och att vi efter koordinattransformationer får ett fel på bara några enstaka meter. Hur långt avstånd skall vi då tillåta mellan fartkamerapositioner för att anse att två positioner från olika källor skall vara samma fartkamera? Å ena sidan bör detta avstånd vara större än felet på våra inmätta kameror, å andra sidan bör avståndet vara mindre än minsta avståndet mellan två kameror. Så hur litet avstånd kan det vara mellan två kameror? Jag skulle tro att det avståndet är något större än bredden på den smalaste väg som har två kameror på varsin sida om vägen som mäter motsatt riktning.

Däri består problemet att automatiskt identifiera nya kameror, precisionen på våra positioner är nog sämre än minsta möjliga avstånd mellan två kameror.

m v h Henrik

[henca]

Jag tror mig kunna fixa batch-konvertering. Då det är mera avancerade beräkningar som skall göras duger det nog inte med enkla script i bash och awk. Sannolikt kommer jag i stället förlita mig på php eller i värsta fall C. Men hur stor nytta har vi av konverterade koordinater när övriga data saknas? Framför allt hastigheten vill vi ju ha till varje punkt. Även vägnamn vill vi helst ha för att kunna fortsätta med samma namnkonvention.

m v h Henrik

[henca]

Ja, man kan ju tycka att det borde räcka med en 3×3-matris för att konvertera mellan olika koordinatsystem, men matten blir ofta mer avancerad än så. Då kart-länkarna till hitta.se byggde på RT90 var jag tvungen att konvertera WGS84-koordinaterna i databasen till RT90 för att stödja kartlänkar till hitta.se. Den delen av phpoi blev såhär:

<?php

/******************************************************************
    Copyright 2008 Henrik Carlqvist

    This file is part of phpoi.

    Phpoi is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    Phpoi is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with phpoi.  If not, see <http://www.gnu.org/licenses/>.
******************************************************************/

function RT90($lat, $long)
{
   $sin_lat = sin(deg2rad($lat));
   $cos_lat = cos(deg2rad($lat));
   $sin_lon = sin(deg2rad($long));
   $cos_lon = cos(deg2rad($long));
   $a_wgs = 6378137.0;
   $e2p_wgs = 0.00669437999013;
   
   $Bessel_Ellipsoid_a   = 6377397.1542;
   $Bessel_Ellipsoid_e2p = 0.006674372232;

   $Gd_dR0 = 424.3;
   $Gd_dR1 = -80.5;
   $Gd_dR2 = 613.1;
   $Gd_dC0 = 0.000021315;
   $Gd_dC1 = -0.0000096313;
   $Gd_dC2 = 0.000025136;

  /* Transform WGS84 Lat/Long/H to geo-centric coordinates */
   $N_wgs = $a_wgs/sqrt( 1.0 - $e2p_wgs * $sin_lat*$sin_lat );
   $X_wgs = $N_wgs*$cos_lat*$cos_lon;
   $Y_wgs = $N_wgs*$cos_lat*$sin_lon;
   $Z_wgs = $N_wgs*(1.0 - $e2p_wgs)*$sin_lat;

   /* Bursa-Wolf transformation */
   $X_gd = ( $X_wgs-$Gd_dC2*$Y_wgs+$Gd_dC1*$Z_wgs ) - $Gd_dR0;
   $Y_gd = ( $Y_wgs+$Gd_dC2*$X_wgs-$Gd_dC0*$Z_wgs ) - $Gd_dR1;
   $Z_gd = ( $Z_wgs-$Gd_dC1*$X_wgs+$Gd_dC0*$Y_wgs ) - $Gd_dR2;
  
   /* Transform geo-centric coordinates to Lat/Long/H in Generic Datum */
   $P = hypot($X_gd, $Y_gd);
   $Q = sqrt(1.0 - $Bessel_Ellipsoid_e2p);
   $R = $Bessel_Ellipsoid_a*$Bessel_Ellipsoid_e2p;
   $Eta = atan( $Z_gd/($P*$Q) );
   $lat_gd = atan( ($Z_gd + $R*pow(sin($Eta),3.0)/$Q) /
                   ($P - $R*pow(cos($Eta),3.0) ) );
   $long_gd = atan2( $Y_gd, $X_gd );
   $N_gd = $Bessel_Ellipsoid_a /
       sqrt( 1.0 - $Bessel_Ellipsoid_e2p*pow(sin($lat_gd), 2));
   $h_gd = $P/cos($lat_gd) - $N_gd;

   $a1 = 6.674372206E-3;        /* Series expansion coefficients */
   $a2 = 3.7073149E-5;
   $a3 = 2.569374E-7;
   $a4 = 1.9482E-9;
   $b1 = 8.3522527E-4;
   $b2 = 7.56302E-7;
   $b3 = 1.193E-9;
   $Y0 = 1500000.0;              /* Mean Meridian Y Coordinate */
   $RTH_L0 = 0.2759064963;       /* Mean Meridian in RT90 (rad) */
   $EarthRadius = 6366742.5194;  /* [m]      Mean Earth Radius */

   $sin_lat = sin($lat_gd);
   $sin2_lat = pow($sin_lat, 2);

   $Lat_gc = $lat_gd -
       $sin_lat*cos($lat_gd) *
       ($a1+$sin2_lat*($a2+$sin2_lat*($a3+$sin2_lat*$a4)));
   $Xsi = atan(tan($Lat_gc) / cos($long_gd-$RTH_L0));
   $Eta = atanh(cos($Lat_gc) * sin($long_gd-$RTH_L0));
   
   $res["X"] = $EarthRadius*($Xsi+$b1*sin(2.0*$Xsi)*cosh(2.0*$Eta) +
                             $b2*sin(4.0*$Xsi)*cosh(4.0*$Eta) +
                             $b3*sin(6.0*$Xsi)*cosh(6.0*$Eta));
   $res["Y"] = $Y0 + $EarthRadius*($Eta+$b1*cos(2.0*$Xsi)*sinh(2.0*$Eta) +
                                   $b2*cos(4.0*$Xsi) * sinh(4.0*$Eta) +
                                   $b3*cos(6.0*$Xsi) * sinh(6.0*$Eta));
   return $res;
} /* RT90 */

?>

Det går säkert åt någonting med motsvarande komplexitetsnivå för att konvertera från Sweref99 till WGS84.

m v h Henrik

[henca]

[quote quote=213]Till Henca ok laddade ner detta program 1.0.4. Här måste jag erkänna att min kunskap brister *.tar.gz är ju ett format som jag kopplar ihop med androider och programvara.

Misstänker att man måste ha något mer, men frågan är om jag ska ge mig in i detta [/quote]

Programmet är ingenting för ”vanliga” GPS-användare. Programmet är bara för dem som vill sätta upp en egen web-server där man kan tanka hem GPS-filer.

m v h Henrik

[henca]

[quote quote=212]Om vi utgår från att vägverkets positionsbestämmelser är mest rätt. Gör en programsnutt som jämför koordinaterna matcha och hämta därefter all data som är kopplat till denna position och gör en ny fil…[/quote]

Problemet med att skapa en ny fil är att vägverkets fil inte innehåller alla data som vi vill ha i den nya filen. Att koordinaterna är på fel format är ett mindre problem. Större problem är de data som saknas, det saknas vägnamn och framförallt saknas hastighet.

m v h Henrik

[henca]

[quote quote=201]har lite idéer om hur det eventuellt skulle kunna gå till att få dit vägnamnet och gällande hastighet.
Återkommer när jag har haft tid o leka lite.[/quote]

Om du lyckas lösa detta så vore det en höjdare! Då skulle man kunna jämföra med de kameror som vi redan har och lägga till de kameror som saknas. Tyvärr kommer det ändock nog krävas en del manuellt arbete då många kameror mätts in manuellt och inte matats in med exakt samma koordinater. På något vis gäller det att identifiera vilka kameror som är nya och vilka som bara skiljer lite i position.

m v h Henrik

[Författare]

Inlägg