def to_openquake_site(self, missing_vs30=None):
"""
Returns the site as an instance of the :class:
openquake.hazardlib.site.Site
"""
if self.vs30:
vs30 = self.vs30
vs30_measured = self.vs30_measured
else:
vs30 = missing_vs30
vs30_measured = False
if self.z1pt0:
z1pt0 = self.z1pt0
else:
z1pt0 = vs30_to_z1pt0_as08(vs30)
if self.z2pt5:
z2pt5 = self.z2pt5
else:
z2pt5 = z1pt0_to_z2pt5(z1pt0)
location = Point(self.longitude,
self.latitude,
-self.altitude / 1000.) # Elevation from m to km
return Site(location,
vs30,
vs30_measured,
z1pt0,
z2pt5,
backarc=self.backarc,
id=self.id)
def to_openquake_site(self, missing_vs30=None):
"""
Returns the site as an instance of the :class:
openquake.hazardlib.site.Site
"""
if self.vs30:
vs30 = self.vs30
vs30_measured = self.vs30_measured
else:
vs30 = missing_vs30
vs30_measured = False
if self.z1pt0:
z1pt0 = self.z1pt0
else:
z1pt0 = vs30_to_z1pt0_as08(vs30)
if self.z2pt5:
z2pt5 = self.z2pt5
else:
z2pt5 = z1pt0_to_z2pt5(z1pt0)
location = Point(self.longitude,
self.latitude,
-self.altitude / 1000.) # Elevation from m to km
oq_site = Site(location,
vs30,
z1pt0,
z2pt5,
vs30measured=vs30_measured,
backarc=self.backarc)
setattr(oq_site, "id", self.id)
return oq_site
def _parse_site_data(self, metadata):
"""
Returns the site data as an instance of the :class:
smtk.sm_database.RecordSite
"""
site = RecordSite(metadata["Station Sequence Number"],
metadata["Station ID No."],
metadata["Station ID No."],
get_float(metadata["Station Longitude"]),
get_float(metadata["Station Latitude"]),
0.0, # Elevation data not given
get_float(metadata["Preferred Vs30 (m/s)"]),
network_code=metadata["Owner"])
site.nehrp = metadata["Preferred NEHRP Based on Vs30"]
site.vs30_measured_type = metadata["Measured/Inferred Class"]
if site.vs30_measured_type in ["0", "5"]:
site.vs30_measured = True
else:
site.vs30_measured = False
site.vs30_uncertainty = get_float(
metadata["Sigma of Vs30 (in natural log Units)"])
site.z1pt0 = get_float(metadata["Z1 (m)"])
site.z1pt5 = get_float(metadata["Z1.5 (m)"])
site.z2pt5 = get_float(metadata["Z2.5 (m)"])
# Implement default values for z1pt0 and z2pt5
if site.z1pt0 is None:
site.z1pt0 = rcfg.vs30_to_z1pt0_as08(site.vs30)
if site.z2pt5 is None:
site.z2pt5 = rcfg.z1pt0_to_z2pt5(site.z1pt0)
site.arc_location = metadata["Forearc/Backarc for subduction events"]
site.instrument_type = metadata["Type of Recording"]
return site
def to_openquake_site(self, missing_vs30=None):
"""
Returns the site as an instance of the :class:
openquake.hazardlib.site.Site
"""
if self.vs30:
vs30 = self.vs30
vs30_measured = self.vs30_measured
else:
vs30 = missing_vs30
vs30_measured = False
if self.z1pt0:
z1pt0 = self.z1pt0
else:
z1pt0 = vs30_to_z1pt0_as08(vs30)
if self.z2pt5:
z2pt5 = self.z2pt5
else:
z2pt5 = z1pt0_to_z2pt5(z1pt0)
return Site(Point(self.longitude, self.latitude),
vs30,
vs30_measured,
z1pt0,
z2pt5,
self.id)
def _parse_site_data(self, metadata):
"""
Returns the site data as an instance of the :class:
smtk.sm_database.RecordSite
"""
site = RecordSite(
self._get_site_id(metadata["Station ID"]),
metadata["Station ID"],
metadata["Station Code"],
get_float(metadata["Station Longitude (deg positive E)"]),
get_float(metadata["Station Latitude (deg positive N)"]),
get_float(metadata["Station Elevation (m)"]),
get_float(metadata["Preferred Vs30 (m/s)"]),
site_class=metadata[
'Site Class (Hard Rock; Rock; Stiff Soil; Soft Soil)'],
network_code=None, # not provided
)
# network_code=metadata["Owner"])
site.nehrp = metadata["Preferred NEHRP Based on Vs30"]
site.vs30_measured_type = metadata["Measured(1)/Inferred(2) Class"]
if metadata["Measured(1)/Inferred(2) Class"] == "1":
site.vs30_measured = True
else:
site.vs30_measured = False
site.vs30_uncertainty = get_float(
metadata["Sigma of Vs30 (in natural log Units)"])
site.z1pt0 = get_float(metadata["Z1 (m)"])
site.z1pt5 = None
# site.z1pt5 = get_float(metadata["Z1.5 (m)"])
site.z2pt5 = get_float(metadata["Z2.5 (m)"])
# Implement default values for z1pt0 and z2pt5
if site.z1pt0 is None:
site.z1pt0 = rcfg.vs30_to_z1pt0_as08(site.vs30)
if site.z2pt5 is None:
site.z2pt5 = rcfg.z1pt0_to_z2pt5(site.z1pt0)
else:
# Need to convert z2pt5 from m to km
site.z2pt5 = site.z2pt5 / 1000.0
if "Backarc" in metadata["Forearc/Backarc for subduction events"]:
site.backarc = True
site.instrument_type = metadata["Digital (D)/Analog (A) Recording"]
return site
def _parse_site_data(self, metadata):
"""
Returns the site data as an instance of the :class:
smtk.sm_database.RecordSite
"""
site = RecordSite(
self._get_site_id(metadata["Station ID"]),
metadata["Station ID"],
metadata["Station Code"],
get_float(metadata["Station Longitude (deg positive E)"]),
get_float(metadata["Station Latitude (deg positive N)"]),
get_float(metadata["Station Elevation (m)"]),
get_float(metadata["Preferred Vs30 (m/s)"]),
site_class=metadata['Site Class (Hard Rock; Rock; Stiff Soil; Soft Soil)'],
network_code=None, # not provided
)
# network_code=metadata["Owner"])
site.nehrp = metadata["Preferred NEHRP Based on Vs30"]
site.vs30_measured_type = metadata["Measured(1)/Inferred(2) Class"]
if metadata["Measured(1)/Inferred(2) Class"] == "1":
site.vs30_measured = True
else:
site.vs30_measured = False
site.vs30_uncertainty = get_float(
metadata["Sigma of Vs30 (in natural log Units)"])
site.z1pt0 = get_float(metadata["Z1 (m)"])
site.z1pt5 = None
# site.z1pt5 = get_float(metadata["Z1.5 (m)"])
site.z2pt5 = get_float(metadata["Z2.5 (m)"])
# Implement default values for z1pt0 and z2pt5
if site.z1pt0 is None:
site.z1pt0 = rcfg.vs30_to_z1pt0_as08(site.vs30)
if site.z2pt5 is None:
site.z2pt5 = rcfg.z1pt0_to_z2pt5(site.z1pt0)
else:
# Need to convert z2pt5 from m to km
site.z2pt5 = site.z2pt5/1000.0
if "Backarc" in metadata["Forearc/Backarc for subduction events"]:
site.backarc = True
site.instrument_type = metadata["Digital (D)/Analog (A) Recording"]
return site