summaryrefslogtreecommitdiff
path: root/src/modules/Terrain.js
blob: 21908ac0774c5e8a79b4b36ae6247caec4e12b2d (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
'use strict';

import Voronoi from './3rdparty/rhill-voronoi-core.js';

import { useAverage } from './Util.js';
import { AABB, QuadTree } from './Geometry.js';


/* from here on up, we always assume that points live in the range [(0,0), (1,1)) */

function lloydRelax(point_set, density) {
	/* setup quadtree and averages */
	let tree = new QuadTree(1,1);
	let averages = {};
	for (let i=0; i<point_set.length; i++) {
		const point = point_set[i];
		point.index = i;
		tree.insert(point);

		let [avg, append] = useAverage();
		const cent_x = { avg, append };
		[avg, append] = useAverage();
		const cent_y = { avg, append };
		averages[i] = { cent_x, cent_y };
	}

	/* compute average centroids */
	for (let x=0; x<1; x += 1/density) {
		for (let y=0; y<1; y += 1/density) {
			const point = { x, y };
			const closest = tree.closest(point);
			const { cent_x, cent_y } = averages[closest.index];
			cent_x.append(point.x);
			cent_y.append(point.y);
		}
	}

	/* return centroid points */
	const result = [];
	for (let i=0; i<point_set.length; i++) {
		const point = { x: averages[i].cent_x.avg(), y: averages[i].cent_y.avg() };
		result.push(point);
	}
	return result;
}


class Terrain {
	constructor() {
		const N_SEED_POINTS = 2**12;
		const N_RELAX_ITERATIONS = 1;
		const RELAX_DENSITY = 400;
		const randomPoint = () => ({x: Math.random(), y: Math.random()});

		this.min_height = 0;
		this.max_height = 0;

		let seed_points = [];
		for (let i=0; i<N_SEED_POINTS; i++) seed_points.push(randomPoint());

		for (let i=0; i<N_RELAX_ITERATIONS; i++)
			lloydRelax(seed_points, RELAX_DENSITY);

		const v = new Voronoi();
		this.graph = v.compute(seed_points, {xl: 0, xr: 1, yt: 0, yb: 1});

		this.tree = new QuadTree(1,1);
		for (let v of this.graph.vertices) {
			v.height = 0;
			this.tree.insert(v);
		}
	}


	applyBrush(x, y, f, strength, radius) {
		const region = new AABB(x-radius, y-radius, 2*radius, 2*radius);
		const points = this.tree.root.getPointsInRegion(region);
		
		const dist2 = (a, b) => (a.x - b.x)**2 + (a.y - b.y)**2;

		const sigma = radius/3;
		const beta = 1/(2*sigma*sigma);
		const center = { x, y };
		const power = pt => Math.exp(-beta * dist2(pt, center));
		
		for (let pt of points) f(pt, strength * power(pt));
	}


	renderGrid(ct) {
		ct.save();
		ct.lineWidth = 0.001;
		for (let edge of this.graph.edges) {
			ct.fillStyle = `hsl(${edge.va.height}, 100%, 50%)`;
			ct.beginPath();
			ct.arc(edge.va.x, edge.va.y, 0.005, 0, 2*Math.PI);
			ct.closePath();
			ct.fill();

			ct.beginPath();
			ct.moveTo(edge.va.x, edge.va.y);
			ct.lineTo(edge.vb.x, edge.vb.y);
			ct.closePath();
			ct.stroke();
		}
		ct.restore();
	}
}

export { lloydRelax };
export default Terrain;