How to convert design matrix to cell means matrix for unbalanced two-way ANOVA?

Question

My goal is to create a cell means model matrix from scratch in r. Please help show how this can be achieved.

I prefer the manual way over per-defined