https://github.com/kadyb/geomorph_clustering
High-resolution geomorphological clustering of Poland
https://github.com/kadyb/geomorph_clustering
clustering geomorphology geomorphometry poland r spatial
Last synced: 7 months ago
JSON representation
High-resolution geomorphological clustering of Poland
- Host: GitHub
- URL: https://github.com/kadyb/geomorph_clustering
- Owner: kadyb
- License: mit
- Created: 2021-08-18T19:27:54.000Z (about 4 years ago)
- Default Branch: main
- Last Pushed: 2022-07-22T12:05:05.000Z (about 3 years ago)
- Last Synced: 2024-01-26T09:18:33.611Z (over 1 year ago)
- Topics: clustering, geomorphology, geomorphometry, poland, r, spatial
- Language: R
- Homepage:
- Size: 3.39 MB
- Stars: 3
- Watchers: 1
- Forks: 0
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
- License: LICENSE
- Citation: CITATION.cff
Awesome Lists containing this project
README
# Gomorphological clustering
This repository contains the code for “[*Toward geomorphometry of plains - Country-level unsupervised classification of low-relief areas (Poland)*](https://www.sciencedirect.com/science/article/pii/S0169555X22002665)” article.
## Reproduction
1. Open the `geomorph_clustering.Rproj` project file in [RStudio](https://rstudio.com/).
2. Run `01_extract.R` to generate a random sample from the rasters.
3. Run `02_clustering.R` to train the data transformer and Gaussian mixture models (GMM).
Finally, the model with the highest BIC value is selected.
4. Run `03A_predict_lowres.R` to perform spatial clustering on low resolution (downsampled) rasters.
The result will be maps with clusters (geomorphological units) and uncertainty.
5. Run `03B_predict_highres.R` to perform high resolution spatial clustering.
This operation is performed for **all** data in blocks of the specified size.
## Dataset
The final high-resolution maps with geomorphological units and uncertainty, Gaussian mixture model and data transformation model, and low-resolution rasters (500 m) with geomorphometric variables are available in the Zenodo repository: https://zenodo.org/record/6415362.
Note that in the project, rasters with 30 m resolution were used as input data.
The "*input_rasters*" folder in the archive should be renamed to "*rasters*" to work with the shared code.