Difference between revisions of "Basin Query Notes"

From SCECpedia
Jump to navigationJump to search
(Created page with "I have created all of the plots discussed (and a few more). I did this by creating a new OpenSHA class, ZFilePlot. That class is posted here: https://github.com/opensha/opens...")
 
 
(4 intermediate revisions by the same user not shown)
Line 1: Line 1:
−
I have created all of the plots discussed (and a few more). I did this by creating a new OpenSHA class, ZFilePlot. That class is posted here:
+
== Example Comparison Plots OpenSHA versus UCVM basin_query: ==
 +
*[http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare/plots/BQ_CCA_Z2.5_Second.png basin_query Z2.5 second crossing]
 +
*[http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare/plots/OpenSHA_Z2.5_CCA_Second.png OpenSHA Z2.5 second crossing]
 +
 
 +
== Log-ratio of basin_query to OpenSHA: ==
 +
*[http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare/plots/BQ_CCA_Z2.5_Second_vs_OpenSHA_Z2.5_CCA_Second_map.png BQ compared to OpenSHA]
 +
 
 +
== OpenSHA-produced Z1.0 and Z2.5 Plots are posted here ==
 +
*[http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare/plots/ OpenSHA Z Plots]
 +
 
 +
== Scatter plot of basin_query vs OpenSHA: ==
 +
*[http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare/plots/BQ_CCA_Z2.5_Second_vs_OpenSHA_Z2.5_CCA_Second_scatter.png Scatter Plot]
 +
 
 +
Those plots do show a difference in the registration of the model in Mei's files. This is probably related to the issue Christine described.
 +
 
 +
A zip file with the code, all data files, and all plots is posted here as that may be more convenient than manually downloading everything:
 +
 
 +
http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare.zip
 +
 
 +
== How Plots were Produced ==
 +
Kevin created all of the plots discussed (and a few more). I did this by creating a new OpenSHA class, ZFilePlot. That class is posted here:
  
 
https://github.com/opensha/opensha-cybershake/blob/master/src/scratch/kevin/cybershake/ZFilePlot.java
 
https://github.com/opensha/opensha-cybershake/blob/master/src/scratch/kevin/cybershake/ZFilePlot.java
Line 17: Line 37:
 
The forward slashes in the label arguments are for spaces within a single argument. While I hope this tool will be useful to Mei in the short-term, I don't believe that it should be bundled with UCVM (it relies on an opensha web service that may change/move at any time), and I still believe that a python/matplotlib plotting script should be created within UCVM.
 
The forward slashes in the label arguments are for spaces within a single argument. While I hope this tool will be useful to Mei in the short-term, I don't believe that it should be bundled with UCVM (it relies on an opensha web service that may change/move at any time), and I still believe that a python/matplotlib plotting script should be created within UCVM.
  
−
Plots are posted here (I'll let Phil and Mei decide which plots are desired for the wiki, as uploading all of them manually is rather onerous):
 
−
 
−
http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare/plots/
 
−
 
−
For example:
 
−
 
−
basin_query Z2.5 second crossing: http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare/plots/BQ_CCA_Z2.5_Second.png
 
−
OpenSHA Z2.5 second crossing: http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare/plots/OpenSHA_Z2.5_CCA_Second.png
 
−
Log-ratio of basin_query to OpenSHA: http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare/plots/BQ_CCA_Z2.5_Second_vs_OpenSHA_Z2.5_CCA_Second_map.png
 
−
Scatter plot of basin_query vs OpenSHA: http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare/plots/BQ_CCA_Z2.5_Second_vs_OpenSHA_Z2.5_CCA_Second_scatter.png
 
−
 
−
Those plots do show a difference in the registration of the model in Mei's files. This is probably related to the issue Christine described.
 
−
 
−
A zip file with the code, all data files, and all plots is posted here as that may be more convenient than manually downloading everything:
 
−
 
−
http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare.zip
 
 
== Related Entries ==
 
== Related Entries ==
 
*[[UCVM]]
 
*[[UCVM]]
 +
*[[UCVM_Review]]
 
*[[USR]]
 
*[[USR]]

Latest revision as of 00:20, 7 December 2017

Example Comparison Plots OpenSHA versus UCVM basin_query:

Log-ratio of basin_query to OpenSHA:

OpenSHA-produced Z1.0 and Z2.5 Plots are posted here

Scatter plot of basin_query vs OpenSHA:

Those plots do show a difference in the registration of the model in Mei's files. This is probably related to the issue Christine described.

A zip file with the code, all data files, and all plots is posted here as that may be more convenient than manually downloading everything:

http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare.zip

How Plots were Produced

Kevin created all of the plots discussed (and a few more). I did this by creating a new OpenSHA class, ZFilePlot. That class is posted here:

https://github.com/opensha/opensha-cybershake/blob/master/src/scratch/kevin/cybershake/ZFilePlot.java

A build of this tool, along with the bash script and all data I used to run it is posted here:

http://opensha.usc.edu/ftp/kmilner/cvm_data/2017_z_file_compare/

The usage of the tool (requiring java 8 and after downloading the jar file) is:

java -cp opensha-cybershake-all.jar scratch.kevin.cybershake.ZFilePlot <min-lat> <min-lon> <nx> <ny> <grid-spacing> <max-z> <data-file> <label> [<compare-file> <compare-label>]

The "max-z" argument scales the color pallet to this value (I use 10 for Z2.5 and 2 for Z1.0). All plots are in km (the files themselves are in meters, but I scale them before plotting). The last 2 arguments are optional, for creating ratio/scatter plots. For example, to compare basin_query Z2.5 to OpenSHA Z2.5:

java -cp opensha-cybershake-all.jar scratch.kevin.cybershake.ZFilePlot 33.35 -123 1551 1201 0.005 10 basin_query_mpi_cca_z2.5_n20_second.result BQ\ CCA\ Z2.5\ Second opensha_cca_second_z2.5.bin OpenSHA\ Z2.5\ CCA\ Second

The forward slashes in the label arguments are for spaces within a single argument. While I hope this tool will be useful to Mei in the short-term, I don't believe that it should be bundled with UCVM (it relies on an opensha web service that may change/move at any time), and I still believe that a python/matplotlib plotting script should be created within UCVM.

Related Entries