Recently I wanted to plot the pixel values of two images against each other. I though it would be good to combine extracting the pixel values (using RIOS) and plotting the data (using matplotlib) in a single script.
The script I wrote (two_band_scatter_plot.py) is available to download from the RSGIS Scripts repository
There are three main steps:
- RIOS is used to iterate through the image in small blocks, for each block only pixels that have data in each band are retained.
- The numpy.histogram2d is used to produce a 2D histogram.
- The histogram is plotted using matplotlib.
The 2D histogram was based on code from http://oceanpython.org/2013/02/25/2d-histogram/. As there are a lot of points, displaying as a 2D histogram makes the information easier to see.
The simplest usage of the script is:
python two_band_scatter_plot.py \ --imageA imageA.kea \ --imageB imageB.kea \ --plot outplot.pdf
Which will plot the first band of each image. To change the bands used, labels for the plot and scaling additional flags can be passed in.
python two_band_scatter_plot.py \ --imageA imageA.kea \ --imageB imageB.kea \ --bandA 1 --bandB 1 \ --plot outplot.pdf \ --plotMin 0 --plotMax 0.5 \ --labelA "13157 - 004 (HH)" \ --labelB "13157 - 006 (HH)" \
The script assumes both images are the same projection and resolution, and RIOS takes care of the spatial information. It would be possible to adapt so RIOS resamples the data, see the RIOS Wiki for more details.