Problem: My study area is a simple rectangular polygon in which the sampled points are spread inside. The extent of this polygon is obviously greater than the extent of the bounding box of input points. Now, the Geostatistical Wizard does not honor the (extent) environment settings. Also, Kriging is not given as an independent tool in the GA toolbox which can be supplied with the extent environment settings. As much as I know, only the wizard can do together (a) data transformations, (b) semivariogram modelling, and (c) neighbourhood parameterization, but it won't produce the result in the entire extent of the polygon!
Alternative adopted: Added four points at the four corners of the study area polygon. I updated their z-values with the mean of the sampled points. This causes the wizard to produce output in the exact extent of the study area polygon. The interpolation in the new 'extra region' would have used data from the sampled points and these 4 extra points.
Discussion: Is the above method a good alternative to the environment settings problem? I believe that a lot of users would come across this problem because a lot of times the output extent would be bigger than the sampling extent. Regarding the use of spatial analyst kriging tool, there's the only the option of simple and ordinary kriging and the tool is less flexible than the GA one. Is there a workflow of tools available to perform the above kriging job? Please enrich my knowledge.