Menu
▾
▴
[Plplot-cvs] SF.net SVN: plplot:[8708] trunk/examples/tcl
|
From: <arj...@us...> - 2008-08-24 10:05:29
|
Revision: 8708
http://plplot.svn.sourceforge.net/plplot/?rev=8708&view=rev
Author: arjenmarkus
Date: 2008-08-24 10:05:39 +0000 (Sun, 24 Aug 2008)
Log Message:
-----------
Added first implementation of example x21 - API not complete yet, so stops
at the call to plgriddata
Added Paths:
-----------
trunk/examples/tcl/x21
trunk/examples/tcl/x21.tcl
Added: trunk/examples/tcl/x21
===================================================================
--- trunk/examples/tcl/x21 (rev 0)
+++ trunk/examples/tcl/x21 2008-08-24 10:05:39 UTC (rev 8708)
@@ -0,0 +1,19 @@
+#!/bin/sh
+#--------------------------------*- Tcl -*------------------------------------#
+# $Id$
+#
+# Maurice LeBrun
+# 12/24/02
+#
+# A front-end to x06.tcl for running directly from the command line, locating
+# pltcl via PATH.
+# Handles all usual plplot command arguments. See "pltcl -h" for info.
+#-----------------------------------------------------------------------------#
+#\
+exec pltcl -f "$0" ${1+"$@"}
+
+source x21.tcl
+
+plinit
+x21
+plend
Property changes on: trunk/examples/tcl/x21
___________________________________________________________________
Added: svn:executable
+ *
Added: svn:keywords
+ Author Date Id Revision
Added: svn:eol-style
+ native
Added: trunk/examples/tcl/x21.tcl
===================================================================
--- trunk/examples/tcl/x21.tcl (rev 0)
+++ trunk/examples/tcl/x21.tcl 2008-08-24 10:05:39 UTC (rev 8708)
@@ -0,0 +1,290 @@
+# $Id$
+# Grid data demo.
+#
+# Copyright (C) 2004 Joao Cardoso
+# Copyright (C) 2008 Andrew Ross
+#
+# This file is part of PLplot.
+#
+# PLplot is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Library Public License as
+# published by the Free Software Foundation; either version 2 of the
+# License, or (at your option) any later version.
+#
+# PLplot is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+# GNU Library General Public License for more details.
+#
+# You should have received a copy of the GNU Library General Public
+# License along with PLplot; if not, write to the Free Software
+# Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
+#
+# Note:
+# This example uses NaNs (not-a-number) and these special "numbers"
+# are only properly supported by Tcl 8.5 and later (in the sense
+# that a NaN does not stop the program with an error message).
+
+proc x21 {{w loopback}} {
+
+ if { [catch {package require Tcl 8.5}] } {
+ puts "This example require Tcl 8.5 or later: use of NaNs"
+ return
+ }
+
+ set PI [expr {4.0*atan(1.0)}]
+
+ set GRID_CSA 1
+ set GRID_DTLI 2
+ set GRID_NNI 3
+ set GRID_NNIDW 4
+ set GRID_NNLI 5
+ set GRID_NNAIDW 6
+
+ set DRAW_LINEXY 3
+ set MAG_COLOR 4
+ set BASE_CONT 8
+
+ set pts 500
+ set xp 25
+ set yp 20
+ set nl 16
+ set knn_order 20
+ set threshold 1.001
+ set wmin -1e3
+ set randn 0
+ set rosen 0
+
+ matrix x f $pts
+ matrix y f $pts
+ matrix z f $pts
+ matrix clev f $pts
+ matrix xg f $xp
+ matrix yg f $yp
+ matrix zg f $xp $yp
+
+ set title {"Cubic Spline Approximation"
+ "Delaunay Linear Interpolation"
+ "Natural Neighbors Interpolation"
+ "KNN Inv. Distance Weighted"
+ "3NN Linear Interpolation"
+ "4NN Around Inv. Dist. Weighted"}
+
+ matrix opt f 6 = { 0.0 0.0 0.0 0.0 0.0 0.0 }
+
+ set xmin -0.2
+ set ymin -0.2
+ set xmax 0.6
+ set ymax 0.6
+
+ opt 3 = $wmin
+ opt 4 = [expr {double($knn_order)}]
+ opt 5 = $threshold
+
+ for {set i 0} {$i < $pts} {incr i} {
+ set xt [expr {($xmax-$xmin)*[plrandd]}]
+ set yt [expr {($ymax-$ymin)*[plrandd]}]
+ if {$randn == 0} {
+ x $i = [expr {$xt + $xmin}]
+ y $i = [expr {$yt + $ymin}]
+ } else {
+ x $i = [expr {sqrt(-2.*log($xt)) * cos(2.*$PI*$yt) + $xmin}]
+ y $i = [expr {sqrt(-2.*log($xt)) * sin(2.*$PI*$yt) + $ymin}]
+ }
+ if {$rosen == 0} {
+ set xx [x $i]
+ set yy [x $i]
+ set r [expr {sqrt($xx*$xx + $yy*$yy)}]
+ z $i = [expr {exp(-$r*$r)*cos(2.*$PI*$r)}]
+ } else {
+ set xx [x $i]
+ set yy [x $i]
+ z $i = [expr {log((1.-$xx)**2 + 100.*($yy-$xx**2)**2)}]
+ }
+ }
+
+ set zmin [z 0]
+ set zmax [z 0]
+ for {set i 1} {$i < $pts} {incr i} {
+ set zmax [max $zmax [z $i]]
+ set zmin [min $zmin [z $i]]
+ }
+
+ for {set i 0} {$i < $xp} {incr i} {
+ xg $i = [expr {$xmin + ($xmax-$xmin)*($i-1.)/($xp-1.)}]
+ }
+ for {set i 0} {$i < $yp} {incr i} {
+ yg $i = [expr {$ymin + ($ymax-$ymin)*($i-1.)/($yp-1.)}]
+ }
+
+ $w cmd plcol0 1
+ $w cmd plenv $xmin $xmax $ymin $ymax 2 0
+ $w cmd plcol0 15
+ $w cmd pllab "X" "Y" "The original data sampling"
+ $w cmd plcol0 2
+ $w cmd plpoin $pts x y 5
+ $w cmd pladv 0
+
+ $w cmd plssub 3 2
+
+ for {set k 0} {$k < 2} {incr k} {
+ $w cmd pladv 0
+ for {set alg 1} {$alg <= 6} {incr alg} {
+
+ $w cmd plgriddata x y z xg yg zg $alg [opt $alg]
+
+# - CSA can generate NaNs (only interpolates? #).
+# - DTLI and NNI can generate NaNs for points outside the convex hull
+# of the data points.
+# - NNLI can generate NaNs if a sufficiently thick triangle is not found
+#
+# PLplot should be NaN/Inf aware, but changing it now is quite a job...
+# so, instead of not plotting the NaN regions, a weighted average over
+# the neighbors is done.
+#
+
+ if {($alg == $GRID_CSA) || ($alg == $GRID_DTLI) ||
+ ($alg == $GRID_NNLI) || ($alg == $GRID_NNI} {
+
+ for {set i 0} {$i < $xp} {incr i} {
+ for {set j 0} {$j < $yp} {incr j} {
+ if { [zg $i $j] == NaN } {
+# average (IDW) over the 8 neighbors
+
+ zg $i $j = 0.
+ set dist 0.
+
+ set ii [expr {$i-1}]
+ while {($ii < $i+1) && ($ii < $xp)} {
+ set jj [expr {$j-1}]
+ while {($jj == $j+1) && ($jj < $yp)} {
+ if {($ii >= 0) && (jj >= 1) &&
+ ![isnan [zg $ii $jj]] } {
+ if {abs($ii-$i) + abs($jj-$j) == 1)} {
+ set d 1.
+ } else {
+ set d 1.4142
+ }
+ zg $i $j = [expr {[zg $i $j] + [zg $ii $jj]/($d*$d)}]
+ set dist [expr {$dist + $d}]
+ }
+ incr jj
+ }
+ incr ii
+ }
+ if {$dist != 0.} {
+ zg $i $j = [expr {[zg $i $j] / $dist}]
+ } else {
+ zg $i $j = $zmin
+ }
+ }
+ }
+ }
+ }
+
+ a2mnmx zg xp yp lzmin lzmax xp
+
+ set lzmin = [min $lzmin $zmin]
+ set lzmax = [max $lzmax $zmax]
+
+ set lzmin = [expr {lzmin - 0.01}]
+ set lzmax = [expr {lzmax + 0.01}]
+
+ $w cmd plcol0 1
+ $w cmd pladv $alg
+
+ if {$k == 0} {
+
+ for {set i 0} {$i < $nl} {incr i} {
+ clev $i = [expr {$lzmin + ($lzmax-$lzmin)/($nl-1.)*($i-1.)}]
+ }
+ $w cmd plenv0 xmin xmax ymin ymax 2 0
+ $w cmd plcol0 15
+ $w cmd pllab "X" "Y" [lindex $title $alg]
+ $w cmd plshades zg $defined $xmin $xmax $ymin \
+ $ymax clev 1 0 1
+ $w cmd plcol0 2
+ } else {
+
+ for {set i 0} {$i < $nl} {incr i} {
+ clev $i = [expr {$lzmin + ($lzmax-$lzmin)/($nl-1.)*($i-1.)}]
+ }
+ cmap1_init
+ $w cmd plvpor 0. 1. 0. 0.9
+ $w cmd plwind -1.1 0.75 -0.65 1.20
+#
+# For the comparison to be fair, all plots should have the
+# same z values, but to get the max/min of the data generated
+# by all algorithms would imply two passes. Keep it simple.
+#
+# plw3d(1., 1., 1., xmin, xmax, ymin, ymax, zmin, zmax, 30, -60);
+#
+
+ $w cmd plw3d 1. 1. 1. $xmin $xmax $ymin $ymax \
+ $lzmin $lzmax 30. -40.
+ $w cmd plbox3 "bntu" "X" 0. 0 \
+ "bntu" "Y" 0. 0 \
+ "bcdfntu" "Z" 0.5 0
+ $w cmd plcol0 15
+ $w cmd pllab "" "" [lindex $title $alg]
+ $w cmd plot3dc xg yg zg ior ior $DRAW_LINEXY
+ $MAG_COLOR $BASE_CONT clev
+ }
+ }
+ }
+
+ $w cmd plend
+}
+
+
+#----------------------------------------------------------------------------
+# proc max and min
+proc max {x y} {
+ expr {$x > $y? $x : $y}
+}
+proc min {x y} {
+ expr {$x > $y? $y : $x}
+}
+
+
+#----------------------------------------------------------------------------
+# proc cmap1_init
+# Set up the colour map
+proc cmap1_init {} {
+
+ matrix i i 2
+ matrix h i 2
+ matrix l i 2
+ matrix s i 2
+
+ i 1 = 0.
+ i 2 = 1.
+
+ h 1 = 240.
+ h 2 = 0.
+
+ l 1 = 0.6
+ l 2 = 0.6
+
+ s 1 = 0.8
+ s 2 = 0.8
+
+ $w cmd plscmap1n 256
+ $w cmd plscmap1l .false. i h l s
+}
+
+#----------------------------------------------------------------------------
+# proc a2mnmx
+# Minimum and the maximum elements of a 2-d array.
+
+proc a2mnmx {f nx ny fmin fmax xdim} {
+
+ set fmax [f 0 0]
+ set fmin $fmax
+ for {set j 0} {$j < $ny} {incr j} {
+ for {set i 0} {$i < $nx} {incr i} {
+ set fmax [max $fmax [f $i $j]]
+ set fmin [min $fmin [f $i $j]]
+ }
+ }
+}
Property changes on: trunk/examples/tcl/x21.tcl
___________________________________________________________________
Added: svn:keywords
+ Author Date Id Revision
Added: svn:eol-style
+ native
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
|
