A computer implemented method for the detection of features such as faults or channels in seismic images. First, edges are detected in a smoothed seismic image (106). To detect a fault line, an image intensity of the edges is projected in multiple spatial directions, for example by performing a Radon transform (118). The directions of maximum intensity are used to define a fault line (124c). To detect channels, smooth curves are identified within the detected edges (810). Sets of parallel smooth curves (812a) are then identified and used to define channels (812).