![](\"../images/tutorial08_cover.png\")
"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Set File Paths and download Tutorial Data\n",
"\n",
"If you downloaded the latest `GemGIS` version from the Github repository, append the path so that the package can be imported successfully. Otherwise, it is recommended to install `GemGIS` via `pip install gemgis` and import `GemGIS` using `import gemgis as gg`. In addition, the file path to the folder where the data is being stored is set. The tutorial data is downloaded using Pooch (https://www.fatiando.org/pooch/latest/index.html) and stored in the specified folder. Use `pip install pooch` if Pooch is not installed on your system yet."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"ExecuteTime": {
"end_time": "2021-03-17T11:05:43.666136Z",
"start_time": "2021-03-17T11:05:41.365748Z"
}
},
"outputs": [],
"source": [
"import gemgis as gg\n",
"\n",
"file_path ='data/08_sampling_interfaces_orientations_from_raster/'"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"ExecuteTime": {
"end_time": "2021-03-17T11:05:43.744631Z",
"start_time": "2021-03-17T11:05:43.715945Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Downloading file '08_sampling_interfaces_orientations_from_raster.zip' from 'https://rwth-aachen.sciebo.de/s/AfXRsZywYDbUF34/download?path=%2F08_sampling_interfaces_orientations_from_raster.zip' to 'C:\\Users\\ale93371\\Documents\\gemgis\\docs\\getting_started\\tutorial\\data\\08_sampling_interfaces_orientations_from_raster'.\n"
]
}
],
"source": [
"gg.download_gemgis_data.download_tutorial_data(filename=\"08_sampling_interfaces_orientations_from_raster.zip\", dirpath=file_path)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Loading raster data"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"ExecuteTime": {
"end_time": "2021-03-14T12:06:13.711835Z",
"start_time": "2021-03-14T12:06:11.518170Z"
}
},
"outputs": [],
"source": [
"import rasterio\n",
"\n",
"raster = rasterio.open(file_path + 'raster.tif')"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"ExecuteTime": {
"end_time": "2021-03-14T12:06:22.043105Z",
"start_time": "2021-03-14T12:06:22.032581Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"BoundingBox(left=-4.440892098500626e-16, bottom=1.3642420526593924e-12, right=971.002897605575, top=1068.1031873661339)"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"raster.bounds"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Plotting the raster"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"ExecuteTime": {
"end_time": "2021-03-14T12:06:13.964096Z",
"start_time": "2021-03-14T12:06:13.746704Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"| \n", " | X | \n", "Y | \n", "Z | \n", "geometry | \n", "formation | \n", "
|---|---|---|---|---|---|
| 0 | \n", "100.00 | \n", "100.00 | \n", "376.72 | \n", "POINT (100.00000 100.00000) | \n", "Test | \n", "
| 1 | \n", "200.00 | \n", "200.00 | \n", "435.22 | \n", "POINT (200.00000 200.00000) | \n", "Test | \n", "
| 2 | \n", "300.00 | \n", "300.00 | \n", "505.29 | \n", "POINT (300.00000 300.00000) | \n", "Test | \n", "
| 3 | \n", "400.00 | \n", "400.00 | \n", "540.68 | \n", "POINT (400.00000 400.00000) | \n", "Test | \n", "
| 4 | \n", "500.00 | \n", "500.00 | \n", "561.65 | \n", "POINT (500.00000 500.00000) | \n", "Test | \n", "
| \n", " | X | \n", "Y | \n", "Z | \n", "geometry | \n", "formation | \n", "
|---|---|---|---|---|---|
| 0 | \n", "100.00 | \n", "100.00 | \n", "378.80 | \n", "POINT (100.00000 100.00000) | \n", "Test | \n", "
| 1 | \n", "200.00 | \n", "200.00 | \n", "435.22 | \n", "POINT (200.00000 200.00000) | \n", "Test | \n", "
| 2 | \n", "300.00 | \n", "300.00 | \n", "505.29 | \n", "POINT (300.00000 300.00000) | \n", "Test | \n", "
| 3 | \n", "400.00 | \n", "400.00 | \n", "541.74 | \n", "POINT (400.00000 400.00000) | \n", "Test | \n", "
| 4 | \n", "500.00 | \n", "500.00 | \n", "562.02 | \n", "POINT (500.00000 500.00000) | \n", "Test | \n", "
| \n", " | X | \n", "Y | \n", "Z | \n", "geometry | \n", "formation | \n", "
|---|---|---|---|---|---|
| 0 | \n", "199.86 | \n", "496.09 | \n", "520.81 | \n", "POINT (199.85519 496.09330) | \n", "Test | \n", "
| 1 | \n", "511.85 | \n", "763.18 | \n", "643.24 | \n", "POINT (511.85172 763.17686) | \n", "Test | \n", "
| 2 | \n", "643.35 | \n", "886.38 | \n", "620.91 | \n", "POINT (643.35122 886.37867) | \n", "Test | \n", "
| 3 | \n", "965.65 | \n", "1045.58 | \n", "620.59 | \n", "POINT (965.65475 1045.58279) | \n", "Test | \n", "
| 4 | \n", "172.10 | \n", "864.05 | \n", "602.61 | \n", "POINT (172.10151 864.05005) | \n", "Test | \n", "
| 5 | \n", "744.41 | \n", "52.90 | \n", "276.00 | \n", "POINT (744.40762 52.90476) | \n", "Test | \n", "
| 6 | \n", "739.94 | \n", "823.36 | \n", "533.44 | \n", "POINT (739.93577 823.35867) | \n", "Test | \n", "
| 7 | \n", "942.70 | \n", "406.58 | \n", "414.82 | \n", "POINT (942.69846 406.58069) | \n", "Test | \n", "
| 8 | \n", "420.50 | \n", "87.25 | \n", "436.52 | \n", "POINT (420.50070 87.24755) | \n", "Test | \n", "
| 9 | \n", "339.68 | \n", "264.21 | \n", "503.23 | \n", "POINT (339.68451 264.21229) | \n", "Test | \n", "
| \n", " | X | \n", "Y | \n", "Z | \n", "geometry | \n", "formation | \n", "
|---|---|---|---|---|---|
| 0 | \n", "950.27 | \n", "424.37 | \n", "426.03 | \n", "POINT (950.26670 424.37299) | \n", "Test | \n", "
| 1 | \n", "205.91 | \n", "83.67 | \n", "410.35 | \n", "POINT (205.90842 83.67224) | \n", "Test | \n", "
| 2 | \n", "833.19 | \n", "424.42 | \n", "388.64 | \n", "POINT (833.18890 424.42058) | \n", "Test | \n", "
| 3 | \n", "442.79 | \n", "68.90 | \n", "427.47 | \n", "POINT (442.79487 68.90223) | \n", "Test | \n", "
| 4 | \n", "333.34 | \n", "642.61 | \n", "620.60 | \n", "POINT (333.33986 642.60911) | \n", "Test | \n", "
| 5 | \n", "905.91 | \n", "585.65 | \n", "484.48 | \n", "POINT (905.90705 585.64586) | \n", "Test | \n", "
| 6 | \n", "359.39 | \n", "1025.74 | \n", "742.35 | \n", "POINT (359.39233 1025.74463) | \n", "Test | \n", "
| 7 | \n", "599.42 | \n", "936.93 | \n", "660.12 | \n", "POINT (599.41731 936.92645) | \n", "Test | \n", "
| 8 | \n", "14.28 | \n", "532.99 | \n", "417.31 | \n", "POINT (14.28131 532.98620) | \n", "Test | \n", "
| 9 | \n", "177.50 | \n", "985.25 | \n", "643.90 | \n", "POINT (177.49500 985.24662) | \n", "Test | \n", "
| \n", " | X | \n", "Y | \n", "Z | \n", "geometry | \n", "formation | \n", "dip | \n", "azimuth | \n", "polarity | \n", "
|---|---|---|---|---|---|---|---|---|
| 0 | \n", "100.00 | \n", "100.00 | \n", "376.72 | \n", "POINT (100.00000 100.00000) | \n", "Test | \n", "23.77 | \n", "255.93 | \n", "1 | \n", "
| 1 | \n", "200.00 | \n", "200.00 | \n", "435.22 | \n", "POINT (200.00000 200.00000) | \n", "Test | \n", "30.46 | \n", "239.55 | \n", "1 | \n", "
| 2 | \n", "300.00 | \n", "300.00 | \n", "505.29 | \n", "POINT (300.00000 300.00000) | \n", "Test | \n", "17.34 | \n", "222.68 | \n", "1 | \n", "
| 3 | \n", "400.00 | \n", "400.00 | \n", "540.68 | \n", "POINT (400.00000 400.00000) | \n", "Test | \n", "15.51 | \n", "178.84 | \n", "1 | \n", "
| 4 | \n", "500.00 | \n", "500.00 | \n", "561.65 | \n", "POINT (500.00000 500.00000) | \n", "Test | \n", "19.26 | \n", "145.55 | \n", "1 | \n", "
| \n", " | X | \n", "Y | \n", "Z | \n", "geometry | \n", "formation | \n", "dip | \n", "azimuth | \n", "polarity | \n", "
|---|---|---|---|---|---|---|---|---|
| 0 | \n", "100.00 | \n", "100.00 | \n", "378.80 | \n", "POINT (100.00000 100.00000) | \n", "Test | \n", "23.75 | \n", "255.93 | \n", "1 | \n", "
| 1 | \n", "200.00 | \n", "200.00 | \n", "435.22 | \n", "POINT (200.00000 200.00000) | \n", "Test | \n", "30.44 | \n", "239.55 | \n", "1 | \n", "
| 2 | \n", "300.00 | \n", "300.00 | \n", "505.29 | \n", "POINT (300.00000 300.00000) | \n", "Test | \n", "17.32 | \n", "222.68 | \n", "1 | \n", "
| 3 | \n", "400.00 | \n", "400.00 | \n", "541.74 | \n", "POINT (400.00000 400.00000) | \n", "Test | \n", "15.50 | \n", "178.84 | \n", "1 | \n", "
| 4 | \n", "500.00 | \n", "500.00 | \n", "562.02 | \n", "POINT (500.00000 500.00000) | \n", "Test | \n", "19.24 | \n", "145.55 | \n", "1 | \n", "
| \n", " | X | \n", "Y | \n", "Z | \n", "geometry | \n", "formation | \n", "dip | \n", "azimuth | \n", "polarity | \n", "
|---|---|---|---|---|---|---|---|---|
| 0 | \n", "770.72 | \n", "454.44 | \n", "411.17 | \n", "POINT (770.72317 454.43623) | \n", "Test | \n", "35.50 | \n", "126.05 | \n", "1 | \n", "
| 1 | \n", "98.85 | \n", "104.73 | \n", "377.13 | \n", "POINT (98.85090 104.72719) | \n", "Test | \n", "24.00 | \n", "256.35 | \n", "1 | \n", "
| 2 | \n", "439.88 | \n", "23.92 | \n", "417.44 | \n", "POINT (439.88303 23.91500) | \n", "Test | \n", "15.43 | \n", "138.91 | \n", "1 | \n", "
| 3 | \n", "444.39 | \n", "434.85 | \n", "548.88 | \n", "POINT (444.39082 434.84986) | \n", "Test | \n", "16.92 | \n", "165.67 | \n", "1 | \n", "
| 4 | \n", "642.46 | \n", "857.77 | \n", "612.54 | \n", "POINT (642.45819 857.77092) | \n", "Test | \n", "33.90 | \n", "115.23 | \n", "1 | \n", "
| 5 | \n", "668.68 | \n", "101.45 | \n", "321.65 | \n", "POINT (668.68137 101.45438) | \n", "Test | \n", "34.24 | \n", "122.96 | \n", "1 | \n", "
| 6 | \n", "833.16 | \n", "737.23 | \n", "491.34 | \n", "POINT (833.16382 737.23242) | \n", "Test | \n", "3.94 | \n", "242.61 | \n", "1 | \n", "
| 7 | \n", "137.75 | \n", "604.27 | \n", "506.55 | \n", "POINT (137.75113 604.27271) | \n", "Test | \n", "39.04 | \n", "252.36 | \n", "1 | \n", "
| 8 | \n", "249.23 | \n", "610.05 | \n", "590.12 | \n", "POINT (249.22976 610.04522) | \n", "Test | \n", "30.48 | \n", "224.83 | \n", "1 | \n", "
| 9 | \n", "799.49 | \n", "1011.20 | \n", "581.62 | \n", "POINT (799.49348 1011.19910) | \n", "Test | \n", "11.49 | \n", "170.42 | \n", "1 | \n", "
| \n", " | X | \n", "Y | \n", "Z | \n", "geometry | \n", "formation | \n", "dip | \n", "azimuth | \n", "polarity | \n", "
|---|---|---|---|---|---|---|---|---|
| 0 | \n", "433.83 | \n", "627.61 | \n", "606.19 | \n", "POINT (433.82765 627.60890) | \n", "Test | \n", "19.29 | \n", "144.83 | \n", "1 | \n", "
| 1 | \n", "732.25 | \n", "1007.74 | \n", "594.51 | \n", "POINT (732.25432 1007.73581) | \n", "Test | \n", "21.30 | \n", "109.13 | \n", "1 | \n", "
| 2 | \n", "849.26 | \n", "539.88 | \n", "444.47 | \n", "POINT (849.25852 539.87992) | \n", "Test | \n", "29.99 | \n", "196.23 | \n", "1 | \n", "
| 3 | \n", "506.54 | \n", "297.89 | \n", "499.94 | \n", "POINT (506.54294 297.89069) | \n", "Test | \n", "22.06 | \n", "134.76 | \n", "1 | \n", "
| 4 | \n", "482.80 | \n", "959.46 | \n", "715.72 | \n", "POINT (482.79546 959.45977) | \n", "Test | \n", "21.73 | \n", "129.97 | \n", "1 | \n", "
| 5 | \n", "914.62 | \n", "832.23 | \n", "543.36 | \n", "POINT (914.61626 832.23124) | \n", "Test | \n", "45.25 | \n", "231.65 | \n", "1 | \n", "
| 6 | \n", "16.88 | \n", "805.99 | \n", "447.03 | \n", "POINT (16.88303 805.98712) | \n", "Test | \n", "37.63 | \n", "260.69 | \n", "1 | \n", "
| 7 | \n", "97.01 | \n", "896.66 | \n", "534.58 | \n", "POINT (97.00884 896.65803) | \n", "Test | \n", "48.06 | \n", "257.25 | \n", "1 | \n", "
| 8 | \n", "528.77 | \n", "325.33 | \n", "500.61 | \n", "POINT (528.77290 325.33217) | \n", "Test | \n", "20.33 | \n", "138.42 | \n", "1 | \n", "
| 9 | \n", "635.98 | \n", "435.47 | \n", "500.45 | \n", "POINT (635.98363 435.47230) | \n", "Test | \n", "29.20 | \n", "118.66 | \n", "1 | \n", "