tests/app/broadcast_areas/custom_polygons.py

BRISTOL = [
    [51.4371, -2.6216],
    [51.4371, -2.5750],
    [51.4668, -2.5750],
    [51.4668, -2.6216],
    [51.4371, -2.6216],
]

SKYE = [
    [57.1004, -6.8280],
    [57.1004, -5.7733],
    [57.7334, -5.7733],
    [57.7334, -6.8280],
    [57.1004, -6.8280],
]
Estimate number of phones in an arbitrary polygon We want to know how many phones are in a user-supplied polygon, so we can show the impact of a broadcast, in the same way that we do when users pick areas from our library. We already know how many phones are in each electoral ward. But there are challenges with an arbitrary polygon: - where it does overlap a ward, the overlap could be partial - it could overlap more than one ward - finding out which wards it overlaps by brute force (looping through all the wards and seeing which ones intersect with our polygon) would be way to slow to do in real time Instead we can use a data structure called an R-tree[1] to build an index which provides a much, much faster way of looking up which polygons overlap another. We can build this tree in advance and save it somewhere, which means there’s a lot of computation we don’t need to do in real time. The R-tree returns a set of objects (ward IDs) which we can go and look up in our library of electoral wards. These wards will be the ones that might have some overlap with our custom polygon. Once we have this small set of wards which might overlap our ward, we can look at the size of the area of overlap (relative to the size of the whole ward) and multiply that by the known count of phones in that ward to get an approximation of the count of phones in the overlap area. Summing these approximations give an estimate for the whole area of the custom polygon. 1. https://en.wikipedia.org/wiki/R-tree 2021-03-18 23:02:32 +00:00			`BRISTOL = [`
			`[51.4371, -2.6216],`
			`[51.4371, -2.5750],`
			`[51.4668, -2.5750],`
			`[51.4668, -2.6216],`
			`[51.4371, -2.6216],`
			`]`

			`SKYE = [`
			`[57.1004, -6.8280],`
			`[57.1004, -5.7733],`
			`[57.7334, -5.7733],`
			`[57.7334, -6.8280],`
			`[57.1004, -6.8280],`
			`]`