We adopt general relativistic ray-tracing (GRRT) schemes to study images of Kerr-MOG black holes surrounded by geometrically thick magnetized equilibrium tori, which belong to steady-state solutions of thick accretion disks within the framework of general relativistic magnetohydrodynamics (GRMHD). The black hole possesses an extra dimensionless MOG parameter described its deviation from usual Kerr one. Our results show that the presence of the MOG parameter leads to smaller disks in size, but enhances the total flux density and peak brightness in their images. Combining with observation data of black hole M87* from the Event Horizon Telescope (EHT), we make a constraint on parameters of the Kerr-MOG black hole and find that the presence of the MOG parameter broadens the allowable range of black hole spin.