search

bcgov / ols-devkit

Developer Tool Kit using the BC Address Geocoder, BC Route Planner and the Geomark Web Service
Apache License 2.0
6 stars 5 forks source link

Write a script to compute admin area of each DRA locality #187

Open mraross opened 3 years ago

mraross commented 3 years ago

Create a python script called somthing like findDraLocalityAdminAreas.py

The script should read the following inputs:

  1. DRA locality coverage containing polygons for each dra locality.
  2. Admin Area coverage containing polygons for a given set of admin areas (e.g., Community Health Service Areas).
  3. A parameter called adminAreaType which contains the official code of the set of admin areas in the coverage (e.g., "CHSA")

The script should perform a spatial overlay of each DRA locality polygon against the admin area coverage to determine the list of intersecting admin areas.

Script should output the following for each dra locality:

  1. localityName
  2. adminAreaType
  3. adminAreaCode

If there is more than one admin area intersecting a locality, set adminAreaCode to a special character like asterisk ("") to represent "more than one". For example, if the dra locality polygon for Victoria intersects five community health service areas, set Victoria's adminAreaCode to "", If the dra locality polygon for Ahousat intersects a single CHSA called "NBAH", set Ahousat's adminAreaCode to "NBAH" If there is no intersecting admin area for a given dra locality polygon, set adminAreaCode to the empty string ("")