estimates suburbs to a square

suburbs.py

import json

with open('wa.json') as f:
	data = json.load(f)
out = {}
for feature in data['features']:
	out[feature["properties"]["wa_local_2"].lower()] = {}
	out[feature["properties"]["wa_local_2"].lower()]["lat"] = [min([i[1] for i in feature["geometry"]["coordinates"][0][0]]), max([i[1] for i in feature["geometry"]["coordinates"][0][0]])]
	out[feature["properties"]["wa_local_2"].lower()]["lon"] = [min([i[0] for i in feature["geometry"]["coordinates"][0][0]]), max([i[0] for i in feature["geometry"]["coordinates"][0][0]])]

a = open("suburbs.json", "w")
json.dump(out, a)

wa.json

downloaded from here -> https://data.gov.au/dataset/ds-dga-6a0ec945-c880-4882-8a81-4dbcb85e74e5/details