Same as read_shapefile but directly transformed to a grid.
The whole thing is cached so that the second call will will be much faster.