sedinet_predict1image.py

## Written by Daniel Buscombe,
## MARDA Science
## daniel@mardascience.com

##> Release v1.3 (July 2020)

###===================================================
# import libraries
import sys, getopt, json, os
from numpy import any as npany

USE_GPU = True

if USE_GPU == True:
   ##use the first available GPU
   os.environ['CUDA_VISIBLE_DEVICES'] = '0' #'1'
else:
   ## to use the CPU (not recommended):
   os.environ['CUDA_VISIBLE_DEVICES'] = '-1'

from sedinet_eval import *

#==============================================================
if __name__ == '__main__':

    argv = sys.argv[1:]
    try:
        opts, args = getopt.getopt(argv,"h:c:i:w:1:2:3:4:")
    except getopt.GetoptError:
        print('python sedinet_predict1image.py -c configfile.json -i path/to/image.ext {-w weightsfile.hdf5} OR {-1 weightsfile_batch1.hdf5 -2 weightsfile_batch2.hdf5 -3 weightsfile_batch3.hdf5 -4 weightsfile_batch4.hdf5}')
        sys.exit(2)
    for opt, arg in opts:
        if opt == '-h':
            print('Example usage (single batch / weights file): python sedinet_predict1image.py \
                   -c config/config_mattole.json \
                   -i images/mattole_images/all/DSCN3521c.JPG \
                   -w mattole/res/mattole_simo_batch7_im512_512_2vars_pinball_aug.hdf5')
            print('Example usage (multiple batches / weights files): python sedinet_predict1image.py -c config/config_9percentiles.json \
                   -1 grain_size_global/res/global_9prcs_simo_batch7_im768_9vars_pinball_noaug.hdf5 \
                   -2 grain_size_global/res/global_9prcs_simo_batch12_im768_9vars_pinball_noaug.hdf5 \
                   -3 grain_size_global/res/global_9prcs_simo_batch14_im768_9vars_pinball_noaug.hdf5')
            sys.exit()
        elif opt in ("-c"):
            configfile = arg
        elif opt in ("-w"):
            weights_path = arg
        elif opt in ("-i"):
            image = arg
        elif opt in ("-1"):
            weights_path1 = arg
        elif opt in ("-2"):
            weights_path2 = arg
        elif opt in ("-3"):
            weights_path3 = arg
        elif opt in ("-4"):
            weights_path4 = arg

    if 'image' not in locals():
        if not os.path.isfile(image):
           print("image folder path does not exist ... exiting")
           sys.exit()

    if 'weights_path1' not in locals():
        if not os.path.isfile(os.getcwd()+os.sep+weights_path):
            if not os.path.isfile(weights_path):
               print("Weights path does not exist ... exiting")
               sys.exit()
    else:
        weights_path = []
        if not os.path.isfile(os.getcwd()+os.sep+weights_path1):
            if not os.path.isfile(weights_path1):
               print("Weights path 1 does not exist ... exiting")
               sys.exit()
            else:
               weights_path.append(weights_path1)
        else:
           weights_path.append(weights_path1)

    if 'weights_path2' in locals():
        if not os.path.isfile(os.getcwd()+os.sep+weights_path2):
            if not os.path.isfile(weights_path2):
               print("Weights path 2 does not exist ... exiting")
               sys.exit()
            else:
               weights_path.append(weights_path2)
        else:
           weights_path.append(weights_path2)

    if 'weights_path3' in locals():
        if not os.path.isfile(os.getcwd()+os.sep+weights_path3):
            if not os.path.isfile(weights_path3):
               print("Weights path 3 does not exist ... exiting")
               sys.exit()
            else:
               weights_path.append(weights_path3)
        else:
           weights_path.append(weights_path3)

    if 'weights_path4' in locals():
        if not os.path.isfile(os.getcwd()+os.sep+weights_path4):
            if not os.path.isfile(weights_path4):
               print("Weights path 4 does not exist ... exiting")
               sys.exit()
            else:
               weights_path.append(weights_path4)
        else:
           weights_path.append(weights_path4)

    try:
       # load the user configs
       with open(os.getcwd()+os.sep+configfile) as f:
          config = json.load(f)
    except:
       # load the user configs
       with open(configfile) as f:
          config = json.load(f)

    ###===================================================
    #csvfile containing image names and class values
    csvfile = config["csvfile"]
    #csvfile containing image names and class values
    res_folder = config["res_folder"]
    #folder containing csv file and that will contain model outputs
    name = config["name"]
    #name prefix for output files
    #convert imagery to greyscale or not
    dropout = config["dropout"]
    #dropout factor
    scale = config["scale"] #do scaling on variable
    greyscale = config['greyscale']

    try:
       numclass = config['numclass']
    except:
       numclass = 0

    try:
       greyscale = config['greyscale']
    except:
       greyscale = 'true'

    #output variables
    vars = [k for k in config.keys() if not npany([k.startswith('base'), k.startswith('MAX_LR'),
            k.startswith('MIN_LR'), k.startswith('DO_AUG'), k.startswith('SHALLOW'),
            k.startswith('res_folder'), k.startswith('train_csvfile'), k.startswith('csvfile'),
            k.startswith('test_csvfile'), k.startswith('name'),
            k.startswith('greyscale'), k.startswith('aux_in'),
            k.startswith('dropout'), k.startswith('N'),
            k.startswith('scale'), k.startswith('numclass')])]
    vars = sorted(vars)

    #this relates to 'mimo' and 'miso' modes that are planned for the future but not currently implemented
    auxin = [k for k in config.keys() if k.startswith('aux_in')]

    if len(auxin) > 0:
       auxin = config[auxin[0]]   ##at least for now, just one 'auxilliary' (numerical/categorical) input in addition to imagery
       if len(vars) ==1:
          mode = 'miso'
       elif len(vars) >1:
          mode = 'mimo'
    else:
       if len(vars) ==1:
          mode = 'siso'
       elif len(vars) >1:
          mode = 'simo'

    print("Mode: %s" % (mode))
    ###==================================================

    if type(BATCH_SIZE) is list and type(weights_path) is not list:
       print("Please specify one weights file per batch size in the list ... exiting")
       sys.exit()

    if (mode=='siso' or mode=='simo'):
       if numclass>0:
          result = estimate_categorical_1image(vars, image, res_folder, dropout,
                                   numclass, greyscale, name, mode, weights_path)
       else:
          result = estimate_siso_simo_1image(vars, image, greyscale,
                             dropout, numclass, scale, name, mode,
                             res_folder, BATCH_SIZE, weights_path)


    print("=====================================")
    print("==============RESULTS ==============")
    if len(vars)>1:
        counter = 0
        for var in vars:
            print(var+": %f" % (result[counter]))
            counter +=1
    else:
        print(vars[0]+": %f" % (result))