area_mean.tcl

This example shows how to build surface by approximating areas with constant real numbers and how to set additional condition (mean value) on surface in area. In this example 6 different areas were used. This example looks like area.tcl example, but one area (area6) with mean condition was added.

See also:

area_mean, area, area_add, completer

Here is the image of the resulting surface (green colour for areas 1,3 and 5, red colour for areas 2 and 4, magenta for area 6):

Pic 1. Surface with mean value = 30, built from areas

This result was obtained with "examples/surfit/curvs/area_mean.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_mean" 

# set tolerance for cg solver
set tol 5.01187e-006 

##
## 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 "area5.txt" "area5"  
 
# load area from text file 
area_read "area6.txt" "area6"  
 
##
## construct grid 
##
grid_get -12 12 0.2 -12 12 0.2 
 
##
## create gridding rules
##

# resulting surface in area = value... 
area 10 "area1"  

# resulting surface in area = value... 
area 30 "area3"  

# resulting surface in area = value... 
area undef "area5"  

# resulting surface mean value in area = value... 
area_mean 40 "area6"  

# resulting surface mean value = value... 
mean 30 

# resulting surface should tend to be constant or plane 
completer 1 1.6  

# add "resulting surface in area = value..." with weight 
area_add 20 0.07 "area2"  

# add "surface resulting in area = value..." with weight 
area_add 40 0.07 "area4"  

##
## run gridding algorithm
##
surfit 

##
## save results 
##

# unload grid from memory
grid_unload 

# save surface to surfit datafile 
surf_save "area_mean.dat" "map_area_mean" 

puts [concat "Surface mean value in area6 : " [surf_area_mean "area6" "map_area_mean"]]
puts [concat "Surface mean value : " [surf_mean "map_area_mean"]]

---

surfit: open source gridding and contouring software.