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, <EOS>