In this study, we present an advanced nanofabrication approach to produce gradient-index photonic crystal structures based on nanoporous anodic alumina. An apodization strategy is for the first time applied to a sinusoidal pulse anodisation process in order to engineer the photonic stop band of nanoporous anodic alumina (NAA) in depth. Four apodization functions are explored, including linear positive, linear negative, logarithmic positive and logarithmic negative, with the aim of finely tuning the characteristic photonic stop band of these photonic crystal structures. We systematically analyse the effect of the amplitude difference (from 0.105 to 0.840 mA cm-2), the pore widening time (from 0 to 6 min), the anodisation period (from 650 to 950 s) and the anodisation time (from 15 to 30 h) on the quality and the position of the characteristic photonic stop band and the interferometric colour of these photonic crystal structures using the aforementioned apodization functions. Our results reveal that a logarithmic negative apodisation function is the most optimal approach to obtain unprecedented well-resolved and narrow photonic stop bands across the UV-visible-NIR spectrum of NAA-based gradient-index photonic crystals. Our study establishes a fully comprehensive rationale towards the development of unique NAA-based photonic crystal structures with finely engineered optical properties for advanced photonic devices such as ultra-sensitive optical sensors, selective optical filters and all-optical platforms for quantum computing.