Add plotting script for correction files

[adeane-ga]

This script is designed to plot the correction files (NumPy), against the original interferogram and the corrected interferogram (GeoTIFFs) to aid user's interpretation of corrections and help debugging with PyRate.

Here is an example where from left to right is original IFG (interferogram_dir/), correction applied (orbit_error_dir/), and resulting corrected IFG (temp_mlooked_dir/):

Help Output

Example Usage:

python3 ~/PyRate/utils/plot_correction_files.py ./interferogram_dir ./orbit_error_dir ./temp_mlooked_dir ./plot_save_dir 1 10

the last two arguments will indicate that it will plot from the 1st to the 10th IFG (10 separate PNGs), this can be changed to any range within the maximum number of IFGs in user PyRate directory. Using 1 1 will plot the 1st only

Notes

• This script is designed to be generic and as simple as possible, so it is easy for any given user to change it on the fly in there own local repository (rather than a script that deals with every possible combination of processing).
• The script is also naïve, in that it relies on the user to decide how valid the comparison is, depending on user's own processing (which corrections were applied etc.).
• I have only tested this script on the Orbit and APS corrections, it theoretically should work on the DEM error correction too but I have not tested it.
• Currently, the files we save out as APS "correction" files, don't seem to contain the full correction, work will be done to investigate this.
• Users can plot a single IFG or the full time series based on the last two arguments (but users need to have awareness of how many IFGs are in their directory).
• If plotting full time series, consider running as a QSUB job (if on NCI), otherwise plotting subsets will not take up too much memory.