1. An image reconstruction method of transmitting and receiving signals to and from a measurement region and determining a quantity of a physical property of the measurement region through numerical calculation using a measuring apparatus, the method comprising: a transmit-receive data acquisition step, by a transmit-receive data acquisition part of the measuring apparatus, of transmitting and receiving signals to and from the measurement region; a physical property value calculation step, by a physical property value calculation part of the measuring apparatus, of calculating a physical property value from an actual receive signal received from the measurement region; and a physical property value-display value conversion step, by a physical property value-display value conversion part of the measuring apparatus, of converting the physical property value into a display value, wherein the physical property value calculation step includes a step of calculating, by a processor in communication with the transmit-receive data acquisition part, the physical property value calculation part, and the physical property value-display value conversion part, an estimate of a deviation quantity of an assumed physical property value, the estimate derived from the transmit-receive data acquired in the transmit-receive data, a distribution of the assumed physical property value, and an equation formed in consideration of and a wave propagation, the deviation quantity of the physical property value is calculated, by the processor, as an evaluation quantity J which is any one of, a linear sum of any two or more of, a product of any two or more of, or an exponential function including as an exponent any one of: an equation residual quantity that is a residual being a difference between an operator term and an external force term of an equation of motion, a non-uniformity detection equation residual quantity that is a residual of an equation of detecting non-uniformity of the physical property value according to the matching degree between solutions of the equation of motion under two types of boundary conditions, and a conditional equation residual quantity that is a residual of a constraint condition; and a sum of the assumed physical property value and the deviation quantity is outputted, by the processor, as a calculated physical property value.