Normal Texturing + SDF Meshing + Embeddable friendly (#123)

* embeddable friendly

* New Circle, Grid, and repeat pattern nodes. Fixes. New Quad Sphere.

* drawing/texturing/normals cleanup, new draw spheres node

* examples updated

* golfball example, dewarping quadsphere UVs, abstracting UVs concepts

* SDF nodes

* Marching cubes optimizations

* Multi point line SDF

Author: EliCDavis

README.md

@@ -75,6 +75,10 @@ Packages that have spawned from polyform's undertaking and have since been refac
 - [sfm](https://github.com/EliCDavis/sfm) - Utilities for interacting with reconstruction data from different SFM programs
 - [bitlib](https://github.com/EliCDavis/bitlib) - Utilities for reading and writing binary data
 
+## Contributing
+
+Learn how to [create your own nodes](./docs/guides/CreatingNodes/README.md) for others to use.
+
 ## Procedural Generation Examples
 
 You can at the different projects under the [examples](/examples/) folder for different examples on how to procedurally generate meshes.

cmd/polyform/main.go

@@ -18,6 +18,7 @@ import (
 	"github.com/EliCDavis/polyform/drawing/texturing"
 	_ "github.com/EliCDavis/polyform/drawing/texturing"
 	_ "github.com/EliCDavis/polyform/drawing/texturing/normals"
+	_ "github.com/EliCDavis/polyform/drawing/texturing/pattern"
 
 	_ "github.com/EliCDavis/polyform/formats/colmap"
 	_ "github.com/EliCDavis/polyform/formats/gltf"
@@ -42,6 +43,7 @@ import (
 	_ "github.com/EliCDavis/polyform/math/geometry"
 	_ "github.com/EliCDavis/polyform/math/noise"
 	_ "github.com/EliCDavis/polyform/math/quaternion"
+	_ "github.com/EliCDavis/polyform/math/sdf"
 	_ "github.com/EliCDavis/polyform/math/trig"
 	_ "github.com/EliCDavis/polyform/math/trs"
 	_ "github.com/EliCDavis/polyform/math/unit"
@@ -50,6 +52,7 @@ import (
 
 	_ "github.com/EliCDavis/polyform/modeling"
 	_ "github.com/EliCDavis/polyform/modeling/extrude"
+	_ "github.com/EliCDavis/polyform/modeling/marching"
 	_ "github.com/EliCDavis/polyform/modeling/meshops"
 	_ "github.com/EliCDavis/polyform/modeling/meshops/gausops"
 	_ "github.com/EliCDavis/polyform/modeling/primitives"

docs/resources/README.md

@@ -163,3 +163,7 @@ Resources either directly contributing to the code, or are just interesting find
 - Morton Encoding
   - [_"Morton encoding/decoding through bit interleaving: Implementations"_ Jeroen Baert's Blog](https://www.forceflow.be/2013/10/07/morton-encodingdecoding-through-bit-interleaving-implementations/)
     - [Code](https://github.com/Forceflow/libmorton)
+- Quad Spheres
+  - [_"Wraparound square tile maps on a sphere"_ by Redblob](https://www.redblobgames.com/x/1938-square-tiling-of-sphere/)
+  - [_"Mapping a Cube to a Sphere"_ by Harry van Langen](https://hvlanalysis.blogspot.com/2023/05/mapping-cube-to-sphere.html?m=1)
+  - [_"Survey of Cube Mapping Methods in Interactive Computer Graphics"_ by M. Lambers](https://marlam.de/publications/cubemaps/lambers2019cubemaps.pdf)

drawing/coloring/operations.go

@@ -236,3 +236,11 @@ func Interpolate(a, b color.Color, t float64) color.Color {
 		A: uint8(aVal),
 	}
 }
+
+func Greyscale(c color.Color) float64 {
+	r, g, b, _ := c.RGBA()
+	rF := float64(r >> 8)
+	gF := float64(g >> 8)
+	bF := float64(b >> 8)
+	return (rF + gF + bF) / (255 * 3)
+}

drawing/texturing/blur.go

@@ -118,8 +118,8 @@ func convolve1DGaussian[T any](space vector.Space[T], src Texture[T], dst Textur
 func RadialGaussianBlur[T any](space vector.Space[T], src Texture[T], radius int, sigma float64) Texture[T] {
 	kernel := gaussianKernel(radius, sigma)
 
-	tmp := NewTexture[T](src.width, src.height)
-	out := NewTexture[T](src.width, src.height)
+	tmp := Empty[T](src.width, src.height)
+	out := Empty[T](src.width, src.height)
 
 	// Horizontal then vertical
 	convolve1DGaussian(space, src, tmp, kernel, true)

drawing/texturing/color.go

@@ -9,7 +9,7 @@ func DivideColor(a, b Texture[coloring.Color]) (*Texture[coloring.Color], error)
 		return nil, ErrMismatchDimensions
 	}
 
-	result := NewTexture[coloring.Color](a.width, a.height)
+	result := Empty[coloring.Color](a.width, a.height)
 	a.ScanParallel(func(x, y int, v coloring.Color) {
 		other := b.Get(x, y)
 		result.Set(x, y, coloring.Color{
@@ -25,7 +25,7 @@ func DivideColor(a, b Texture[coloring.Color]) (*Texture[coloring.Color], error)
 
 func MaxColor(a Texture[coloring.Color], other float64) Texture[coloring.Color] {
 
-	result := NewTexture[coloring.Color](a.width, a.height)
+	result := Empty[coloring.Color](a.width, a.height)
 	a.ScanParallel(func(x, y int, v coloring.Color) {
 		result.Set(x, y, coloring.Color{
 			R: max(v.R, other),
@@ -40,7 +40,7 @@ func MaxColor(a Texture[coloring.Color], other float64) Texture[coloring.Color]
 
 func ClampColor(a Texture[coloring.Color], minimum, maximum float64) Texture[coloring.Color] {
 
-	result := NewTexture[coloring.Color](a.width, a.height)
+	result := Empty[coloring.Color](a.width, a.height)
 	a.ScanParallel(func(x, y int, v coloring.Color) {
 		result.Set(x, y, coloring.Color{
 			R: min(max(v.R, minimum), maximum),

drawing/texturing/errors.go

@@ -1,5 +1,17 @@
 package texturing
 
-import "errors"
+import (
+	"errors"
+	"fmt"
+
+	"github.com/EliCDavis/vector/vector2"
+)
 
 var ErrMismatchDimensions = errors.New("mismatch texture resolutions")
+
+type InvalidDimension vector2.Vector[int]
+
+func (id InvalidDimension) Error() string {
+	v := vector2.Vector[int](id)
+	return fmt.Sprintf("invalid texture dimensions %dx%d", v.X(), v.Y())
+}

drawing/texturing/float.go

@@ -7,7 +7,7 @@ import (
 )
 
 func OneMinus(in Texture[float64]) Texture[float64] {
-	result := NewTexture[float64](in.width, in.height)
+	result := Empty[float64](in.width, in.height)
 	for y := range in.height {
 		for x := range in.width {
 			result.Set(x, y, 1-in.Get(x, y))
@@ -47,7 +47,7 @@ func (n MultiplyFloat1Node) Result(out *nodes.StructOutput[Texture[float64]]) {
 		return
 	}
 
-	result := NewTexture[float64](textures[0].Width(), textures[0].Height())
+	result := Empty[float64](textures[0].Width(), textures[0].Height())
 	for y := range result.Height() {
 		for x := range result.Width() {
 			v := textures[0].Get(x, y)

drawing/texturing/gradient.go

@@ -18,7 +18,7 @@ type RadialGradient[T any] struct {
 }
 
 func (rg RadialGradient[T]) Texture() Texture[T] {
-	tex := NewTexture[T](rg.Width, rg.Height)
+	tex := Empty[T](rg.Width, rg.Height)
 	center := vector2.New(rg.Width, rg.Height).
 		ToFloat64().
 		Scale(0.5)
@@ -45,7 +45,7 @@ func (n RadialGradientNode[T]) LinearGradient(out *nodes.StructOutput[Texture[T]
 	height := nodes.TryGetOutputValue(out, n.Height, 1)
 
 	if n.Gradient == nil {
-		out.Set(NewTexture[T](width, height))
+		out.Set(Empty[T](width, height))
 		return
 	}
 
@@ -69,7 +69,7 @@ type LinearGradient[T any] struct {
 }
 
 func (lg LinearGradient[T]) Texture() Texture[T] {
-	tex := NewTexture[T](lg.Width, lg.Height)
+	tex := Empty[T](lg.Width, lg.Height)
 
 	boxCenter := vector3.New(lg.Width, lg.Height, 0).
 		ToFloat64().
@@ -109,7 +109,7 @@ func (n LinearGradientNode[T]) LinearGradient(out *nodes.StructOutput[Texture[T]
 	height := nodes.TryGetOutputValue(out, n.Height, 1)
 
 	if n.Gradient == nil {
-		out.Set(NewTexture[T](width, height))
+		out.Set(Empty[T](width, height))
 		return
 	}
 
@@ -124,7 +124,7 @@ func (n LinearGradientNode[T]) LinearGradient(out *nodes.StructOutput[Texture[T]
 }
 
 func ApplyGradient[T any](time Texture[float64], gradient coloring.Gradient[T]) Texture[T] {
-	result := NewTexture[T](time.width, time.height)
+	result := Empty[T](time.width, time.height)
 	for y := range result.height {
 		for x := range result.width {
 			result.Set(x, y, gradient.Sample(time.Get(x, y)))
@@ -145,7 +145,7 @@ func (n ApplyGradientNode[T]) Texture(out *nodes.StructOutput[Texture[T]]) {
 
 	time := nodes.GetOutputValue(out, n.Time)
 	if n.Gradient == nil {
-		out.Set(NewTexture[T](time.Width(), time.Height()))
+		out.Set(Empty[T](time.Width(), time.Height()))
 		return
 	}
 

drawing/texturing/noise.go

@@ -34,7 +34,7 @@ func texture(
 		out.CaptureError(fmt.Errorf("invalid height dimension: %d", height))
 		return
 	}
-	t := NewTexture[float64](width, height)
+	t := Empty[float64](width, height)
 
 	scale := nodes.TryGetOutputValue(out, Scale, vector2.One[float64]())
 	offset := nodes.TryGetOutputValue(out, Offset, 0)