Full 3D CMB bispectrum recovered from the Planck foreground-cleaned maps, including SMICA (left), NILC (centre) and SEVEM (right), using the hybrid Fourier mode coefficients illustrated in Fig. 9, These are plotted in three-dimensions with multipole coordinates { ℓ₁,ℓ₂,ℓ₃ } on the tetrahedral domain shown in Fig. 1 out to ℓ_max = 2000. Several density contours are plotted with red positive and blue negative. The bispectra extracted from the different foreground-separated maps appear to be almost indistinguishable.