Merge pull request #241 from ArtesiaWater/epsg-28992-issue

fix Projection issue in backgroundmaps #240
gwmod · Aug 25, 2023 · c73f524 · c73f524
2 parents 8ecca8a + 33d6478
commit c73f524
Show file tree

Hide file tree

Showing 3 changed files with 11 additions and 3 deletions.
diff --git a/nlmod/dims/base.py b/nlmod/dims/base.py
@@ -7,6 +7,7 @@
 from scipy.spatial import cKDTree
 
 from .. import util
+from ..epsg28992 import EPSG_28992
 from . import resample
 from .layers import fill_nan_top_botm_kh_kv, set_idomain
 
@@ -490,7 +491,7 @@ def get_ds(
     botm=None,
     kh=10.0,
     kv=1.0,
-    crs=28992,
+    crs=EPSG_28992,
     xorigin=0.0,
     yorigin=0.0,
     angrot=0.0,

diff --git a/nlmod/epsg28992.py b/nlmod/epsg28992.py
@@ -0,0 +1,7 @@
+# NOTE: this is the correct epsg:28992 definition for plotting backgroundmaps in RD
+EPSG_28992 = (
+    "+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 "
+    "+x_0=155000 +y_0=463000 +ellps=bessel "
+    "+towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725 +units=m "
+    "+no_defs"
+)
diff --git a/nlmod/plot/plotutil.py b/nlmod/plot/plotutil.py
@@ -4,7 +4,7 @@
 from matplotlib.ticker import FuncFormatter, MultipleLocator
 
 from ..dims.resample import get_affine_mod_to_world
-
+from ..epsg28992 import EPSG_28992
 
 def get_patches(ds, rotated=False):
     """Get the matplotlib patches for a vertex grid."""
@@ -88,7 +88,7 @@ def get_map(
     fmt="{:.0f}",
     sharex=False,
     sharey=True,
-    crs=28992,
+    crs=EPSG_28992,
     background=False,
     alpha=0.5,
     tight_layout=True,