Kaardisoft ja koordinaatide teisendamine

[felch]

Päevast!
Tean et see pole just otseselt elektroonikafoorumi teema kuis siinvaria rubriigis ehk väga ei sega ja ehk keegi teab vastust.
Nimelt on selline tore vabavaraline GIS-soft nagu MapWindow http://www.mapwindow.org/
Saab teha oma programme, kasutades active-x kaardikomponenti. Just see, mida vaja. Aga...GPS väljastab koordinaate kraadides ja minutites. Kaardile neid märkida ei saa. Vähemasti mitte otseselt. On olemas funktsioonid, mis projitseerivad punkti koordinaatidest ekraani pixeliteks ja vastupidi. Seda soovitataksegi kasutada. Paraku annab ümberprojekteerimine tulemuseks ristkoordinaadid (Lambert-EST süsteemis), mitte aga kraade ja minuteid (või kraade kümnendkujul).
Kuidas konverteerida Lambert-EST "tavalisteks" koordinaatideks? Kindlasti on see elementaarne kuid mitte minusugusele GIS-i alal algajale
Tähh!

[fre]

Vs: Kaardisoft ja koordinaatide teisendamine

Kas meil levinud ja täielikult kümnendsüsteemis töötav
L-EST 92 ei sobi kasutada? Kas see soft sellist toetab?
Selles viimane (üheliste) koht peaks minuteada täpselt
1 meetrit tähendama. Arvan ei GPS peaks seda ka
väljastama hakkama (näitama küll hakkab) kui ta on
sellele seadistatud.
Seadistamise juhend on mul juba ammu, aga nüüd
leidsin selle ka siit viimaselt leheküljelt:

http://www.mnt.ee/atp/failid/02_Aruanne_II.pdf

Geopeituse tegelased ehk teavad täpsemalt...

[felch]

Vs: Kaardisoft ja koordinaatide teisendamine

Tänapäeval peaks kehtima L-EST 97 mis ei pidavat 92-st eriti erinema. Minu ülesanne on teisendada laius/pikkuskraadid just sellesse süsteemi (ja vastupidi). Digikaardilt saan vastava funktsiooniga kätte nimelt L-EST koordinaadid (need X ja Y). GPS on Falcomi JP-7, pisike aparaati ehitatav moodul. Tema täielikku käsustikku ei oma, Falcomi leht on paras peitusemäng... Ehk siis peab ise teisendama. Paraku teisenduse tegemiseks valemit ei ole juhtunud leidma. Selles ongi kogu probleem.

[pupsik46]

Vs: Kaardisoft ja koordinaatide teisendamine

Olen teisaldanud kraadid-minutid L-EST'i - assembleris
Otsin homme valemid üles.

[felch]

Vs: Kaardisoft ja koordinaatide teisendamine

Antud juhul kirjutan C-s (võehh rõvedust, hakkan vist vanaks jääma... ). Aga algoritmi saab Sinu koodist ikkagi. Suurimad tänud!

PS. Kui teisendus on piisavalt lihtne et teda asmis realiseerida, tekib kohe tahtmine see otse spektromeetrisse toppida. PICile...

[fre]

Vs: Kaardisoft ja koordinaatide teisendamine

Siin on palju valemeid, kas neist aga abi on... :

http://www.posc.org/Epicentre.2_2/DataModel/ExamplesofUsage/eu_cs35.html

Google andis ka võrdlevat lugemist (kui huvitab):

http://www.geo.ut.ee/~raivo/ESTCOORD.HTML

[Matjuhha]

Vs: Kaardisoft ja koordinaatide teisendamine

Mingi katkend...
Teisendab L-Est alusellipsoidile (ETRS miski) ja tagasi.
Stiili üle paluks mitte joriseda

#define kesk_mer 24
#define alf 0.85417585805
#define quu 4020205.479
#define kaa 11473298.03838
#define false_n 6375000
#define false_e 500000
#define pi 3.14159265358979323846264
//#define AAA 0.003356551
//#define BEE 0.000006571
//#define CEE 0.000000017
//#define DEE 0
#define EEE sqrt(0.00669438)
#define AAA pow(EEE,2.0)/2.0+5.0*pow(EEE,4.0)/24.0+pow(EEE,6.0)/12.0+13.0*pow(EEE,8.0)/360.0
#define BEE 7.0*pow(EEE,4.0)/48.0+29.0*pow(EEE,6.0)/240.0+811.0*pow(EEE,8.0)/11520.0
#define CEE 7.0*pow(EEE,6.0)/120.0+81.0*pow(EEE,8.0)/1120.0
#define DEE 4279.0*pow(EEE,8.0)/161280.0

void Lest_Geo(double * pik, double * lai, double n, double e)
{
double a=AAA;
double b=BEE;
double c=CEE;
double d=DEE;
*pik = atan((e - false_e) / (quu - (n - false_n)));
*pik = *pik / alf;
*pik = *pik * 180 / pi + kesk_mer;
*lai = pow(e - false_e, 2.0);
*lai = *lai + pow(quu - (n - false_n), 2.0);
*lai = pow(*lai, 0.5) / kaa;
*lai = pow(*lai, 1.0 / alf);
*lai = pi / 2.0 - 2.0 * atan(*lai);
*lai = *lai + a * sin(2.0 * *lai) +
b * sin(4.0 * *lai) +
c * sin(6.0 * *lai) +
d * sin(8.0 * *lai);
*lai = *lai * 180 / pi;
}
void Geo_Lest(double * n, double * e, double pik, double lai)
{
double p,b;
pik = (pik - kesk_mer) / 180.0 * pi;
lai = lai / 180.0 * pi;
b = pik * alf;
p = asin(EEE * sin(lai));
p = pow(tan(pi / 4.0 + p / 2.0), EEE);
p = tan(pi / 4.0 + lai /2.0) / p;
p = kaa / pow(p, alf);
*n = quu - p * cos(b) + false_n;
*e = p * sin(b) + false_e;
}

[felch]

Vs: Kaardisoft ja koordinaatide teisendamine

Tänud kaasamõtlejatele, kood töötab kenasti. Aga kui juba kaardiga jamamisest juttu oli - kuidas teie lahendaksite olukorra, kui mõõdetavad koordinaadid kõiguvad ühest kaardiruudust teise? Laeb 2 kaarti korraga? Või annab neid kuidagi mergeda?

[Timo]

Vs: Kaardisoft ja koordinaatide teisendamine

Lisan siia mõned lingid mu enda lingikogust, mis mind aidanud koordinaatide konvertimise vajadusel:
http://www.fizzymagic.net/Geocaching...c/GeoCalc.html
http://boulter.com/gps/
http://www.rvdumps.com/gpscalc.htm
http://www.jeeep.com/details/coord/