![](\"../images/tutorial02_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-17T10:52:58.444997Z",
"start_time": "2021-03-17T10:52:56.150000Z"
}
},
"outputs": [],
"source": [
"import gemgis as gg\n",
"\n",
"file_path ='data/02_extract_xyz/'"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"ExecuteTime": {
"end_time": "2021-03-17T10:52:58.525121Z",
"start_time": "2021-03-17T10:52:58.495075Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Downloading file '02_extract_xyz.zip' from 'https://rwth-aachen.sciebo.de/s/AfXRsZywYDbUF34/download?path=%2F02_extract_xyz.zip' to 'C:\\Users\\ale93371\\Documents\\gemgis\\docs\\getting_started\\tutorial\\data\\02_extract_xyz'.\n"
]
}
],
"source": [
"gg.download_gemgis_data.download_tutorial_data(filename=\"02_extract_xyz.zip\", dirpath=file_path)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Point Data\n",
"\n",
"The point data stored as shape file will be loaded as GeoDataFrame. The raster will be loaded using rasterio. "
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"ExecuteTime": {
"end_time": "2020-12-29T18:22:50.980061Z",
"start_time": "2020-12-29T18:22:44.811443Z"
}
},
"outputs": [],
"source": [
"import rasterio\n",
"from rasterio.plot import show\n",
"import geopandas as gpd\n",
"\n",
"gdf = gpd.read_file(file_path + 'interfaces_points.shp')\n",
"\n",
"dem = rasterio.open(file_path + 'raster.tif')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The GeoDataFrame consists of ``id``, ``formation`` and ``geometry`` columns."
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"ExecuteTime": {
"end_time": "2020-12-29T18:22:51.010195Z",
"start_time": "2020-12-29T18:22:50.983063Z"
}
},
"outputs": [
{
"data": {
"text/html": [
"| \n", " | id | \n", "formation | \n", "geometry | \n", "
|---|---|---|---|
| 0 | \n", "None | \n", "Ton | \n", "POINT (19.15013 293.31349) | \n", "
| 1 | \n", "None | \n", "Ton | \n", "POINT (61.93437 381.45933) | \n", "
| 2 | \n", "None | \n", "Ton | \n", "POINT (109.35786 480.94557) | \n", "
| 3 | \n", "None | \n", "Ton | \n", "POINT (157.81230 615.99943) | \n", "
| 4 | \n", "None | \n", "Ton | \n", "POINT (191.31803 719.09398) | \n", "
| \n", " | formation | \n", "geometry | \n", "X | \n", "Y | \n", "Z | \n", "
|---|---|---|---|---|---|
| 0 | \n", "Ton | \n", "POINT (19.15013 293.31349) | \n", "19.15 | \n", "293.31 | \n", "364.99 | \n", "
| 1 | \n", "Ton | \n", "POINT (61.93437 381.45933) | \n", "61.93 | \n", "381.46 | \n", "400.34 | \n", "
| 2 | \n", "Ton | \n", "POINT (109.35786 480.94557) | \n", "109.36 | \n", "480.95 | \n", "459.55 | \n", "
| 3 | \n", "Ton | \n", "POINT (157.81230 615.99943) | \n", "157.81 | \n", "616.00 | \n", "525.69 | \n", "
| 4 | \n", "Ton | \n", "POINT (191.31803 719.09398) | \n", "191.32 | \n", "719.09 | \n", "597.63 | \n", "
| \n", " | id | \n", "formation | \n", "geometry | \n", "
|---|---|---|---|
| 0 | \n", "None | \n", "Sand1 | \n", "LINESTRING (0.25633 264.86215, 10.59347 276.73... | \n", "
| 1 | \n", "None | \n", "Ton | \n", "LINESTRING (0.18819 495.78721, 8.84067 504.141... | \n", "
| 2 | \n", "None | \n", "Ton | \n", "LINESTRING (970.67663 833.05262, 959.37243 800... | \n", "
| \n", " | formation | \n", "geometry | \n", "X | \n", "Y | \n", "Z | \n", "
|---|---|---|---|---|---|
| 0 | \n", "Sand1 | \n", "POINT (0.25633 264.86215) | \n", "0.26 | \n", "264.86 | \n", "353.97 | \n", "
| 1 | \n", "Sand1 | \n", "POINT (10.59347 276.73371) | \n", "10.59 | \n", "276.73 | \n", "359.04 | \n", "
| 2 | \n", "Sand1 | \n", "POINT (17.13494 289.08982) | \n", "17.13 | \n", "289.09 | \n", "364.28 | \n", "
| 3 | \n", "Sand1 | \n", "POINT (19.15013 293.31349) | \n", "19.15 | \n", "293.31 | \n", "364.99 | \n", "
| 4 | \n", "Sand1 | \n", "POINT (27.79512 310.57169) | \n", "27.80 | \n", "310.57 | \n", "372.81 | \n", "
| \n", " | id | \n", "formation | \n", "geometry | \n", "
|---|---|---|---|
| 0 | \n", "None | \n", "Sand1 | \n", "POLYGON ((0.25633 264.86215, 10.59347 276.7337... | \n", "
| 1 | \n", "None | \n", "Ton | \n", "POLYGON ((0.25633 264.86215, 0.18819 495.78721... | \n", "
| 2 | \n", "None | \n", "Sand2 | \n", "POLYGON ((0.18819 495.78721, 0.24897 1068.7595... | \n", "
| 3 | \n", "None | \n", "Sand2 | \n", "POLYGON ((511.67477 1068.85246, 971.69794 1068... | \n", "
| \n", " | formation | \n", "geometry | \n", "X | \n", "Y | \n", "Z | \n", "
|---|---|---|---|---|---|
| 0 | \n", "Sand1 | \n", "POINT (10.59347 276.73371) | \n", "10.59 | \n", "276.73 | \n", "359.04 | \n", "
| 1 | \n", "Sand1 | \n", "POINT (17.13494 289.08982) | \n", "17.13 | \n", "289.09 | \n", "364.28 | \n", "
| 2 | \n", "Sand1 | \n", "POINT (19.15013 293.31349) | \n", "19.15 | \n", "293.31 | \n", "364.99 | \n", "
| 3 | \n", "Sand1 | \n", "POINT (27.79512 310.57169) | \n", "27.80 | \n", "310.57 | \n", "372.81 | \n", "
| 4 | \n", "Sand1 | \n", "POINT (34.41735 324.13919) | \n", "34.42 | \n", "324.14 | \n", "377.43 | \n", "
| \n", " | formation | \n", "geometry | \n", "X | \n", "Y | \n", "Z | \n", "
|---|---|---|---|---|---|
| 0 | \n", "Sand1 | \n", "POINT (0.25633 264.86215) | \n", "0.26 | \n", "264.86 | \n", "353.97 | \n", "
| 1 | \n", "Sand1 | \n", "POINT (10.59347 276.73371) | \n", "10.59 | \n", "276.73 | \n", "359.04 | \n", "
| 2 | \n", "Sand1 | \n", "POINT (17.13494 289.08982) | \n", "17.13 | \n", "289.09 | \n", "364.28 | \n", "
| 3 | \n", "Sand1 | \n", "POINT (19.15013 293.31349) | \n", "19.15 | \n", "293.31 | \n", "364.99 | \n", "
| 4 | \n", "Sand1 | \n", "POINT (27.79512 310.57169) | \n", "27.80 | \n", "310.57 | \n", "372.81 | \n", "
| \n", " | \n", " | id | \n", "formation | \n", "geometry | \n", "points | \n", "X | \n", "Y | \n", "Z | \n", "
|---|---|---|---|---|---|---|---|---|
| 0 | \n", "0 | \n", "None | \n", "Sand1 | \n", "POINT (0.25633 264.86215) | \n", "[0.256327195431048, 264.86214748436396] | \n", "0.26 | \n", "264.86 | \n", "353.97 | \n", "
| 0 | \n", "None | \n", "Sand1 | \n", "POINT (10.59347 276.73371) | \n", "[10.59346813871597, 276.73370778641777] | \n", "10.59 | \n", "276.73 | \n", "359.04 | \n", "|
| 0 | \n", "None | \n", "Sand1 | \n", "POINT (17.13494 289.08982) | \n", "[17.134940141888464, 289.089821570188] | \n", "17.13 | \n", "289.09 | \n", "364.28 | \n", "|
| 0 | \n", "None | \n", "Sand1 | \n", "POINT (19.15013 293.31349) | \n", "[19.150128045807676, 293.313485355882] | \n", "19.15 | \n", "293.31 | \n", "364.99 | \n", "|
| 0 | \n", "None | \n", "Sand1 | \n", "POINT (27.79512 310.57169) | \n", "[27.79511673965105, 310.571692592952] | \n", "27.80 | \n", "310.57 | \n", "372.81 | \n", "
| \n", " | level_0 | \n", "level_1 | \n", "id | \n", "formation | \n", "geometry | \n", "points | \n", "X | \n", "Y | \n", "Z | \n", "
|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "0 | \n", "0 | \n", "None | \n", "Sand1 | \n", "POINT (0.25633 264.86215) | \n", "[0.256327195431048, 264.86214748436396] | \n", "0.26 | \n", "264.86 | \n", "353.97 | \n", "
| 1 | \n", "0 | \n", "0 | \n", "None | \n", "Sand1 | \n", "POINT (10.59347 276.73371) | \n", "[10.59346813871597, 276.73370778641777] | \n", "10.59 | \n", "276.73 | \n", "359.04 | \n", "
| 2 | \n", "0 | \n", "0 | \n", "None | \n", "Sand1 | \n", "POINT (17.13494 289.08982) | \n", "[17.134940141888464, 289.089821570188] | \n", "17.13 | \n", "289.09 | \n", "364.28 | \n", "
| 3 | \n", "0 | \n", "0 | \n", "None | \n", "Sand1 | \n", "POINT (19.15013 293.31349) | \n", "[19.150128045807676, 293.313485355882] | \n", "19.15 | \n", "293.31 | \n", "364.99 | \n", "
| 4 | \n", "0 | \n", "0 | \n", "None | \n", "Sand1 | \n", "POINT (27.79512 310.57169) | \n", "[27.79511673965105, 310.571692592952] | \n", "27.80 | \n", "310.57 | \n", "372.81 | \n", "