Description:

This code written by Brian Bue takes orthocorreced methane matched filter results and the locations of plumes manually identified from a spreadsheet and generates a number of data products. Most notable are Integrated Methane Enhancements (IME) and plume lengths (fetch) calculations saved to .txt files, GeoTIFFs and GeoJSONs that are used as part of the Methane Source Finder project. There are a number of input variables that can be adjusted, including -p (concentration threshold in ppmm), -m (merge distance in m, used to isolate the plume), and -f (maximum fetch in meters). The file path to the matched filter results are provided within the code.


To run:

/shared/anaconda2/bin/python2.7 /lustre/athorpe/ch4/compute_ime_040318_mf/compute_ime.py /lustre/athorpe/ch4/compute_ime_040318_mf/input_spreadsheet/FA16_Master_Source_List_IME_GasStorage.xlsx "FA16 plumes IME" -o /lustre/athorpe/ch4/compute_ime_040318_mf/output/FA16_Gas_Storage_1000ppm_150fetch_5_10_20merge --calids v1n2 v2m4 -p 1000 -m 10 20 50 -f 150




Some addtional description from Brian:

bbue@gibson:srcfinder$ python compute_ime/compute_ime.py --help
usage: compute_ime.py [-h] [-m MERGEDISTS [MERGEDISTS ...]] [-c] [-o OUTDIR]
                      [-f FETCHMAX] [-p PPMMTHR] [--nodetfilt] [-i inputfile]
                      [-l listfile] [--randomize] [--filter filter]
                      xlsfile sheetname

CMF IME Computation + Product Generation

positional arguments:
  xlsfile               Source master list spreadsheet
  sheetname             Source master list spreadsheet sheet name

optional arguments:
  -h, --help            show this help message and exit
  -m MERGEDISTS [MERGEDISTS ...], --mergedists MERGEDISTS [MERGEDISTS ...]
                        List of one or more distances (in meters) to merge
                        neighboring blobs (default=[10, 20, 50])
  -c, --clobber         Clobber (overwrite) existing outputs
  -o OUTDIR, --outdir OUTDIR
                        Save filtered detection image to OUTDIR
  -f FETCHMAX, --fetchmax FETCHMAX
                        Maximum fetch value (default=unlimited)
  -p PPMMTHR, --ppmmthr PPMMTHR
                        Minimum ppmm threshold (default=500.00)
  --nodetfilt           Do not filter out small or weakly connected detections
  -i inputfile, --inputfile inputfile
                        CMF image (assumed listed in spreadsheet)
  -l listfile, --listfile listfile
                        List of CMF images to process (subset of images in
                        provided spreadsheet)
  --randomize           Randomize spreadsheet processing order (by source id)
  --filter filter       Compute IME for spreadsheet rows matching filters of
                        the form "spreadsheet_column:filter_string (e.g.,
                        "flightline:ang201605")"