gpsgek

Waarom niet? Ik kom in Delft uit

Ik heb niet zozeer iets aan de berekening veranderd, maar wel de invoer in het script gezet en de uitvoer wordt nu op de gebruikelijke manier GG MM.MMM weergegeven....

Prof Lupardi, ik ben ermee aan het stoeien gegaan maar loop een klein beetje vast. Voor de laatste transformatie is een hoogte nodig. Ik heb daarvoor 0 ingevoerd, maar krijg aan het einde van de berekening wèl een hoogte in het WGS formaat Bovendien klopt de positie in WGS ook niet helemaal. Ik heb de code zo aangepast: % rd2bessel % % Omrekenen geo coordinaten (tov Bessel) naar RD-grid coordinaten % Formules van Geodetic Datum of Australia website % 1 december 2001 %function geo = rd2bessel(rd); %INVOER rd(1)=84000; rd(2)=448000; RDx = rd(1); RDy = rd(2); deg = pi/180; % omrekenfactor naar radialen % Bessel 1841 datum: a = 6377397.155; f = 1/299.1528128; % RD grid parameters: lab0 = 5.3876333*deg; % origin Efalse = 155000; % false easting K0 = 0.9999079; % scale factor Nfalse = -5316592; % false northing % omrekenen: e2 = f*(2-f); e = sqrt(e2); a0 = 1 - 1/4*e^2 - 3/64*e^4 - 5/256*e^6; a2 = 3/8*(e^2 + 1/4*e^4 + 15/128*e^6); a4 = 15/256*(e^4 + 3/4*e^6); a6 = 35/3072*e^6; m = (RDy-Nfalse)/K0; n = f/(2-f); G = a*(1-n)*(1-n^2)*(1+9/4*n^2 + 225/64*n^4)*pi/180; sig = m/G*pi/180; % phic = phi-accent. phic = sig +((3*n/2)-(27*n^3/32))*sin(2*sig) + ((21*n^2/16)-(55*n^4/32))*sin(4*sig)... +(151*n^3/96)*sin(6*sig) + (1097*n^4/512)*sin(8*sig); rho = a*(1-e^2)/(1-e^2*(sin(phic))^2)^(3/2); nu = a/sqrt(1-e^2*(sin(phic))^2); psi = nu/rho; % Grid to Geographical: t = tan(phic); Ec = RDx - Efalse; x = Ec/(K0*nu); term1 = (t/(K0*rho))*(x*Ec/2); term2 = (t/(K0*rho))*(Ec*x^3/24)*(-4*psi^2+9*psi*(1-t^2)+12*t^2); term3 = (t/(K0*rho))*((Ec*x^5)/720)*(8*psi^4*(11-24*t^2)-12*psi^3*(21-71*t^2)+15*(psi^2*15-98*t^2+15*t^4)+180*psi*(5*t^2-3*t^4)+360*t^4); term4 = (t/(K0*rho))*(Ec*x^7/40320)*(1385+3633*t^2+4095*t^4+1575*t^6); phi = phic - term1 + term2 - term3 + term4; term1 = x*sec(phi); term2 = (x^3/6)*sec(phi)*(psi+2*t^2); term3 = (x^5/120)*sec(phi)*(-4*psi^3*(1-6*t^2)+psi^2*(9-68*t^2)+72*psi*t^2+24*t^4); term4 = (x^7/5040)*sec(phi)*(61+662*t^2+1320*t^4+720*t^6); omega = term1-term2+term3-term4; lab = lab0 + omega; geo(1) = phi/deg; % lengte [graden] geo(2) = lab/deg; % breedte [graden] % %------------------------------------------------------------------- % molodensky datum omzetting van Bessel1841 naar WGS84 % Formules volgens website Peter Dana. % Gecontroleerd met rekenvoorbeeld en GPS. % In/uitvoer in graden en meters % verbeterd nav check met EPS: fout in berekening Rm gehaald. % 8 feb 2001 % 1 dec 2001: hoogte niet *deg! %function [to_geo]=bessel2wgs(from_geo); %INVOER from_geo(1)=geo(1); from_geo(2)=geo(2); from_geo(3)=0; deg = pi/180; % omrekenfactor naar radialen from_phi = from_geo(1)*deg; % from_hoogte boven evenaar from_lab = from_geo(2)*deg; % from_oosterlengte from_h = from_geo(3); % from_hoogte boven zeespiegel % parameters voor WGS 84: to_a = 6378137.000; % semi-major axis of ellipsoid to_f = 1/298.257223563; % 1/flattening % parameters voor Bessel 1841: from_a = 6377397.155; from_f = 1/299.1528128; % correctiefactoren uitrekenen: da = to_a - from_a; df = to_f - from_f; bda = 1-from_f; from_es = 2*from_f - from_f^2; dx = 593; dy = 26; dz = 478; % Nu coordinaten omrekenen van Bessel 1841 naar WGS84 sphi = sin(from_phi); cphi = cos(from_phi); slab = sin(from_lab); clab = cos(from_lab); Rn = from_a/sqrt(1.0 - from_es*(sin(from_phi))^2); % oud/fout: Rm = from_a*(1-from_es)/(1.0-from_es*(sin(from_phi)^2)^(3/2)); Rm = from_a*(1-from_es)/((1.0-from_es*(sin(from_phi)^2))^(3/2)); term1 = (((-dx*sphi*clab-dy*sphi*slab) + dz*cphi) + da*Rn*from_es*sphi*cphi/from_a); term2 = df*(Rm/bda + Rn*bda)*sphi*cphi; dphi = (term1 + term2)/(Rm+from_h); dlab = (-dx*slab + dy*clab)/((Rn+from_h)*cphi); dh = dx*cphi*clab + dy*cphi*slab + dz*sphi - da*from_a/Rn + df*bda*Rn*sphi*sphi; to_phi = from_phi+dphi; to_lab = from_lab + dlab; to_h = from_h + dh; to_phi = to_phi/deg; to_lab = to_lab/deg; to_h = from_h + dh; LatitudeGraden = floor(to_phi) LongitudeGraden = floor(to_lab) LatitudeMinuten = (to_phi - LatitudeGraden) * 60.0 LongitudeMinuten = (to_lab - LongitudeGraden) * 60.0 %to_geo = [to_phi,to_lab,to_h] % %-------------------------------------------

Das helemaal mooi Gpsgek is enthousiast Matlab gebruiker en heeft bovenstaande code ook daarin toegepast

Zelf heb ik het niet nodig, maar wellicht is het handig het voor de volledigheid aan het draadje toe te voegen. Voor mensen die later iets soortgelijks zoeken!

Bedankt Bruine Beren, was precies wat ik zocht! Ja hoor, dat gaat prima... Antwoord gegeven: Draadje kan gesloten worden. Ruud Genoeg te doen voor Kruimel... Desondanks was het draadje zeker nuttig

Bedankt Bruine Beren, was precies wat ik zocht! Ja hoor, dat gaat prima...

Niet zozeer een newbie-vraag, geloof ik, maar toch. Kruimel, als je een betere plaats weet mag je het topic verplaatsen... Wie kan mij helpen aan een - relatief eenvoudige - formule waarmee RD wordt omgerekend naar WGS84 of andersom? Geen tooltjes dus om dit te doen, die heb ik zat. Ik wil graag de formule implementeren in een rekenprogramma (nee, geen excel...).

Gebruik onder de opties knop eens de instelling standaard, dat zal het ongetwijfeld zijn. Groetjes Jan-Willem. Thanks, dat was het. Soms kan een oplossing toch zo eenvoudig zijn... Misschien toch per ongeluk ergens op geklikt. Bepaalde PC's hebben daar last van. Op de TU is de site ook altijd in eerste instantie met een boomstructuur...

Ik gebruikte ze al voor augustus op de PDA. Toen vond ik al dat het handiger was als OziCE. Toen ik eind augustus de 60CSx heb aangeschaft ben ik de kaarten ook gelijk op de GPS gaan gebruiken. Eigenlijk is er maar een woord te zeggen over de kaarten: FANTASTISCH!

Volgens mij is de layout al even zo... Overigens wil ik ook wel cash vinden...

Gelukkig hoef jij het blad niet in je eentje vol te schrijven Want, als honderd man ieder voor zich maar 1 interessant onderwerp kunnen verzinnen, zou het bij elkaar nog wel een fors blad kunnen worden. Ik ga in ieder geval alvast op zoek naar een verstevigd kruipruimteluik onder mijn digitale brievenbus. Das ook weer zo... Nou, ik ben ook benieuwd...

Voordeel voor de makers is dat er dan ook weinig inhoud aan gegeven hoeft te worden! Ik vraag me altijd af wat je in zo'n blad zou willen schrijven. Vergelijkingen van GPSsen vinden we immers op GC.nl al, caches "bespreken" lijkt me niet leuk voor de legger, in munten heeft lang niet iedere cacher interesse, voor cache-leg-technieken bestaan al beschrijvingen en roddelbladen zijn er al genoeg dus dat zal het ook niet zijn... EDIT: het enige wat ik zou kunnen verzinnen is een artikeltje over "Hoe schrijf ik een goede log." Daaraan kan - wat mij betreft - in 98 van de 100 found-logs behoorlijk wat verbeteren...

Zucht.... Ipods...

Dit was een klein jaartje geleden bij de Rottemeren. Ik liep/fietste mijn derde cache. En wat voor een... Het vroor maar 10 graden... Uiteindelijk crashte ik met een lekke fietsband ongeveer 12 km van huis zonder plakspullen... Oei wat had ik er de pest in.

Gelukkig liegen ze niet, hebben ze wel gelijk OziExplorer en BeelineGPS hebben mijn PPC-cache-tijd veraangenaamd. Je zou ook Garmin Que kunnen overwegen.

Is die voorste straal bezopen ofzo?

OEMPF

Vooruit dan, voor een dagje

Dan zie je dat je helemaal niets ziet . Even HTML inschakelen, dat werkt vast beter Ook even eventuele tags code of /code weghalen

Bedankt HH!

Toen ik dit bericht in de bus naar huis las, viel er een diepe stilte. Bekenden en familie wil ik heel veel sterkte toewensen. Reintje, we'll never forget you...

Ik kan me slecht voorstellen dat ze er ècht gelukkig mee zijn...

In het kader van ècht gek: TeamAlamo. De eerste 18 tonner...

Hmmm, in dit jaargetijde ga ik toch meer voor de Sinterklaaskaarten