A method of generating diffraction images based on wave equations includes generating a source wavefield and a receiver wavefield. Based on the source wavefield, a first source wavefield propagating in a first direction and a second source wavefield propagating in a second direction are generated. Based on the receiver wavefield, a first receiver wavefield propagating in the first direction and a second receiver wavefield propagating in the second direction are generated. A first seismic image is generated based on the first source wavefield and the first receiver wavefield. A second seismic image is generated based on the second source wavefield and the second receiver wavefield. A final seismic image is generated based on the first seismic image and the second seismic image.