Hi,
To compute the msd with images, provide the images directly from the gsd trajectory. The script provided was very close, something like the following will provide the right input for the images argument:
```
for i, frame in enumerate(traj):
positions[i] = frame.particles.position
images[i] = frame.particles.image
msd = freud.msd.MSD(box, mode='window')
msd.compute(positions=positions, images=images)
```
Best,
Tommy