if row.crime_type == "BURGLARY":
times.append(row.datetime)
x, y = proj(*row.point)
xcs.append(x)
ycs.append(y)
return open_cp.data.TimedPoints.from_coords(times, xcs, ycs)
def load_geometry():
"""Load the geometry for Chicago; we'll use Southside, as ever..."""
open_cp.sources.chicago.set_data_directory(datadir)
return open_cp.sources.chicago.get_side(side)
with scripted.Data(load_points,
load_geometry,
start=datetime.datetime(2016, 1, 1)) as state:
time_range = scripted.TimeRange(datetime.datetime(2016, 10, 1),
datetime.datetime(2017, 1, 1),
datetime.timedelta(days=1))
state.add_prediction(scripted.NaiveProvider, time_range)
provider = predictors.ExpDecayGridProvider(state.grid, state.timed_points,
time_range.first)
state.add_prediction(provider, time_range)
for bw in [0.05, 0.15, 0.3, 1]:
provider = predictors.NonParamGridProvider(state.grid,
state.timed_points, bw,
with bz2.open("points2.pic.bz") as f:
return pickle.load(f)
def load_geometry():
with bz2.open("geo.pic.bz") as f:
return pickle.load(f)
import datetime, itertools
grid = open_cp.data.Grid(xsize=150, ysize=150, xoffset=0, yoffset=0)
grid = open_cp.geometry.mask_grid_by_intersection(load_geometry(), grid)
with scripted.Data(load_points,
load_geometry,
start=datetime.datetime(2007, 1, 1),
grid=grid) as state:
time_range = scripted.TimeRange(datetime.datetime(2007, 10, 1),
datetime.datetime(2008, 1, 1),
datetime.timedelta(days=1))
state.add_prediction(scripted.NaiveProvider, time_range)
for bw in range(50, 301, 10):
weight = open_cp.retrohotspot.Quartic(bw)
state.add_prediction(scripted.RetroHotspotProvider(weight), time_range)
for bw in range(50, 301, 10):
weight = open_cp.retrohotspot.Quartic(bw)
state.add_prediction(scripted.RetroHotspotCtsProvider(weight),
def load_geometry():
frame = gpd.read_file("SouthSide")
return frame.geometry[0]
def stringToDatetime(s):
return datetime.datetime(*(s.split("-")))
earliest_time = "2016-01-01"
start_time = "2016-10-01"
end_time = "2017-01-01"
# Perform the predictions; see `scripted_intro.md`
with scripted.Data(load_points,
load_geometry,
start=stringToDatetime(earliest_time)) as state:
time_range = scripted.TimeRange(stringToDatetime(start_time),
stringToDatetime(end_time),
datetime.timedelta(days=1))
state.add_prediction(scripted.NaiveProvider, time_range)
state.score(scripted.HitRateEvaluator)
state.score(scripted.HitCountEvaluator)
state.process(scripted.HitRateSave("ratesTest.csv"))
state.process(scripted.HitCountSave("countsTest.csv"))