Compressive sensing of noisy multispectral images is considered in this letter. Multispectral images in remote sensing applications are multichannel and inherently noisy. An approach using Bregman split method for optimization in both spatial and transform domains is proposed. The performance of the proposed algorithm is evaluated by comparing with other approaches. It is shown that the proposed algorithm performs favorably compared with other approaches with noisy multispectral images in experiments.