We estimate the rate of aftershocks triggered by a heterogeneous stress change, using the rate-and-state model of Dieterich. We show that an exponential stress distribution Pτ(τ) ~ exp(-τ/τ0) gives an Omori law decay of aftershocks with time ~1/tp, with an exponent p = 1 - Aσn/τ0, where A is a parameter of the rate-and-state friction law and σn is the normal stress. Omori exponent p thus decreases if the stress "heterogeneity" τ0 decreases. We also invert the stress distribution Pτ(τ) from the seismicity rate R(t), assuming that the stress does not change with time. We apply this method to a synthetic stress map, using the (modified) scale invariant "k2" slip model (Herrero and Bernard). We generate synthetic aftershock catalogs from this stress change. The seismicity rate on the rupture area shows a huge increase at short times, even if the stress decreases on average. Aftershocks are clustered in the regions of low slip, but the spatial distribution is more diffuse than for a simple slip dislocation. Because the stress field is very heterogeneous, there are many patches of positive stress changes everywhere on the fault. This stochastic slip model gives a Gaussian stress distribution but nevertheless produces an aftershock rate which is very close to Omori's law, with an effective p ≤ 1, which increases slowly with time. We obtain a good estimation of the stress distribution for realistic catalogs when we constrain the shape of the distribution. However, there are probably other factors which also affect the temporal decay of aftershocks with time. In particular, heterogeneity of Aσn can also modify the parameters p and c of Omori's law. Finally, we show that stress shadows are very difficult to observe in a heterogeneous stress context.