Natural Constraints Post Processing

import geolime as geo
from pyproj import CRS, Transformer
import numpy as np
import osmnx as ox
import json

geo.Project().set_crs(CRS('EPSG:20350'))

bm = geo.read_file("../data/block_model.geo")

bm_fe_map = geo.plot_2d(bm, 'Fe_kriged', 'mean')

bm_fe_map