Summary Imaging dense clusters of seismicity is crucial to many problems in seismology: to delineate complex systems of faults, provide constraints on the causes of volcanic and cryogenic swarms, and to shed light on possible means to prevent damaging induced seismicity in mining, geothermal and oil and gas extraction activities. Current imaging methods rely upon high-resolution relative location techniques, commonly requiring arrival-time picks for seismic phases. This paper examines an alternative approach, based upon concepts drawn from differential geometry, that images directly from waveform data. It relies upon the common assumption of spatial continuity of seismic wavefield observations, which implies that a differentiable map exists between the source region to be imaged and waveform observations considered as elements of a vector space. The map creates an image of event clusters on a Riemannian manifold embedded in that vector space. The image can be visualized by projecting the observations into a tangent space of the manifold and is a distorted rendering of cluster geometry. However, the distortion can be predicted and removed if a model for wavefield propagation is available. This visualization approach is applicable to clusters of uniform events with highly similar waveforms, such as are commonly acquired with correlation detectors or other pattern matching techniques. To assess its performance, it is applied to the closely-related reciprocal problem of imaging the (known) geometry of an array from observations by the array of several regional events. Differences between the original problem and its reciprocal analog are noted and controlled for in the analysis. Chief among the differences is the necessity for aligning the waveforms in the original problem, which, to maintain consistency with the original problem, is solved in the reciprocal problem by a generalization of the VanDecar-Crosson algorithm. The VanDecar-Crosson algorithm exhibits a bias, shown through an analysis of the situation when the observed wavefields are adequately modeled as plane waves. In that circumstance, the bias can be predicted and removed. In a test using a portion of a large-N array, this imaging approach is shown to successfully reconstruct the array geometry. The method is applicable directly to infinitesimal array apertures, but is extended to a larger aperture by partitioning the image into local, effectively infinitesimal overlapping subsets. These are inverted, then assembled into a global picture of the array geometry using constraints provided by the overlapped regions. Although demonstrated in a reciprocal array context, the method appears viable for imaging clusters of events with highly similar source mechanisms and time histories.