Monte Carlo Based Non-Linear Mixture Model of Earth Observation Satellite Imagery Pixel Data

Monte Carlo based non-linear mixel (mixed pixel) model of visible to near infrared radiometer of earth observation satellite imagery is proposed. Through comparative studies with actual real earth observation satellite imagery data between conventional linear mixel model and the proposed non-linear mixel model, it is found that the proposed mixel model represents the pixels in concern much precisely rather than the conventional linear mixel model.


INTRODUCTION
The pixels in earth observed images which are acquired with Visible to Near Infrared: VNIR sensors onboard remote sensing satellites are, essentially mixed pixels (mixels) which consists of several ground cover materials [1].Some mixel model is required for analysis such as un-mixing of the mixel in concern [2], [3].Typical mixel is linear mixing model which is represented by linear combination of several ground cover materials with mixing ratio for each material [4].It is not always true that the linear mixel model is appropriate [5].Due to the influences from multiple reflections between the atmosphere and ground, multiple scattering in the atmosphere on the observed radiance from the ground surface, pixel mixture model is essentially non-linear rather than linear.These influence is interpreted as adjacency effect [6].[7].
Method for representation of non-linear mixel model is not so easy.In particular, there is not sophisticated multi reflection model between ground materials.The representation method for non-linear mixel model is based on Monte Carlo Ray Tracing: MCRT model [8].It is rather easy to designate surface slopes on the ground and multi reflection among trees for MCRT model.The proposed MCRT based non-linear mixel model is applied to real earth observation satellite imagery data of Advanced Spaceborn Thermal Emission and Reflection Radiometer / Visible and Near Infrared Radiometer: ASTER/VNIR onboard on Terra satellite.A comparison of radiance between the conventional linear mixel model and the proposed non-linear mixel model is conducted.As a result, validity of the proposed model is confirmed.
The following section describes the proposed non-linear mixel model based on MCRT followed by some experiments for validation of the proposed model.Then, finally, conclusions with some discussions are described.

A. Monte CarloRay Tracing Simulation
In order to show a validity of the proposed non-linear mixel model, MCRT simulation study and field experimental study is conducted.MCRT allows simulation of polarization characteristics of sea surface with designated parameters of the atmospheric conditions and sea surface and sea water conditions.Illustrative view of MCRT is shown in Fig. 1.Photon from the sun is input from the top of the atmosphere (the top of the simulation cell).Travel length of the photon is calculated with optical depth of the atmospheric molecule and that of aerosol.There are two components in the atmosphere; molecule and aerosol particles while three are also two components, water and particles; suspended solid and phytoplankton in the ocean.When the photon meets molecule or aerosol (the meeting probability with molecule and aerosol depends on their optical depth), then the photon scattered in accordance with scattering properties of molecule and aerosol.The scattering property is called as phase function 1 .In the visible to near infrared wavelength region, the scattering by molecule is followed by Rayleigh scattering law [10] while that by aerosol is followed by Mie scattering law [10].Example of phase function of Mie scattering is shown in Fig. 2 (a) while that of Rayleigh scattering is shown in Fig. 2 (b).www.ijacsa.thesai.orgIn the atmosphere, there are absorption due to water vapor, ozone and aerosols together with scattering due to the atmospheric molecules, aerosols.Atmospheric Optical Depth: AOD (optical thickness) in total, Optical Depth: OD due to water vapor (H 2 O), ozone (O 3 ), molecules (MOL), aerosols (AER), and real observed OD (OBS) are plotted in Fig. 3 as an example.
For simplifying the calculations of the atmospheric influences, it is assumed that the atmosphere containing only molecules and aerosols.As shown in Fig. 3, this assumption is not so bad.Thus the travel length of the photon at once, L is expressed with equation (1).
where Z max, τ, RND(i) are maximum length, altitude of the atmosphere, optical depth, and i-th random number, respectively.In this equation, τ is optical depth of molecule or aerosol.The photon meets molecule when the random number is greater than τ.Meanwhile, if the random number is less than τ, then the photon meats aerosol.The photon is scattered at the molecule or aerosol to the direction which is determined with the aforementioned phase function and with the rest of the travel length of the photon.

B. Ground Surface with Slopes
When the photon reaches on the ground, the photon reflects at the ground surface to the direction which is determined by random number.Lambertian surface [11] is assumed.Therefore, reflectance is constant for all the directions.The reflected photon travels with the rest of travel length.Two adjacent slopes of Lambertian surfaces are assumed on the ground as shown in Fig. 4. Slope angles for both are α、β while their reflectance are Γ A and Γ B where N + is the number of photons which are gathered by VNIR, N total denotes the number of photons input to the simulation cell.Also I denote extraterrestrial irradiance at the top of the atmosphere.www.ijacsa.thesai.org

A. Validity of the Monte CarloRay Tracing Simulation
In order to confirm that the developed MCRT is valid, a comparative study is conducted between radiative transfer code of Gauss Seidel method and the MCRT derived TOA radiance.Because the Gauss Seidel method allows calculation of TOA radiance with flat surface of ground, 0.2 of reflectance of flat surface is assumed in the comparison.Also, 0.02 and 0.03 of optical depths are assumed for aerosol and molecule.
The size of simulation cell is determined as 50 km by 50 km by 50 km.Solar zenith angle is set at 30 degree while solar azimuth is set at 120 degree.700,000 of photons are input to the simulation cell.TOA radiance derived from the Gauss Seidel method is 0.565 (mW/m 2 /sr/μm) while that from the MCRT is 0.579 (mW/m 2 /sr/μm) at the 500nm of wavelength.For both cases, IFOV of the VNIR radiometer is assumed to be 2 π , all of the photons output from the top of the atmosphere are counted.Therefore, the developed MCRT seems valid enough.

B. TOA Radiance for the Different Combination of Optical
Depths of Aerosol and Molecule and for the Ground with the Different Slopes TOA radiance at 500 nm of wavelength for the different combination of optical depths of aerosol and molecule which ranges from 0.01 to 0.04 and for the ground with the different slopes, 0 and 20 degree are calculated.
Again, IFOV of the VNIR radiometer is assumed to be 2π, all of the photons output from the top of the atmosphere are counted.The reflectance for both slopes are same as 0.      Other than these, topographic map of three test sites which is corresponding to the Google map is shown in Fig. 6 (c) while the extracted portion of each test site on ASTER/VNIR image is shown in Fig. 6 (d), (e) and (f), respectively.These digital elevation models for three test sites are taken into account in the MCRT simulations.www.ijacsa.thesai.orgAlso, solar zenith angle of 58 degree and solar azimuth angle of 17 degree are taken into account in the simulations.From the atmospheric optical depth measurement data with sun photometer, optical depth of total atmosphere is calculated.Furthermore, molecule optical depth τ R is calculated with equation ( 4) as a function of atmospheric pressure P which is measured on the ground.(4) where P 0 denotes standard atmospheric pressure on the ground (1013 hPa) while λ denotes wavelength.Then aerosol optical depth is calculated from total atmospheric optical depth by subtracting molecule optical depth.www.ijacsa.thesai.org Comparative study is conducted between ASTER/VNIR derived radiance of Band 2 (Green band) and the radiance which derived from the conventional linear mixel model and the proposed non-linear mixel model.
Table 3 shows the calculated radiance in unit of μW/m 2 /sr/μm and the radiance difference between ASTER/VNIR and the estimated with the conventional and the proposed mixel models.It is found that the estimated radiance with the proposed non-linear mixel model is much closer rather than that with the conventional linear mixel model.

IV. CONCLUSION
Monte Carlo based non-linear mixel (mixed pixel) model of visible to near infrared radiometer of earth observation satellite imagery is proposed.Through comparative studies between ASTER/VNIR derived radiance and the conventional linear mixel model derived radiance as well as the proposed non-linear mixel model derived radiance, it is found that the estimated radiance with the proposed non-linear mixel model is much closer to ASTER/VNIR derived radiance (around 6%) rather than that with the conventional linear mixel model.
One of the disadvantages of the proposed non-linear mixel model based on MCRT is time consumable computations.Acceleration is highly required.

Figure 1
Figure 1 Illustrative view of MCRT for the atmosphere and sea water (a)Mie scattering (aerosol) (b)Rayleigh scattering (Molecule) Figure 2 Phase functions for Mie and Rayleigh scattering

Figure 3
Figure 3 Example of observed atmospheric optical depth in total and the best fit curves of optical depth due to water vapor, ozone, molecules, and aerosols calculated with MODTRAN of atmospheric radiative transfer software code..

Figure 4
Figure 4 Two adjacent slopes of Lambertian surfaces which are assumed on the ground C. Top of the Atmosphere: TOA Radiance Calculation If the photon reaches on the wall of the simulation cell, the photon disappears at the wall and it appears from the corresponding position on the opposite side wall.Then it travels with the rest of travel length.Eventually, the photons which are reached at the top of the atmosphere are gathered with the Instantaneous Field of View: IFOV of the Visible to Near Infrared Radiometer: VNIR onboard satellite.At sensor radiance, I + with direction and IFOV of μ, μ 0 can be calculated with equation (3)

C.
Validity of the Proposed Non-Linear Mixel Model with Real VNIR Data The proposed non-linear mixel model based on MCRT is validated with real earth observation satellite imagery data of ASTER/VNIR onboard Terra satellite which is acquired at 11:09 Japanese Standard Time: JST on December 15 2004.IFOV of ASTER/VNIR is 15m with 60km of swath width.Whole scene of ASTER/VNIR is shown in Fig.5 (a) while Fig.5 (b) shows a portion of the scene.
Figure 5 ASTER/VNIR image used for experiment Three test sites, Area #1, 2, 3 are extracted from the scene.Attribute information of these sites are listed in Table2. .

Fig. 6 (
Fig.6 (a) shows three test sites on ASTER/VNIR image while Fig.6 (b) shows three test sites on Google map.Other than these, topographic map of three test sites which is corresponding to the Google map is shown in Fig.6 (c) while the extracted portion of each test site on ASTER/VNIR image is shown in Fig.6 (d), (e) and (f), respectively.These digital elevation models for three test sites are taken into account in the MCRT simulations.
5. The results are shown in Table 1.In the table, τ aer , τ mol are optical depths of aerosol and molecule, respectively.

TABLE I .
TOP OF THE ATMOSPHERE: TOA RADIANCE FOR THE COMBINATIONS OF ATMOSPHERIC CONDITIONS

TABLE II .
ATTRIBUTIONS FOR THE TEST SITE WITH SLOPES