Keyword: Aerosol optical depth, structure function, SPOT

Abstract:
In applying the contrast method, such s the structure function method, to estimate the atmospheric aerosol depth with satellite data, a uniform land cover area in image in chosen as a test are to avoid the probable errors induced by the poor structure function patterns . unfortunately, the uniform are is not easy to pick up or does not exist in some complex land use regions, such as in Taiwan. In order to pursue the potential application of the structure function method, especially for complex terrain area, our study extends the original single-directional structure ;function to multi-directional, and introduces an " optimal number" in to our procedure to improve the accuracy; of aerosol depth estimation. The comparison between the estimated and sun photometer-observed aerosol depths shows that the accuracy is improved significantly by our improvements.

1. Introduction
Basically, satellite remote sensed data are affected mainly by the scattering effects of atmospheric molecular and aerosols in the visible band, Fraser et al., 1984.some atmospheric correcting schemes have been proposed in applying to the satellite visible band data, such as Landsat and SPOT images in the past studies ( Griggs, 1975; Mekler et al., 1977; Tanre et al., Rao et al., 1989; holben B.n. et.al., 1990 and Liu et al., 1997). Some studies showed that the structure function method can be used to estimate accurately the atmospheric aerosol depth, however obvious errors could be induced in same cases by the structure function patterns. Further analysis indicates that athe abnormal patterns probably are caused by the change of satellite observation geometry, the change or the complexity of land cover. Therefore , the aim of this study is to try to provide a solution reduce the probable errors in applying the structure method and improve the accuracy of optical depth estimation.

2.Methodology
Assuming the surface observed by satellite sensor is Lambertian , the apparent reflectance , ?*, observed by satellite is ( Tanre et al., 1988).

where m_s= cosq_s , q_s is the solar zenith angle, msub>v= cosq_v, q_v is the observed zenith angle, p^a is the atmospheric reflectance, f is the relative azimuth anlgle between the sun and the satellite, t is the optical depth, p is the surface reflectance, T is the transmittance from the sun to surface, s is the atmospheric albideo, <p> is the mean surface reflectance, and t_d is the diffused transmittance from surface to satellite.

The nulti-scattering effect of the surface and atmosphere is small and can be neglected. In other words, <?> s=0, Assuming the <?> values to be the same in local area, the apparent reflectance differences between two neighboring pixels , (I,j) ? (I, j+d), in distance d can be written as

??*_(I,j) = ?? _(I,j) T (m_s) exp[-t/mv]            (2)

where I, j are the row and column index of image, respectively. If we assume the ? _(I,j) value being constant in time, the ?*_(I,j) will be the function oft , which depends upon atmosphere condition. Tanre et al. (1988) has defined a structure function parameter, M, as



Where N is the total pixel number in the test are. In their original method, a single direction difference is calculated in the structure function. However , in our study, we consider a multi-directional( I, j and cross (c ) directions shown in Fig .1) structure function, which can be expressed as

and the structure function M* (d) derived by satellite observation can is written as

M*² (d) =M² (d) T² (m_s) exp [-2t/m_v]             (5)

Where M is the real surface structure function .
d=0          1      2      3    ...

	?j
?i	?c

Figure 1 the directions of I, j and c of the multi-directional structure function.
If assuming the land cover remaining unchanged from t1 to t2 observation time i.e. M² (d,t₁) = M²(d,t₂) , the relationship of observed structure functions between t1 to t2 is,

Equation ( 6) shows that if one of the aerosol optical depth in t₁ or t₂ is known, the other optical depth can be derived from satellite observation.

Analysis shows that the estimated optical depth could be different for different d values. In this general cases, the mean value of optical depth from d=1 to d=10 is used, but significant errors could be observed in some cases when the correlation is low for different distances, so the " optimal number " is adopted to remove the abnormal structure function areas 9 Liu et al., 1997). The introduction of " optimal number can remove the poor structure functions and improve the accuracy of the optical depth estimation .