Fish length frequency histograms from research surveys are of prime importance for identifying habitats of different life stages, as well as for stock assessment. However, no method has thus far been available for mapping these histograms as spatially varying curves. Here, a procedure is applied to map spatially connected curves, and detail is given on how it can be applied to map the length frequency histograms. At each sample location, a fish length frequency histogram is given as a vector of non-independent values. The histogram is first modelled as a polynomial expansion on the basis of orthogonal polynomials. Then, the polynomial coefficients are mapped by co-kriging, after fitting a model of co-regionalization. The length frequency map is finally derived by linearly combining maps of polynomial coefficients. An estimation variance associated with the map is also derived. Maps of anchovy length distributions are produced by applying the method to midwater trawl length data from the PELGAS acoustic surveys in the Bay of Biscay. This novel approach extends the application of kriging techniques to curves or functions, opening new perspectives for mapping more complex information than just the values of fish density.