<jupyter_start><jupyter_text>Test for MohammedThis container was started with**sudo docker run -d -p 433:8888 --name=sar -v /home/mort/imagery/mohammed/Data:/home/imagery mort/sardocker**<jupyter_code>%matplotlib inline<jupyter_output><empty_output><jupyter_text>Here are the RadarSat-2 quadpol coherency matrix image directories as created from the Sentinel-1 Toolbox:<jupyter_code>ls /home/imagery<jupyter_output><empty_output><jupyter_text>To combine the matrix bands into a single GeoTiff image, we run the python script ingestrs2quad.py:<jupyter_code>run /home/ingestrs2quad /home/imagery/RS2_OK82571_PK721079_DK650144_FQ17W_20160403_230258_HH_VV_HV_VH_SLC/
run /home/ingestrs2quad /home/imagery/RS2_OK82571_PK721080_DK650145_FQ17W_20160427_230257_HH_VV_HV_VH_SLC/
run /home/ingestrs2quad /home/imagery/RS2_OK82571_PK721081_DK650146_FQ17W_20160614_230256_HH_VV_HV_VH_SLC/<jupyter_output><empty_output><jupyter_text>Here is an RGB display of the three diagonal matrix elements of the above image (bands 1,6 and 9):<jupyter_code>run /home/dispms -f /home/imagery/RS2_OK82571_PK721081_DK650146_FQ17W_20160614_230256_HH_VV_HV_VH_SLC/polSAR.tif \
-p [1,6,9]<jupyter_output><empty_output><jupyter_text>To estimate the equivalent number of looks, run the python script enlml.py:<jupyter_code>run /home/enlml /home/imagery/RS2_OK82571_PK721081_DK650146_FQ17W_20160614_230256_HH_VV_HV_VH_SLC/polSAR.tif<jupyter_output><empty_output><jupyter_text>So the ENL would appear to be about 5.To run the change sequential change detection on the three images, run the bash script sar_seq_rs2quad.sh. It gathers the three images together and calls the python script sar_seq.py which does the change detection. By choosing a spatial subset (in this case 400x400), the images are clipped and co-registered to the first image. This might be unnecessary if the images are well registered anyway.If you have a multicore processor you can eneable parallel computation by openeing a terminal window in the container (new terminal) and runningipcluster start -n 4<jupyter_code>!/home/sar_seq_rs2quad.sh 20160403 20160427 20160614 [50,50,400,400] 5 0.01<jupyter_output><empty_output><jupyter_text>Here is the change map for the most recent changes:<jupyter_code>run /home/dispms \
-f /home/imagery/RS2_OK82571_PK721079_DK650144_FQ17W_20160403_230258_HH_VV_HV_VH_SLC/sarseq(20160403-1-20160614)_cmap.tif -c<jupyter_output><empty_output>