mgdevereux · March 29, 2022 23:46 · Mar 23, 2015
diff --git a/contrast_s.m b/contrast_s.m
@@ -0,0 +1 @@
+s  =  (cos(s * pi - pi) + 1) / 2;
diff --git a/noise3D.m b/noise3D.m
@@ -0,0 +1,12 @@
+function s = noise3D (m,k)
+  s = zeros([m,m,m]);
+  w = m;
+  i = 0;
+  while w > 3
+    i = i + 1;
+    d = smooth3(randn([m,m,m]), 'gaussian',k,i);
+    s = s + i * d(1:m, 1:m, 1:m);
+    w = w - ceil(w/2 - 1);
+  end
+  s = (s - min(min(s(:,:)))) ./ (max(max(s(:,:))) - min(min(s(:,:))));
+end
diff --git a/norm_s.m b/norm_s.m
@@ -0,0 +1 @@
+s = (s - min(min(s(:,:)))) ./ (max(max(s(:,:))) - min(min(s(:,:))));
diff --git a/perlin.m b/perlin.m
@@ -0,0 +1,11 @@
+function s = perlin (m)
+  s = zeros(m);    # output image
+  w = m;           # width of current layer
+  i = 0;           # iterations
+  while w > 3
+    i = i + 1;
+    d = interp2(randn(w), i-1, "spline");
+    s = s + i * d(1:m, 1:m);
+    w -= ceil(w/2 - 1);
+  end
+end
diff --git a/perlin2D.m b/perlin2D.m
@@ -0,0 +1,12 @@
+function s = perlin2D (m)
+  s = zeros([m,m]);     % Prepare output image (size: m x m)
+  w = m;
+  i = 0;
+  while w > 3
+    i = i + 1;
+    d = interp2(randn([m,m]), i-1, 'spline');
+    s = s + i * d(1:m, 1:m);
+    w = w - ceil(w/2 - 1);
+  end
+  s = (s - min(min(s(:,:)))) ./ (max(max(s(:,:))) - min(min(s(:,:))));
+end
diff --git a/perlin3D.m b/perlin3D.m
@@ -0,0 +1,12 @@
+function s = perlin3D (m)
+  s = zeros([m,m,m]);     % Prepare output image (size: m x m x m)
+  w = m;
+  i = 0;
+  while w > 3
+    i = i + 1;
+    d = interp3(randn([m,m,m]), i-1, 'spline');
+    s = s + i * d(1:m, 1:m, 1:m);
+    w = w - ceil(w/2 - 1);
+  end
+  s = (s - min(min(s(:,:)))) ./ (max(max(s(:,:))) - min(min(s(:,:))));
+end
Original file line number	Diff line number	Diff line change
		@@ -0,0 +1 @@
		s = (s - min(min(s(:,:)))) ./ (max(max(s(:,:))) - min(min(s(:,:))));