Accurate postmortem interval estimation is vital in the investigation of homicides, suicides, and accidental deaths. It is key in narrowing suspect lists, improving crime-solving efficiency, and offering solace to bereaved families. The intra-puparial period, comprising about half of a fly’s developmental cycle, presents challenges for morphological age estimation. External changes are limited to color shifts and the appearance of respiratory horns on the puparium only within several hours after pupariation, while detailed internal development analysis often requires invasive methods like removing the puparium, which can be damaging. Additionally, these techniques usually depend on a forensic entomologist’s expertise, which lead to subjective biases. This study employed attenuated total reflection-fourier transform infrared spectroscopy, a rapid, non-destructive method for analyzing proteins, chitosan, and chitin in puparia. Data showed a consistent reduction in the concentration of the amide I band within the puparium during the intra-puparial development at five constant temperatures (19 °C, 22 °C, 25 °C, 28 °C and 31 °C). This trend in the spectral data effectively distinguishes pupae at various stages of intra-puparial development, facilitating precise age estimation, which is critical for the estimation of the minimum postmortem interval (PMImin). Finally, this work combined the total reflection-fourier transform infrared spectroscopy with chemometric analysis and successfully developed a partial least squares discriminant analysis model and a random forest model, with accuracies of 88 % and 81 %, respectively. These models enable the non-invasive age estimation of P. regina in its intra-puparial period, a stage traditionally difficult to assess morphologically, thus laying the groundwork for PMImin estimation using fly pupae.