Install modified rasterstats from https://github.com/AsgerPetersen/python-rasterstats/tree/boundless to use this script as is.
Or just modify the bool in this line https://github.com/perrygeo/python-rasterstats/blob/master/src/rasterstats/io.py#L319
Script output:
❯ time python test.py
Using boundless: True
python test.py 72.11s user 9.82s system 99% cpu 1:21.96 total