Skip to content
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.

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

Level 3 Islands #13

Open
alrichardbollans opened this issue Nov 28, 2024 · 0 comments
Open

Level 3 Islands #13

alrichardbollans opened this issue Nov 28, 2024 · 0 comments

Comments

@alrichardbollans
Copy link

I'm interested in the proportion of Level-3 regions that are considered islands. Firstly, I've determined there are 369 level 3 regions in total using:

    import cartopy.io.shapereader as shpreader
    tdwg3_shp = shpreader.Reader(
        os.path.join('inputs', 'wgsrpd-master', 'level3', 'level3.shp'))

    level3_codes = []
    for c in tdwg3_shp.records():
        tdwg_code = c.attributes['LEVEL3_COD']
        level3_codes.append(tdwg_code)
    level3_codes_set = set(level3_codes)
    print(len(level3_codes_set))

but now I'm wondering if there's any reasonable way to determine which of these regions can be considered islands or groups of islands?

Obviously if 'Is.' is in the LEVEL3_NAM, this is probably a good start, but there are lots that for which this isn't the case e.g. Kazan-retto or Maldives

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant