We introduce a new representation of cortical regions via distribution functions of their features. The distribution functions are estimated non-parametrically from the data and are observed to be non Gaussian. Cortical pattern matching is enabled by using the information-based Jensen-Shannon divergence as a measure between features. Our approach explicitly avoids pairwise registrations between brains, but instead focuses on modeling and discriminating between the cortical structural patterns. We demonstrate our approach on 120 subject brains from an Alzheimer's dataset, and present applications to clustering, classification, and dimension reduction.