Rewritten using new method based on numpy image processing
Add land area