Use FSRMatrix sparse storage for all projection-related operations.
Allow Operator keywords in fromfits.