We present a processing method, based on the multivariate curve resolution approach (MCR), to denoise 2D solid-state NMR spectra, yielding a substantial S/N ratio increase while preserving the lineshapes and relative signal intensities. These spectral features are particularly important in the quantification of silicon species, where sensitivity is limited by the low natural abundance of the 29Si nuclei and by the dilution of the intrinsic protons of silica, but can be of interest also when dealing with other intermediate-to-low receptivity nuclei. This method also offers the possibility of coprocessing multiple 2D spectra that have the signals at the same frequencies but with different intensities (e.g.: as a result of a variation in the mixing time). The processing can be carried out on the time-domain data, thus preserving the possibility of applying further processing to the data. As a demonstration, we have applied Cadzow denoising on the MCR-processed FIDs, achieving a further increase in the S/N ratio and more effective denoising also on the transients at longer indirect evolution times. We have applied the combined denoising on a set of experimental data from a lysozyme-silica composite.