Prof. Y. Lupardi

Helaas, helaas. De cache Hold Your Horses is niet meer sinds vandaag. Er kwam bericht van SBB dat door werkzaamheden de cache niet meer via het (paarden) pad te bereiken is met het verzoek hem te verwijderen. Na 7 jaar is deze dus exit. Jammer. De cache heeft gedurende zijn leven bijna 450 bezoekers gehad.

Ik ben zo'n oud persoon dat ik het WWW heb zien komen. Eerst was er het fido-netwerk. En ik was bij xs4all die als eerste (in die prille begintijd nog met 8 analoge telefoonmodems) internet aanbood voor niet-universitairen. Daarvoor zat ik in de schrijverij/uitgeverij en net zoals het internet, heb je dan te maken met lieden waar geen land mee te bezeilen valt en een paar die ronduit gek zijn en gevaarlijk. Nog altijd hangt mijn huis vol met rookmelders en brandblussers. En dat is geen overbodige luxe gebleken (toendertijd). Dus toen geocaching kwam bedacht ik: wat ga ik doen? En toen heb ik mij de naam prof. Y. Lupardi aangemeten. Want ik wil altijd graag uitleggen en onderwijzen. Aangezien ik nogal eens wat moest uitvogelen heb ik als tweede naam daarachter Piet Vogelaar zitten. Piet heeft een kaasboerderij. Piet zijn kaasboerderij Dus ja, ik heb moedwillig een barriere opgeworpen tussen mij en de internet activiteiten. Vrijheid van meningsuiting. Ja maar hoe zit het met veiligheid van meningsuiting?

Als je gebruik maakt van de digitale topografische kaarten (rasterkaarten) 1:25k van onze topodienst in Ozi Explorer dan zit je max. zo'n 10 meter mis is mijn ervaring. En een GPS'r in een stad wekt door slecht zicht op satellieten ook niet zo nauwkeurig als je zou wensen. Ook met uitmeten op papieren kaarten 1:25k kom je een heel eind.

Toen geocaching begon in ons land, kwamen er wat deelnemers uit de EO-hoek die caches een goede manier vonden om evangelisatie te doen. De weg naar een cache gelardeerd met bijbelse (Nieuwe Testament) vragen en de cache vol met evangelisatie literatuur. Er was toen een discussie of dergelijke caches zouden vallen onder 'commercieel'. Ik geloof dat het initiatief is overgewaaid. Maar met zoveel caches in ons landje kan het best zijn dat er evangelisatie- dan wel missie-caches zijn. Mij staat ook een astrologie-cache bij: je moest dan te weten komen of Pluto (of een andere) momenteel in een of ander huis stond dan wel in een sterrebeeld oid. Ook zou er een aardstralen-cache zijn, te vinden via ley-lines of whatever.

Ja dat is bekend. En graden staan garant voor verwarring. Wat je trouwens ook van minuten kan zeggen. Hoeveel minuten moet ik wandelen om 5 minuten in oostelijke richting te gaan?

De temperatuur die de meeste GPSR's laten zien (mijn GPS12 b.v.) is die van het kristal dat de oscillatie-frequentie maakt. De kristallen die gebruikt worden zijn speciaal zo gemaakt dat ze op de korte termijn heel precies lopen. Dus niet 'hikken' bijvoorbeeld of een beetje veel faze-verschuiving hebben. De GPSr haalt de juiste tijd (absoluut) van de klokken van de GPS satellieten. Dus houdt die bij hoeveel het kristal mis loopt. Maar dat mislopen is temperatuurafhankelijk. Jouw GPSr bouwt een tabel op van miswijzing per temperatuur. Op die manier geeft de kristaloscillator een heel kleine bijdrage in de positie-fout.

Het heeft heel lang geduurd (een maand mbv. Murphy) maar ze zijn nu ontvangen hier: de nieuwe Auto-TB's. Vandaag doe ik alle achterstallige bestellingen op de post.

Er was snel actie. The link is hersteld en de page heeft nu een button 'Activate Geocoin" Maar ja, een Travel Patch is een geocoin volgens de logica. Wel een beetje verwarrend. Ik heb gevraagd om het wat duidelijker te maken. http://www.cachezone.de/jump/ dit is de link waar ik het over heb.

Op het moment is de link via www.cachezone.com voor het activeren van de Travel Patches van de Cache Zone Shop in Duitsland kapot. Met wat zoeken zullen de meeste er wel uit komen maar hier is dan de directe deep-link: http://www.cachezone.de/activation/ Excuses voor het ongemak. Ik hoop dat het gauw wordt opgelost.

Hier duikt dan de gewenste professor op. Met een Mededeling voor het Nut van het Algemeen. Lees dit: goodiesshop nieuws Er is weinig aan de situatie veranderd sinds mijn laatste actie in deze. En als een firma zijn webpage laat verslonzen dan is dat geen goed teken.

Dank voor het LOF De shop begon klein, dus ik begon met weinig investering uit eigen zak. Beetje bij beetje kon ik vanwege winst gaan investeren in voorraad. Het aantal klanten groeide dus er kwam steeds meer geld onder de mensen te zitten en ook de voorraad moest groter worden. Die financiering vanuit de winst is tot op heden voldoende gebleken en zodoende kan ik op rekening blijven verkopen. No problemo.

Er waren grote moeilijkheden om de reflectiefolie te krijgen voor de fire-needles. Maar nu heb ik weer de gele folie en de fire-needles worden met monnikenwerk in elkaar gezet. Er zijn inmiddels voldoende aantallen om iedereen blij te maken.

Er zijn weer wat nano-caches hier aangekomen. Ik heb ze in zwart en in groen, maar ik moet heel goed kijken om verschil te zien (en ik ben niet groen/zwart kleurenblind) Zoals gezegd, ze gaan met een reserveset logrolletjes de deur uit. Liefhebbers: gewoon even de shoppage op en in het opmerkingenveld de wensen noteren.

Dat wordt even rekenen. De nano's kosten in de USA 4,50 dollar. Groundspeak De logrolletjes uit Duitsland (4 stuks) zijn 3 euro. logrolls En dan komen er verzendkosten bij.

Er was zonet een draadje over iemand die wanhopig op zoek was naar zo'n kleine nano-cache houder. En ja, die heb ik niet als coulant artikel in de shoplijst staan. Wat niet wegneemt dat ik ze wel op bestelling kan leveren. Momenteel heb ik een paar blauwe en zwarte nano's in de pijplijn vanuit de USA zitten. Ik verkoop ze altijd (geen uitzonderingen!) met een reserve set logrolletjes uit Duitsland (die ook momenteel in de pijplijn zitten).

Even ter verklaring van het waarom: TomTom en ook Garmin met Mapsource gebruiken vector-kaarten. Ozi gebruikt raster-kaarten. Het vergt erg veel rekenkracht om alle pixels te herberekenen. En die rekenkracht is niet beschikbaar.

Gezien berichten van de organisatoren aan mij gericht, behoor ik ook tot het gezelschap. Ik zit wat onregelmatig, maar als het lukt zal ik ook aanwezig zijn.

Als een cacher een foto van een ander bij gc.com plaatst en dat komt uit en de auteurrechthebbende maakt bezwaar dan ben jij en niet gc.com het haasje. Ooit had Google met een site de bepaling dat alles wat deelnemers erop zetten van Google werd. Dat hebben ze bezuurd want iemand plaatste een foto van een brood-fotograaf en die eistte en kreeg heel veel vergoeding!,. Ook omdat Google zich de foto wederrechtelijk had toegeëigend!

Even rekenen: 10 $+ 200(gram) x 1, $c = 12$= 9,6euro afgerond + BTW+invoer geeft 11,50 euro. Verzend is 2,50 dus dat wordt 14,00 euro . En ik heb nog 1 exemplaar liggen.

Zeer zeker, alles van Groundspeak kan ik bij ze bestellen. Maar ik ga pas volgende maand opnieuw bestellen bij ze. Geef maar precies op (met link svp ter voorkomen van misverstanden) wat gewenst is. Let op: er komt wel wat bij vanwege de luchtvracht van 1 dollarcent per gram en BTW en invoer. Gewoon op de bestelpage van de goodiesshop in vak 'opmerkingen' zeggen wat gewenst is.

Ik heb er niet over rondgetoeterd maar de Auto-TB's waren afgelopen maand uitverkocht geraakt en de levering vanuit de USA verliep moeilijk en traag (productie en communicatie problemen tgv. vakanties) Nu zijn ze er weer, dus liefhebbers: kom maar op! Tevens is er nu een application-manual (een papiertje) oftewel gebruiksaanwijzing in het Engels bij iedere levering. Nog een puntje over fire-needles: soms worden die in stammen van bomen geprikt zodat de lichtbundel recht op de top schijnt. Dat werk niet goed. Ze moeten echt zijwaards aangestraald worden. Met de Trail Tack punaises is het juist: recht aanstralen. Nog iets: tijdens het nestseizoen zijn er diverse vogeltjes die glinsterdingen verzamelen, zelfs al zijn onze fire-needles prikkelig, dat houdt ze niet tegen kennelijk....

De affaire van de Heldense valkuilen is nog niet afgelopen. Niet dat er nieuwe kuilen zijn ontdekt, maar er zijn nieuwe verdachten. Tweede verdachte valkuilen gepakt Hoeveel man heb je nodig om een kuil te graven? Als de KPN bij mij in de straat graaft dan zijn ze ook met zijn drieën: eentje op de graafmachine; eentje bij de kuil om te waarschuwen als ze een kabel tegen komen een eentje die staat te opzichten.

Het onderwerp van deze draad was mij even ontgaan. Maar ik ben vandaag opgeport om mijn professorale licht over deze zaak te laten schijnen. Eerst over de input files. twf bevat linker bovenhoek en pixelgrootte (plus eventuele rotatie) in metrische coordinaten. Hetzelfde met een GeoTIF: ook metrisch. Maar met de complicatie dat de notatie niet standaard is. Divers GeoTIF bij mij laten Ozi Exploreer klagen bij import: geen geodata in de GeoTIF terwijl vele andere progjes het wel kunnen vinden. Bijvoorbeeld GeoTIF Examiner (Mentor software) laat zien wat er is. Met ECW loopt het beter: dat is een fabrikant-eigen format. De geoinfo in een ECW bekijk ik met ECWSpy of met ECW Header Editor van Earth Resource Mapping. De info is wat een twf ook bevat plus de datum (naam) en de projectie (naam). Maar als een importerend programma beide niet kent, dan moet je externe info aandragen. De projectie Onze kaarten zijn met een De DubbelProjectie volgens Schreiber gemaakt. Dat is een heel obscure die voor ons prima werkt omdat ons landje nagenoeg in een vierkant past. Gelukkig kunnen we deze projectie benaderen met iets anders. De fout die optreed is niet merkbaar op kaarten. Het lukt prima met de wel bekende Transverse Mercator. Wel moet je dan de juiste parameters hebben. Dat is uitgezocht en die zijn: Lat origin= 0.000000000 N Central meridian 5.387633333 Scale Factor 0.999907900 False Easting 155000.00, False Northing -5316592.00 Als software niet met de 0.00 N overweg kan (zero divide!)dan kan je de lat org. zetten op het nulpunt van ons raster in Amersfoort. Dan krijg je wel een andere False Northing. Welke? dat staat in de TRANSNAP documentatie. Als het programma de Bessel ellipsoïde kent dan is dit voldoende om terug te kunnen rekenen van meters naar graden (op de Bessel ellips). Ik heb hier de zaak in MatLab geprogrammerd. Dat ziet er dan zo uit: % rd2bessel % % Omrekenen geo coordinaten (tov Bessel) naar RD-grid coordinaten % Formules van Geodetic Datum of Australia website % 1 december 2001 function geo = rd2bessel(rd); 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] Heb je belangstelling vor het omgekeerde (laat ik het voor de volledigheid maar geven): % bessel2rd % % Omrekenen geo coordinaten (tov Bessel) naar RD-grid coordinaten % Formules van Geodetic Datum of Australia website % Gecontroleerd met GPS % 7 juni 2000 function coordinaten = bessel2rd(geo); deg = pi/180; % omrekenfactor naar radialen phi = geo(1)*deg; % lengte [graden] lab = geo(2)*deg; % breedte [graden] % 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 = a*(a0*phi - a2*sin(2*phi) + a4*sin(4*phi) - a6*sin(6*phi)); 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; rho = a*(1-e^2)/(1-e^2*(sin(phi))^2)^(3/2); nu = a/sqrt(1-e^2*(sin(phi))^2); psi = nu/rho; % Geographical to grid: t = tan(phi); omega = lab - lab0; term1 = 1/6*omega^2*(cos(phi))^2*(psi-t^2); term2 = 1/120*omega^4*(cos(phi))^4*(4*psi^3*(1-6*t^2)+psi^2*(1+8*t^2)-psi*2*t^2+t^4); term3 = 1/5040*omega^6*(cos(phi))^6*(61-479*t^2+179*t^4-t^6); E = (K0*nu*omega*cos(phi))*(1+term1+term2+term3) + Efalse; term1 = 1/2*omega^2*nu*sin(phi)*cos(phi); term2 = 1/24*omega^4*nu*sin(phi)*(cos(phi))^3*(4*psi^2+psi-t^2); term3 = 1/720*omega^6*nu*sin(phi)*(cos(phi))^5*(8*psi^4*(11-24*t^2)-28*psi^3*(1-6*t^2)+psi^2*(1-32*t^2)-psi*(2*t^2)+t^4); term4 = 1/40320*omega^8*nu*sin(phi)*(cos(phi))^7*(1385-3111*t^2+543*t^4-t^6); N = K0*(m + term1 + term2 + term3 + term4) + Nfalse; coordinaten = [E,N]; Heb je phi en lambda (in graden) op de Bessel Ellips dan zijn er standaard manieren om die positie over te zetten naar een andere ellipsoïde. In ons geval die van WGS84. Voor ons is de 5-parameter transformatie voldoende nauwkeurig. Zijn de parameters van beide ellipsen bekend dan dit door te rekenen. Hier gebeurt dit zo: % 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); 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; to_geo = [to_phi,to_lab,to_h]; Je ziet hoe gauw je een foutje maakt in de programmering van deze formules! In omgekeerde richting is het: % molodensky datum omzetting van WGS84 naar Bessel 1841 % Formules volgens website Peter Dana. % Gecontroleerd met rekenvoorbeeld en GPS (echter later nog fout in berekening Rm gehaald) % In/uitvoer in graden en meters % 1 dec 2001: hoogte niet *deg omrekenen! function [to_geo]=wgs2bessel(from_geo); 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: from_a = 6378137.000; % semi-major axis of ellipsoid from_f = 1/298.257223563; % 1/flattening % parameters voor Bessel 1841: to_a = 6377397.155; to_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 WGS 84 naar Bessel 1841: 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); % fout en oud: 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_phi = to_phi/deg; to_lab = to_lab/deg; to_h = from_h + dh; to_geo = [to_phi,to_lab,to_h]; Zodat mijn transfer functie simpel is: % omzetten RD x-y coordinaten (in meters) naar graden (tov WGS84) % 1 december 2001 function geo = rd2wgs(rd); bessel = rd2bessel(rd); bessel(3) = 0; %hoogte nul meter geo = bessel2wgs(bessel); Al deze bovenstaande formules zijn in de praktijk getest en tot op heden correct bevonden. Er is een andere manier natuurlijk om van meters RD naar graden WGS84 te gaan door middel van een directe polynome aanpassing. Het voordeel hiervan is dat er minder rekenkracht nodig is voor een polynoom. Mocht de omrekening te veel rekenkracht slurpen in de omzetting RD-> Bessel-graden (veel sinussen en zo) dan is er een reeksontwikkeleing mogelijk. Die staat beschreven op pag 31 en verder van het document Coördinaattransformaties en kaartprojecties van de meetkundige dienst.

Als dit de juiste griezel is, dan blijf ik met het raadsel zitten van de andere eerder opgepakte verdachte: Waar was die kuil wel voor bedoeld?

Uit de USA hebik liggen: 2 card games 1 frog pop klein potloden 1 truckers hat 2 groene petten verder camo -tape 1 nano cache met reserve logrolletjes