A two-phase sampling design has been applied to forest inventory. First, a
large number of first phase sample plots were defined with a square grid in
a geographic coordinate system for two study areas of about 1800 and 4500
ha. The first phase sample plots were supplied by auxiliary data of Landsat
TM and IRS-1C with principal component transformation for stratification a
nd drawing the second phase sample (field sample). Proportional allocation
was used to draw the second phase sample. The number of field sample plots
in the two study areas was 300 and 380.
The local estimates of five continuous forest stand variables, mean diamete
r, mean height, age, basal area, and stem volume, were calculated for each
of the first phase sample plots. This was done separately by using one auxi
liary data source at a time together with the held sample information. Howe
ver, if the first phase sample plot for which the stand variables were to b
e estimated was also a field sample plot, the information of that field sam
ple plot was eliminated according to the cross validation principle. This w
as because it was then possible to calculate mean square errors of estimate
s related to a specific auxiliary data source.
The procedure produced as many estimates for each first phase sample plot a
nd forest stand variable as was the number of auxiliary data sources, i.e.
seven estimates: These were based on Landsat TM, IRS-1C, digitized aerial p
hotos, ocular stereoscopic interpretation from aerial photographs, data fro
m old forest inventory made by compartments, Landsat TM95-TM89 difference i
mage and IRS96-TM95 difference image. The final estimates were calculated a
s weighted averages where the weights were inversely proportional to mean s
quare errors. The alternative estimates were calculated by applying simple
rules based on knowledge and the outliers were defined. The study shows tha
t this kind of system for finding outliers for elimination and a weighting
procedure improves the accuracy of stand variable estimation.