In linearized gravity, two linearized metrics are considered gauge-equivalent, hµν ∼ hµν + Kµν [v], when they differ by the image of the Killing operator, Kµν [v] = ∇µvν + ∇ν vµ. A universal (or complete) compatibility operator for K is a differential operator K1 such that K1 • K = 0 and any other operator annihilating K must factor through K1. The components of K1 can be interpreted as a complete (or generating) set of local gauge-invariant observables in linearized gravity. By appealing to known results in the formal theory of overdetermined PDEs and basic notions from homological algebra, we solve the problem of constructing the Killing compatibility operator K1 on an arbitrary background geometry, as well as of extending it to a full compatibility complex Ki (i ≥ 1), meaning that for each Ki the operator Ki+1 is its universal compatibility operator. Our solution is practical enough that we apply it explicitly in two examples, giving the first construction of full compatibility complexes for the Killing operator on these geometries. The first example consists of the cosmological FLRW spacetimes, in any dimension. The second consists of a generalization of the Schwarzschild-Tangherlini black hole spacetimes, also in any dimension. The generalization allows an arbitrary cosmological constant and the replacement of spherical symmetry by planar or pseudo-spherical symmetry.