function [pca_ds, pca_params] = cosmo_map_pca(ds, varargin)
cosmo map2surface hdr
cosmo mask dim intersect hdr