area_surf.tcl

This example shows how to build surface by approximating other surface in areas. In this example 5 different areas were used. Suppose that for four of them (area1, area2, area3 and area4) we know surface values approximately, and for one of them (area6) we know them exactly. So, we want to build surface, which approximates as exactly as possible values for area6 and "not so exactly" for areas 1-4.

Here is the image of the resulting surface (red colour for areas 1-4, green colour for area6):

Pic 1. Surface built from regions of other surface

This result was obtained with "examples/surfit/curvs/area_surf.tcl" script.

#
# This script was generated with funner v2.1beta (www.gridding.info)

# 
# load plugins
# 
load libsurfit[info sharedlibextension]

# remove all previous data and gridding rules
clear_data 

# set name of surface
set map_name "map_area_surf" 

# set solver
set_solver "cg" 

# set tolerance for solver
set tol 1e-005 

##
## load initial data 
##
 
# load area from text file 
area_read "area1.txt" "area1"  
 
# load area from text file 
area_read "area2.txt" "area2"  
 
# load area from text file 
area_read "area3.txt" "area3"  
 
# load area from text file 
area_read "area4.txt" "area4"  
 
# load area from text file 
area_read "area6.txt" "area6"  
 
# load surface from surfit datafile 
surf_load "func.func" "test_func"  
 
##
## construct grid 
##
grid_get -12 12 0.2 -12 12 0.2 
 
##
## create gridding rules
##

# resulting surface in area = surface values... 
area_surf "test_func" "area6"  

# resulting surface should tend to be constant or plane 
completer 

# add "resulting surface in area = surface values..." with weight 
area_surf_add "test_func" 0.07 "area1"  

# add "resulting surface in area = surface values..." with weight 
area_surf_add "test_func" 0.1 "area2"  

# add "resulting surface in area = surface values..." with weight 
area_surf_add "test_func" 0.1 "area3"  

# add "resulting surface in area = surface values..." with weight 
area_surf_add "test_func" 0.05 "area4"  

##
## run gridding algorithm
##
surfit 

##
## save results 
##

# unload grid from memory
grid_unload 

# save surface to surfit datafile 
surf_save "area_surf.dat" "map_area_surf" 

---

surfit: open source gridding and contouring software.