An image processing device that is used in a radiography apparatus performing tomosynthesis imaging by irradiating a breast that is placed on an imaging table and is compressed by a compression member with radiation emitted from a radiation source at each of a plurality of irradiation positions (80 1 to 80 13, 80 H, 80 L) having different irradiation angles to capture projection images of the compressed breast at each of the irradiation positions (80 1 to 80 13, 80 H, 80 L), the image processing device comprising: at least one processor that is configured to: acquire overall information corresponding to the thickness and the area of the compressed breast and determine an overall imaging irradiation angle (AR a) on the basis of the acquired overall imaging information, wherein the overall imaging irradiation angle range is a maximum irradiation angle range in which the tomographic image including the entire compressed breast, using the projection images obtained at each of the plurality of irradiation positions, can be obtained; acquire a projection image group including a plurality of projection images (84 1 to 84 13) obtained by the tomosynthesis imaging in an irradiation angle range wider than the overall imaging irradiation angle range (AR a); generate a plurality of first tomographic images (86 1) including a part of the compressed breast, using a first group of projection images (84 1 to 84 13) obtained by the tomosynthesis imaging in a first irradiation angle range (AR 1) wider than the overall imaging irradiation angle range (AR a) among the projection images (84 1 to 84 13) included in the projection image group; and generate a plurality of second tomographic images (86 2) including the entire compressed breast, using a second group of projection images (84 5 to 84 9) obtained by the tomosynthesis imaging in a second irradiation angle range (AR 2) that is equal to or narrower than the overall imaging irradiation angle range (AR a) among the projection images (84 1 to 84 13) included in the projection image group. acquire overall information corresponding to the thickness and the area of the compressed breast and determine an overall imaging irradiation angle (AR a) on the basis of the acquired overall imaging information, wherein the overall imaging irradiation angle range is a maximum irradiation angle range in which the tomographic image including the entire compressed breast, using the projection images obtained at each of the plurality of irradiation positions, can be obtained; acquire a projection image group including a plurality of projection images (84 1 to 84 13) obtained by the tomosynthesis imaging in an irradiation angle range wider than the overall imaging irradiation angle range (AR a); generate a plurality of first tomographic images (86 1) including a part of the compressed breast, using a first group of projection images (84 1 to 84 13) obtained by the tomosynthesis imaging in a first irradiation angle range (AR 1) wider than the overall imaging irradiation angle range (AR a) among the projection images (84 1 to 84 13) included in the projection image group; and generate a plurality of second tomographic images (86 2) including the entire compressed breast,