Difference between revisions of "UCVMC How to process bin data"
From SCECpedia
Jump to navigationJump to search (→Plot) |
|||
| (6 intermediate revisions by the same user not shown) | |||
| Line 1: | Line 1: | ||
| − | + | Binary data can be generated from UCVMC's plotting tool in python or by | |
| + | one of the UCVMC app that is written to work with the UCVMC api. | ||
| − | == | + | == example plot == |
| + | |||
| + | This is horizontal Vs slice plot from UCVMC's plot_horizontal_slice.py | ||
{| | {| | ||
| Line 9: | Line 12: | ||
| − | == | + | == metadata file == |
| + | |||
| + | Metadata file contains many fields. Some of the fields are specific to the type | ||
| + | of the plot that were made. These are the fields for a horizontal slice plot. | ||
| − | + | lon_list | |
| + | lat_list | ||
| + | lat1 | ||
| + | lon1 | ||
| + | lat2 | ||
| + | lon2 | ||
| + | datafile | ||
| + | num_x | ||
| + | num_y | ||
| + | min | ||
| + | max | ||
| + | color | ||
| + | data_type | ||
| + | outfile | ||
| + | cvm_selected | ||
| + | datapoints | ||
| + | depth | ||
| − | + | (lon1,lat1) is the lower left corner of the plot and (lon2,lat2) is the upper right | |
| + | corner of the plot. lon_list and lat_list are the tick values of the axes. | ||
[http://hypocenter.usc.edu/research/UCVM/cvmh_no_ely_gtl/cvmh_nogtl_vs_0_map_meta.json metadata in json] | [http://hypocenter.usc.edu/research/UCVM/cvmh_no_ely_gtl/cvmh_nogtl_vs_0_map_meta.json metadata in json] | ||
| + | |||
| + | == binary data file == | ||
| + | |||
| + | The data is written out as an array of float32. It can be imported with numpy | ||
| + | in python, | ||
| + | |||
| + | <pre> | ||
| + | import numpy as np | ||
| + | fh = open(rawfile, 'r') | ||
| + | floats = np.fromfile(fh, dtype=np.float32) | ||
| + | </pre> | ||
| + | |||
| + | and fill into the 2D array as, | ||
| + | |||
| + | <pre> | ||
| + | datapoints = np.arange(num_x * num_y,dtype=np.float32).reshape(num_y, num_x) | ||
| + | i=0 | ||
| + | j=0 | ||
| + | for f in floats: | ||
| + | datapoints[i][j] = f | ||
| + | j = j + 1 | ||
| + | if j >= num_x: | ||
| + | j = 0 | ||
| + | i = i + 1 | ||
| + | </pre> | ||
| + | |||
| + | |||
| + | [http://hypocenter.usc.edu/research/UCVM/cvmh_no_ely_gtl/cvmh_nogtl_vs_0_map_data.bin binary data] | ||
| + | |||
== Related Entries == | == Related Entries == | ||
*[[CVM-H 15.1 Maps]] | *[[CVM-H 15.1 Maps]] | ||
*[[UCVM]] | *[[UCVM]] | ||
Latest revision as of 19:21, 5 July 2018
Binary data can be generated from UCVMC's plotting tool in python or by one of the UCVMC app that is written to work with the UCVMC api.
example plot
This is horizontal Vs slice plot from UCVMC's plot_horizontal_slice.py
metadata file
Metadata file contains many fields. Some of the fields are specific to the type of the plot that were made. These are the fields for a horizontal slice plot.
lon_list lat_list lat1 lon1 lat2 lon2 datafile num_x num_y min max color data_type outfile cvm_selected datapoints depth
(lon1,lat1) is the lower left corner of the plot and (lon2,lat2) is the upper right corner of the plot. lon_list and lat_list are the tick values of the axes.
binary data file
The data is written out as an array of float32. It can be imported with numpy in python,
import numpy as np
fh = open(rawfile, 'r')
floats = np.fromfile(fh, dtype=np.float32)
and fill into the 2D array as,
datapoints = np.arange(num_x * num_y,dtype=np.float32).reshape(num_y, num_x)
i=0
j=0
for f in floats:
datapoints[i][j] = f
j = j + 1
if j >= num_x:
j = 0
i = i + 1
