In a method of chromagenic illuminant estimation pixels from mutually-corresponding images with different filtering are compared, a fraction of the brightest pixels being selected for a subsequent chromagenic estimation. The pixels may be at corresponding locations or they may correspond in that their mean brightness is in the same rank order. In one method, in which, in a first preprocessing stage, for a database of m lights Ei(λ) and n surfaces Sj(λ) there is calculated Ti˜QFQ+ where Q1F and QF represent the matrices of unfiltered and filtered sensor responses to the n surfaces under the i th light and + denotes an inverse, and in a second operation stage, given P surfaces in an image and 3×P matrices Q and QF, from these matrices there are chosen the r % brightest pixels giving the matrices Q′ and Q′F, and the scene illuminant Pest is estimated where formula (I) and (II).