function [Pox, invCov] = PixelStates(rock, ice)

R = length(rock);
I = length(ice);

for i = 1:R
   OR(i) = ObservedState(rock(i,:));
end

for i = 1:I
   OI(i) = ObservedState(ice(i,:));
end

Nox = zeros(13,2);
for i = 1:R
   Nox(OR(i),1) = Nox(OR(i),1) + 1;
end
Pox(:,1) = Nox(:,1)./sum(Nox(:,1));

for i = 1:I
   Nox(OI(i),2) = Nox(OI(i),2) + 1;
end
Pox(:,2) = Nox(:,2)./sum(Nox(:,2));

invCov = inv(cov([rock; ice]));
