Python command line tools

This module contains several command line tools that might be useful for spatial data analysis and modelling. The development of these tools aims for the ease of use. Most of the tools are built upon the following packages such that installing them are necessary:

• numpy
• pandas
• geopandas
• rasterio
• rasterstats

To install the required package, open the command prompt or Anaconda and use pip install [package name]. Alternatively, conda install [package name] can also be used in the Anaconda.

Zonal Statistics

zonal_statistics.py is a script to extract zonal statistics from raster input based on the list of features provided. Copy the script to the working directory to use it. Type python zonal_statistics.py --help to know how to use this script.

Dynamic World Wrapper

dworld_wrapper.py can be used to acquire Dynamic World land use/land cover (LULC) class for a specific area. Tiles of annually aggregated Dynamic World data (from 2016-2021) are available at the WorldPop drive. The spatial resolution of the rasters is 100 m. The dworld_wrapper.py script helps to merge and clip the tiles.

Usage: python dworld_wrapper.py -y year -b band -c boundary -o output.tif

• Available years: 2016, 2017, 2018, 2019, 2020, 2021
• Available bands: water, trees, grass, flooded_vegetation, crops, shrub_and_scrub, built, bare, snow_and_ice
• Clipping boundary: can either be a shapefile or coordinate boundary (xmin, xmax, ymin, ymax)

Example:

• python dworld_wrapper.py -y 2020 -b crops,built -c 100,105,0,5 -o sumatra.tif will produce a raster with two bands (crops, built) clipped on the rectangular boundary between 100-105 longitude and 0-5 latitude. The output will be binary integer, 1 represents pixel categorised as a certain (crops, built) class.
• python dworld_wrapper.py -y 2021 -b built -c gadm41_IND_0.shp -o india.tif will produce a raster of artificially built area in India (defined by gadm41_IDN_0.shp.