
DIRECTORY G2dBasic, G3dBasic, RefTab, Rope;

G3dVoxel: CEDAR DEFINITIONS

~ BEGIN
ROPE:						TYPE ~ Rope.ROPE;

Table:						TYPE ~ RefTab.Ref;
Nat3:						TYPE ~ G2dBasic.Nat3;
Type:						TYPE ~ {hash, array};

Triple:					TYPE ~ G2dBasic.Triple;				-- a point in space
RealSequence:			TYPE ~ G2dBasic.RealSequence;		-- sequence of x
RealSequenceRep:		TYPE ~ G2dBasic.RealSequenceRep;
Box:						TYPE ~ G3dBasic.Box;					-- bounding box
Error: SIGNAL [reason: ATOM];
Refs:						TYPE ~ RECORD [SEQUENCE maxLength: CARDINAL OF REF];
Refs2d:					TYPE ~ RECORD [SEQUENCE maxLength: CARDINAL OF REF Refs];
Refs3d:					TYPE ~ RECORD [SEQUENCE maxLength: CARDINAL OF REF Refs2d];

Voxels:					TYPE ~ REF VoxelsRep;
VoxelsRep:				TYPE ~ RECORD [
name:								ROPE м NIL,					-- name of Voxels
type:								Type м hash,				-- type of Voxels 
hashTable:						RefTab.Ref м NIL,			-- hash table
size:								Nat3 м [0, 0, 0],				-- size of array
array:								REF Refs3d м NIL,			-- array
clientData:						REF ANY м NIL				-- for miscellaneous data
];

CreateVoxels: PROC [sparse: BOOL м TRUE, size: Nat3 м [0, 0, 0], name: ROPE м NIL]
RETURNS [Voxels];

Store: PROC [voxels: Voxels, address: Nat3, voxel: REF];

RetrieveFromVoxels: PROC [voxels: Voxels, address: Nat3] RETURNS [REF];
Reals2d:					TYPE ~ REF Reals2dRep;				-- sequence of [x][y]
Reals3d:					TYPE ~ REF Reals3dRep;				-- sequence of [x][y][z]

Reals2dRep:				TYPE ~ RECORD [
length:								CARDINAL м 0,
element:								SEQUENCE maxLength: CARDINAL OF RealSequence];
Reals3dRep:				TYPE ~ RECORD [
length:								CARDINAL м 0,
element:								SEQUENCE maxLength: CARDINAL OF Reals2d];

RealProc:					TYPE ~ PROC [ref: REF ANY] RETURNS [value: REAL];

Grid:						TYPE ~ REF GridRep;
GridRep:					TYPE ~ RECORD [
res:								Nat3 м [],					-- resolution of data
evenlySpaced:					BOOL м TRUE,				-- if spacing is constant
scale:								Triple м [],					-- scale factor for indexing
box:								Box м [],				-- domain of the grid
x, y, z:							RealSequence м NIL,		-- indices along the three axes
values:							Reals3d м NIL,				-- functional values
clientData:						REF м NIL					-- for miscellaneous data
];

CreateGrid: PROC [
values: Reals3d,					-- values at the grid points
box: Box м [],					-- domain of the grid, presuming equal spacing
x, y, z: RealSequence м NIL,	-- domain of the grid, presuming unequal spacing
clientData: REF м NIL]
RETURNS [Grid];

SetGridBox: PROC [grid: Grid, box: Box];
AllocateGridValues: PROC [res: Nat3] RETURNS[Reals3d];
PointFromIJK: PROC [grid: Grid, i, j, k: NAT] RETURNS [Triple];
RetrieveFromGrid: PROC [grid: Grid, q: Triple] RETURNS [REAL];
ReadGridFromFile: PROC [fileName: ROPE] RETURNS [Grid];
WriteGridToFile: PROC [grid: Grid, fileName: ROPE, comment: ROPE м NIL];
GridFromVoxels: PROC [voxels: Voxels, realProc: RealProc] RETURNS [Grid];
END.
   Ў  
G3dVoxel.mesa
Copyright ╙ 1985, 1992 by Xerox Corporation.  All rights reserved.
Bloomenthal, July 15, 1992 11:07 pm PDT
Type Declarations
Errors
Storage and Retrieval of Voxel Data
Initialize the voxels data structure.  If sparse, a hash table is used.
Otherwise a three-dimensional array dimensioned size.i by size.j by size.k is created.
! Error[$sizeTooGreat].
Store a ref into voxels.
! Error[$addressOutOfRange]
! Error[$noVoxels]
Retrieve a ref from voxels.
! Error[$addressOutOfRange]
! Error[$noVoxels]
Storage and Retrieval of Function Values Within A Grid
Create a grid.  The indices x, y, and z must be in ascending order.
Can raise Error[$badIndices].
Scale and offset grid.x, grid.y and grid.z sequences so that the grid occupies box.

Allocate the refs.

Return the point that maps to grid[i][j][k].

Retrieve a value from the grid.  If q is out of grid's range, then Error[$outOfRange].

Read the grid data from a file, presumed to be formatted as:
NumberOfIndices~ <xRes> <yRes> <zRes>
XLocations~
sequence of reals (if only min and max specified, presumes even spacing)
YLocations~
sequence of reals (if only min and max specified, presumes even spacing)
ZLocations~
sequence of reals (if only min and max specified, presumes even spacing)
Values~
sequence of reals, looping most rapidly through zRes and most slowly through xRes.
! Error[$NoSuchFile].
! Error[$UnspecifiedDomain].

Write the grid data to the named file.

Create a grid from the voxel array; a voxel value is determined by realProc.

Э╩у  ХNewlineDelimiter
Ц"cedarcode" styleЩ ЩJЪЬ	╧eЬ6ЩBJЩ'JШ JЪ╧k	Ь"Ш+JШ ЧЪ╨blЬЮЬЮШJШ ЧJЪЬЮШheadЪ╧lЩЪЮЬЮЬЮЬШJШ ЧJЪЬЮ	ЬШJЪЬЮЬШ JЪЬЮЬШ JШ JЪЬ
ЮЬ╧cШ:JЪЬЮЬбШ?JЪЬЮЬШ2JЪЬЮЬбШ2ЧЪаЩJЪ╧nЬЮЬ
ЮЬШЧЪа#Щ#JЪЬЮ	ЬЮЬЮЬЮЬЮЬЮЬШ?JЪЬЮЬЮЬЮЬЮЬЮЬЮЬШEЪЬЮЬЮЬЮЬЮЬЮЬЮЬШGJШ ЧJЪЬЮЬЮЬШ!ЪЬЮЬЮЬШJЪЬЮЬЮЬбШ.JЪЬбШ/JЪЬЮЬбШ1JЪЬ"бШ2JЪЬЮЬ
ЮЬбШ*JЪЬЮЬЮЬЮЬбШ;JШJШ ЧЪвЬЮЬ
ЮЬЮЬ ЮЬЮЬШRJЪЮЬ
ШJЩGЩVJЩЧJШ ЧЪвЬЮЬ(ЮЬШ8ЩJЩJЩЧJШ ЧЪвЬЮЬ!ЮЬЮЬШGЩJЩJЩЧЧЧЪа6Щ6JЪЬЮЬЮЬбШ<ЪЬЮЬЮЬбШ?JШ ЧЪЬЮЬЮЬШJЪЬЮЬШJЪЬЮЬЮЬЮЬШ>ЧЪЬЮЬЮЬШJЪЬЮЬШJЪЬЮЬЮЬЮЬ
Ш9JШ ЧЪЬЮЬЮЬЮЬЮЬЮЬ	ЮЬШ?JШ ЧJЪЬ	ЮЬЮЬ	ШЪЬЮЬЮЬШJЪЬбШ0JЪЬЮЬЮЬбШ;JЪЬбШ;JЪЬбШ.JЪЬЮЬбШCJЪЬЮЬбШ4JЪЬЮЬЮЬбШ8JШJШ ЧЪв
ЬЮЬШJЪЬбШ1JЪЬб.ШAJЪЬЮЬб0ШMJЪЬЮЬЮЬШJЪЮЬШJЪЬCЩCJЪЬЩJШ ЧЪв
ЬЮЬШ(JЩSJЩ ЧЪвЬЮЬЮЬ
Ш6JЩJЩ ЧЪвЬЮЬЮЬЮЬ
Ш?JЩ,JЩ ЧЪвЬЮЬЮЬЮЬШ>JЩVJЩ ЧЪвЬЮЬЮЬЮЬШ7Щ<JЩ%ЪЬЩJЩHЧЩJЩHЧЩJЩHЧЩJЩRЧJЪЬЩЪЬЩJЩ ЧЧЧЪ
вЬЮЬЮЬЮЬЮЬШHJЩ&JЩ ЧЪвЬЮЬ&ЮЬШIJЩLЧJЩ ЧJЪЮЬШЧ ЕЧ       ∙  