extent schalen about nlmod HOT 4 CLOSED

Volgens mij gaan we er nu vanuit dat regis-data gedefinieerd is op op ronde 100 getallen, dus op x = 200100, 200200, 200300, enz. Volgens mij is dit een foutje in de netcdf van regis. In deze netcdf zit ook een variabele genaamd x_bounds en y_bounds, wat de randen van de cellen voorstellen. De dimensies x en y moeten dan zijn:

url=r'http://www.dinodata.nl:80/opendap/REGIS/REGIS.nc'
regis = xr.open_dataset(url, decode_times=False)
regis['x'] = regis.x_bounds.mean('bounds')
regis['y'] = regis.y_bounds.mean('bounds')

Dit zorgt ervoor dat de coördinaten dan in het midden van de cellen staan, dus x = 200150, 200250, 200350, enz. Dit lijkt me ook een logischer indeling. Oftewel: ik denk dat we nu een foutje maken door alles 50 meter te veel naar het westen en 50 meter te veel naar het zuiden te verplaatsen. Niet dat dat op de schaal van regis tot veel problemen leidt, maar we moeten het wel goed hebben.

from nlmod.

Je hebt gelijk daar is nog geen rekening mee gehouden. Ik zal de get_regis_dataset en de fit_extent_to_regis functies aanpassen met behulp van de code die je stuurt.

Ik denk dat hetzelfde geldt voor de geotop dataset. Daarin zijn de x- en y-coordinaten nu ook gehele hondertallen. Helaas zit bij die dataset geen bounds variabele om dat te checken.

from nlmod.

Het is best lastig te achterhalen hoe dit voor geotop zit. In deze pdf staat het volgende:

Voor de geïnterpreteerde boormonsterprofielen geven de coördinaten de ligging van de boorlocatie aan maaiveld aan. Voor het lagenmodel geldt de conventie dat de ligging van een rastercelwordt beschreven door de coördinaten van de linkeronderhoek (“lower left corner”). Voor het voxelmodel geldt de conventie dat de horizontale ligging van een voxel wordt beschreven door de coördinaten van het middelpunt (“cell center”) van de voxel.

Ik ga er nu even vanuit dat de ligging van een rastercel in de geotop netcdf wordt beschreven door de coördinaten van de linkeronderhoek (net als in REGIS). Ik zal ook het inlezen van de geotop netcdf aanpassen.

from nlmod.

Je kan geotop en regis in QGIS openen toch? Wellicht dat je daarmee de juiste ligging kan bepalen?

from nlmod.