We have developed a novel method for random noise attenuation in seismic data by applying regularized nonstationary autoregression (RNA) in the frequency-space () domain. The method adaptively predicts the signal with spatial changes in dip or amplitude using RNA. The key idea is to overcome the assumption of linearity and stationarity of the signal in conventional domain prediction technique. The conventional domain prediction technique uses short temporal and spatial analysis windows to cope with the nonstationary of the seismic data. The new method does not require windowing strategies in spatial direction. We implement the algorithm by an iterated scheme using the conjugate-gradient method. We constrain the coefficients of nonstationary autoregression (NA) to be smooth along space and frequency in the domain. The shaping regularization in least-square inversion controls the smoothness of the coefficients of RNA. There are two key parameters in the proposed method: filter length and radius of shaping operator. Tests on synthetic and field data examples showed that, compared with domain and time-space domain prediction methods, RNA can be more effective in suppressing random noise and preserving the signals, especially for complex geological structure.