Assume that $X\sim \mathcal N(\sigma_1,\mu_1)$ and $Y\sim \mathcal N(\sigma_2,\mu_2)$.
I want to estimate $\frac{\mu_1+\mu_2}{2}$ after observing $X,Y$.
In my setting, $\sigma_1,\sigma_2$ are known and we want to estimate the average of the means (e.g., what is the MLE of it?).
In the special case where we know that $\mu=\mu_1=\mu_2$ (i.e., we estimate the same quantity from observations with different variances), it is known that the MLE is:
$$ X\cdot \frac{\sigma_2^2}{\sigma_1^2+\sigma_2^2} + Y\cdot \frac{\sigma_1^2}{\sigma_1^2+\sigma_2^2}. $$
How does this generalize to arbitrary $\mu_1,\mu_2$?