def test_csv2building_polygons(self):
"""test_csv2building_polygons
"""
path = get_pathname_from_package('anuga.shallow_water')
testfile = os.path.join(path, 'tests', 'data', 'polygon_values_example.csv')
polygons, values = load_csv_as_building_polygons(testfile,
floor_height=3)
assert len(polygons) == 7, 'Must have seven polygons'
assert len(values) == 7, 'Must have seven values'
# Known floor values
floors = {'1': 6, '2': 0, '3': 3, '4': 6, '5': 0, '8': 3, '9': 3}
for id in ['1', '2', '3', '4', '5' ,'8' ,'9']:
assert id in polygons.keys()
assert id in values.keys()
assert float(values[id]) == float(floors[id])
filename=project.combined_elevation+'.pts',
use_cache=True,
verbose=True,
alpha=project.alpha)
Stage = domain.quantities['stage']
Elev = domain.quantities['elevation']
domain.set_quantity('stage', 0.8, use_cache=False, verbose=True)
if project.use_buildings:
# Add buildings from file
log.critical('Reading building polygons')
building_polygons, building_heights = load_csv_as_building_polygons(project.building_polygon, floor_height=3.)
#clipping_polygons=project.building_area_polygons)
log.critical('Creating %d building polygons' % len(building_polygons))
def create_polygon_function(building_polygons, geo_reference=None):
L = []
for i, key in enumerate(building_polygons):
if i%100==0: log.critical(i)
poly = building_polygons[key]
elev = building_heights[key]
L.append((poly, elev))
buildings = anuga.Polygon_function(L, default=0.0,
geo_reference=geo_reference)
return buildings
domain.set_quantity('elevation',
filename=project.combined_elevation + '.pts',
use_cache=True,
verbose=True,
alpha=project.alpha)
Stage = domain.quantities['stage']
Elev = domain.quantities['elevation']
domain.set_quantity('stage', 0.8, use_cache=False, verbose=True)
if project.use_buildings:
# Add buildings from file
log.critical('Reading building polygons')
building_polygons, building_heights = load_csv_as_building_polygons(
project.building_polygon, floor_height=3.)
#clipping_polygons=project.building_area_polygons)
log.critical('Creating %d building polygons' % len(building_polygons))
def create_polygon_function(building_polygons, geo_reference=None):
L = []
for i, key in enumerate(building_polygons):
if i % 100 == 0: log.critical(i)
poly = building_polygons[key]
elev = building_heights[key]
L.append((poly, elev))
buildings = anuga.Polygon_function(L,
default=0.0,
geo_reference=geo_reference)
