1 A few checks of the habitatmap_stdized data source

result <- read_habitatmap_stdized("../../data")

Do all polygons from the sf object have a unique ID? (Beware, this calculation takes approx. 2-3 minutes; it has been disabled)

result$habitatmap_polygons %>% 
  count(polygon_id) %>% 
  filter(n > 1) %>% 
  nrow == 0

Is the number of unique IDs in the dataframe the same as the number of polygons in the sf object?

result$habitatmap_patches %>% 
  distinct(polygon_id) %>% 
  nrow == nrow(result$habitatmap_polygons)
[1] TRUE

Do all polygon IDs from the dataframe coincide with those of the sf object?

result$habitatmap_patches %>% 
  inner_join(result$habitatmap_polygons) %>% 
  nrow == nrow(result$habitatmap_patches)
Column `polygon_id` joining factor and character vector, coercing into character vector
[1] TRUE

Number of polygons with phab = 0:

result$habitatmap_patches %>% 
  filter(phab == 0) %>% 
  nrow()

[1] 755
LS0tCnRpdGxlOiAiSGFuZGxpbmcgdGhlIGhhYml0YXRtYXAiCmRhdGU6ICdgciBwYXN0ZSgiVmVyc2lvbiIsbHVicmlkYXRlOjpub3coKSlgJwpvdXRwdXQ6CiAgaHRtbF9ub3RlYm9vazoKICAgIG51bWJlcl9zZWN0aW9uczogeWVzCiAgICBjb2RlX2ZvbGRpbmc6IHNob3cKICAgIGluY2x1ZGVzOgogICAgICBpbl9oZWFkZXI6IC4uL2hlYWRlci5odG1sCiAgICB0b2M6IHllcwogICAgdG9jX2Zsb2F0OgogICAgICBjb2xsYXBzZWQ6IG5vCiAgICAgIHNtb290aF9zY3JvbGw6IG5vCi0tLQoKYGBge3Igc2V0dXAsIG1lc3NhZ2U9RkFMU0UsIGVjaG89RkFMU0V9Cm9wdGlvbnMoc3RyaW5nc0FzRmFjdG9ycyA9IEZBTFNFKQojIGxpYnJhcnkoc3ApCmxpYnJhcnkoc2YpCmxpYnJhcnkocmFzdGVyKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShuMmtoYWIpCiMgIyBsaWJyYXJ5KHBsb3RseSkKIyBsaWJyYXJ5KHJhc3RlclZpcykKIyBsaWJyYXJ5KHN0YXJzKQojIGxpYnJhcnkodW5pdHMpCiMgbGlicmFyeSh0bWFwKQpsaWJyYXJ5KGtuaXRyKQpvcHRzX2NodW5rJHNldCgKICBlY2hvID0gVFJVRSwKICBkcGkgPSAzMDAKKQpgYGAKCgojIEEgZmV3IGNoZWNrcyBvZiB0aGUgaGFiaXRhdG1hcF9zdGRpemVkIGRhdGEgc291cmNlCgpgYGB7cn0KcmVzdWx0IDwtIHJlYWRfaGFiaXRhdG1hcF9zdGRpemVkKCIuLi8uLi9kYXRhIikKYGBgCgoKRG8gYWxsIHBvbHlnb25zIGZyb20gdGhlIHNmIG9iamVjdCBoYXZlIGEgdW5pcXVlIElEPyAoQmV3YXJlLCB0aGlzIGNhbGN1bGF0aW9uIHRha2VzIGFwcHJveC4gMi0zIG1pbnV0ZXM7IGl0IGhhcyBiZWVuIGRpc2FibGVkKQoKYGBge3IgZXZhbD1GQUxTRX0KcmVzdWx0JGhhYml0YXRtYXBfcG9seWdvbnMgJT4lIAogIGNvdW50KHBvbHlnb25faWQpICU+JSAKICBmaWx0ZXIobiA+IDEpICU+JSAKICBucm93ID09IDAKYGBgCgpJcyB0aGUgbnVtYmVyIG9mIHVuaXF1ZSBJRHMgaW4gdGhlIGRhdGFmcmFtZSB0aGUgc2FtZSBhcyB0aGUgbnVtYmVyIG9mIHBvbHlnb25zIGluIHRoZSBzZiBvYmplY3Q/CgpgYGB7cn0KcmVzdWx0JGhhYml0YXRtYXBfcGF0Y2hlcyAlPiUgCiAgZGlzdGluY3QocG9seWdvbl9pZCkgJT4lIAogIG5yb3cgPT0gbnJvdyhyZXN1bHQkaGFiaXRhdG1hcF9wb2x5Z29ucykKYGBgCgpEbyBhbGwgcG9seWdvbiBJRHMgZnJvbSB0aGUgZGF0YWZyYW1lIGNvaW5jaWRlIHdpdGggdGhvc2Ugb2YgdGhlIHNmIG9iamVjdD8KCmBgYHtyIG1lc3NhZ2U9RkFMU0V9CnJlc3VsdCRoYWJpdGF0bWFwX3BhdGNoZXMgJT4lIAogIGlubmVyX2pvaW4ocmVzdWx0JGhhYml0YXRtYXBfcG9seWdvbnMpICU+JSAKICBucm93ID09IG5yb3cocmVzdWx0JGhhYml0YXRtYXBfcGF0Y2hlcykKYGBgCgpOdW1iZXIgb2YgcG9seWdvbnMgd2l0aCBwaGFiID0gMDoKCmBgYHtyfQpyZXN1bHQkaGFiaXRhdG1hcF9wYXRjaGVzICU+JSAKICBmaWx0ZXIocGhhYiA9PSAwKSAlPiUgCiAgbnJvdygpCmBgYAoK