For a three channel remote-sensing data the following is an example of input ascii file
Eight models can be choosen by user after prompted by the program. They are:
The user has three options for check: multx1, multx2, and multx3. These options check multi-collinearity of the input data. multx1 calculates $R sup 2$ for channel x1 being replaced by coverage y. multx2 for channel x2 if input data are composed of two remote-sensing channels. And multx3 for channel x3 if input data are composed of at least three channels.
The program provides three methods to conduct the nonlinear regression calculation. These methods can be invoked by the user as options of calculat. The first one is rational fraction method. If no any option of calculat is given by the user, the program will adopt the first method. If calculat = rvi_linear the linearization method is taken to conduct the nonlinear relaxation index regression. If the user set calculat = rvi_nonlinear the general nonlinear regression method is used. If calculat = all, the program will conduct nonlinear regression first using the rational fraction method and then using the obtained coefficients as initial values to conduct the second method calculation and finally using the third method to improve the results. For nonlinear regression calculation there may exist multi-minimums. We can not be in reliance on the existing usual algorithms which can find one minimum only. The computer will not scan the possible minimums in order to save computer time. The user should judge and select a least minimum during the iteration. The necessary parameters which will facilitate the judgement for each iteration are displayed in the screen and also output to the output file.
For plots of relavant features of the regression calculation the user can set option for plot. There are four options for the plot. If plot = radiance-coverage computer will send radiance data against coverage data for further plotting radiance-coverage curves or figures to show the scatter of input data. This data will stored in ascii files named "curve.radiance_coverage_x1y", "curve.radiance_coverage_x2y" and "curve.radiance_coverage_x3y". If plot = vegetation-soil three ascii data files will be generated by the program for further figuers. Their names are "curve.veget_soil_x1x2", "curve.veget_soil_x1x3" and "curve.veget_soil_x2x3". If the user set plot = adequacy the program will generate residual data files for plotting figures of residuals against predictors and dependent variable to show the adequacy. These files named curve_residual_x1", "curve_residual_x2", "curve_residual_x3" and "curve_residual_y". If plot = all the program will generate all these above mentioned data files for further plotting of different figures.
The program not only can calculate regression models but also can make prediction to new remote sensing data using the obtained model. The user can set option of predict = same to calculate regression model using part of one imagery data and predict vegetation coverage for another part of remote sensing data in the same imagery. If predict = other the program will use one imagery to get regression model and predict vegetation coverage for other imagery. If the user did not give any option for predict the program calculates model based on the whole input data and not conduct any prediction. The program still generate a set of "calculated vegetation coverage" using the obtained model and the same input data in order for the user to check the utility, goodness and confidencial status of the regression. r.rational.regression will be run non-interactively if the user specifies program arguments on the command line, using the form:
But after run, the computer will prompt the user to select model number.
Alternately, the user can simply type:
on the command line without program arguments. In this case, the user will be prompted for parameter values using the standard GRASS user interface described in the manual entry for parser.
Michael Shapiro, U.S. Army Construction Engineering Research Laboratory.
Last changed: $Date: 2002/01/25 05:45:34 $