Difference between revisions of "Final CCA06 Crossings Result"

From SCECpedia
Jump to navigationJump to search
 
(17 intermediate revisions by the same user not shown)
Line 1: Line 1:
 +
As the result of this crossing review study, it is apparent that it is not possible to
 +
have a single crossing algorithm that is applicable to all the velocity models and
 +
a reasonable crossing algorithm requires domain scientists' post processing evaluation.
 +
Therefore, we are retaining the 5 algorithms that were produced from the [[Basin Query Result]] and
 +
also create a new output file that contains raw crossing information that can be used
 +
in later time.
 +
 +
== Raw crossing result file ==
  
 
An option was added in basin_query_mpi(-o) to create an ascii file of top 3 crossings plus the last crossing per grid point within ROI region. This result file can be post processed to determine the best
 
An option was added in basin_query_mpi(-o) to create an ascii file of top 3 crossings plus the last crossing per grid point within ROI region. This result file can be post processed to determine the best
Line 24: Line 32:
  
  
 +
Following is the final set of CCA06 crossing result with updated projection configuration from [[CVM_Projection_Issue_Discussion]].
  
The final set of CCA06 crossing result with updated projection configuration from *[[CVM_Projection_Issue_Discussion]].
+
The target crossing to be used for CCA06 is the 2nd row, the firstOrSecond set.
 
 
The target crossing to be used for CCA06 is the 2nd row, the firstOrSecond set.  
 
 
 
  
 +
== Basin Query Algorithm Result for CCA 06 without ELY GTL ==
 +
{| class="wikitable sortable" border="2"
 +
|-
 +
! scope="col" style="width: 30%;"| Algorithm
 +
! scope="col" style="width: 35%;"| Threshold = 1000m/s
 +
! scope="col" style="width: 35%;"| Threshold = 2500m/s
 +
|-
 +
! scope="row" | First Crossing           
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_nGTL.first data]
 +
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_nGTL_first_map.png png]
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_nGTL.first  data]
 +
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_nGTL_first_map.png png]
 +
|-
 +
! scope="row" | First Or Second Crossing           
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_nGTL.firstOrSecond data]
 +
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_nGTL_firstOrSecond_map.png png]
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_nGTL.firstOrSecond  data]
 +
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_nGTL_firstOrSecond_map.png png]
 +
|-
 +
! scope="row" | Last Crossing           
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_nGTL.last data]
 +
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_nGTL_last_map.png png]
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_nGTL.last  data]
 +
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_nGTL_last_map.png png]|
 +
|-
 +
! scope="row" | Second Only  Crossing       
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_nGTL.secondOnly data]
 +
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_nGTL_secondOnly_map.png png]
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_nGTL.secondOnly  data]
 +
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_nGTL_secondOnly_map.png png]
 +
|-
 +
! scope="row" | Last but
 +
not First or Second Crossing         
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_nGTL.threeLast data]
 +
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_nGTL_threeLast_map.png png]
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_nGTL.threeLast  data]
 +
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_nGTL_threeLast_map.png png]
 +
|-
 +
! scope="row" | Raw crossing result     
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_nGTL.result data]
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_nGTL.result  data]
 +
|-
 +
! scope="row" | Metadata result     
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_nGTL.meta.json data]
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_nGTL.meta.json  data]
 +
|}
  
== Basin Query Algorithm Result for CCA 06 ==
+
== Basin Query Algorithm Result for CCA 06 with ELY GTL ==
 
{| class="wikitable sortable" border="2"
 
{| class="wikitable sortable" border="2"
 
|-
 
|-
Line 39: Line 91:
 
|-
 
|-
 
! scope="row" | First Crossing             
 
! scope="row" | First Crossing             
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca_z1_final.first data]
+
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0.first data]
[http://hypocenter.usc.edu/research/ucvmc_result/cca_z1_final_first_map.png png]
+
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_first_map.png png]
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_final.first  data]
+
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5.first  data]
[http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_final_first_map.png png]
+
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_first_map.png png]
 
|-
 
|-
 
! scope="row" | First Or Second Crossing             
 
! scope="row" | First Or Second Crossing             
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca_z1_final.firstOrSecond data]
+
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0.firstOrSecond data]
[http://hypocenter.usc.edu/research/ucvmc_result/cca_z1_final_firstOrSecond_map.png png]
+
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_firstOrSecond_map.png png]
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_final.firstOrSecond  data]
+
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5.firstOrSecond  data]
[http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_final_firstOrSecond_map.png png]
+
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_firstOrSecond_map.png png]
 
|-
 
|-
 
! scope="row" | Last Crossing             
 
! scope="row" | Last Crossing             
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca_z1_final.last data]
+
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0.last data]
[http://hypocenter.usc.edu/research/ucvmc_result/cca_z1_final_last_map.png png]
+
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_last_map.png png]
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_final.last  data]
+
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5.last  data]
[http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_final_last_map.png png]|
+
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_last_map.png png]|
 
|-
 
|-
 
! scope="row" | Second Only  Crossing         
 
! scope="row" | Second Only  Crossing         
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca_z1_final.secondOnly data]
+
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0.secondOnly data]
[http://hypocenter.usc.edu/research/ucvmc_result/cca_z1_final_secondOnly_map.png png]
+
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_secondOnly_map.png png]
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_final.secondOnly  data]
+
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5.secondOnly  data]
[http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_final_secondOnly_map.png png]
+
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_secondOnly_map.png png]
 
|-
 
|-
 
! scope="row" | Last but
 
! scope="row" | Last but
 
not First or Second Crossing           
 
not First or Second Crossing           
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca_z1_final.threeLast data]
+
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0.threeLast data]
[http://hypocenter.usc.edu/research/ucvmc_result/cca_z1_final_threeLast_map.png png]
+
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0_threeLast_map.png png]
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_final.threeLast  data]
+
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5.threeLast  data]
[http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_final_threeLast_map.png png]
+
[http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5_threeLast_map.png png]
 +
|-
 +
! scope="row" | Raw crossing result     
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0.result data]
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5.result  data]
 
|-
 
|-
 +
! scope="row" | Metadata result     
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z1.0.meta.json data]
 +
|| [http://hypocenter.usc.edu/research/ucvmc_result/cca-final/cca_z2.5.meta.json  data]
 
|}
 
|}
 +
  
 
The basin query command that was used to generate the result:
 
The basin query command that was used to generate the result:
  
mpirun -np 20 ./basin_query_mpi  
+
  mpirun -np 20 ./basin_query_mpi  
 
   -b ${TEST_TOP_PATH}/RESULT/cca_z2.5.first,
 
   -b ${TEST_TOP_PATH}/RESULT/cca_z2.5.first,
 
   ${TEST_TOP_PATH}/RESULT/cca_z2.5.firstOrSecond,
 
   ${TEST_TOP_PATH}/RESULT/cca_z2.5.firstOrSecond,
Line 85: Line 145:
 
and,
 
and,
  
mpirun -np 20 ./basin_query_mpi  
+
  mpirun -np 20 ./basin_query_mpi  
 
   -b ${TEST_TOP_PATH}/RESULT/cca_z1.0.first,
 
   -b ${TEST_TOP_PATH}/RESULT/cca_z1.0.first,
 
   ${TEST_TOP_PATH}/RESULT/cca_z1.0.firstOrSecond,
 
   ${TEST_TOP_PATH}/RESULT/cca_z1.0.firstOrSecond,
Line 97: Line 157:
 
and the plotting commands that were used are:
 
and the plotting commands that were used are:
  
./plot_z25_map.py -f cca_z2.5.firstOrSecond  
+
  ./plot_z25_map.py -f cca_z2.5.firstOrSecond  
 
   -o cca_z2.5_firstOrSecond_map.png  
 
   -o cca_z2.5_firstOrSecond_map.png  
 
   -e "(first or second crossing)" -b 33.349998,-123.000000  
 
   -e "(first or second crossing)" -b 33.349998,-123.000000  
Line 104: Line 164:
 
and,
 
and,
  
./plot_z10_map.py -f cca_z1.firstOrSecond
+
  ./plot_z10_map.py -f cca_z1.firstOrSecond
 
   -o cca_z1_firstOrSecond_map.png -e "(first or second crossing)"
 
   -o cca_z1_firstOrSecond_map.png -e "(first or second crossing)"
 
   -b 33.349998,-123.000000 -u 39.355000,-115.245003 -c cca  
 
   -b 33.349998,-123.000000 -u 39.355000,-115.245003 -c cca  
 
   -s 0.005000 -x 1551 -y 1201
 
   -s 0.005000 -x 1551 -y 1201
  
 +
 +
 +
CCA06 no GTL, section of comparable area to Kevin's zoomed in area plot :
 +
 +
{|
 +
| [[File:zoom_first.png|thumb|400px|first crossing]]
 +
| [[File:zoom.firstOrSecond.png|thumb|400px|first or second crossing]]
 +
|}
 +
 +
Scatter plot :
 +
 +
{|
 +
| [[File:compare.png|thumb|400px|first crossing]]
 +
|}
  
 
== Related Entries ==
 
== Related Entries ==

Latest revision as of 06:28, 12 June 2018

As the result of this crossing review study, it is apparent that it is not possible to have a single crossing algorithm that is applicable to all the velocity models and a reasonable crossing algorithm requires domain scientists' post processing evaluation. Therefore, we are retaining the 5 algorithms that were produced from the Basin Query Result and also create a new output file that contains raw crossing information that can be used in later time.

Raw crossing result file

An option was added in basin_query_mpi(-o) to create an ascii file of top 3 crossings plus the last crossing per grid point within ROI region. This result file can be post processed to determine the best strategy to use in selecting the suitable crossing to use for the model.

The format of the content:

lon lat count crossing1 crossing2 crossing3 crossing4

examples,

 *no crossing
     lon  lat  0  -1  -1  -1  -1
 *1 crossing (c1 = first crossing = last crossing)
     lon  lat  1  c1  -1  -1  -1
 *2 crossings (c1 = first, c2 = 2nd = last)
     lon  lat  2  c1  c2  -1  -1
 *3 crossings (c1 = first, c2 = 2nd, c3 = third = last)
     lon  lat  3  c1  c2  c3  -1
 *4 crossings (c1 = first, c2 = 2nd, c3 = third, c4 = fourth = last)
     lon  lat  4  c1  c2  c3  c4
 *5 crossings (c1 = first, c2 = 2nd, c3 = third, c5 = fifth = last)
     lon  lat  5  c1  c2  c3  c5
 *6 crossings (c1 = first, c2 = 2nd, c3 = third, c6 = sixth = last)
     lon  lat  6  c1  c2  c3  c6


Following is the final set of CCA06 crossing result with updated projection configuration from CVM_Projection_Issue_Discussion.

The target crossing to be used for CCA06 is the 2nd row, the firstOrSecond set.

Basin Query Algorithm Result for CCA 06 without ELY GTL

Algorithm Threshold = 1000m/s Threshold = 2500m/s
First Crossing data

png

data

png

First Or Second Crossing data

png

data

png

Last Crossing data

png

data

png|

Second Only Crossing data

png

data

png

Last but

not First or Second Crossing

data

png

data

png

Raw crossing result data data
Metadata result data data

Basin Query Algorithm Result for CCA 06 with ELY GTL

Algorithm Threshold = 1000m/s Threshold = 2500m/s
First Crossing data

png

data

png

First Or Second Crossing data

png

data

png

Last Crossing data

png

data

png|

Second Only Crossing data

png

data

png

Last but

not First or Second Crossing

data

png

data

png

Raw crossing result data data
Metadata result data data


The basin query command that was used to generate the result:

 mpirun -np 20 ./basin_query_mpi 
 -b ${TEST_TOP_PATH}/RESULT/cca_z2.5.first,
 ${TEST_TOP_PATH}/RESULT/cca_z2.5.firstOrSecond,
 ${TEST_TOP_PATH}/RESULT/cca_z2.5.last,
 ${TEST_TOP_PATH}/RESULT/cca_z2.5.secondOnly,
 ${TEST_TOP_PATH}/RESULT/cca_z2.5.threeLast 
 -o ${TEST_TOP_PATH}/RESULT/cca_z2.5.result,
 ${TEST_TOP_PATH}/RESULT/cca_z2.5.json 
 -f ucvm.conf -m cca -i 20 -v 2500 -l 33.35,-123.0 -s 0.005 -x 1551 -y 1201

and,

 mpirun -np 20 ./basin_query_mpi 
 -b ${TEST_TOP_PATH}/RESULT/cca_z1.0.first,
 ${TEST_TOP_PATH}/RESULT/cca_z1.0.firstOrSecond,
 ${TEST_TOP_PATH}/RESULT/cca_z1.0.last,
 ${TEST_TOP_PATH}/RESULT/cca_z1.0.secondOnly,
 ${TEST_TOP_PATH}/RESULT/cca_z1.0.threeLast 
 -o ${TEST_TOP_PATH}/RESULT/cca_z1.0.result 
 -f ucvm.conf -m cca -i 20 -v 1000 -l 33.35,-123.0 -s 0.005 -x 1551 -y 1201 


and the plotting commands that were used are:

 ./plot_z25_map.py -f cca_z2.5.firstOrSecond 
 -o cca_z2.5_firstOrSecond_map.png 
 -e "(first or second crossing)" -b 33.349998,-123.000000 
 -u 39.355000,-115.245003 -c cca -s 0.005000 -x 1551 -y 1201

and,

 ./plot_z10_map.py -f cca_z1.firstOrSecond
 -o cca_z1_firstOrSecond_map.png -e "(first or second crossing)"
 -b 33.349998,-123.000000 -u 39.355000,-115.245003 -c cca 
 -s 0.005000 -x 1551 -y 1201


CCA06 no GTL, section of comparable area to Kevin's zoomed in area plot :

first crossing
first or second crossing

Scatter plot :

first crossing

Related Entries