Data-Driven PSF Modeling using B-Splines