It has been postulated that particular patterns of macular pigment (MP) distribution may be associated with the risk for eye diseases such as age-related macular degeneration (AMD). This work investigates the potential of Zernike polynomials (ZP) to characterise the level and distribution of MP, and their suitability as a representation for analysis of the effects of age and AMD on MP patterns. As the case study, MP distribution maps computed using an experimental method based on fundus reflectance (MRIA) were obtained for ninety volunteers representing three groups: under-fifty without AMD, fifty and over without AMD, and fifty and over with AMD. ZP with 105 coefficients were fitted to the maps using least-squares optimisation and found to represent MP maps accurately (RMSE<10-1). One-way MANOVA analysis carried out on ZP representations showed that the three subject groups have significantly different means (Wilk's Lambda 0.125, p<0.0001). Linear discriminant analysis with leave-one-out scheme resulted in accuracy, sensitivity and specificity of classification according to, respectively, disease status regardless of age (81% all); disease status in the age-matched groups (87%, 88%, 86%); age irrespective of disease status (81%, 83%, 73%); and age for subjects without AMD (83%, 88%, 80%). Mean MP distributions computed from ZP coefficients for the three groups showed more elevated and more peaked MP for the healthy under-fifty group; more irregular and more elevated peripheral levels in over-fifty AMD group than in over-fifty non-AMD group; and moderate radial asymmetry in non-AMD over-50 group. The results suggest that ZP coefficients are capable of accurately representing MP in a way that captures certain spatial patterns of its distribution. Using the ZP representation MP maps could be classified according to both age and disease status with accuracy significantly greater than chance, with peak elevation, pattern irregularity and radial asymmetry identified as important features.