Recursive image subsampling which yields support areas approaching fractals is described and analyzed using iterated function systems. The subsampling scheme is suitable in, e.g., hierarchical image processing and image coding schemes. For hexagonally sampled images a hierarchical subsampling structure is given which yields hexagon-like regions with fractal borders.