Feature extraction is an important factor in determining the accuracy that can be attained in the classification of multispectral images. The traditional per point classification methods do not use all the available information, since they disregard the spatial relationships that exist among pixels belonging to the same class. In this paper, methods are developed to extract additional image spatial features by means of linear and non-linear local filtering. Feature selection methods are also developed, since it is usually not possible to use all the generated features. The classification stage is performed in a supervised mode using the maximum likelihood criterion. A quantitative analysis of the performance of the spatial features show that an overall increase in precision of classification is achieved, although at the expense of increased rejection levels, particularly on the borders between different fields.