computeFeldRegion function

Retrieving Feldpausch regions