Signals from a plurality of sources are observed by a plurality of sensors disposed in two dimensions, and the observed signals are subject to a short-time Fourier transform, from which frequency domain signals are derived. Using the independent component analysis process on the frequency domain signals, separation matrices are produced, and an inverse matrix of each separation matrix is calculated. The direction of each source is calculated based on a ratio of a pair of elements in each column of the inverse matrix.