shape@ main_shape()

{   

   image2d@ image = image2d("ektepar_ca.png");

   uint w = image.width();

   uint h = image.height();

   double s =0.002;

   double dx=400;

   double dy=400;

   auto surface = loft3d( image,dx,dy,w);

   return  mirror(0,1,0)*

           scale(s/5,s/5,2*s)*

     polyhedron(surface,vec3d(0,0,-300),h);

}

const double pi = 4.0*atan(1.0);

double to_rad(double deg) { return pi*deg/180.0; }  // from degree to radian

shape@ main_shape()

{   

   double h = 30;

   double r = 30;

   double e = 3;

   solid@ cyl = cylinder(h+10,r+e);

   // model space distance between the data points

   double dx=1;

   double dy=1;

   // specify that default segmentation shall be used, i.e.

   // all points in the file in s-dir (=x) and t-dir (=y)

   int ns_seg = -1;

   int nt_seg = -1;

   loft3d@ surface = loft3d("model.data.txt",dx,dy,ns_seg);

   solid@  wolf = translate(-23.1,-0.5,h+5) *

                  rotate_x(to_rad(-90)) *

                  scale(0.154, 0.109, 1)*

                  polyhedron(surface,vec3d(0,0,-30),nt_seg);   

   return union3d(cyl,wolf);

}

void main()

{

   shape@ obj = main_shape();

   obj.write_csg(GetOutputFullPath('.csg'));

}

2d utilities

pos2d    (position in 2d space)

array type [x,y]

spline2d (spline curve in 2d space)

N/A

image2d  (bitmap image support)

'file' in surface command

2d primitives

shape2d  (abstract type)

N/A

circle(double r);

circle(radius)

square(double size, bool center=false);

square(size,center)

rectangle(double dx, double dy, bool center=false);

square([width,height],center)

polygon(array<pos2d@> points);

polygon([points])

polygon(pos2d@ p1,...);

2d booleans

union2d(array<shape2d@> shapes);

union() { ... }  

union2d(shape2d@ shape, …);

difference2d(array<shape2d@> shapes);

difference() { ... }

difference2d(shape2d@ shape, …);

intersection2d(array<shape2d@> shapes);

intersection() { ... }

intersection2d(shape2d@ shape, …);

hull2d(array<shape2d@> shapes);

hull() {...}

hull2d(shape2d@ shape, …);

minkowski2d(array<shape2d@> shapes);

minkowski() {...}

minkowski2d(shape2d@ shape, …);

offset2d(shape2d@ shape, double r);

offset(r=...)

3d utilities

pos3d    (position in 3d space)

array type [x,y,z]

vec3d    (vector in 3d space)

array type [x,y,z]

spline3d (spline curve in 3d space)

N/A

loft3d   (lofted surface)

N/A

3d primitives

solid   (abstract type)

N/A

cone(double h, double r);

cylinder(h,r1,r2);

cube(double size);

cube(size);

cuboid(double dx, double dy, double dz);

cube([width,depth,height]);

cylinder(double h, double r);

cylinder(h,r);

sphere(double r);

sphere(r):

polyhedron(points,faces)

polyhedron(points,faces)

polyhedron(surface,offset)

polyhedron(points,faces)

3d booleans

union3d(array<solid@> shapes);

union() { ... }  

union3d(solid@ shape, …);

difference3d(array<solid@> shapes);

difference() { ... }

difference3d(solid@ shape, …);

intersection3d(array<solid@> shapes);

intersection() { ... }

intersection3d(solid@ shape, …);

hull3d(array<solid@> shapes);

hull() {...}

hull3d(solid@ shape, …);

minkowski3d(array<solid@> shapes);

minkowski() {...}

minkowski3d(solid@ shape, …);

linear_extrude(shape2d@ shape, double height);

linear_extrude(height=..) { ... }

rotate_extrude(shape2d@ shape, double angle);

rotate_extrude(angle=..) { ... }

2d/3d transformations

tmatrix   (abstract type)

multimatrix([ …. ]);

rotate_x(double rx);

rotate([rx,0,0]);

rotate_y(double ry);

rotate([0,ry,0]);

rotate_z(double rz);

rotate([0,0,rz]);

scale(double sx, double sy, double sz=1.0);

scale([x,y,z]);

translate(double dx, double dy, double dz=0.0);

translate([dx,dy,dz]);

mirror(double dx, double dy, double dz);

mirror([dx,dy,dz]);