2D Boundary Conditions with EPSG:27700 gives 3Di Gridbuilder error (TOPDESK 2311 0136)

leendertvanwolfswinkel commented 10 months ago

3Di model (https://management.3di.live/schematisations/8456/53265). This fails on make_gridadmin, with the following feedback: "2D boundary condition 1 does not touch any edge cell". Unfortunately this is not the case and the 2D boundary condition line does touch an edge cell.

It probably depends on the EPSG code of the model. This one is 27700. Once the model was converted to a German projection, the problem was no longer there.

We have had issues with EPSG:27700 before. The 3Di documenation used list this "known issue". This may be helpful when analysing this problem.

" There are some issues related to projections. We are encountered this for our UK-based models (.sqlite). These models will show a shift in the geometries in QGIS. The reprojection from EPGS:27700 (British National Grid) to WGS84 is based on an outdated reprojection in the spatialite. There is a solution for this:

In QGIS open the (table without geometry) table "spatial_ref_sys"
go to row srid 27700
update column proj4text (overwrite existing value) with +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +datum=OSGB36 +units=m +no_defs "

leendertvanwolfswinkel commented 7 months ago

I suspect that perhaps #257 may have fixed this

leendertvanwolfswinkel commented 6 months ago

Test with https://management.staging.3di.live/schematisations/571/6012?organisation=b08433fa47c1401eb9cbd4156034c679

nens / threedigrid-builder

2D Boundary Conditions with EPSG:27700 gives 3Di Gridbuilder error (TOPDESK 2311 0136) #343