Skip to content

Tutorial 2: Change "gdf" to "gpd" in variable name | Show output of ".head()" and ".crs" explicitly #39

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
yvonnefroehlich merged 4 commits into from
Dec 31, 2024
Merged

Tutorial 2: Change "gdf" to "gpd" in variable name | Show output of ".head()" and ".crs" explicitly #39

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
efbb705
Use 'gdf' instead of 'gpd' to avoid confusions | Split cells to show …
yvonnefroehlich Dec 26, 2024
8eb360a
Fix typo
yvonnefroehlich Dec 26, 2024
f22680d
Merge remote-tracking branch 'origin/main' into split-cells-gdf-crs
yvonnefroehlich Dec 31, 2024
72dbdd0
Add a bit more explanation
yvonnefroehlich Dec 31, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion book/tut01_firstfigure.ipynb
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -522,7 +522,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.12.6"
"version": "3.12.7"
},
"panel-cell-order": [
"eb4b43cc-dc76-45d1-a370-7266bd943910",
Expand Down
97 changes: 79 additions & 18 deletions book/tut02_spe_pd_gpd.ipynb
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -258,33 +258,94 @@
"metadata": {},
"outputs": [],
"source": [
"gpd_rivers_org = gpd.read_file(\n",
"gdf_rivers_org = gpd.read_file(\n",
" \"https://www.eea.europa.eu/data-and-maps/data/wise-large-rivers-and-large-lakes/\"\n",
" + \"zipped-shapefile-with-wise-large-rivers-vector-line/zipped-shapefile-with-wise-large-rivers-vector-line/\"\n",
" + \"at_download/file/wise_large_rivers.zip\"\n",
")\n",
"# gpd_rivers_org = pd.read_file(\"wise_large_rivers.zip\")\n",
"gpd_rivers_org.head()"
"# gdf_rivers_org = gpd.read_file(\"wise_large_rivers.zip\")"
]
},
{
"cell_type": "markdown",
"id": "f8db1757-f80d-4cc3-b650-83109eecc390",
"metadata": {},
"source": [
"Have a look at the values in the geometry column. The coordinates are currently not given in the geographic coordinate reference system (longitude/latitude) and have to be converted. This can be done directly with `GeoPandas`."
"Have a look at the data and especially at the values in the geometry column:"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "cc67a58c-45f3-4819-9df4-04e5aa9d06f7",
"id": "926955ee-3859-468b-9b41-0398ed6c7b99",
"metadata": {},
"outputs": [],
"source": [
"gpd_rivers_org.crs\n",
"gpd_rivers = gpd_rivers_org.to_crs(\"EPSG:4326\")\n",
"gpd_rivers.head()"
"gdf_rivers_org.head()"
]
},
{
"cell_type": "markdown",
"id": "b868d284-23ff-4607-a950-de47da71f7cd",
"metadata": {},
"source": [
"The coordinates are currently not given in the geographic coordinate reference system (longitude/latitude):"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "f033930d-fe1c-4883-a27c-6a768bb29598",
"metadata": {},
"outputs": [],
"source": [
"gdf_rivers_org.crs"
]
},
{
"cell_type": "markdown",
"id": "59edb4bf-9466-41a9-9cf6-f2ae3d7a5a0d",
"metadata": {},
"source": [
"Thus, they have to be converted which can be done directly with `GeoPandas`."
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "c6a07d4d-ec70-4624-836b-41f713218f11",
"metadata": {},
"outputs": [],
"source": [
"gdf_rivers = gdf_rivers_org.to_crs(\"EPSG:4326\")"
]
},
{
"cell_type": "markdown",
"id": "cd9b9f2b-d15b-4336-b2ee-4370bb6ffb25",
"metadata": {},
"source": [
"Again have a look at the coordinate system and the data:"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "06124f49-4dc7-4518-8ceb-b85257a30920",
"metadata": {},
"outputs": [],
"source": [
"gdf_rivers.crs"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "537e141d-6f64-4b99-9e30-75e1c217d3c5",
"metadata": {},
"outputs": [],
"source": [
"gdf_rivers.head()"
]
},
{
Expand Down Expand Up @@ -314,7 +375,7 @@
" frame=True,\n",
")\n",
"\n",
"fig.plot(data=gpd_rivers, pen=\"0.5p,steelblue,solid\")\n",
"fig.plot(data=gdf_rivers, pen=\"0.5p,steelblue,solid\")\n",
"\n",
"fig.show(dpi=img_dpi)"
]
Expand All @@ -341,8 +402,8 @@
"# Split the dataset into two subsets of shorter and longer rivers\n",
"# Feel free to play around with the limit\n",
"len_limit = 700000 # in meters\n",
"gpd_rivers_short = gpd_rivers[gpd_rivers[\"Shape_Leng\"] < len_limit]\n",
"gpd_rivers_long = gpd_rivers[gpd_rivers[\"Shape_Leng\"] > len_limit]\n",
"gdf_rivers_short = gdf_rivers[gdf_rivers[\"Shape_Leng\"] < len_limit]\n",
"gdf_rivers_long = gdf_rivers[gdf_rivers[\"Shape_Leng\"] > len_limit]\n",
"\n",
"\n",
"fig = pygmt.Figure()\n",
Expand All @@ -357,11 +418,11 @@
"\n",
"# Plot the subsets differently and specify the text for the legend entries\n",
"fig.plot(\n",
" data=gpd_rivers_short,\n",
" data=gdf_rivers_short,\n",
" pen=\"0.5p,orange\",\n",
" label=f\"shorter {len_limit} m+Hriver length+f9p\",\n",
")\n",
"fig.plot(data=gpd_rivers_long, pen=\"0.5p,darkred\", label=f\"longer {len_limit} m\")\n",
"fig.plot(data=gdf_rivers_long, pen=\"0.5p,darkred\", label=f\"longer {len_limit} m\")\n",
"\n",
"# Place the legend at the Top Left corner with an offset of 0.1 centimeters from the map frame\n",
"# Add a box with semi-transparent (@30) white fill (+g) and a 0.1-points thick gray outline (+p)\n",
Expand Down Expand Up @@ -520,14 +581,14 @@
")\n",
"\n",
"pygmt.makecpt(\n",
" cmap=\"SCM/oslo\", series=[gpd_rivers.Shape_Leng.min(), 1500000], reverse=True\n",
" cmap=\"SCM/oslo\", series=[gdf_rivers.Shape_Leng.min(), 1500000], reverse=True\n",
")\n",
"fig.colorbar(frame=[\"x+lriver length\", \"y+lm\"], position=\"+ef0.2c\")\n",
"\n",
"for i_river in range(len(gpd_rivers)):\n",
"for i_river in range(len(gdf_rivers)):\n",
" fig.plot(\n",
" data=gpd_rivers[gpd_rivers.index == i_river],\n",
" zvalue=gpd_rivers.loc[i_river, \"Shape_Leng\"],\n",
" data=gdf_rivers[gdf_rivers.index == i_river],\n",
" zvalue=gdf_rivers.loc[i_river, \"Shape_Leng\"],\n",
" pen=\"0.5p\",\n",
" cmap=True,\n",
" )\n",
Expand All @@ -552,7 +613,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.12.6"
"version": "3.12.7"
}
},
"nbformat": 4,
Expand Down