cosmo map pca hdr

function [pca_ds, pca_params] = cosmo_map_pca(ds, varargin)