AGN contour plots obtained from the core-shift-determining sessions, and the AGN jet directions determined with arcfits. For the jet direction determinations, image models obtained at higher observing frequencies are preferred, unless extended jet features are not identified by arcfits. The flux density on the m-th (m = 1,2,3,…) contour equals to $3\cdot \text{rms}\cdot {(\sqrt{2})}^{m-1}$ , where rms is 0.18 mJy beam⁻¹, 7.3 mJy beam⁻¹, 0.50 mJy beam⁻¹ and 0.65 mJy beam⁻¹ for J2221, J2218, J2219 and J2226, respectively. Δα and Δδ are, respectively, relative right ascension (left to the east) and relative declination in the unit of image pixel. The reference position of each image is provided over the image. In each panel, the gray ellipse shows the size of the synthesized beam; the AGN jet direction is marked with the orange arrow, while the uncertainty of the direction is illustrated with the two dashed lines.