Modelling spatial heteroskedasticity by volatility modulated moving averages