def loadAreaSourceFromSHP(filename_in,
mmax_data,
background_layer,
filename_out=None,
layer2file=False,
ui_mode=True):
"""Load area source layer from Shapefile, independent of QGis UI."""
temp_area_source_layer = QgsVectorLayer(filename_in, "Area Sources", "ogr")
# PostGIS SRID 4326 is allocated for WGS84
crs = QgsCoordinateReferenceSystem(
4326, QgsCoordinateReferenceSystem.PostgisCrsId)
layer = utils.shp2memory(temp_area_source_layer, "Area Sources")
layer.setCrs(crs)
# check if all features are okay
_checkAreaSourceLayer(layer, ui_mode=ui_mode)
# assign Mmax from Meletti data set
assignMmaxfromMelettiDataset(layer, mmax_data)
# assign attributes from background zones
assignAttributesFromBackgroundZones(layer,
background_layer,
MCDIST_ATTRIBUTES,
ui_mode=ui_mode)
# write memory layer to disk (as a Shapefile)
if layer2file is True:
utils.writeLayerToShapefile(layer, filename_out, crs, ui_mode=ui_mode)
return layer
def loadAreaSourceFromSHP(filename_in, mmax_data, background_layer,
filename_out=None, layer2file=False, ui_mode=True):
"""Load area source layer from Shapefile, independent of QGis UI."""
temp_area_source_layer = QgsVectorLayer(filename_in, "Area Sources", "ogr")
# PostGIS SRID 4326 is allocated for WGS84
crs = QgsCoordinateReferenceSystem(4326,
QgsCoordinateReferenceSystem.PostgisCrsId)
layer = utils.shp2memory(temp_area_source_layer, "Area Sources")
layer.setCrs(crs)
# check if all features are okay
_checkAreaSourceLayer(layer, ui_mode=ui_mode)
# assign Mmax from Meletti data set
assignMmaxfromMelettiDataset(layer, mmax_data)
# assign attributes from background zones
assignAttributesFromBackgroundZones(layer, background_layer,
MCDIST_ATTRIBUTES, ui_mode=ui_mode)
# write memory layer to disk (as a Shapefile)
if layer2file is True:
utils.writeLayerToShapefile(layer, filename_out, crs, ui_mode=ui_mode)
return layer
def loadFaultSourceFromSHP(filename_in, filename_out=None, layer2file=False):
"""Load fault source layer from Shapefile, independent of QGis UI."""
temp_fault_source_layer = QgsVectorLayer(filename_in, "Fault Sources",
"ogr")
# PostGIS SRID 4326 is allocated for WGS84
crs = QgsCoordinateReferenceSystem(
4326, QgsCoordinateReferenceSystem.PostgisCrsId)
layer = utils.shp2memory(temp_fault_source_layer, "Fault Sources")
layer.setCrs(crs)
# write memory layer to disk (as a Shapefile)
if layer2file is True:
utils.writeLayerToShapefile(layer, filename_out, crs)
return layer
def loadFaultSourceFromSHP(filename_in, filename_out=None, layer2file=False):
"""Load fault source layer from Shapefile, independent of QGis UI."""
temp_fault_source_layer = QgsVectorLayer(filename_in, "Fault Sources",
"ogr")
# PostGIS SRID 4326 is allocated for WGS84
crs = QgsCoordinateReferenceSystem(4326,
QgsCoordinateReferenceSystem.PostgisCrsId)
layer = utils.shp2memory(temp_fault_source_layer, "Fault Sources")
layer.setCrs(crs)
# write memory layer to disk (as a Shapefile)
if layer2file is True:
utils.writeLayerToShapefile(layer, filename_out, crs)
return layer
