More precise beam hardening correction is performed at low computational loads to objects having a wider range of effective atomic number (Z eff ) elements, contributing to the presentation of more quantitative X-ray images. X-ray loss amount (μt) of two X-ray energy BINs selected from two or more X-ray energy BINs is normalized, and at least one standardized X-ray loss amount is obtained for each pixel area, and the standardized X-ray loss amount and elemental From reference information indicating a theoretical correspondence with the effective atomic number, at least one effective atomic number is estimated for each pixel region, and the estimated at least one effective atomic number (Z High , Z Low ) and beam hardening correction are performed. When doing so, the coincidence state of at least two atomic numbers among the effective atomic numbers Zm specified in advance is determined.보다 고정밀한 빔하드닝 보정을, 보다 광범위한 실효 원자 번호(Zeff)의 원소를 가지는 대상까지 저연산 부하로 실시해, 보다 정량적인 X선 화상의 제시에 기여한다. 2개 이상의 X선 에너지 BIN 중 선택되는 2개의 X선 에너지 BIN의 X선 감약량(μt)을 규격화하여, 적어도 하나의 규격화 X선 감약량을 화소 영역 마다 구하고, 규격화 X선 감약량과 원소의 실효 원자 번호와의 이론적인 대응 관계를 나타내는 참조 정보로부터, 적어도 하나의 실효 원자 번호를 화소 영역 마다 추정하여, 추정되는 적어도 하나의 실효 원자 번호(ZHigh, ZLow)와 빔하드닝 보정을 실시할 때에 미리 지정되는 실효 원자 번호(Zm) 중 적어도 2개의 원자 번호의 일치 상태를 판정한다.