Abstract. We propose a two-step approach for the analysis of functional magnetic resonance images, in the context of natural stimuli. In the first step, elements of functional brain activity emerge, based on spatial independence assumptions. The second step exploits temporal covariation between the elements and given features of the natural stimuli to identify functional networks. The networks can have complex activation patterns related to common task goals.