IMAGE INFORMATION PROCESSING DEVICE, IMAGE INFORMATION PROCESSING SYSTEM, IMAGE INFORMATION PROCESSING METHOD, AND RECORDING MEDIUM STORING IMAGE INFORMATION PROCESSING PROGRAM

20170344845 · 2017-11-30

Assignee

Inventors

Cpc classification

International classification

Abstract

In an image representing an observed ground surface area, the device calculates an optical-path-radiance with a high degree of accuracy. An image information processing device 3 includes: a storage unit 31 that associates and stores observation images representing the results of observing electromagnetic waves of a plurality of different wavelength bands reflected from a ground surface, information representing the wavelength bands, and information representing the observation environment; a first intermediate-optical-path-radiance calculation unit that, for each of the wavelength band, makes the radiance of a dark pixel meeting a radiance standard from among pixels composing an associated observation image an intermediate-optical-path-radiance; an irradiance calculation unit that calculates the irradiance from sunlight for each of the wavelength band on the basis of the information representing the observation environment; and a final-optical-path-radiance calculation unit that calculates a final-optical-path-radiance for each of the wavelength band on the basis of the irradiance and intermediate-optical-path-radiance.

Claims

1. An image information processing device comprising: a storage unit configured to store information, in which observation images representing results of observing electromagnetic waves of a plurality of different wavelength bands reflected from a ground surface, information representing the wavelength bands, and information representing an observation environment, which the electromagnetic waves are observed, are associated; a first intermediate-optical-path-radiance calculation unit configured to calculate, for each of the wavelength band, radiance of a dark pixel meeting a radiance standard from among pixels composing the observation image, which is associated with the wavelength band, as a first intermediate-optical-path-radiance; an irradiance calculation unit configured to calculate, for each of the wavelength band, irradiance from sunlight on the basis of the information representing the observation environment; and a final-optical-path-radiance calculation unit configured to calculate, for each of the wavelength band, a final-optical-path-radiance on the basis of the irradiance and the first intermediate-optical-path-radiance.

2. The image information processing device according to claim 1, wherein the final-optical-path-radiance calculation unit includes: a cost minimization unit configured to sequentially generate atmosphere information on a state of the atmosphere by updating a value(s) of a variable(s) representing the atmosphere information; a second intermediate-optical-path-radiance calculation unit configured to calculate a second intermediate-optical-path-radiance on the basis of the atmosphere information and the irradiance at every time when the cost minimization unit generates the atmosphere information; and a cost calculation unit configured to calculate a value of a cost indicating a degree of difference between the second intermediate-optical-path-radiance and actual-optical-path-radiance which is a result of observing the electromagnetic wave, on the basis of the first and the second intermediate-optical-path-radiance, wherein the cost minimization unit calculates the second intermediate-optical-path-radiance, when the value of the cost meets a predetermined condition, as the final-optical-path-radiance.

3. The image information processing device according to claim 1, wherein the final-optical-path-radiance calculation unit, comprising: a cost minimization unit configured to sequentially generate atmosphere information on a state of the atmosphere by updating the value(s) of the variable(s) representing the atmosphere information; a second intermediate-optical-path-radiance calculation unit configured to calculate a second intermediate-optical-path-radiance on the basis of the atmosphere information and the irradiance at every time when the cost minimization unit generates the atmosphere information; a wavelength selection unit configured to select one or more wavelength bands out of a plurality of wavelength bands on the basis of the first and the second intermediate-optical-path-radiance; and a wavelength cost calculation unit configured to calculate a value of a cost indicating a degree of difference between the second intermediate-optical-path-radiance, and actual-optical-path-radiance which is a result of observing the electromagnetic wave, on the basis of information indicating the first and the second intermediate-optical-path-radiance and the wavelength band which is selected by the wavelength selection unit, and wherein the cost minimization unit calculates the second intermediate-optical-path-radiance, when the value of the cost meets a predetermined condition, as the final-optical-path-radiance.

4. The image information processing device according to claim 3, wherein the wavelength selection unit selects the wavelength band in which the first intermediate-optical-path-radiance is higher than the second intermediate-optical-path-radiance, and wherein the wavelength cost calculation unit calculates the value of the cost giving more weight to the wavelength band which is not selected by the wavelength selection unit than the wavelength band which is selected by the wavelength selection unit.

5. The image information processing device according to claim 2, wherein the cost minimization unit calculates the second intermediate-optical-path-radiance, which is calculated when the value of cost is minimized, as the final-optical-path-radiance.

6. The image information processing device according to claim 2, wherein the cost minimization unit outputs the atmosphere information, which is generated when the value of the cost meets the predetermined condition, to an outside.

7. The image information processing device according to claim 1, further comprising: an optical-path-radiance correction unit configured to generate, for each the wavelength band, an optical-path-radiance correction image by removing a component which is indicated by the final-optical-path-radiance from the observation image.

8. An image information processing system, comprising: an image information processing device including; a storage unit configured to store information, in which observation images representing results of observing electromagnetic waves of a plurality of different wavelength bands reflected from a ground surface, information representing the wavelength bands, and information representing an observation environment, which the electromagnetic waves are observed, are associated; a first intermediate-optical-path-radiance calculation unit configured to calculate, for each of the wavelength band, radiance of a dark pixel meeting a radiance standard from among pixels composing the observation image, which is associated with the wavelength band, as a first intermediate-optical-path-radiance; an irradiance calculation unit configured to calculate, for each of the wavelength band, irradiance from sunlight on the basis of the information representing the observation environment; and a final-optical-path-radiance calculation unit configured to calculate, for each of the wavelength band, a final-optical-path-radiance on the basis of the irradiance and the first intermediate-optical-path-radiance; and an image supply device that is connected with the image information processing device so as to be able to communicate, and associates the observation image with the wavelength band in which the observation image is observed, and the information representing the observation environment, and inputs the observation image into the image information processing device.

9. An image information processing method comprising: in a case that observation images representing results of observing electromagnetic waves of a plurality of different wavelength bands reflected from a ground surface, information representing the wavelength bands, and information representing an observation environment, which the electromagnetic waves are observed, are associated and stored in a storage unit, by an information processing device: calculating, for each of the wavelength band, radiance of a dark pixel meeting a radiance standard from among pixels composing the observation image which is associated with the wavelength band, as a first intermediate-optical-path-radiance; calculating, for each of the wavelength band, irradiance from sunlight on the basis of the information representing the observation environment; and calculating, for each of the wavelength band, a final-optical-path-radiance on the basis of the irradiance and the first intermediate-optical-path-radiance.

10. A non-transitory computer readable recording medium storing an image information processing program for making a computer, which is capable of accessing a storage unit in which observation images representing results of observing electromagnetic waves of a plurality of different wavelength bands reflected from a ground surface, information representing the wavelength bands, and information representing an observation environment, which the electromagnetic waves are observed, are associated and stored, carry out following processing: an intermediate-optical-path-radiance calculation processing of calculating, for each of the wavelength band, radiance of a dark pixel meeting a radiance standard from among pixels composing the observation image which is associated with the wavelength band, as a first intermediate-optical-path-radiance; an irradiance calculation processing of calculating, for each of the wavelength band, irradiance from sunlight on the basis of the information representing the observation environment; and a final-optical-path-radiance calculation processing of calculating, for each of the wavelength band, a final-optical-path-radiance on the basis of the irradiance and the first intermediate-optical-path-radiance.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0020] FIG. 1 is a block diagram showing a configuration of an image information processing system according to a first example embodiment of the invention of the present application.

[0021] FIG. 2 is a flowchart showing an operation of the image information processing system according to the first example embodiment of the invention of the present application.

[0022] FIG. 3 is a block diagram showing a configuration of an image information processing system according to a second example embodiment of the invention of the present application.

[0023] FIG. 4 is a block diagram showing a configuration of an image information processing system according to a third example embodiment of the invention of the present application.

[0024] FIG. 5 is a diagram showing a relation between observation light and an environment fluctuation component in an image which is photographed at a height.

[0025] FIG. 6 is a block diagram showing a configuration of an image information processing device which uses an art disclosed by PTL 1.

[0026] FIG. 7 is a block diagram showing a configuration of an information processing device which can carry out the image information processing device according to each example embodiment of the invention of the present application.

DESCRIPTION OF EMBODIMENTS

[0027] Hereinafter, an example embodiment of the invention of the present application will be explained in detail with reference to drawings.

First Example Embodiment

[0028] FIG. 1 is a block diagram showing an image information processing system 100 of a first example embodiment of the present invention. The image information processing system 100 according to the present example embodiment includes an image information processing device 1, an image supply device 5 and an output device 6 as a main block. The image information processing device 1 is connected with the image supply device 5 and the output device 6 so as to be able to communicate each other. Here, a direction of arrow, which is described in the block diagram such as FIG. 1 and the like according to the present application is merely exemplified, and the direction does not limit a direction of signal between the blocks.

[0029] The image supply device 5 is, for example, a photographing device which, for each of the plural wavelength bands different each other, observes intensity of an electromagnetic wave (observation light) reflected from an observation target, and outputs the observation result as an observation image. Or, the image supply device 5 may be a storage device such as a hard disk or the like, or an information processing device such as a server device or the like which stores the observation image corresponding to the observation result outputted by the photographing device. In the following explanation, N (N is an integer larger than 1) is number of wavelength bands which are used when the image supply device 5 (photographing device) observes the observation object.

[0030] The image supply device 5 supplies the image information processing device 1 with the observation images whose number is N and each of which is associated with the corresponding wavelength band. The Image supply device 5 may supply the image information processing device 1 with N observation images, and the maximum value and the minimum value of the wavelength band which is associated with each observation image. Moreover, the image supply device 5 supplies the image information processing device 1 with information which expresses an observation environment existing at a time of observation. The information which expresses the observation environment includes, for example, altitude of the sun at a time of observation, an observation date, observation time, latitude and longitude of an observation target area, topography of the observation target area, quantity of steam and quantity of aerosol in the observation target area, or the like.

[0031] The Image supply device 5 observes distribution of intensity of observation light which is emitted by the observation target to reach the image supply device 5 through a medium such as the atmosphere, dust or the like. In each example embodiment of the invention of the present application, the observation light includes visible light, infrared rays and ultraviolet rays. However, the image supply device 5 may observe distribution of intensity of an electromagnetic wave other than the visible light, the infrared rays and the ultraviolet rays.

[0032] As will be mentioned later, the observation light which reaches the image supply device 5 includes a component which reaches through the medium in addition to a component which reaches by reflection at the observation target. The observation target is, for example, the ground surface. The image supply device 5 is installed, for example, on an airplane or a man-made satellite. Then, by making a plurality of wavelength bands, which are different each other, the observation target, the image supply device 5 observes reflection light, which is reflected from the ground surface, in the air as a part of the observation light. The image supply device 5 is, for example, a camera which photographs through a band-pass filter. The image supply device 5 may carry out observation by using N band-pass filters each of which selectively transmits any one of N different wavelength bands. Here, in place of observing the ground surface in the air, the image supply device 5 may observe the distant ground surface from the ground surface or a near place above the ground surface. The wavelength bandwidth may not be uniform. The image supply device 5 outputs the observation result as N observation images.

[0033] Each of N observation images is an image which is generated by observing any one of N different wavelength bands and which expresses distribution of radiance related to the observation target. In each observation image, the radiance value of each pixel expresses intensity of the observation light which reaches from a direction according to the pixel. The above-mentioned N observation images are images which are obtained by observing the same observation target. In this case, N images are, for example, results which are obtained by the image supplying device 5's observing the observation target by using N wavelength bands. Moreover, in the case that a value of the optical path radiance is the same or almost the same in all observation areas of each image, the above-mentioned N images may be images which are obtained by observing the different observation targets. In this case, a manager who manages, for example, the image supply device 5 or the image information processing device 1 may select N images.

[0034] The output device 6 is, for example, a display unit. Moreover, the output device 6 may be, for example, an object identification device or the like which extracts information on a material or a state of the photographing object from the observation image.

[0035] The image information processing device 1 includes a storage unit 11, a first intermediate-optical-path-radiance calculation unit 12, an irradiance calculation unit 13, a final-optical-path-radiance calculation unit 14, an optical-path-radiance correction unit 15, and an output unit 19. The storage unit 11 is a storage device such as a magnetic disk, an electronic memory or the like. There is a case that the first intermediate-optical-path-radiance calculation unit 12, the irradiance calculation unit 13, the final-optical-path-radiance calculation unit 14, the optical-path-radiance correction unit 15, and the output unit 19 are realized by using electronic circuits, and there is another case that they are realized by using a computer program and a processor which works according to the computer program.

[0036] The storage unit 11 reads N observation images, which are associated with the observed wavelength bands respectively, from the image supply device 5. The storage unit 11 associates the read N observation images with the wavelength bands by which the observation images are observed, and stores the read N observation images. The wavelength band which is associated with the observation image may be, for example, a combination of the minimum value and the maximum value of the wavelength which indicate the wavelength band. The wavelength band which is associated with the observation image may be an identifier assigned to the combination of the minimum value and the maximum value of the wavelength which indicate the wavelength band. The storage unit 11 associates and stores the information representing the observation environment with the observation image. The information which expresses the observation environment includes, for example, altitude of the sun at a time of observation, an observation date, observation time, latitude and longitude of an observation target area, topography of the observation target area, quantity of steam, quantity of aerosol or the like at a time of observation.

[0037] The first intermediate-optical-path-radiance calculation unit 12 extracts a radiance value of the pixel, which meets a radiance standard, out of all pixels composing the observation images which are stored by the storage unit 11 and which are associated with N wavelength bands respectively. Then, the first intermediate-optical-path-radiance calculation unit 12 makes each of the extracted N radiance values the first intermediate-optical-path-radiance of each wavelength band. For example, the first intermediate-optical-path-radiance calculation unit 12 extracts the minimum radiance value as the radiance value which meets the radiance standard. Moreover, for example, in the case that the radiance values of all pixels are lined up in an order of smallness, the first intermediate-optical-path-radiance calculation unit 12 may extract the radiance value, whose rank in the line is coincident with number obtained by multiplying number of all pixels by a predetermined ratio, as the radiance value which meets the radiance standard.

[0038] In the following explanation, a pixel, which is corresponding to the radiance value extracted by the first intermediate-optical-path-radiance calculation unit 12, is referred to as ‘dark pixel’. Moreover, an assumption that the surface reflectance of a dark area of the ground surface, which is measured as the radiance value of the dark pixel, can be considered to be 0 is referred to as ‘dark pixel assumption’. Moreover, i'th (i is an integer of any one of 1 to N) wavelength band out of N wavelength bands, which are associated with N observation images respectively, is denoted as γ.sub.i. Moreover, the extracted radiance value of the dark pixel in the wavelength band γ.sub.i out of N wavelength bands associated with the observation images is referred to as the first intermediate optical path radiance β.sup.D(γ.sub.i).

[0039] The irradiance calculation unit 13 calculates the irradiance from the sunlight in the wavelength band γ.sub.i out of N wavelength bands on the basis of the information which expresses the observation environment and which is stored by the storage unit 11. In the following explanation, the irradiance from the sunlight, which is calculated by the irradiance calculation unit 13, in the wavelength band γ.sub.i out of N wavelength bands is denoted as I(γ.sub.i). For example, the irradiance calculation unit 13 calculates the irradiance I(γ.sub.i) by using a physical model shown by NPL 1 which expresses the irradiance. By using the physical model, the irradiance calculation unit 13 simulates direct sunlight and light, which is scattered by the atmosphere, by using parameters representing the zenith angle of sun, which is calculated on the basis of a location and time, and a state of the atmosphere. By virtue of the above, the irradiance calculation unit 13 can calculate, for example, the spectral irradiance of the sunlight, which is observed at a sunny time, for each wavelength.

[0040] In the case of using the physical model, firstly, the irradiance calculation unit 13 sets a predetermined value of the parameter which has a small amount of influence on the simulation result. Then, the irradiance calculation unit 13 calculates the spectral irradiance, which is irradiance per a unit wavelength of the sunlight at a sunny time, by using the predetermined value, and the zenith angle of sun and a value representing the state of atmosphere out of the information representing the observation environment. Then, the irradiance calculation unit 13 calculates a total amount of the spectral irradiance of the sunlight at the sunny time from the minimum value of the wavelength up to the maximum value of the wavelength out of the wavelength band γ.sub.i. By carrying out the above, the irradiance calculation unit 13 calculates the irradiance I(γ.sub.i) for each wavelength band γ.sub.i out of N wavelength bands.

[0041] Moreover, the irradiance calculation unit 13 may stores the standard irradiance from the sunlight for each wavelength band γ.sub.i out of N wavelength bands, and may use the standard irradiance in place of the irradiance I(γ.sub.i). As shown, for example, by NPL2, the irradiance calculation unit 13 calculates a total amount of the spectral irradiance of the sunlight which is observed at a position where the sunlight enters into the inner space from the outer space, from the minimum value of the wavelength up to the maximum value of the wavelength out of the wavelength band γ.sub.i. By carrying out the above, the irradiance calculation unit 13 may calculate and store the irradiance for each wavelength band γ.sub.i out of N wavelength bands.

[0042] The final-optical-path-radiance calculation unit 14 includes a cost calculation unit 141, a second intermediate-optical-path-radiance calculation unit 142, and a cost minimization unit 143.

[0043] The second intermediate-optical-path-radiance calculation unit 142 calculates a second intermediate-optical-path-radiance by using the irradiance I(γ.sub.i) which is calculated by the irradiance calculation unit 13, and atmosphere information which is outputted by the cost minimization unit 143. In the following explanation, a parameter out of M parameters representing the atmosphere information outputted by the cost minimization unit 143 is denoted as X.sub.j (j is an integer of any one of 1 to M and is a subscript for identifying the parameter.). The atmosphere information X.sub.j expresses information on a particle, which is included by the atmosphere, such as a molecule, aerosol, sand, soot, or the like. The atmosphere information X.sub.j is, for example, an angstrom index, atmosphere turbidity, atmosphere transmittance, an optical thickness of atmosphere molecule, an optical thickness of aerosol, visibility or a value which is calculated by using these. Moreover, the second intermediate-optical-path-radiance, which is calculated by the second intermediate-optical-path-radiance calculation unit 142, in the wavelength band γ.sub.i out of N wavelength bands is denoted as β.sup.M(γ.sub.i). Each of Formulas expressed as Formula. 2 to Formula. 4 is an example of Formula representing a calculation method for calculating the second intermediate-optical-path-radiance β.sup.M(γ.sub.i) in the wavelength band γ.sub.i. Each of Formulas expressed as Formula. 2 to Formula. 4 is Formula based on a rule of thumb. The second intermediate-optical-path-radiance calculation unit 142 calculates the second intermediate-optical-path-radiance β.sup.M(γ.sub.i), for example, by using any one of Formulas of Formula. 2 to Formula. 4.


β.sup.M(γ.sub.i)=X.sub.1I(γ.sub.i)γ.sub.i.sup.−X.sup.2   [Formula. 2]

[0044] A unit (dimension) of β.sup.M(γ.sub.i) in Formula. 2 is W (watt)/steradian×(m (meter)).sup.−3. A unit of I(γ.sub.i) is W/m.sup.−3. A unit of γ.sub.i is m. X.sub.1 and X.sub.2 in Formula. 2 are parameters which express the atmosphere information having the following means respectively and which are determined based on experiences. Moreover, X.sub.1 and X.sub.2 have dimensions so that a part other than the term I(γ.sub.i) in a right side of Formula. 2 may have no dimension.

[0045] X.sub.1: a parameter which expresses a degree of atmosphere turbidity such as a coefficient of atmosphere turbidity or the like

[0046] X.sub.2: a parameter, which expresses an inclination (a large value when there are many small size particles) in a size of particle included by the atmosphere, such as the angstrom index or the like.


β.sup.M(γ.sub.i)=X.sub.3I(γ.sub.i)(1−e.sup.−X.sup.1.sup.γsi.sup.−X.sup.2)   [Formula. 3]

[0047] Units of β.sup.M(γ.sub.i), I(γ.sub.i) and y.sub.i in Formula. 3 are the same as the units in Formula. 2. X.sub.1 to X.sub.3 in Formula. 3 are parameters which express the atmosphere information having the following means respectively and which are determined based on experiences. Moreover, X.sub.1 to X.sub.3 have dimensions so that a part other than the term I(γ.sub.i) in a right side of Formula. 3 may have no dimension.

[0048] X.sub.1 and X.sub.2: the same as the parameters in Formula. 1 respectively

[0049] X.sub.3: a parameter (depending on a relation between a position of the sun and a position of a sensor) which expresses a ratio of scattered light


β.sup.M(γ.sub.i)=X.sub.3I(γ.sub.i)(1−e.sup.−(X.sup.1.sup.γsi.sup.−X.sup.2.sup.+X.sup.4.sup.τ(λ)))   [Formula. 4]

[0050] Units of β.sup.M(γ.sub.i), I(γ.sub.i) and γ.sub.i in Formula. 4 are the same as the units in Formula. 2. and Formula. 3. X.sub.1 to X.sub.4 in Formula. 4 are parameters which express the atmosphere information having the following means respectively and which are determined based on experiences. Moreover, X.sub.1 to X.sub.4 have dimensions so that a part other than the term I(γ.sub.i) in a right side of Formula. 4 may have no dimension.

[0051] X.sub.1 to X.sub.3: the same as the parameters in Formula.3 respectively

[0052] X.sub.4: a parameter which expresses a degree of influence of the atmosphere transmittance on the optical-path-radiance

[0053] The cost calculation unit 141 calculates a cost C by using the first intermediate-optical-path-radiance β.sup.D(γ.sub.i) which is outputted by the first intermediate-optical-path-radiance calculation unit 12, and the second intermediate-optical-path-radiance β.sup.M(γ.sub.i) which is calculated by the second intermediate-optical-path-radiance calculation unit 142. Here, the cost C is a value which indicates a degree of difference between the second intermediate-optical-path-radiance, and actual-optical-path-radiance at a time of observation. Formula. 5 is an example of Formula which expresses a cost C calculation method. The cost calculation unit 141 calculates the cost C, for example, by using a mathematical Formula shown in Formula. 5.


C=Σ.sub.i=1.sup.N(β.sup.D(γ.sub.i)−β.sup.M(γ.sub.i)).sup.2   [Formula. 5]

[0054] The cost minimization unit 143 sequentially updates values of M parameters to generate the atmosphere information X.sub.j. When a value of the cost C calculated by the cost calculation unit 141 meets a predetermined condition, the cost minimization unit 143 outputs the final-optical-path-radiance and estimated atmosphere information. The final-optical-path-radiance outputted by the cost minimization unit 143 is a value, which makes the value of the cost C that is calculated by the cost calculation unit 141 minimum, out of the values of the second intermediate-optical-path-radiance which are calculated by the second intermediate-optical-path-radiance calculation unit 142. Moreover, the estimated atmosphere information is atmosphere information which is obtained when the cost calculated by the cost calculation unit 141 is minimized. The final-optical-path-radiance is a calculation value of the final-optical-path-radiance which is calculated at a time of observation. Moreover, the estimated atmosphere information is a value of the atmosphere information which is finally estimated at a time of observation.

[0055] In the following explanation, the final-optical-path-radiance, which is calculated by the cost minimization unit 143, in the wavelength band γ.sub.i out of N wavelength bands is denoted as β.sup.E(γ.sub.i). Moreover, the estimated atmosphere information is denoted as X.sup.E.sub.J. Moreover, as intermediate parameters which are used in a process for calculating the final-optical-path-radiance β.sup.E(γ.sub.i), the cost minimization unit 143 stores the minimum cost C.sup.m which is obtained at the present point of time, the atmosphere information X.sup.m.sub.j which is obtained when the value of the cost is the minimum cost C.sup.m, and the atmosphere information X.sup.b.sub.j which is obtained at a previous time when sequentially updating M values of parameters. As initial values of the minimum cost C.sup.m, the atmosphere information X.sup.m.sub.j at the time of the minimum cost, and the atmosphere information X.sub.j, the cost minimization unit 143 stores predetermined values.

[0056] The cost minimization unit 143 outputs the atmosphere information X.sub.j by sequentially updating the values of M parameters with, for example, the following method. Firstly, the cost minimization unit 143 outputs the stored initial value of the atmosphere information X.sub.j. At every time when inputting the cost C, the cost minimization unit 143 carries out the following process, and outputs the atmosphere information X.sub.j.

[0057] The cost minimization unit 143 compares the cost C and the minimum cost C.sup.m. In the case that the cost C is smaller than the minimum cost C.sup.m, the cost minimization unit 143 substitutes the cost C for the minimum cost C.sup.m, and substitutes the atmosphere information X.sub.j for the atmosphere information X.sup.m.sub.j which is obtained at the time of the minimum cost. The cost minimization unit 143 stores the updated minimum cost C.sup.m and the updated atmosphere information X.sup.m.sub.j which is obtained at the time of the minimum cost. The cost minimization unit 143 substitutes a value which is obtained by adding a predetermined value to the parameter of the atmosphere information X.sup.b.sub.j, which is obtained at a previous time, for the parameter of the atmosphere information X.sub.j, and outputs the substitution value. The cost minimization unit 143 substitutes the updated value of the atmosphere information X.sub.j for the atmosphere information X.sup.b.sub.j which is obtained at a previous time, and stores the substitution value.

[0058] Moreover, the cost minimization unit 143 outputs the final optical path radiance β.sup.E(γ.sub.i) and the estimated atmosphere information X.sup.E.sub.J when the value of the cost C meets the predetermined condition with the following method. In the case that the value of the cost C is smaller than a predetermined value, the cost minimization unit 143 uses a mathematical Formula which is the same as the mathematical formula used when the second intermediate-optical-path-radiance calculation unit 142 calculates the second intermediate-optical-path-radiance β.sup.M(γ.sub.i). For example, in the case that the second intermediate-optical-path-radiance calculation unit 142 calculates the second intermediate-optical-path-radiance β.sup.M(γ.sub.i) by using Formula. 2, the cost minimization unit 143 works as follows.

[0059] The cost minimization unit 143 substitutes the value of atmosphere information X.sup.b.sub.j, which is obtained at a previous time, for the atmosphere information X.sub.j. The cost minimization unit 143 calculates the second intermediate-optical-path-radiance β.sup.M(γ.sub.i) by using Formula. 2. The cost minimization unit 143 substitutes the calculated value of the second intermediate-optical-path-radiance β.sup.M(γ.sub.i) for the final-optical-path-radiance (γ.sub.i). The cost minimization unit 143 substitutes the value of the atmosphere information X.sup.m.sub.j, which is obtained at the time of the minimum cost, for the estimated atmosphere information X.sup.E.sub.J. Then, the cost minimization unit 143 outputs the final optical path radiance β.sup.E(γ.sub.i) and the estimated atmosphere information X.sup.E.sub.J. Similarly, in the case that the second intermediate-optical-path-radiance calculation unit 142 calculates the second intermediate-optical-path-radiance β.sup.M(γ.sub.i) by using Formula. 3 or Formula. 4, also the cost minimization unit 143 calculates the second intermediate-optical-path-radiance β.sup.M(γ.sub.i) by using Formula. 3 or Formula. 4.

[0060] Moreover, in the case of calculating the second intermediate-optical-path-radiance β.sup.M(γ.sub.i), for example, by using Formula. 2, the cost minimization unit 143 works as follows. That is, the cost minimization unit 143 sequentially updates the value of the atmosphere information X.sub.j, and outputs the updated value of the atmosphere information X.sub.j as shown in the following. Firstly, the cost minimization unit 143 outputs the initial value of the atmosphere information X.sub.j. Hereinafter, at every time when inputting the cost C, the cost minimization unit 143 calculates the value of the atmosphere information X.sub.j by using Formula. 6 (matrix formula), and outputs the value of the atmosphere information X.sub.j.

[00001] [ X 1 X 2 ] = [ X 1 b X 2 b ] - α [ - 2 .Math. .Math. i = 1 N .Math. I ( γ i ) .Math. γ i - X 2 b 2 .Math. X 1 b .Math. .Math. i = 1 N .Math. I ( γ i ) .Math. ( log e .Math. .Math. γ i ) .Math. γ i - X 2 b ] [ Formula . .Math. .Math. 6 ]

[0061] The cost minimization unit 143 substitutes the value of the atmosphere information X.sub.j for the atmosphere information X.sup.b.sub.j which is obtained at a previous time, and substitutes the value of cost C for a cost C.sup.b which is obtained at a previous time. The cost minimization unit 143 stores the atmosphere information X.sup.b.sub.j and the cost C.sup.b whose values are changed and which are obtained at the previous time.

[0062] Moreover, in the case that the cost minimization unit 143 calculates, for example, the second intermediate-optical-path-radiance β.sup.M(γ.sub.i) by using Formula. 2, by working as follows, the cost minimization unit 143 judges that the cost C meets the predetermined condition, and afterward outputs the final-optical-path-radiance β.sup.E(γ.sub.i) and the estimated atmosphere information X.sup.E.sub.J. The cost minimization unit 143 compares a difference between the cost C and the cost C.sup.b, which is obtained at the previous time, with a predetermined value. When the difference between the cost C and the cost C.sup.b, which is obtained at the previous time, is smaller than the predetermined value, the cost minimization unit 143 works as follows.

[0063] The cost minimization unit 143 substitutes the value of the atmosphere information X.sup.b.sub.j, which is obtained at a previous time, for the atmosphere information X.sub.j. The cost minimization unit 143 calculates the second intermediate-optical-path-radiance β.sup.M(γ.sub.i) by using Formula. 2. The cost minimization unit 143 substitutes the calculated value of the second intermediate-optical-path-radiance β.sup.M(γ.sub.i) for the final-optical-path-radiance β.sup.E(γ.sub.i). The cost minimization unit 143 substitutes the value of the atmosphere information X.sup.m.sub.j, which is obtained at the time of the minimum cost, for the estimated atmosphere information X.sup.E.sub.J. Then, the cost minimization unit 143 outputs the final-optical-path-radiance β.sup.E(γ.sub.i) and the estimated atmosphere information X.sup.E.sub.J.

[0064] The optical-path-radiance correction unit 15 corrects a component of the optical-path-radiance, which is included by the observation image, by using the observation image which the storage unit 11 stores, and the final-optical-path-radiance β.sup.E(γ.sub.i), which the cost minimization unit 143 outputs, to generate an optical path radiance correction image. For example, the optical-path-radiance correction unit 15 works as follows.

[0065] The optical-path-radiance correction unit 15 reads the observation image, which is associated with the wavelength band γ.sub.i out of N wavelength bands, from the storage unit 11. The optical-path-radiance correction unit 15 carries out correction of subtracting the final-optical-path-radiance β.sup.E(γ.sub.i), which the cost minimization unit 143 outputs, from the radiance value L(γ.sub.i) for all pixels composing the read observation image. That is, the corrected observation image is an image obtained by removing the component of the optical-path-radiance, which is calculated by the final-optical-path-radiance calculation unit 14, from the observation image which the storage unit 11 stores. The image is defined as the optical-path-radiance correction image. The optical-path-radiance correction unit 15 carries out the above-mentioned correction process to all of the observation images which are associated with N wavelength bands including the wavelength band γ.sub.i.

[0066] The output unit 19 outputs the optical-path-radiance correction image which is generated by the optical-path-radiance correction unit 15, and the estimated atmosphere information X.sup.E.sub.J, which is outputted by the cost minimization unit 143, to the output device 6.

[0067] Next, with reference to a flowchart shown in FIG. 2, an operation (process) of the image information processing system 100 according to the present example embodiment will be explained in detail.

[0068] The storage unit 11 reads N observation images associated with N wavelength bands by which observation is carried out, and the information which expresses the observation environment existing at a time of observation from the image supply device 5. The storage unit 11 associates the read observation image with the wavelength band by which the observation image is observed, and stores the observation image and the wavelength band, which are associated, together with information which expresses an observation environment (Step S101). The first intermediate-optical-path-radiance calculation unit 12 extracts the radiance value of the dark pixel, which meets the radiance standard, out of all pixels of each of the observation images which are associated with the wavelength bands respectively, and outputs the extracted radiance values related to N wavelength bands as the first intermediate-optical-path-radiance (Step S102).

[0069] The irradiance calculation unit 13 calculates the irradiance from the sunlight on the basis of the information which expresses the observation environment (Step S103). At every time when the cost calculation unit 141 calculates a new cost, the cost minimization unit 143 sequentially updates the value of the atmosphere information, and outputs the updated value of the atmosphere information (Step S104). The second intermediate-optical-path-radiance calculation unit 142 calculates the second intermediate-optical-path-radiance on the basis of the irradiance from the sunlight which the irradiance calculation unit 13 calculates, and the atmosphere information which the cost minimization unit 143 outputs (step S105).

[0070] The cost calculation unit 141 calculates the cost on the basis of the first intermediate-optical-path-radiance which the first intermediate-optical-path-radiance calculation unit 12 outputs, and the second intermediate-optical-path-radiance which the second intermediate-optical-path-radiance calculation unit 142 calculates (Step S106).

[0071] In the case that the cost meets the predetermined condition (Yes in Step S107), the cost minimization unit 143 calculates the final-optical-path-radiance and the estimated atmosphere information (Step S108). In the case that the cost does not meet the predetermined condition (No in Step S107), the process returns to step S104.

[0072] The optical-path-radiance correction unit 15 removes the component of the optical-path-radiance, which is included by the observation image, by using the observation image and the final-optical-path-radiance to generate the optical-path-radiance correction image (Step S109). The output unit 19 outputs the optical-path-radiance correction image, which the optical-path-radiance correction unit 15 generates, to the output device 6 (step S110), and then the whole process has been finished.

[0073] The image information processing system 100 according to the present example embodiment can calculate the optical-path-radiance with a high degree of accuracy even in the wavelength band in which the surface reflectance of the dark area in the image representing the observed ground surface area cannot be considered to be 0. The reason is that the final-optical-path-radiance calculation unit 14 calculates the cost, which indicates the degree of difference between the second intermediate-optical-path-radiance and the actual-optical-path-radiance, by using the first intermediate-optical-path-radiance which is calculated by the first intermediate-optical-path-radiance calculation unit 12, and the second intermediate-optical-path-radiance which is calculated by the second intermediate-optical-path-radiance calculation unit 142, and outputs the second intermediate-optical-path-radiance, which is calculated at the time when the cost becomes minimum, as the final-optical-path-radiance.

[0074] An effect which is obtained by the image information processing system 100 according to the present example embodiment will be explained in the following in detail. In the case of the general image information processing device 4 shown in FIG. 6, it is impossible to calculate the optical-path-radiance with a high degree of accuracy in the wavelength band in which the surface reflectance of the area, that is corresponding to the pixel having small radiance, out of the ground surface area observed as the image cannot be considered to be 0. That is, by the dark pixel assumption, it is difficult to calculate the actual-optical-path-radiance at a time of observation with a high degree of accuracy.

[0075] In contrast, in the case of the image information processing device 1 according to the present example embodiment, the first intermediate-optical-path-radiance calculation unit 12 extracts the radiance value of the dark pixel, which meets the radiance standard, from the observation image which is associated with the wavelength band, and outputs the extracted radiance value as the first intermediate-optical-path-radiance. The second intermediate-optical-path-radiance calculation unit 142 repeatedly calculates the second intermediate-optical-path-radiance with updating the irradiance from the sunlight which is calculated by the irradiance calculation unit 13, and the atmosphere information which is calculated by the cost minimization unit 143. Then, as the final-optical-path-radiance, the final-optical-path-radiance calculation unit 14 outputs the second intermediate-optical-path-radiance obtained when the cost, which is calculated on the basis of the first intermediate-optical-path-radiance and the second intermediate-optical-path-radiance, is minimized. Accordingly, the image information processing device 1 according to the present example embodiment can calculate the optical-path-radiance with a high degree of accuracy even in the wavelength band in which the surface reflectance of the dark area in the image representing the observed ground surface area cannot be considered to be 0.

[0076] Moreover, the final-optical-path-radiance calculation unit 14 of the image information processing device 1 according to the present application example embodiment outputs the estimated atmosphere information together with the final-optical-path-radiance through the output unit 19. The estimated atmosphere information is information which indicates the state of the atmosphere existing at a time of observation. Therefore, the image information processing device 1 according to the present example embodiment can obtain the atmosphere information which can be used effectively as weather information.

Second Example Embodiment

[0077] FIG. 3 is a block diagram which conceptually shows a configuration of an image information processing system 200 according to a second example embodiment. In the following description, explanation is omitted by assigning a code, which is the same as the code used in the first example embodiment, to a component which works identically to the component of the first example embodiment.

[0078] The image information processing system 200 according to the present example embodiment includes an image information processing device 2, the image supply device 5 and the output device 6 as a main component.

[0079] The image information processing device 2 includes the storage unit 11, the first intermediate-optical-path-radiance calculation unit 12, the irradiance calculation unit 13, a final-optical-path-radiance calculation unit 24, the optical-path-radiance correction unit 15, and the output unit 19. The final-optical-path-radiance calculation unit 24 may be realized by an electronic circuit in some cases, and may be realized by a computer program and a processor, which works according to the computer program, in the other cases.

[0080] The final-optical-path-radiance calculation unit 24 includes a wavelength selection unit 2411, a wavelength cost calculation unit 2412, the second intermediate-optical-path-radiance calculation unit 142 and the cost minimization unit 143.

[0081] The final-optical-path-radiance calculation unit 24 calculates the final-optical-path-radiance β.sup.E(γ.sub.i) (i is an integer of any one of 1 to N) and the estimated atmosphere information X.sup.E.sub.J by using the first intermediate-optical-path-radiance β.sup.D(γ.sub.i) which the first intermediate-optical-path-radiance calculation unit 12 outputs, and the irradiance I(γ.sub.i) from the sunlight which the irradiance calculation unit 13 calculates. The final-optical-path-radiance calculation unit 24 selects the wavelength band γi out of N wavelength bands on the basis of a predetermined standard, and calculates a value in the selected wavelength band, and a value in the not-selected wavelength band with different methods respectively, and afterward calculates addition of two values. Then, the final-optical-path-radiance calculation unit 24 makes the result of addition a cost. Moreover, the final-optical-path-radiance calculation unit 24 outputs the second intermediate-optical-path-radiance, which is obtained when the cost meets a predetermined condition, as the final-optical-path-radiance. The final-optical-path-radiance calculation unit 24 outputs the atmosphere information, which is obtained when the cost is minimized, as the estimated atmosphere information X.sup.E.sub.J.

[0082] The wavelength selection unit 2411 selects 1 or more wavelength bands from N observation wavelength bands by using the first intermediate-optical-path-radiance which the first intermediate-optical-path-radiance calculation unit 12 outputs, and the second intermediate-optical-path-radiance which the second intermediate-optical-path-radiance calculation unit 142 calculates, and outputs the selected wavelength band. Number of the wavelength bands which the wavelength selection unit 2411 selects is P which is any one of 1 to (N−1). In the following explanation, the wavelength band which the wavelength selection unit 2411 selects is denoted as a selected wavelength band γ.sub.j (j is an integer of any one of s.sub.1 to s.sub.P).

[0083] For example, the wavelength selection unit 2411 compares the first intermediate-optical-path-radiance β.sup.D(γ.sub.i) and the second intermediate-optical-path-radiance β.sup.M(γ.sub.i) in turn for each wavelength band γ.sub.i out of N wavelength bands. In the case that the first intermediate-optical-path-radiance β.sup.D(γ.sub.i) is higher (larger) than the second intermediate-optical-path-radiance β.sup.M(γ.sub.i), the wavelength selection unit 2411 selects the wavelength band γ.sub.i and stores the wavelength band γ.sub.i as the selected wavelength band γ.sub.j. At every time when the wavelength selection unit 2411 stores the selected wavelength band γ.sub.j, the wavelength selection unit 2411 updates a value of j, and stores γs.sub.1 to γs.sub.P in turn.

[0084] The wavelength cost calculation unit 2412 calculates the cost C by using the first intermediate-optical-path-radiance β.sup.D(γ.sub.i) which the first intermediate-optical-path-radiance calculation unit 12 outputs, the second intermediate-optical-path-radiance β.sup.M(γ.sub.i) which the second intermediate-optical-path-radiance calculation unit 142 calculates, and the wavelength band γ.sub.i which the wavelength selection unit 2411 selects. For example, the wavelength cost calculation unit 2412 stores a predetermined value of a selected-wavelength cost coefficient A and a predetermined value of a not-selected-wavelength cost coefficient B. The selected-wavelength cost coefficient A is a coefficient representing a degree of influence of the first intermediate-optical-path-radiance on calculation of the final-optical-path-radiance in the wavelength band in which the dark pixel assumption is not met. The not-selected-wavelength cost coefficient B is a coefficient representing a degree of influence of the first intermediate-optical-path-radiance on calculation of the final-optical-path-radiance in the wavelength band in which the dark pixel assumption is met. The predetermined value of the selected-wavelength cost coefficient A is set to be smaller than the predetermined value of the not-selected-wavelength cost coefficient B. That is, the wavelength cost calculation unit 2412 gives weight, which is higher (larger) than weight of the wavelength band selected by the wavelength selection unit 2411, to the wavelength band not selected by the wavelength selection unit 2411, and calculates the cost C.

[0085] The wavelength cost calculation unit 2412 makes the wavelength band other than the selected wavelength band γ.sub.j out of N wavelength bands a not-selected wavelength band γ.sub.k (k is an integer of any one of t.sub.1 to t.sub.(N−P)). The wavelength cost calculation unit 2412 calculates the cost C by using Formula. 7. Here, a mark ‘.Math.’ in Formula. 7 indicates an operator of multiplication.

[00002] C = .Math. j = s 1 s P .Math. A .Math. ( β D ( γ j ) - β M ( γ j ) ) 2 + .Math. k = t 1 t ( N - P ) .Math. B .Math. ( β D ( γ k ) - β M ( γ j ) ) 2 [ Formula . .Math. 7 ]

[0086] Similarly to the first example embodiment, the image information processing device 200 according to the present example embodiment can calculate the optical-path-radiance with a high degree of accuracy even in the wavelength band in which the surface reflectance of the dark area in the image representing the observed ground surface area cannot be considered to be 0. The reason is that the final-optical-path-radiance calculation unit 24 calculates the cost, which indicates the degree of difference between the second-intermediate-optical-path-radiance and the actual-optical-path-radiance, by using the first intermediate-optical-path-radiance which is calculated by the first intermediate-optical-path-radiance calculation unit 12, and the second intermediate-optical-path-radiance which is calculated by the second intermediate-optical-path-radiance calculation unit 142, and outputs the second intermediate-optical-path-radiance, which is obtained at the time when the cost is minimized, as the final-optical-path-radiance.

[0087] The predetermined value of the selected-wavelength cost coefficient A, which is stored by the wavelength cost calculation unit 2412, is set to be smaller than the predetermined value of the not-selected-wavelength cost coefficient B. Therefore, the image information processing device 2 according to the present example embodiment reduces influence of the estimation error, that relates to estimation of the first intermediate-optical-path-radiance and that is generated in the wavelength band in which the dark pixel assumption is not met, on calculation of the final-optical-path-radiance. As a result, the image information processing device 2 according to the present example embodiment can calculate the optical-path-radiance with a higher degree of accuracy in comparison with the first example embodiment.

Third Example Embodiment

[0088] FIG. 4 is a block diagram which conceptually shows a configuration of an image information processing device 3 according to a third example embodiment.

[0089] The image information processing device 3 according to the present example embodiment includes a storage unit 31, a first intermediate-optical-path-radiance calculation unit 32, an irradiance calculation unit 33 and a final-optical-path-radiance calculation unit 34.

[0090] In the storage unit 31, observation images representing results of observing electromagnetic waves of a plurality of different wavelength bands reflected from the ground surface, information representing the wavelength bands, and information representing an observation environment existing at a time when the electromagnetic waves are observed are associated and stored.

[0091] For each of the wavelength band, the first intermediate-optical-path-radiance calculation unit 32 calculates radiance of the dark pixel, which meets the radiance standard, out of pixels composing the associated observation image which is associated with the wavelength band, as a first intermediate-optical-path-radiance.

[0092] The irradiance calculation unit 33 calculates, for each of the wavelength band, the irradiance from the sunlight on the basis of the information representing the observation environment.

[0093] The final-optical-path-radiance calculation unit 34 calculates, for each of the wavelength band, a final-optical-path-radiance on the basis of the irradiance and the first intermediate-optical-path-radiance.

[0094] The image information processing device 3 according to the present example embodiment can calculate the optical-path-radiance with a high degree of accuracy even in the wavelength band in which the surface reflectance of the dark area in the image representing the observed ground surface area cannot be considered to be 0. The reason is that the final-optical-path-radiance calculation unit 34 calculates the final-optical-path-radiance, whose value is near to a value of the actual-optical-path-radiance, for each wavelength band by using the first intermediate-optical-path-radiance which is calculated by the first intermediate-optical-path-radiance calculation unit 32, and the irradiance from the sunlight which is calculated by the irradiance calculation unit 33.

[0095] <Example of Hardware Configuration>

[0096] Each unit shown in FIG. 1, FIG. 3 and FIG. 4 according to each example embodiment can be realized by dedicated HW (HardWare) (electronic circuit).

[0097] Moreover, at least, each of the first intermediate-optical-path-radiance calculation units 12 and 32, the irradiance calculation units 13 and 33, and the final-optical-path-radiance calculation units 14, 24 and 34 can be considered to be a function (process) unit (software module) of a software program. However, division into each unit shown in each drawing is a configuration determined temporarily for convenience of explanation, and various configurations can be conceived at a time of implementation. An example of a hardware environment in this case will be explained in the following with reference to FIG. 7.

[0098] FIG. 7 is a diagram for exemplifying and explaining a configuration of an information processing device 900 (computer) which can carry out the image information processing device according to the example embodiment of the present invention. That is, FIG. 7 is the configuration of the computer (information processing device) which can realize the image information processing devices shown in FIG. 1, FIG. 3 and FIG. 4, and shows the hardware environment which can realize each function in the above-mentioned example embodiment.

[0099] The information processing device 900 shown in FIG. 7 includes components shown in the following. [0100] CPU (Central_Processing_Unit) 901 [0101] ROM (Read_Only_Memory) 902 [0102] RAM (Random_Access_Memory) 903 [0103] Hard disk (storage device) 904 [0104] Communication interface 905 communicating with an external device [0105] Bus 906 (communication line), [0106] Reader-writer 908 which can read and write data stored by a recording medium 907 such as CD-ROM (Compact Disc Read Only Memory) or the like [0107] Input-output interface 909
The information processing device 900 is a general computer in which these components are connected each other through the bus 906.

[0108] According to the present invention which has been explained by using the above-mentioned example embodiments as examples, a computer program, which can realize the following functions, is supplied to the information processing device 900 shown in FIG. 7. The functions are functions of the first intermediate optical path radiance calculation units 12 and 32, the irradiance calculation units 13 and 33, and the final intermediate optical path radiance calculation units 14, 24 and 34 in the block diagrams (FIG. 1, FIG. 3 and FIG. 4) which are referred to in explanation of the example embodiment, or functions of the flowchart (FIG. 2). The present invention is achieved by reading the computer program, and supplying CPU 901, which is the hardware, with the computer program, and interpreting the computer program, and then carrying out the computer program. Also, the computer program which is supplied to the device may be stored in a volatile storage device (RAM 903) which is capable of reading and writing, or a non-volatile storage device such as the hard disk 904 or the like.

[0109] Moreover, in the above-mentioned case, a procedure which is general at the present time is applicable to a method for supplying the hardware with the computer program. As the procedure, for example, a method of installation in the device through the various recording media 907 such as CD-ROM or the like, or a method of down-loading through a communication circuit such as The Internet or the like is exemplified. Moreover, in this case, it is conceivable that the present invention is made from codes which compose the computer program, or from the recording medium 907 which stores the codes.

[0110] While the invention has been particularly shown and described with reference to example embodiments thereof, the invention is not limited to these embodiments. It will be understood by those of ordinary skill in the art that various changes in form and details may be made therein without departing from the spirit and scope of the present invention as defined by the claims.

[0111] This application is based upon and claims the benefit of priority from Japanese patent application No. 2014-257039, filed on Dec. 19, 2014, the disclosure of which is incorporated herein in its entirety by reference.

INDUSTRIAL APPLICABILITY

[0112] The present invention is available for farming support and the resources exploration which are based on the result of measuring the ground surface from the height.

REFERENCE SIGNS LIST

[0113] 1 image information processing device [0114] 11 storage unit [0115] 12 first intermediate-optical-path-radiance calculation unit [0116] 13 irradiance calculation unit [0117] 14 final-optical-path-radiance calculation unit [0118] 141 cost calculation unit [0119] 142 second intermediate-optical-path-radiance calculation unit [0120] 143 cost minimization unit [0121] 15 optical-path-radiance correction unit [0122] 19 output unit [0123] 5 image supply device [0124] 6 output device [0125] 100 image information processing system [0126] 2 image information processing device [0127] 24 final-optical-path-radiance calculation unit [0128] 2411 wavelength selection unit [0129] 2412 wavelength cost calculation unit [0130] 200 image information processing system [0131] 3 image information processing device [0132] 31 storage unit [0133] 32 first intermediate-optical-path-radiance calculation unit [0134] 33 irradiance calculation unit [0135] 34 final-optical-path-radiance calculation unit [0136] 4 image information processing device [0137] 41 image reading unit [0138] 42 dark pixel identification unit [0139] 43 sensor parameter storage unit [0140] 44 atmosphere propagation correction unit [0141] 900 information processing device [0142] 901 CPU [0143] 902 ROM [0144] 903 RAM [0145] 904 hard disk (storage unit) [0146] 905 communication interface [0147] 906 bus [0148] 907 recording medium [0149] 908 reader-writer [0150] 909 Input-output interface