Python SimpleFaultSurface.get_surface_vertexes Examples

Example #1

File: source_model_converter.py Project: gbonorino/nrml_converters

def set_simple_fault_3D_geometry(w, src):
    """
    Set simple fault surface coordinates as shapefile geometry.

    :parameter w:
        Writer
    :parameter src:
        NRML source object
    """
    # Create an oq object for the source
    hc = HC(1., 10., 0.1, 50.)
    converter = sss.NrmlHazardlibConverter(hc)
    src_oq = converter._nrml_source_to_hazardlib(src)
    lon, lat = \
        SimpleFaultSurface.get_surface_vertexes(src_oq.fault_trace,
                                                src_oq.upper_seismogenic_depth,
                                                src_oq.lower_seismogenic_depth,
                                                src_oq.dip)
    # Reorder the vertexes
    print lon, lat
    lons = numpy.concatenate([lon[::2], numpy.flipud(lon[1::2])])
    lats = numpy.concatenate([lat[::2], numpy.flipud(lat[1::2])])
    depths = numpy.concatenate([numpy.ones_like(lon[::2]) *
                                src_oq.upper_seismogenic_depth,
                                numpy.ones_like(lon[::2]) *
                                src_oq.lower_seismogenic_depth])
    # Create the 3D polygon
    w.poly(parts=[[[tlon, tlat, tdep] for tlon, tlat, tdep
                  in zip(list(lons), list(lats), list(depths))]])

Example #2

File: source_model_converter.py Project: preinh/nrml_converters

def set_simple_fault_3D_geometry(w, src):
    """
    Set simple fault surface coordinates as shapefile geometry.

    :parameter w:
        Writer
    :parameter src:
        NRML source object
    """
    # Create an oq object for the source
    converter = sss.NrmlHazardlibConverter(
        investigation_time=50.0, rupture_mesh_spacing=1.0, width_of_mfd_bin=0.1, area_source_discretization=10.0
    )
    src_oq = converter(src)
    lon, lat = SimpleFaultSurface.get_surface_vertexes(
        src_oq.fault_trace, src_oq.upper_seismogenic_depth, src_oq.lower_seismogenic_depth, src_oq.dip
    )
    # Reorder the vertexes
    lons = numpy.concatenate([lon[::2], numpy.flipud(lon[1::2])])
    lats = numpy.concatenate([lat[::2], numpy.flipud(lat[1::2])])
    depths = numpy.concatenate(
        [
            numpy.ones_like(lon[::2]) * src_oq.upper_seismogenic_depth,
            numpy.ones_like(lon[::2]) * src_oq.lower_seismogenic_depth,
        ]
    )
    # Create the 3D polygon
    w.poly(parts=[[[tlon, tlat, tdep] for tlon, tlat, tdep in zip(list(lons), list(lats), list(depths))]])