Prof. Y. Lupardi

Met alleen je handheld GPS kan je alsvolgt te werk gaan: Houdt de mapdatum op WGS 84 maar zet de weergave van de posities op UTM. Dan kan je waypoints met de hand inkloppen met hun metrische UTM coördinaten. Wil je zien wat dat in graden is dan zet je de weergave op graden (en/of minuten en/of seconden) Het metrische RD-raster van ons landje is slechts geldig voor ons land, niet in Oostenrijk. En ons RD raster is niet gebaseerd op mapdatum WGS84 maar op Bessel 1841. In Oostenrijk hebben ze hun eigen Austrian Grid (zit waarschijnlijk ook in je GPSr) maar staat mogelijk niet op jouw kaart aangegeven. UTM is over de hele wereld bruikbaar als je maar rekening houdt met de zone waar je in zit.

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?

Er is de laatste tijd iets raars aan het gebeuren met mijn caches. Zo kreeg ik eerst bericht dat men 'Eindpunt' niet kon vinden. In plaats daarvan werd een speelgoed autobusje van Connexion aangetroffen. Ik moet dit nog onderzoeken. Eerder in de maand was er opeens bericht dat mijn cache 'Drassig' na vele jaren verdwenen was. Gisteren kreeg ik bericht dat mijn vervangende cache (op de oude plaats) gewoon gevonden wordt en dat men een brug verder en 100 meter de originele cachedoos terugvond verstopt op analoge manier (bij boomwortels naast het pad maar zonder geel lintje). Is er een nieuwe geocacher aan de gang die het allemaal nog niet goed begrijpt?

Ik weet van gebruikers met een ARM-CPU in de PDA die OziCE gebruiken met kaarten in ecw format die heel groot zijn. Zo''n 600 Mb. Of zulke grote kaarten uberhaupt gemaakt kunnen worden in het format ozfx, dat weet ik niet. Ik dacht van wel. Mij is niet bekend dat realtime navigatie moeilijk doet op grote kaarten; alleen het deel dat benodigd is voor schermweergave wordt in het snelle geheugen geladen. En over het algemeen zijn geheugenkaarten/systeem/ozi snel genoeg om op tijd nieuw benodigde kaartblokken te zoeken, vinden en in te laden. Het ecw format vergt wat meer rekenkracht (en een passende CPU) dan het ozf format.

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

Nee helaas. Ik had gewenst om goedkope waterbestendige logboekjes te kunnen leveren maar helaas ze zijn niet meer te vinden voor een prijs die redelijk is. En het heeft ook lang geduurd (door de prijs?) voordat ik al de groene boekjes kwijt was. Dus ik ga er niet meer achteraan.

Vroeger mat men hoogte verschillen met een waterpas (lange slang). Het bleek dat als de hoogte bouten nu ongeveer op kniehoogte zitten, dat dat het handigste is voor de meetopstelling. Vaak staat er op de ronde bout NAP. Maar niet op allen. Sommige kunnen gezet zijn door het polder/waterschap.

De jongste stadia van de teek leven van kleine zoogdieren. Muizen en dat soort grut. Kleine jonge teekjes hoeven niet erg hoog te klimmen. Worden ze groter dan komen b.v. konijnenoren in aanmerking. En katten en honden. Hiervoor klimmen ze ook wat hoger. Het eindstadium doet het met herten en zwijnen en zo en klimt dus hoog. Dus hoger dan een meter hoeft er niet geklommen te worden.

Niet alleen bij veerponten in warme streken is er teken-risico, maar ook in ons landje kleeft aan het gebruik van het openbaar vervoer een zeker risico. Was er kort geleden niet in het draadje over de komkommertijd een eng bericht over een vrouw die in de Amsterdamse tram besprongen werd door.... teken! ??? TEKEN UIT AMSTERDAM

De uitgaven van de topo 25k van het Kadaster (topodienst) van vorig jaar bevatten rasterkaarten vanaf 1999 tot 2006. Vanaf het moment dat ze via Garmin de vector-topo gingen maken tot de uitgave ging zeker meer dan een jaar voorbij. Die moet je bij de veroudering mee tellen. Met topografische gegevens moet je erop rekenen dat je soms 4 jaar achter loopt. Er zijn landen waar de zaak 10 jaar achter loopt!!! (Balkan)

Hoe kiezen ze een Miss World? Of nog beter een Miss Universe? Zoveel formaten, smaken en inzichten. Zoveel verschillen. etc. En toch hebben we voortdurend overal Miss verkiezingen. Ik stel voor de methode van de Miss-verkiezing ook voor caches te gebruiken.

Het is mij al twee keer overkomen bij cache "Biertje" dat ze het mooie legergroene verbanddoosje hebben meegenomen en de inhoud in een plastic zak gestopt. Nu zijn mijn mooie doosjes op (weet iemand ze nog ergens te koop) en gebruik ik simpele tupperware van de Wibra. Jammer maar ja, als je een cache plaats is die overgeleverd aan zijn eigen lot. Je kan er weinig aan doen verder.

Ooit vroeger had je kaarten gebasserd kwa posities op de meridiaan van Parijs. Moderne Franse kaarten van hun topo-dienst IGN gebruiken NTF en vier verschillende grids (zones) Ik heb het ooit in mijn oude gPS12 gezet: > - POSITION FRMT: USER GRID > - MAP DATUM: USER (Clark 1880 ellipsoid, IGN modified) > DX: -168 mt > DY: -60 mt > DZ: +320 mt > DA: -112 mt > DF: -0.5473884 > - original meridian is Paris and not Greenwich > - offset Paris is 2° 20' 14.025" east of Greenwich. > - angles are in grades : 100 grades = 90 degrees. Maar de Gele Michelin moet iets anders hebben zegt men: > Projection : Bonne (Conical) > Ellipsoid : Plessis > Origin : 2°20'14" E (Paris), 48°51'36" N DX, DY, DZ in relation with > WGS 84 > a = axis = 6376523 m > f = flattening = 0.0032400188611 > DX = 1118 m > DY = 23 m > DZ = 66 m > > These data are for map Nr 51 à 89, 230 à 246 > > N43° 15' 08" > E05° 33' 25" Dus verwarring alom. Je zou het eens aan een inboorling kunnen vragen

Voor de echte vreselijk fanatieke cache-makers heb ik momenteel 3 (drie) van die miezerig kleine nano-caches (zwart). Ik verkoop ze met een setje reservelogboekjes. Er zijn heden ook weer een paar groene en zwarte petten beschikbaar.

De letters eu vormen geen ligatuur in druktechnische zin maar hebben wel een verbod om losgesplitst te worden (afbreekregels). Nieuwerwets mag je niet afbreken als ni-euwerwets... of nie-uwerwets.. of nieu-werwets. Iets dat een ligatuur vormt mag je niet afbreken; iets dat niet afgebroken mag worden is daarmee nog geen ligatuur.

Typografisch gezien is de ij 1 letter. Dat heet een ligatuur. Zo kan je bij het afbreken van een woord nooit deze splitsen van elkaar. Vermi-jding gaat niet. Vermij-ding wel. Evenzo is er effect bij gebruik van kapitalen. Men schrijft IJmuiden. En 'IJsberen dreigen te verzuipen'. Als je een voor West Europa gemaakt letertype gebruikt dan zit daar de ij in net zoals andere samentrekkingen (a met e, s met s) Onze ij heeft dan typografisch minder ruimte dan de losse i + j . Ook heeft het originele Nederlandse typemachine toetsenbord een aparte 'ij' die de standaard ruimte opvult. Een i met een j erachter neemt 2 standardruimten in.

ik heb begrepen dat de plaatselijke vogelbescherming nu de kinderen begeleid. Barny Friese kinderen krijgen bescherming tegen buizerd "It is a jungle out there!!" De gevaarlijke situatie is opgelost. Alsof het een auto-toertje door een safari park betreft: raampjes dicht vanwege de Wilde Beesten. Dit lees ik bij Weer agressieve buizerd in Flevoland

Zo te lezen in de specs is het mogelijk om Ozi CE te draaien erop en dan heb je heel wat kaarten ter beschikking zonder afhankelijk te zijn van een of ander specifiek fileformat of programma. Dat er een ARM CPU in zit is heel welkom: dan zijn ook ecw-format kaarten met Ozi Ce te gebruiken. Helemaal ideaal. Heeft iemand het al geprobeerd?

Ik kreeg het verzoek dit onder de aandacht te brengen. Na de Auto-Travel Bug nu de GPSr Travel Bug? Garmin Colorado 300 - Geocaching.com Trackable Special Limited Edition Het was aan de goodiesshop gevraagd om mee te helpen aan een groepsbestelling met ietwat korting, maar de 'policy' van de goodiesshop is om zich verre te houden van de verkoop van alles wat met hardware/electronica zoals GPSr's van doen heeft.