WipNav-WipCom: ricerca formato autovelox

[mira308sw]

ancora non lo so

[mira308sw]

proseguo la mia analisi del file LWD, ma sono sempre lontano dalla soluzione

sto analizzando i dati relativi alle isole (Canarie, Baleari, Malta, Sardegna, Islanda) in quanto sono piu' facili da localizzare sia nel file LWD sia nel file per Google Earth.
per quel che riguarda le isole Canarie sono ormai convinto che il blocco di dati di 125 byte che vi avevo elencato nei messaggi precedenti vada cosi' interpretato:


B7 00 7D 58 BC D5 = radar 1 -17.85222, 28.65361, 50
E0 90 C9 1D 87 75 BA DA = radar 2 -16,42722, 28,47269, 100
23 D9 D2 FD 83 DC = radar 3 -16,37444, 28,49108, 120
A4 9A CC 1B FF D6 = radar 4 -16,30750, 28,47933, 60
84 CB 65 FF D6 = radar 5 -16,30722, 28,47614, 60
DB E4 98 5E C6 97 DC = radar 6 -16,36278, 28,38219, 120
FD C3 F6 2E B8 C4 DC = radar 7 -16,58306, 28,05800, 120
10 42 28 5D 81 92 D8 = radar 8 -15,43764, 28,06375, 80
A1 EC D5 24 FF D8 = radar 9 -15,42369, 28,04731, 80
31 23 6B F4 DD D8 = radar 10 -15,41611, 28,07944, 80
10 23 6D 9B D8 = radar 11 -15,41611, 28,08667, 80
10 19 FC D7 D8 = radar 12 -15,41222, 28,09583, 80
38 85 A5 5F AA D8 = radar 13 -15,41806, 28,11511, 80
29 45 95 5A E4 D8 = radar 14 -15,42722, 28,12514, 80
E3 71 31 75 CA DA DA = radar 15 -15,38744, 27,95053, 100
ED B1 B9 97 2F C1 DA = radar 16 -15,55000, 27,77500, 100
00 D8 F9 6A 2F C0 A3 D9 = radar 17 -13,86833, 28,47483, 90
1D 0A AD DF D9 = radar 18 -13,86889, 28,48100, 90
76 8D 99 3A FB 82 DA = radar 19 -13,57833, 28,95789, 100
8E 32 BE DA = radar 20 -13,57833, 28,95806, 100
00

le ultime tre colonne sono Longitudine, Latitudine, Velocita' presi dal file per Google Earth.
la divisione dei dati in quel modo e' del tutto empirica e l'ho fatta solo perche' ho notato che ci sono ripetuti tanti byte che iniziano con 'D' e il numero seguente e' facilmente riconducibile alla velocita' del radar / 10.
quindi credo che quella su elencata sia la giusta suddivisione dei dati per ogni radar presente.
se e' vero che l'ultimo byte esprime la velocita' del radar, i byte precedenti devono contenere le coordinate (secondo me in metri X,Y) del radar, ma chi sa' in quale formato/compressione.

la stessa divisione l'ho fatta anche per i radar delle isole Baleari, e il risultato e' corretto anche per loro.

per l' Islanda invece si nota una particolarita':
8E 30 D2 68 FF D9
85 24 2D D7 39 FF D9
47 B2 C4 54 8A FF D1
41 1D 56 17 FF D1
52 B0 55 FF D1
DA 10 69 0B FF D1
A2 C5 02 9A FF D1
22 6A 14 FA FF D1
BD 9C 35 B7 08 FF D1 
E7 61 89 2C DF FF D9
B3 8C 9F 69 FF D9
il byte precedente a quello con la velocita' e' sempre FF, quindi questo byte non fa parte delle coordinate X,Y ma "forse" racchiude la direzione del radar, che in questo caso (cioe' il valore FF) non e' specificato.

dal momento che la lunghezza in byte di ogni singolo radar risulta in questo modo variabile, secondo me anche il primo byte non fa parte delle coordinate X,Y ma deve indicare in qualche modo quanti byte/bit sono lunghi i dati del radar.
se queste supposizioni sono corrette c'e' un gran bel mistero su come sono espresse le coordinate X,Y dei radar in quanto ci sono radar che usano solo 1 byte (radar 20) e radar che ne usano 5 (radar 17)

per ora e' tutto

[navtom]

Bravo pour ce travail

[macb]

Boujour,

Désolé je ne parle pas italien. Alors j'essaye en français et en anglais.

D'abord bravo pour ce travail. First thanks for the good work !

Une petite contribution : A small contribution :
lorsqu'il y a une ligne courte, les coordonnées du radar sont très proches de celle du précédent. 3 exemples :
When the line is short, the coordinates are close to the previous one. 3 examples :
A4 9A CC 1B FF D6 = radar 4 -16,30750, 28,47933, 60
84 CB 65 FF D6 = radar 5 -16,30722, 28,47614, 60

00 D8 F9 6A 2F C0 A3 D9 = radar 17 -13,86833, 28,47483, 90
1D 0A AD DF D9 = radar 18 -13,86889, 28,48100, 90

76 8D 99 3A FB 82 DA = radar 19 -13,57833, 28,95789, 100
8E 32 BE DA = radar 20 -13,57833, 28,95806, 100

Donc les premiers octets représentent peut-être le déplacement par rapport au précédent et pas des coordonnées absolues.
So the first bytes may be the movement from the previous radar and not absolute coordinates.

Ce peut être un vecteur (x,y) ou un angle et une longueur. Le dernier cas est intéressant car le déplacement est très faible et il y a très peu d'octets. Voir aussi la série de 3 assez proches :
It can be a vector (x,y) or an angle and a length. The last case is interesting because the movement is very small and there are very few bytes. Also interesting the series of 3 relatively close :
31 23 6B F4 DD D8 = radar 10 -15,41611, 28,07944, 80
10 23 6D 9B D8 = radar 11 -15,41611, 28,08667, 80
10 19 FC D7 D8 = radar 12 -15,41222, 28,09583, 80


Hope this helps

[mira308sw]

ho deciso di scrivere in un documento (che trovate in allegato) quello che sono riuscito a capire fino ad ora del formato del file .LWD.
e' in inglese, non vogliatemene, ma cosi' si puo' allargare il ventaglio dei probabili collaboratori.

ed ora alcune novita':
in alcuni database PSA c'e' un radar (si tratta di un semaforo a Melilla, citta spagnola in Marocco) che risulta essere l'unico della sua area, e quindi si presta a facili confronti.
le sue coordinate sono:
Lon,Lat: -2.94214, 35.28556
non so ancora di preciso come convertire queste coordinate in metri, ma un programmino free mi da questi valori
X,Y: -267562, 3908679
come dicevo questo radar risulta l'unico della sua area in 4 database, di cui di seguito vi elenco le coordinate dell'area in cui e' presente:

X Y W H x y w h coord dir+speed
-1285360 3117554 1700153 1529409 -222765 3882258 212519 191176 AA 75 B5 2B D8 00 D5 71
-1285360 3117554 1934855 1621394 -317933 3117554 967427 810697 E9 30 9E 75 40 80 F8 D5 71
-1285360 3117554 1949294 1700666 -310713 3117554 974647 850333 E9 30 82 41 40 80 F8 D5 71
-1285360 3086111 1956525 1717982 -307098 3086111 978262 858991 E9 30 74 22 BB 80 CB D5 71

la prima area e' piu' piccola, mentre le altre 3 sono quasi identiche, e come gia' sapevamo, i byte che descrivono il radar sono 1 in meno nell'area piu' piccola.
i numeri nelle colonne X,Y,W,H (le prime 4) si riferiscono all'area principale del blocco
mentre i numeri nelle altre 4 colonne x,y,w,h si riferiscono alle coordinate della sub-area dopo le divisioni eseguite secondo la tabella che descrive le suddivisioni del blocco principale.
dopo una breve analisi dei valori delle righe 2 3 4, si puo' capire che i byte segnati in verde contengono la coordinata X del radar, mentre quelli blu la coordinata Y, e sono realtivi alla sub-area.
quindi abbiamo:
X hex Y hex = X dec Y dec
9E75 40F8 = 40565 16632
8241 40F8 = 33345 16632
7422 BBCB = 29730 48075

aggiungendo le coordinate della sub-area:
-317933+40565=-277368 3117554+16632=3134186
-310713+33345=-277368 3117554+16632=3134186
-307098+29730=-277368 3086111+48075=3134186

oltre ad essere lontani dalle coordinate X,Y del radar Melilla, siamo anche oltre alle coordinate x,y della prima riga che sono -222765,3882258 (angolo inferiore sinistro)
ho quindi aggiunto dei numeri a caso prima delle 4 cifre esadecimali individuate prima, ottenendo:
X hex Y hex = X dec Y dec
19E75 C40F8 = 106101 803064
18241 C40F8 = 98881 803064
17422 CBBCB = 95266 834507
che sommati alle coordinate di origine delle rispettive sub-area danno: -211832,3920618

ora ragionando al contrario, togliamo da queste coordinate l'origine della prima sub-area:
-211832-(-222765) 3920618-3882258 = 10933,38360
che in esadecimale si scrivono: 02AB5,095D8
di cui solo alcuni byte sono presenti nei dati del primo radar, quelli verdi e blu.

se tutto questo ragionamento sia giusto non lo so, pero' continuo su questa strada per cercare di scoprire il significato dei byte che mancano.

troppo complicato ?

[baronecausio]

Ciao Mira!
Innanzitutto ti faccio i complimenti ed i ringraziamenti per ciò che stai facendo.
Io non sono un informatico, nè tantomeno un matematico, ma sono un geografo (almeno sulla carta), per cui ti pongo una domanda magari banale.
Non è che bisogna capire le coordinate a quale sistema di proiezioni fanno riferimento?
Nel senso che, magari, individuando il sistema di proiezione (ce ne sono decine) riusciamo ad individuare l'area a cui fanno riferimento le coordinate di cui sopra.

P.S.: magari ho detto solo castronerie... [-o<

[mira308sw]

Non è che bisogna capire le coordinate a quale sistema di proiezioni fanno riferimento?
Nel senso che, magari, individuando il sistema di proiezione (ce ne sono decine) riusciamo ad individuare l'area a cui fanno riferimento le coordinate di cui sopra.

P.S.: magari ho detto solo castronerie... [-o<

hai detto quasi certamente una cosa giusta.
ci metto il quasi perche' non so esattamente cosa sia un "sistema di proiezione"
se vuoi provare a fare un po di conti anche tu, sei il benvenuto!

come ho scritto, ho usato del codice freeware per convertire in UTM, e le coordinate del radar di Melilla
Lon,Lat: -2.94214, 35.28556

mi vengono :
X,Y: -267562, 3908679

invece, se le mie supposizioni sono corrette, il navigatore lo pone a :
X,Y: -211832, 3920618

sapresti identificare quale sistema di proiezione usa il navigatore ?

[baronecausio]

Ciao Mira.
Secondo me potrebbero non fare parte della proiezione UTM, ma alla WGS84, che è un sistema convenzionale terrestre a livello mondiale...
Questa è solo un'ipotesi, perchè comunque la proiezione Universale Trasversa di Mercatore è pressochè usata solo a livello nazionale, e ad oggi neanche più di tanto dato che è stata sostituita dalla WGS84. Ne risulta che, facendo riferimento ad altro datum, le coordinate siano sfalsate rispetto all'altra proiezione.

[mira308sw]

mi puoi indicare un link per documentarmi meglio sul WGS84 ?
e se lo conosci anche un link ad un programma che ne implementi l'algoritmo, in qualsiasi linguaggio.

[baronecausio]

A parte Wikipedia http://en.wikipedia.org/wiki/World_Geodetic_System" onclick="window.open(this.href);return false;
leggi http://www.lordalcol.com/gis/Trasformaz ... 20Maps.pdf" onclick="window.open(this.href);return false; in cui c'è scritto anche a livello informatico come avviene la trasformazione, e il relativo programma di conversione http://www.lordalcol.com/gis/" onclick="window.open(this.href);return false;
Poi http://webgis.regione.sardegna.it/raswe ... .zip#start" onclick="window.open(this.href);return false; oppure http://webgis.regione.sardegna.it/raswe ... uale=punto" onclick="window.open(this.href);return false;
Se hai voglia di leggere in inglese: http://www.epsg.org/guides/docs/G7-2.pdf" onclick="window.open(this.href);return false;
Poi c'è un interessante copia-incolla qui, 2° post: http://www.mtb-forum.it/community/forum ... p?t=155629" onclick="window.open(this.href);return false;
E qui da pag. 15 è spiegato bene: http://users.unimi.it/br1/geodesia.pdf" onclick="window.open(this.href);return false;
In altra maniera, al momento, non so aiutarti...mi spiace... [-o<

[mira308sw]

A parte Wikipedia http://en.wikipedia.org/wiki/World_Geodetic_System" onclick="window.open(this.href);return false;
leggi http://www.lordalcol.com/gis/Trasformaz ... 20Maps.pdf" onclick="window.open(this.href);return false; in cui c'è scritto anche a livello informatico come avviene la trasformazione, e il relativo programma di conversione http://www.lordalcol.com/gis/" onclick="window.open(this.href);return false;
Poi http://webgis.regione.sardegna.it/raswe ... .zip#start" onclick="window.open(this.href);return false; oppure http://webgis.regione.sardegna.it/raswe ... uale=punto" onclick="window.open(this.href);return false;
Se hai voglia di leggere in inglese: http://www.epsg.org/guides/docs/G7-2.pdf" onclick="window.open(this.href);return false;
Poi c'è un interessante copia-incolla qui, 2° post: http://www.mtb-forum.it/community/forum ... p?t=155629" onclick="window.open(this.href);return false;
E qui da pag. 15 è spiegato bene: http://users.unimi.it/br1/geodesia.pdf" onclick="window.open(this.href);return false;
In altra maniera, al momento, non so aiutarti...mi spiace... [-o<

allora mi metto a studiare
per il momento ho provato solo questo: http://www.lordalcol.com/gis/" onclick="window.open(this.href);return false;
ma non mi da i valori che vorrei.

comunque potresti aiutarmi facendo anche tu i conti per il radar di Melilla, ogni volta che trovi un sito che fa conversioni da Lon,Lat in metri, prova a metterci le coordinate del radar di Melilla e poi vedi se il risultato si avvicina ai numeri che ho scritto nel post precedente.

[baronecausio]

Ok! Per quanto sia imbranato coi calcoli...cercherò di darti una mano per quanto mi è possibile!

[baronecausio]

Ciao Mira!
Non mi sono dimenticato di te in questo periodo...
Però purtroppo ho provato anche io vari programmini free, che permettono di trasformare le coordinate, e l'unico che mi da un valore approssimato al tuo si trova qui: http://www.cartografia.regione.lombardia.it/convxy/" onclick="window.open(this.href);return false;
Ma credo siamo lontani dal radar di melilla.
Ho bisogno di capire alcune cose: il valore che hai ottenuto delle coordinate lon-lat del radar, è quello reale o ottenuto da calcoli? E poi il valore che hai di queste coordinate, è già UTM?
Sta cosa mi ingrippa il cervello...

[mira308sw]

Ho bisogno di capire alcune cose: il valore che hai ottenuto delle coordinate lon-lat del radar, è quello reale o ottenuto da calcoli?

lon,lat sono quelle prese dal database radar di SCDB, quindi reali.

E poi il valore che hai di queste coordinate, è già UTM?

le quote X,Y sono invece derivate dal mio ragionamento, non ho idea se siano UTM.

[ghis91]

bonjour,
y a t'il du nouveau pour le décodage du format LWD ?
cordialement