Generate 3D geological model inputs from geographical maps — a high-level implementation and extension of https://github.com/Loop3D/map2loop by Prof. Mark Jessell at UWA. To see an example interactive model built with map2loop and LoopStructural, follow this link:
3D Model from the Hamersley region, Western Australia
You will need some flavour of conda (a python package manager, see here), as well as Python ≥ 3.6
To just use map2loop, issue the following
conda install -c conda-forge -c loop3d map2loop -yIf you want to tinker yourself/contribute, clone the source code with
git clone https://github.com/Loop3D/map2loop-2.gitOr get the source + example notebooks with
git clone https://github.com/Loop3D/map2loop-2.git
git clone --single-branch --branch yohan https://github.com/Loop3D/map2loop2-notebooksNavigate into map2loop-2, and issue the following to install map2loop and its dependencies. Note: The 'develop' flag makes your source changes take effect on saving, so you only need to run this once
python setup.py developFair warning, we recommend conda to almost everyone. With great software development power comes great environment setup inconvenience. You'll need to download and install the docker containerisation software, and the docker and docker-compose CLI.
-
Clone this repo and navigate inside as per above
-
Run the following and click on the Jupyter server forwarded link to access and edit the notebooks
docker-compose up --build
-
To hop into a bash shell in a running container, open a terminal and issue
docker ps
Find the container name or ID and then run
docker exec -it <container_NAMEorID> bash # Probably -> docker exec -it map2loop-2_dev_1 bash
Our notebooks cover use cases in more detail, but here is an example of processing Loop's South Australia remote geospatial data in just 20 lines of Python.
First, lets import map2loop and define a bounding box. You can use GIS software to find one or use Loop's Graphical User Interface for the best experience and complete toolset. Remember what projection your coordinates are in!
from map2loop.project import Project
bbox_3d = {
'minx': 250805.1529856466,
'miny': 6405084.328058686,
'maxx': 336682.921539395,
'maxy': 6458336.085975628,
'base': -3200,
'top': 1200
}
Then, specify: the state, directory for the output, the bounding box and projection from above - and hit go! That's it.
proj = Project(loopdata_state="SA")
proj.update_config(out_dir='sa-remote',
overwrite="true",
bbox_3d=bbox_3d,
proj_crs={'init': 'EPSG:28354'},
# drift_prefix=['T', 'Q', 'water', 'void'],
# quiet='no-figures',
)
proj.run()This is a minimal example and a small part of Loop.
Our documentation and other resources outline how to extend map2loop and port to the LoopStructural modelling engine. We are working to incorporate geophysical tools and best provide visualisation and workflow consolidation in the GUI.
Loop is led by Laurent Ailleres (Monash University) with a team of Work Package leaders from:
- Monash University: Roy Thomson, Lachlan Grose and Robin Armit
- University of Western Australia: Mark Jessell, Jeremie Giraud, Mark Lindsay and Guillaume Pirot
- Geological Survey of Canada: Boyan Brodaric and Eric de Kemp
- Developing with docker on Windows means you won't have GPU passthrough and can’t use a discrete graphics card in the container even if you have one.
- If Jupyter links require a token or password, it may mean port 8888 is already in use. To fix, either make docker map to another port on the host ie -p 8889:8888 or stop any other instances on 8888.
- Sometimes the submodules misbehave. Ensure you specify the recurse-submodules flag with an 's' as opposed to recurse-submodule, and double-check your .git directory is clean, see Loop3D#41