I have 10 day temperature forecast data that is hourly initially and then every 3 hours. I would like to predict the hourly values for the full 10 days. Linear interpolation fails as sunrise and other factors (cloud cover, precipitation, etc.) significantly affect the temperature within that three hour period. For example if the temperature is 40 degrees at 5am and 50 degrees at 8am, the temperature at 6am may be 41 degrees and the temperature at 7am 49 degrees.
There are several methods available to forecast the temperature beyond the end of the period (such as SARIMAX and NeuralProphet), however these also don't work for the "missing" hourly data because they do not take into account the following temperature data point.
A spline is the best solution we have determined so far, but this does not account for things like cloud cover, where the temperature rise over a sunrise may in fact be more linear. PyPOTS and specifically the SAITS algorithm (https://github.com/WenjieDu/SAITS) is promising but I do not see a way to include exogenous variables.