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