Common spatial pattern (CSP) method is widely used in brain machine interface (BMI) applications to extract features from the multichannel neural activity through a set of spatial projections. These spatial projections minimize the Rayleigh quotient (RQ) as the objective function, which is the variance ratio of the classes. The CSP method easily overfits the data when the number of training trials is not sufficiently large and itis sensitive to daily variation of multichannel electrode placement, which limits its applicability for everyday use in BMI systems. To overcome these problems, the amount of channels that is used in projections, should be limited to some adequate number. We introduce a spatially sparse projection (SSP) method that exploits the unconstrained minimization of a new objective function with approximated 1 penalty. Unlike the RQ, this new objective function depends on the magnitude of the sparse filter. The SSP method is employed to classify the multiclass ECoG and two class EEG data sets. We compared our results with a recently introduced sparse CSP solution based on 0 norm. Our method outperforms the standard CSP method and provides comparable results to 0 norm based solution and it is associated with less computational complexity. We also conducted several simulation studies on the effect of noisy channel and intersession variability on the performance of the CSP and sparse filters.