add more exercises to part 4

Author: eemilhaa

kurssimateriaali/4_paikkatietoanalyysi.ipynb

@@ -65,9 +65,7 @@
   {
    "cell_type": "markdown",
    "id": "3cfd0549-f677-43b3-9eda-0182a8a27ef0",
-   "metadata": {
-    "jp-MarkdownHeadingCollapsed": true
-   },
+   "metadata": {},
    "source": [
     "## Ratkaisu"
    ]
@@ -199,7 +197,7 @@
    "id": "1f6c04d5-d59c-4f35-88e2-3b533d7b6678",
    "metadata": {},
    "source": [
-    "Projisoi salamahavainnot samaan koordinaattijärjestelmään, kuin kohta käytettävä kunta-aineisto (epsg=3067). Huomaa, että sijainnit ovat nyt asteita, eli emme voi suoraan asettaa ETRS89 / TM35FIN koordinaattijärjestelmää.\n",
+    "Projisoi salamahavainnot samaan koordinaattijärjestelmään, kuin kohta käytettävä kunta-aineisto (`epsg=3067`). Huomaa, että sijainnit ovat nyt asteita - aineisto on wgs 84 -koordinaatistossa. Emme siis voi suoraan asettaa ETRS89 / TM35FIN koordinaattijärjestelmää.\n",
     "\n",
     "1. Aseta ensin `geo_lightnings`-GeoDataFrame wgs 84 koordinaattijärjestelmään (`set_crs`). Epsg-koodi on 4326.\n",
     "2. Uudelleenprojisoi (`to_crs`) vasta kun koordinaattijärjestelmä on asetettu. Ylikirjoita `geo_lightnings`-muuttuja uudella epsg 3067 -versiolla."
@@ -317,11 +315,37 @@
    "id": "14baae2c-591b-41d9-a7dd-009f0affce1f",
    "metadata": {},
    "source": [
-    "## Liitos sijainnin perusteella\n",
+    "## Harjoitus - liitos sijainnin perusteella\n",
     "\n",
     "Nyt voimme edetä tekemään sijainteihin perustuvan liitoksen salamahavaintojen ja kuntien välillä. Tarkoituksena on, että saisimme jokaiseen salamaan tiedon siitä, minkä kunnan alueelle se iski.\n",
     "\n",
-    "Tämä kannattaa toteuttaa `sjoin` eli spatial join -metodilla. Yhdistetään siis salamahavaintoihin tieto kunnista, perusteena, että havainto on kunnan alueella (\"within\"). Liitoksen tyypiksi otetaan \"left\", sillä haluamme tuloksena vain salamahavainnot."
+    "Tämä kannattaa toteuttaa `sjoin`, eli spatial join, -funktiolla. Yhdistetään siis salamahavaintoihin tieto kunnista, perusteena, että havainto on kunnan alueella.\n",
+    "\n",
+    "1. Käytä geopandasin `sjoin`-funktiota (`gpd.sjoin`)\n",
+    "2. Valitse parametreihin `left_df` ja `right_df` halutut GeoDataFramet, eli: `left_df=geo_lightnings` ja `right_df=municipalities`.\n",
+    "3. Haluamme tiedon sen perusteella, minkä kunnan sisällä havainto on. Valitse siis predikaatiksi, eli `predicate`-parametrin arvoksi, `\"within\"`.\n",
+    "4. Haluamme tuloksena salamahavainnot, joihin on lisättynä tieto kunnasta, jonne salama on iskenyt. Annetaan siis `how`-parametrille arvoksi `\"left\"`: Näin tuloksena on vain salamahavainnot.\n",
+    "5. Kun liitos näyttää toimivan, tallenna tuloksena syntyvä GeoDataFrame muuttujaan `lightnings_with_municipality`."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "ccb44cf2-38a7-4d09-9f6a-2fab93869f3b",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Kirjoita ratkaisu\n"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "e30768ff-f50c-491e-b4c5-3078203b9a5e",
+   "metadata": {
+    "jp-MarkdownHeadingCollapsed": true
+   },
+   "source": [
+    "## Ratkaisu"
    ]
   },
   {
@@ -337,8 +361,9 @@
    },
    "outputs": [],
    "source": [
-    "lightnings_with_municipality = geo_lightnings.sjoin(\n",
-    "    municipalities,\n",
+    "lightnings_with_municipality = gpd.sjoin(\n",
+    "    left_df=geo_lightnings,\n",
+    "    right_df=municipalities,\n",
     "    predicate=\"within\",\n",
     "    how=\"left\",\n",
     ")"
@@ -479,8 +504,9 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "municipalities = municipalities.merge(\n",
-    "    lightning_counts,\n",
+    "municipalities = pd.merge(\n",
+    "    left=municipalities,\n",
+    "    right=lightning_counts,\n",
     "    left_on=\"nimi\",\n",
     "    right_on=\"nimi\",\n",
     "    how=\"left\",\n",
@@ -895,7 +921,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.12.7"
+   "version": "3.12.8"
   }
  },
  "nbformat": 4,