#develop: vrt.transform_points
#develop nansat.crop(): use the original GCPs in ROI. If it is not enough, create and add 100 GCPs.