NetCDF to COG Conversion Tool

Convert Open Datacube NetCDF files to Cloud Optimised GeoTIFFs.

It runs prior to uploading from the NCI file systems to AWS S3.

Use dea-cogger convert to COG:

Installation

conda create -n cogger datacube=1.8.3 python=3.8 gdal boto3 boltons structlog mpi4py colorama requests chardet
conda activate cogger
pip install requests-aws4auth
pip install -e .

Usage

Example commands for calculating which datasets require conversion to COG by downloading a subset of an S3 Inventory and comparing it with an ODC Index.

dea-cogger save-s3-inventory -p ls7_fc_albers -o tmp/ 
dea-cogger generate-work-list --product-name ls7_fc_albers -o tmp

Running Conversion in parallel on Gadi

Example PBS submission script to run in parallel on Gadi.

#!/bin/bash
#PBS -l wd,walltime=5:00:00,mem=190GB,ncpus=48,jobfs=1GB
#PBS -P v10
#PBS -q normal
#PBS -l storage=gdata/v10+gdata/fk4+gdata/rs0+gdata/if87
#PBS -W umask=33
#PBS -N cog_ls8_nbar_albers 

module load dea
module load openmpi/3.1.4

mpirun --tag-output dea-cogger mpi-convert --product-name "{{params.product}}" --output-dir "{{work_dir}}/out/" ls8_nbar_albers_file_list.txt

Example configuration file

    products:
      product_name:
       prefix: WOfS/WOFLs/v2.1.5/combined
       name_template: x_{x}/y_{y}/{time:%Y}/{time:%m}/{time:%d}/LS_WATER_3577_{x}_{y}_{time:%Y%m%d%H%M%S%f}
       predictor: 2
       no_overviews: ["source", "observed"]
       default_resampling: average
       white_list: None
       black_list: None

• product_name: A unique user defined string (required)
• prefix: Define the cogs folder structure and name (required)
• name_template: Define how to decipher the input file names (required)
• default_resampling: Define the resampling method of pyramid view (default: average)
• predictor: Define the predictor in COG convert (default: 2)
• no_overviews: A list of keywords of bands which don't require resampling (optional)
• white_list: A list of keywords of bands to be converted (optional)
• black_list: A list of keywords of bands excluded in cog convert (optional)

Note: no_overviews contains the key words of the band names which one doesn't want to generate overviews. This element cannot be used with other products as this 'cause it will match as source'. For most products, this element is not needed. So far, only fractional cover percentile use this.

Command: save-s3-inventory

Scan through S3 bucket for the specified product and fetch the file path of the uploaded files. Save those file into a pickle file for further processing. Uses a configuration file to define the file naming schema.

Command: generate-work-list

Compares ODC URI's against an S3 bucket and writes the list of datasets for COG conversion into a file.

Uses a configuration file to define the file naming schema.

Command: mpi-convert

Bulk COG Convert netcdf files to COG format using MPI tool. Iterate over the file list and assign MPI worker for processing. Split the input file by the number of workers, each MPI worker completes every nth task. Also, detect and fail early if not using full resources in an MPI job.

Reads the file naming schema from the configuration file.

Command: verify

Verify converted GeoTIFF files are (Geo)TIFF with cloud optimized compatible structure. Mandatory Requirement: validate_cloud_optimized_geotiff.py gdal file.

COG Creation Settings (What to set for predictor and resampling)

Predictor