NEOWISE images of the 3.4 µm emission towards IRAS 20126+4104 at 21 epochs. The number in the bottom right corner of each panel indicates the observation date given in the first column of Table 1. The offsets are relative to position α(J2000)=20^h14^m26.s0364, δ(J2000)=41°13′32.″516. The dashed and dotted patterns outline the polygons used to derive the mid-IR flux densities of the SE (redshifted) lobe and NW (blueshifted) lobe of the bipolar outflow.