A linear relationship y = ax + b was assumed between the
reflectance (x) and the depth (y). The rectangular domain D may be viewed as a
matrix of order m x n, which is comprised of reflectance matrix D
r
and depth matrix D
d . In the above relationship, a and b were
calculated using the least square approximation in the domain D as
=
The depths at the 2(m+n+2) immediate neibouring pixels of the
domain D were computed using the transformation equation
-------------------- (1)
This process is repeated by updating the domain until it merges
with the full region of interest.
A model region (study area) was selected in conjunction with
the satellite imagery and this can be viewed as a 1024 x 1024 pixels. In this
region, the bathymetry data was available only for a small area [91x193 pixels].
A 10 x 10 matrix with significant depth variation was selected as the basic
domain D, where both reflectance and depth values were available. The Digital
Numbers representing reflectance and depth matrices of D are shown in Table 1
and Table 2.
Statistical analysis shows that these two matrices are
correlated with a correlation coefficient of -0.99. The immediate neighbouring
elements of D were computed using the transformation equation (1) and D was
updated to 12 x12 matrix, including the new elements. The reflectance and depth
matrices of order 12x12 corresponding to the updated domain are shown in Table 3
and Table 4 . The observed depth values of this new domain are shown in Table 5.
The comparison of observed and computed depths versus reflectance is given in
Fig. 2. The process was repeated with the updated domain until the whole region
is covered.
Results and Conclusions
The Domain Rippling Scheme is a
useful technique for the estimation of water depths over shallow water region.
The digital satellite image with reflectance values is successfully transformed
to depth image using this technique and is shown in Plate 2. A rigorous
correlation analysis was done for the observed and computed depth values with
various matrix sizes. The correlation analysis results are briefed in Table 6.
It is found that the computed and observed depth values have good correlation,
particularly in the shallow water region (Fig. 2). It can be seen that for the
area under study, the intensity variation with depth is proportionate in shallow
water region upto -3m depth contour. Hence it can be concluded that the present
method can be restricted to shallow water depths. The analysis can be improved
by assuming a nonlinear relationship between the reflectance and the water depth
for deeper water region. Also it is important to consider the parameters like
suspended sediment concentration and bed characteristics in the model which
influence the reflectance.