After that, it's like a GIS: you use the modules 1) to open, save the shapefiles and carry out the treatments with other modules like numpy or scipy, if you want. you can use shapely for manipulation ...
from [Python Script examples for geoprocessing shapefiles without using arcpy](http://gis.stackexchange.com/questions/64146/python-script-examples-for-geoprocessing ...
The aim of this step-by-step procedure is the generation of a burn severity map for the assessment of the areas affected by wildfires. The Normalized Burn Ratio (NBR) is used as it was designed to ...