Zet de coördinaten op de kaart: ‘geocoding’ in Orange

Je hebt een Excel file met daarin de kaart coördinaten (Latitude, Longitude)(WGS-84 format). NB! Als je gegevens hebt in RijksDriehoekscoördinaten (RD-coördinaten) is er code beschikbaar om die in Excel om te zetten naar WGS-84 coördinaten … zie de bijlage onderaan deze pagina.

Maar hoe laat ik die op een geografische kaart zien?

voorbeeldje. Let op de kolommen Latitude en Longitude …

Welnu, met het (gratis) programma Orange kan dat heel simpel: download en installeer het programma Orange op jouw laptop, start het programma en laad de onderstaande file. (Tsja, dus eerst ook die file even downloaden en onthouden waar die staat, doorgaans in de map ‘Downloads’)

Je krijgt na het openen van de bovenstaande file in Orange een grafisch beeld van het programma te zien in de vorm van icoontjes – de functionele onderdelen van het programma – en lijntjes die de data-stromen weergeven. In feite het proces van gegevensverwerking van links naar rechts in één oogopslag!

Dat ziet er dan zò uit als je op de eerste en de laatste icoon hebt geklikt:

In bovenstaand voorbeeld is een voorbeeldbestand geselecteerd, en de uitkomst (de kaart met de coördinaat-punten als resultaat)

Je kunt dit proces (de onderstaande file) dus voor allerhande datafiles die geografische coördinaten bevatten eindeloos hergebruiken zonder ook maar iets aan de file met het proces (de icoontjes, lijntjes en het de er bij gegeven aanwijzingen) (*.ows) te veranderen. Het proces blijft zoals het is, en als je een andere invoer tabel selecteert veranderd de resulterende kaart uiteraard…

Dus het enige dat je natuurlijk wel echt moet doen is de datafile selecteren en op de uitkomst klikken. Simpeler is het gewoon niet te krijgen.

Kortom: hier is de file nog eens die je in Orange moet openen nadat je Orange hebt geïnstalleerd:

en nog even een voorbeeldje van een andere invoer (AWN Afd. 12);

Sukses!


Bijlage: Code voor omzetting RD- naar WGS-coördinaten

( Input = RD[X_Coord], RD[Y_Coord] )


XLong =

VAR dX = (RD[X_COORD] – 155000) * 10 ^ -5

VAR dY = (RD[Y_COORD] – 463000) * 10 ^ -5

VAR SomN = (3235.65389 * dY) + (-32.58297 * dX ^ 2) + (-0.2475 * dY ^ 2) + (-0.84978 * dX ^ 2 * dY) + (-0.0655 * dY ^ 3) + (-0.01709 * dX ^ 2 * dY ^ 2) + (-0.00738 * dX) + (0.0053 * dX ^ 4) + (-0.00039 * dX ^ 2 * dY ^ 3) + (0.00033 * dX ^ 4 * dY) + (-0.00012 * dX * dY)

VAR SomE = (5260.52916 * dX) + (105.94684 * dX * dY) + (2.45656 * dX * dY ^ 2) + (-0.81885 * dX ^ 3) + (0.05594 * dX * dY ^ 3) + (-0.05607 * dX ^ 3 * dY) + (0.01199 * dY) + (-0.00256 * dX ^ 3 * dY ^ 2) + (0.00128 * dX * dY ^ 4) + (0.00022 * dY ^ 2) + (-0.00022 * dX ^ 2) + (0.00026 * dX ^ 5)

VAR Latitude = 52.15517 + (SomN / 3600)

VAR Longitude = 5.387206 + (SomE / 3600)

VAR LatitudeGraden = INT(Latitude)

VAR LongitudeGraden = INT(Longitude)

VAR LatitudeMinuten = (Latitude – LatitudeGraden) * 60

VAR LongitudeMinuten = (Longitude – LongitudeGraden) * 60

RETURN Longitude


YLat =

VAR dX = (rce_wng[X_COORD] – 155000) * 10 ^ -5

VAR dY = (rce_wng[Y_COORD] – 463000) * 10 ^ -5

VAR SomN = (3235.65389 * dY) + (-32.58297 * dX ^ 2) + (-0.2475 * dY ^ 2) + (-0.84978 * dX ^ 2 * dY) + (-0.0655 * dY ^ 3) + (-0.01709 * dX ^ 2 * dY ^ 2) + (-0.00738 * dX) + (0.0053 * dX ^ 4) + (-0.00039 * dX ^ 2 * dY ^ 3) + (0.00033 * dX ^ 4 * dY) + (-0.00012 * dX * dY)

VAR SomE = (5260.52916 * dX) + (105.94684 * dX * dY) + (2.45656 * dX * dY ^ 2) + (-0.81885 * dX ^ 3) + (0.05594 * dX * dY ^ 3) + (-0.05607 * dX ^ 3 * dY) + (0.01199 * dY) + (-0.00256 * dX ^ 3 * dY ^ 2) + (0.00128 * dX * dY ^ 4) + (0.00022 * dY ^ 2) + (-0.00022 * dX ^ 2) + (0.00026 * dX ^ 5)

VAR Latitude = 52.15517 + (SomN / 3600)

VAR Longitude = 5.387206 + (SomE / 3600)

VAR LatitudeGraden = INT(Latitude)

VAR LongitudeGraden = INT(Longitude)

VAR LatitudeMinuten = (Latitude – LatitudeGraden) * 60

VAR LongitudeMinuten = (Longitude – LongitudeGraden) * 60

RETURN Latitude