Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

  1. Open a new project in QGIS

  2. Add the FEWSGrid.txt as an delimited text layer to a empty project (Layer -> Add Layer -> Add Delimited Text Layer). In the popup window select 'First record has field names', e.g.:

  3. Add the shape files (Layer -> Add Layer -> Add Vector Layer). After this the screen should look like this: 

  4. Introduce unique IDs for each sub-area (or "sector") and each receptor
    1. Select a layer and open its attribute table (Right click -> Open attribute table)
    2. 'Toggle editing mode' (pencil icon)
    3. 'Open field calculator'
      1. Check the 'Create a new field' box and choose an 'output field name' (receptorID or areaID depending on the layer you selected)
      2. Set 'Expression' to $rownum
      3. Click 'Ok'
    4. 'Save edits' and once more 'Toggle editing mode'
    5. Repeat a-d for the other layer
  5. Intersect the receptor layer with the area layer (Vector -> Geoprocessing Tools -> Intersect) and save the new layer under a new name (in this example ResBuildings_hotspot.shp)
    The resulting layer will contain all the buildings within the hotspot and associate them with the areaID of the area there in. 
  6.  Create a new polygon shape file from the FEWSGrid layer (Vector -> Geometry tools -> Voronoi Polygons) and save under a new name (in this example FEWSGrid_polygon.shp). 
  7. Intersect this FEWSGrid_polygon layer with the ResBuildings_hotspot layer. (Vector -> Geoprocessing Tools -> Intersect) and save under new name (e.g. ResBuilding_for_BN). Each polygon in the resulting shape layer is associated with a unique combination of gridID, areaID and receptorID. This information is used by the BN Adaptor to understand in which area a receptor is located and which model grid points should be used to calculated its hazard from the hazard.nc. The resulting layer should look similar to this:
  8. Save the columns areaID, receptorID and gridID (in this order!!) to a <receptor>.txt file and verify that the values in receptorID are in ascending order. 
    1. Right click on ResBuilding_for_BN -> Open Attribute Table
    2. 'Toggle editing mode'
    3. 'Delete column': all except areaID, receptorID and gridID
    4. 'Save edits', again 'Toggle editing mode' and close attribute table 
    5. Save layer as Comma Separated File (Right click on ResBuilding_for_BN -> Save As e.g. 'ResBuildings.csv')
    6. Open in excel and, if necessary, reorder the columns to areaID, receptorID and gridID
  9. Open ResBuildings.csv in Matlab for testing
  10. Save as <receptor>.txt file 

...