Move layer creation code to separate function. Add tests.
Add clip_by_polygon function
Clip roads by the extent before impact calculation