  We build a model for the density and integrated Sachs-Wolfe (ISW) profile of
supervoid and supercluster structures. Our model assumes that fluctuations
evolve linearly from an initial Gaussian random field. We find these
assumptions capable of describing N-body simulations and simulated ISW maps
remarkably well on large scales. We construct an ISW map based on locations of
superstructures identified previously in the SDSS Luminous Red Galaxy sample. A
matched filter analysis of the cosmic microwave background confirms a signal at
the $3.2-\sigma$ confidence level and estimates the radius of the underlying
structures to be $55 \pm 28 ~h^{-1}$Mpc. The amplitude of the signal, however,
is $2-\sigma$ higher than $\Lambda$CDM predictions.
