Skip to content

Commit

Permalink
Merge branch '2-vector-cf' into main
Browse files Browse the repository at this point in the history
  • Loading branch information
cforgaci authored Feb 14, 2024
2 parents ccb229e + 462c0ae commit 6ba6916
Show file tree
Hide file tree
Showing 9 changed files with 193 additions and 0 deletions.
1 change: 1 addition & 0 deletions config.yaml
Original file line number Diff line number Diff line change
Expand Up @@ -63,6 +63,7 @@ episodes:
- 02-data-structures.Rmd
- 03-explore-data.Rmd
- 04-intro-to-visualisation.Rmd
- 08-intro-to-geospatial-concepts.Rmd
- 09-open-and-plot-vector-layers.Rmd
- 10-explore-and-plot-by-vector-layer-attributes.Rmd
- 11-plot-multiple-shape-files.Rmd
Expand Down
103 changes: 103 additions & 0 deletions episodes/08-intro-to-geospatial-concepts.Rmd
Original file line number Diff line number Diff line change
@@ -0,0 +1,103 @@
---
title: 'Introduction to Geospatial Concepts'
teaching: 10
exercises: 2
bibliography: references.bib
zotero: true
---

::: questions
- How ...
:::

::: objectives
- Describe ...
:::

## The shape of the Earth

The shape of the Earth is approximately a sphere which is slightly wider than it is tall, and which is called **ellipsoid**. The true shape of the Earth is an irregular ellipsoid, the so-called **geoid** (@fig-earth).

![The shape of the Earth. Source: United Nations Statistics Division and International Cartographic Association (2012a)](https://unstats.un.org/unsd/geoinfo/ungegn/docs/_data_icacourses/_ImagesModules/_selfstudy/S06_images/S06_03_a00.jpg){#fig-earth}

Check warning on line 21 in episodes/08-intro-to-geospatial-concepts.Rmd

View workflow job for this annotation

GitHub Actions / Build Full Site

[image missing alt-text]: https://unstats.un.org/unsd/geoinfo/ungegn/docs/_data_icacourses/_ImagesModules/_selfstudy/S06_images/S06_03_a00.jpg

The most common and basic representation of the position of points on the Earth is the combination of the **geographical latitude and longitude**.

![Geographical latitude and longitude. Source: van der Marel (2014).](fig/latlon.png)

Check warning on line 25 in episodes/08-intro-to-geospatial-concepts.Rmd

View workflow job for this annotation

GitHub Actions / Build Full Site

[image missing alt-text]: fig/latlon.png

**Meridians** are **vertical** circles with constant longitude, called **great circles**, which run from the North Pole to the South Pole. **Parallels** are **horizontal** circles with constant latitude, which are called **small circles**. Only the equator (the largest parallel) is also a great circle.

The black lines in @fig-latlon show the equator and the prime meridian running through Greenwich, with latitude and longitude labels. The red dotted lines show the meridian and parallel running through Karachi, Pakistan (25°45’N, 67°01’E).

## Map projection

**Map projection** is a systematic transformation of the latitudes and longitudes of locations on the surface of an ellipsoid into locations on a plane. It is a transformation of the three-dimensional Earth’s surface into its two-dimensional representation on a sheet of paper or computer screen (see @fig-projection for a comparison with flattening of an orange peel).

![Map projection represented as flattening an orange peel. Source: Data Carpentry (2023)](https://datacarpentry.org/organization-geospatial/fig/orange-peel-earth.jpg){#fig-projection}

Check warning on line 35 in episodes/08-intro-to-geospatial-concepts.Rmd

View workflow job for this annotation

GitHub Actions / Build Full Site

[image missing alt-text]: https://datacarpentry.org/organization-geospatial/fig/orange-peel-earth.jpg

Many different map projections are in use for different purposes. Generally, they can be categorised into the following groups: cylindrical, conic, and azimuthal (see @fig-projections).

![Cylindrical, conic, and azimuthal map projections. Source: Knippers (2009)](https://kartoweb.itc.nl/geometrics/Bitmaps/Intro%201.9a.gif){#fig-projections}

Check warning on line 39 in episodes/08-intro-to-geospatial-concepts.Rmd

View workflow job for this annotation

GitHub Actions / Build Full Site

[image missing alt-text]: https://kartoweb.itc.nl/geometrics/Bitmaps/Intro%201.9a.gif

Each map projection introduces a **distortion** in geometrical elements – **distance**, **angle**, and **area**. Depending on which of these geometrical elements are more relevant for a specific map, we can choose an appropriate map projection. **Conformal projections** are the best for preserving angles between any two curves; **equal area (equivalent) projections** preserve the area or scale; **equal distance (conventional) projections** are the best for preserving distances.

## Coordinate reference systems (CRS)

A **coordinate reference system (CRS)** is a coordinate-based local, regional or global system for locating geographical entities, which uses a specific map projection. It defines how the two-dimensional, projected map relates to real places on the Earth.

All coordinate reference systems are included in a public registry called the **EPSG Geodetic Parameter Dataset (EPSG registry)**, initiated in 1985 by a member of the European Petroleum Survey Group (EPSG). Each CRS has a unique **EPSG code**, which makes it possible to easily identify them among the large number of CRS. This is particularly important for transforming spatial data from one CRS to another.

Some of the most commonly used CRS in the Netherlands are the following:

- **World Geodetic System 1984 (WGS84)** is the best known global reference system (EPSG:4326).
- **European Terrestrial Reference System 1989 (ETRS89)** is the standard coordinate system for Europe (EPSG:4258).
- The most popular projected CRS in the Netherlands is ‘Stelsel van de Rijksdriehoeksmeting (RD)’ registered in EPSG as **Amersfoort / RD New (EPSG:28992)**.

The main parameters of each CRS are the following:

- **Datum** is a model of the shape of the Earth, which specifies how a coordinate system is linked to the Earth, e.g. how to define the origin of the coordinate axis – where (0,0) is. It has angular units (degrees).
- **Projection** is mathematical transformation of the angular measurements on the Earth to linear units (e.g. meters) on a flat surface (paper or a computer screen).
- **Additional parameters**, such as a definition of the centre of the map, are often necessary to create the full CRS.

In this workshop, we use two CRS shown in @tbl-crs.

| | WGS 84 (EPSG:4326) | Amersfoort / RD New (EPSG:28992) |
|-----------------|----------------------------|----------------------------|
| Units | degrees | meters |
| Projection | Geographic (uses latitude and longitude for coordinates) | |
| | Dynamic (relies on a datum which is not plate-fixed) | Static (relies on a datum which is plate-fixed) |
| Celestial body | Earth | |
| Datum | World Geodetic System 1984 ensemble | Amersfoort |
| Method | Lat/long (Geodetic alias) | Oblique Stereographic Alternative |
| Prime meridian | Greenwich | |
| Ellipsoid | | Bessel 1841 |

: Main properties of WGS 84 and Amersfoort / RD New coordinate reference systems {#tbl-crs}

## Map scale

**Map scale** measures the ratio between distance on a map and the corresponding distance on the ground. For example, on a 1:100 000 scale map, 1cm on the map equals 1km (100 000 cm) on the ground. Map scale can be expressed in the following three ways:

- Verbal: 1 centimetre represents 250 meters
- Fraction: 1:25000
- Graphic:

::: callout
# Useful resources

- Campbell, J., Shin, M. E. (2011). Essentials of Geographic Information Systems. Textbooks. 2. <https://digitalcommons.liberty.edu/textbooks/2> (Accessed 22-01-2024)

- Data Carpentry (2023): Introduction to Geospatial Concepts. Coordinate Reference Systems. <https://datacarpentry.org/organization-geospatial/03-crs.html> (Accessed 22-01-2024)

- GeoRepository (2024): EPSG Geodetic Parameter Dataset <https://epsg.org/home.html> (Accessed 22-01-2024)

- Klokan Technologies GmbH (2022) <https://epsg.io/> (Accessed 22-01-2024)

- United Nations Statistics Division and International Cartographic Association (2012b): UNGEGN-ICA webcourse on Toponymy. <https://unstats.un.org/unsd/geoinfo/ungegn/docs/_data_icacourses/2012_Home.html> (Accessed 22-01-2024)
:::

::: keypoints
- Use `.md` files for episodes when you want static content
- Use `.Rmd` files for episodes when you need to generate output
- Run `sandpaper::check_lesson()` to identify any issues with your lesson
- Run `sandpaper::build_lesson()` to preview your lesson locally
:::
Binary file added episodes/fig/CHM_TUD.tiff
Binary file not shown.
Binary file added episodes/fig/latlon.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file added episodes/fig/scalebar.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
8 changes: 8 additions & 0 deletions episodes/references.bib
Original file line number Diff line number Diff line change
@@ -0,0 +1,8 @@

@article{vandermarel2014,
title = {Reference Systems for Surveying and Mapping. Lecture Notes},
author = {van der Marel, Hans},
year = {2014},
date = {2014},
langid = {en}
}
20 changes: 20 additions & 0 deletions renv.lock
Original file line number Diff line number Diff line change
@@ -0,0 +1,20 @@
{
"R": {
"Version": "4.3.0",
"Repositories": [
{
"Name": "CRAN",
"URL": "https://cran.rstudio.com"
}
]
},
"Packages": {
"renv": {
"Package": "renv",
"Version": "1.0.3",
"OS_type": null,
"Repository": "CRAN",
"Source": "Repository"
}
}
}
7 changes: 7 additions & 0 deletions renv/.gitignore
Original file line number Diff line number Diff line change
@@ -0,0 +1,7 @@
library/
local/
cellar/
lock/
python/
sandbox/
staging/
54 changes: 54 additions & 0 deletions renv/profiles/lesson-requirements/renv.lock
Original file line number Diff line number Diff line change
Expand Up @@ -104,6 +104,16 @@
],
"Hash": "cad6cf7f1d5f6e906700b9d3e718c796"
},
"assertthat": {
"Package": "assertthat",
"Version": "0.2.1",
"Source": "Repository",
"Repository": "CRAN",
"Requirements": [
"tools"
],
"Hash": "50c838a310445e954bc13f26f26a6ecf"
},
"backports": {
"Package": "backports",
"Version": "1.4.1",
Expand Down Expand Up @@ -479,6 +489,30 @@
],
"Hash": "bb0eec2fe32e88d9e2836c2f73ea2077"
},
"emo": {
"Package": "emo",
"Version": "0.0.0.9000",
"Source": "GitHub",
"RemoteType": "github",
"RemoteHost": "api.github.com",
"RemoteRepo": "emo",
"RemoteUsername": "hadley",
"RemoteRef": "HEAD",
"RemoteSha": "3f03b11491ce3d6fc5601e210927eff73bf8e350",
"Requirements": [
"R",
"assertthat",
"crayon",
"glue",
"lubridate",
"magrittr",
"purrr",
"rlang",
"stringr",
"utils"
],
"Hash": "10f2391a926b0087a671ed9965ea30a3"
},
"evaluate": {
"Package": "evaluate",
"Version": "0.23",
Expand Down Expand Up @@ -1041,6 +1075,19 @@
],
"Hash": "faa2193fdec94a45b4390aefc1280a62"
},
"prismatic": {
"Package": "prismatic",
"Version": "1.1.1",
"Source": "Repository",
"Repository": "CRAN",
"Requirements": [
"R",
"farver",
"grDevices",
"graphics"
],
"Hash": "faa2193fdec94a45b4390aefc1280a62"
},
"processx": {
"Package": "processx",
"Version": "3.8.3",
Expand Down Expand Up @@ -1918,6 +1965,13 @@
"Repository": "CRAN",
"Hash": "c6abfa47a46d281a7d5159d0a8891e88"
},
"whisker": {
"Package": "whisker",
"Version": "0.4.1",
"Source": "Repository",
"Repository": "CRAN",
"Hash": "c6abfa47a46d281a7d5159d0a8891e88"
},
"withr": {
"Package": "withr",
"Version": "3.0.0",
Expand Down

0 comments on commit 6ba6916

Please sign in to comment.