A method of local feature identification comprises retrieving data representing a lightfield, forming an epipolar volume (V) from the retrieved data, applying a transform to epipolar planes retrieved from the epipolar volume, so as to represent the epipolar planes in a new space, identifying a plurality of epipolar lines, and identifying local features (fp) based on at least some of the epipolar lines.