###################################################################
import numpy as np
from scipy.interpolate import griddata
import matplotlib.pyplot as plt
import sys
import os.path
import gamma_reader
import tree
import covaraps
import argparse
import logmgr
import pyproj
#######Logger for varrespy
logger = logmgr.logger('varres')
######For parsing command line options
def parse():
'''Command line parser for varres. Use -h option to get complete list of options.'''
parser = argparse.ArgumentParser(
description=
'Variable resolution resampler for unwrapped interferograms.')
parser.add_argument(
'-i',
action='store',
required=True,
dest="I_name",
help="ROI-PAC geocoded unwrapped file with .rsc in the same directory",
type=str)
import argparse
#import seaborn as sns
from tempfile import mkdtemp
from gridutils import *
from datautils import * # has varres_data class
import logmgr
import pymc3 as pm
from pymc3.step_methods import smc
import pyproj
import pandas as pd
from scipy.interpolate import RegularGridInterpolator
#######Logger for varrespy
logger = logmgr.logger('pymode')
# set up pymc3
logger.info("Using PyMC3 version " + str(pm.__version__))
def parse():
"""
Command line parser for PyMoDE. Use -h option to get complete list of options.
Returns: Tuple (arguments,[list of varres_data objects])
"""
parser = argparse.ArgumentParser(
description=
'Python Modelling of Deformation Events using unwrapped InSAR interferograms'
)
parser.add_argument(