Minimal python code to run use gdal for layer reading