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 loadBackgroundZoneFromFile(background_mmax_path,
background_completeness_path,
filename_out=None,
layer2file=False,
ui_mode=True):
"""Load layer of background zones with completeness history and
Mmax from ASCII files, independent of QGis UI.
"""
# read data dicts from files (coordinates from mmax file)
background_mmax = readBackgroundMmax(background_mmax_path, ui_mode=ui_mode)
#QMessageBox.information(None, "Mmax", "%s" % background_mmax)
background_completeness = readBackgroundCompleteness(
background_completeness_path, ui_mode=ui_mode)
#QMessageBox.information(None, "Mc", "%s" % background_completeness)
# PostGIS SRID 4326 is allocated for WGS84
crs = QgsCoordinateReferenceSystem(
4326, QgsCoordinateReferenceSystem.PostgisCrsId)
# create layer
layer = QgsVectorLayer("Polygon", "Background Zones", "memory")
layer.setCrs(crs)
pr = layer.dataProvider()
# add attributes
# - ID, Mmax, completeness magnitude history
pr.addAttributes([
QgsField("ID", QVariant.String),
QgsField("mmax", QVariant.Double),
QgsField("mcdist", QVariant.String)
])
# add zones as features
for zone_id, zone in background_mmax.items():
f = QgsFeature()
poly = [QgsPoint(x, y) for (x, y) in zone['coord']]
f.setGeometry(QgsGeometry.fromPolygon([poly]))
f[0] = QVariant(zone_id)
f[1] = QVariant(zone['mmax'])
f[2] = QVariant(background_completeness[zone_id]['mcdist'])
pr.addFeatures([f])
# 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 loadBackgroundZoneFromFile(background_mmax_path,
background_completeness_path, filename_out=None, layer2file=False,
ui_mode=True):
"""Load layer of background zones with completeness history and
Mmax from ASCII files, independent of QGis UI.
"""
# read data dicts from files (coordinates from mmax file)
background_mmax = readBackgroundMmax(background_mmax_path, ui_mode=ui_mode)
#QMessageBox.information(None, "Mmax", "%s" % background_mmax)
background_completeness = readBackgroundCompleteness(
background_completeness_path, ui_mode=ui_mode)
#QMessageBox.information(None, "Mc", "%s" % background_completeness)
# PostGIS SRID 4326 is allocated for WGS84
crs = QgsCoordinateReferenceSystem(4326,
QgsCoordinateReferenceSystem.PostgisCrsId)
# create layer
layer = QgsVectorLayer("Polygon", "Background Zones", "memory")
layer.setCrs(crs)
pr = layer.dataProvider()
# add attributes
# - ID, Mmax, completeness magnitude history
pr.addAttributes([QgsField("ID", QVariant.String),
QgsField("mmax", QVariant.Double),
QgsField("mcdist", QVariant.String)])
# add zones as features
for zone_id, zone in background_mmax.items():
f = QgsFeature()
poly = [QgsPoint(x, y) for (x, y) in zone['coord']]
f.setGeometry(QgsGeometry.fromPolygon([poly]))
f[0] = QVariant(zone_id)
f[1] = QVariant(zone['mmax'])
f[2] = QVariant(background_completeness[zone_id]['mcdist'])
pr.addFeatures([f])
# 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