I think my question breaks down into two parts.
Let's say you have a two qubit system and you can perform projective measurements.
Each round of measurements will consist of results looking like below.
\begin{array}{|c|c|c|c|} \hline & \sigma_x \otimes \sigma_x& \sigma_x \otimes \sigma_y & \sigma_x \otimes \sigma_z& \sigma_y \otimes \sigma_x &\sigma_y \otimes \sigma_y&\sigma_y \otimes\sigma_z&\sigma_z \otimes \sigma_x&\sigma_z \otimes \sigma_y&\sigma_z \otimes \sigma_z \\ \hline round 1 &\uparrow \downarrow &\uparrow \uparrow &\downarrow \uparrow &\uparrow \uparrow&\downarrow \uparrow&\downarrow \downarrow&\uparrow \downarrow&\downarrow \uparrow& \uparrow\uparrow\\ \hline \end{array}
You repeat multiple rounds of this sort of measurements.
So you have experimentally acquired values of $\langle \sigma_i \otimes \sigma_j \rangle$ (expectation values of products of Pauli matrices) with i and j being x,y, and z, not including the identity.
Given what I found in this paper, I think you can construct the density matrix of the two-qubit system if you have expectation values of all possible products of Pauli matrices, including the identity matrix.
So you need to know $\langle \sigma_I \otimes \sigma_I \rangle$, $\langle \sigma_I \otimes \sigma_x \rangle$, $\langle \sigma_I \otimes \sigma_y \rangle$, $\langle \sigma_I \otimes \sigma_z \rangle$, $\langle \sigma_x \otimes \sigma_I \rangle$,$\langle \sigma_x \otimes \sigma_x \rangle$,$\langle \sigma_x \otimes \sigma_y \rangle$, $\langle \sigma_x \otimes \sigma_z \rangle$,$\langle \sigma_y \otimes \sigma_I \rangle$ ,$\langle \sigma_y \otimes \sigma_x \rangle$, $\langle \sigma_y \otimes \sigma_y \rangle$, $\langle \sigma_y \otimes \sigma_z \rangle$, $\langle \sigma_z \otimes \sigma_I \rangle$, $\langle \sigma_z \otimes \sigma_x \rangle$, $\langle \sigma_z \otimes \sigma_y \rangle$, $\langle \sigma_z \otimes \sigma_z \rangle$, for full tomograhy.
Can I just get away with measuring expectation values of the products of Pauli operators without the identity matrix to figure out expectation values of all 16 products of two Pauli operators?
Can I construct the full density matrix of a two-qubit system if I know the expectation values of all 16 products of two Pauli operators?