mcflugen
commented
1 week ago @ChristianMohr I think your suspicion is correct. read_esri_ascii will try to add the data from your ASCII file to your grid's nodes. However, it sounds like your grid is sized such that the data should be placed at cells.
A couple ideas:
- You could use a smaller grid that whose shape matches that of your data file and then load the data onto its nodes.
- You could create a new grid from your data file using
read_esri_asciiand then place that onto the cells of your current grid.
Your issue has reminded me that I had been working on enhancing the ESRI ASCII functionality a while ago but got distracted with something else and never finished the work. I'll get back on that and once I finish, I think it should make what you are trying to do easier.
ChristianMohr
Dear all, I am in the process of creating a Landlab model that combines different components, e.g. the ecohydrology and landslide components. While I had no problems using raster data provided to me, the model now has problems with the data I created myself. I created raster data via QGIS, which I exported as ASCII. All rasters have the same resolution, are correctly aligned and have the same spatial extent. Everything has been checked several times. My raster data includes rasters of plant functional types, vegetation cover, soil density, soil thickness, biomass and a DEM. When I load the data as ASCII files into Landlab, I get the following error message:
ValueError: unable to reshape array to field size (24011 != 24660)
I suspect that this problem stems from where the data is stored, i.e. nodes or cells, which leads to this form error. Do you have any idea how to solve this problem? Does anyone know this problem and have any solutions? Cheers and thanks in advance, Christian