joonsrenderer 1.1-java

data/src/main/java/org/sunflow/image/RegularSpectralCurve.java

@@ -0,0 +1,30 @@
+package org.sunflow.image;
+
+public class RegularSpectralCurve extends SpectralCurve {
+
+    private final float[] spectrum;
+    private final float lambdaMin, lambdaMax;
+    private final float delta, invDelta;
+
+    public RegularSpectralCurve(float[] spectrum, float lambdaMin, float lambdaMax) {
+        this.lambdaMin = lambdaMin;
+        this.lambdaMax = lambdaMax;
+        this.spectrum = spectrum;
+        delta = (lambdaMax - lambdaMin) / (spectrum.length - 1);
+        invDelta = 1 / delta;
+    }
+
+    @Override
+    public float sample(float lambda) {
+        // reject wavelengths outside the valid range
+        if (lambda < lambdaMin || lambda > lambdaMax) {
+            return 0;
+        }
+        // interpolate the two closest samples linearly
+        float x = (lambda - lambdaMin) * invDelta;
+        int b0 = (int) x;
+        int b1 = Math.min(b0 + 1, spectrum.length - 1);
+        float dx = x - b0;
+        return (1 - dx) * spectrum[b0] + dx * spectrum[b1];
+    }
+}

data/src/main/java/org/sunflow/image/SpectralCurve.java

@@ -0,0 +1,118 @@
+package org.sunflow.image;
+
+/**
+ * This class is an abstract interface to sampled or analytic spectral data.
+ */
+public abstract class SpectralCurve {
+
+    /**
+     * This function determines the actual spectral curve data. Note that the
+     * lambda parameter is assumed to be in nanometers.
+     *
+     * @param lambda wavelength to sample in nanometers
+     * @return the value of the spectral curve at this point
+     */
+    public abstract float sample(float lambda);
+    private static final int WAVELENGTH_MIN = 360;
+    private static final int WAVELENGTH_MAX = 830;
+    private static final double[] CIE_XBAR = {0.000129900000, 0.000232100000,
+        0.000414900000, 0.000741600000, 0.001368000000, 0.002236000000,
+        0.004243000000, 0.007650000000, 0.014310000000, 0.023190000000,
+        0.043510000000, 0.077630000000, 0.134380000000, 0.214770000000,
+        0.283900000000, 0.328500000000, 0.348280000000, 0.348060000000,
+        0.336200000000, 0.318700000000, 0.290800000000, 0.251100000000,
+        0.195360000000, 0.142100000000, 0.095640000000, 0.057950010000,
+        0.032010000000, 0.014700000000, 0.004900000000, 0.002400000000,
+        0.009300000000, 0.029100000000, 0.063270000000, 0.109600000000,
+        0.165500000000, 0.225749900000, 0.290400000000, 0.359700000000,
+        0.433449900000, 0.512050100000, 0.594500000000, 0.678400000000,
+        0.762100000000, 0.842500000000, 0.916300000000, 0.978600000000,
+        1.026300000000, 1.056700000000, 1.062200000000, 1.045600000000,
+        1.002600000000, 0.938400000000, 0.854449900000, 0.751400000000,
+        0.642400000000, 0.541900000000, 0.447900000000, 0.360800000000,
+        0.283500000000, 0.218700000000, 0.164900000000, 0.121200000000,
+        0.087400000000, 0.063600000000, 0.046770000000, 0.032900000000,
+        0.022700000000, 0.015840000000, 0.011359160000, 0.008110916000,
+        0.005790346000, 0.004106457000, 0.002899327000, 0.002049190000,
+        0.001439971000, 0.000999949300, 0.000690078600, 0.000476021300,
+        0.000332301100, 0.000234826100, 0.000166150500, 0.000117413000,
+        0.000083075270, 0.000058706520, 0.000041509940, 0.000029353260,
+        0.000020673830, 0.000014559770, 0.000010253980, 0.000007221456,
+        0.000005085868, 0.000003581652, 0.000002522525, 0.000001776509,
+        0.000001251141,};
+    private static final double[] CIE_YBAR = {0.000003917000, 0.000006965000,
+        0.000012390000, 0.000022020000, 0.000039000000, 0.000064000000,
+        0.000120000000, 0.000217000000, 0.000396000000, 0.000640000000,
+        0.001210000000, 0.002180000000, 0.004000000000, 0.007300000000,
+        0.011600000000, 0.016840000000, 0.023000000000, 0.029800000000,
+        0.038000000000, 0.048000000000, 0.060000000000, 0.073900000000,
+        0.090980000000, 0.112600000000, 0.139020000000, 0.169300000000,
+        0.208020000000, 0.258600000000, 0.323000000000, 0.407300000000,
+        0.503000000000, 0.608200000000, 0.710000000000, 0.793200000000,
+        0.862000000000, 0.914850100000, 0.954000000000, 0.980300000000,
+        0.994950100000, 1.000000000000, 0.995000000000, 0.978600000000,
+        0.952000000000, 0.915400000000, 0.870000000000, 0.816300000000,
+        0.757000000000, 0.694900000000, 0.631000000000, 0.566800000000,
+        0.503000000000, 0.441200000000, 0.381000000000, 0.321000000000,
+        0.265000000000, 0.217000000000, 0.175000000000, 0.138200000000,
+        0.107000000000, 0.081600000000, 0.061000000000, 0.044580000000,
+        0.032000000000, 0.023200000000, 0.017000000000, 0.011920000000,
+        0.008210000000, 0.005723000000, 0.004102000000, 0.002929000000,
+        0.002091000000, 0.001484000000, 0.001047000000, 0.000740000000,
+        0.000520000000, 0.000361100000, 0.000249200000, 0.000171900000,
+        0.000120000000, 0.000084800000, 0.000060000000, 0.000042400000,
+        0.000030000000, 0.000021200000, 0.000014990000, 0.000010600000,
+        0.000007465700, 0.000005257800, 0.000003702900, 0.000002607800,
+        0.000001836600, 0.000001293400, 0.000000910930, 0.000000641530,
+        0.000000451810,};
+    private static final double[] CIE_ZBAR = {0.000606100000, 0.001086000000,
+        0.001946000000, 0.003486000000, 0.006450001000, 0.010549990000,
+        0.020050010000, 0.036210000000, 0.067850010000, 0.110200000000,
+        0.207400000000, 0.371300000000, 0.645600000000, 1.039050100000,
+        1.385600000000, 1.622960000000, 1.747060000000, 1.782600000000,
+        1.772110000000, 1.744100000000, 1.669200000000, 1.528100000000,
+        1.287640000000, 1.041900000000, 0.812950100000, 0.616200000000,
+        0.465180000000, 0.353300000000, 0.272000000000, 0.212300000000,
+        0.158200000000, 0.111700000000, 0.078249990000, 0.057250010000,
+        0.042160000000, 0.029840000000, 0.020300000000, 0.013400000000,
+        0.008749999000, 0.005749999000, 0.003900000000, 0.002749999000,
+        0.002100000000, 0.001800000000, 0.001650001000, 0.001400000000,
+        0.001100000000, 0.001000000000, 0.000800000000, 0.000600000000,
+        0.000340000000, 0.000240000000, 0.000190000000, 0.000100000000,
+        0.000049999990, 0.000030000000, 0.000020000000, 0.000010000000,
+        0.000000000000, 0.000000000000, 0.000000000000, 0.000000000000,
+        0.000000000000, 0.000000000000, 0.000000000000, 0.000000000000,
+        0.000000000000, 0.000000000000, 0.000000000000, 0.000000000000,
+        0.000000000000, 0.000000000000, 0.000000000000, 0.000000000000,
+        0.000000000000, 0.000000000000, 0.000000000000, 0.000000000000,
+        0.000000000000, 0.000000000000, 0.000000000000, 0.000000000000,
+        0.000000000000, 0.000000000000, 0.000000000000, 0.000000000000,
+        0.000000000000, 0.000000000000, 0.000000000000, 0.000000000000,
+        0.000000000000, 0.000000000000, 0.000000000000, 0.000000000000,
+        0.000000000000,};
+    private static final int WAVELENGTH_STEP = (WAVELENGTH_MAX - WAVELENGTH_MIN) / (CIE_XBAR.length - 1);
+
+    static {
+        if (WAVELENGTH_STEP * (CIE_XBAR.length - 1) != WAVELENGTH_MAX - WAVELENGTH_MIN) {
+            String err = String.format("Internal error - spectrum static data is inconsistent!\n  * min = %d\n  * max = %d\n  * step = %d\n  * num = %d", WAVELENGTH_MIN, WAVELENGTH_MAX, WAVELENGTH_STEP, CIE_XBAR.length);
+            throw new RuntimeException(err);
+        }
+    }
+
+    /**
+     * Convert this curve to a tristimulus CIE XYZ color by integrating against
+     * the CIE color matching functions.
+     *
+     * @return XYZColor that represents this spectra
+     */
+    public final XYZColor toXYZ() {
+        float X = 0, Y = 0, Z = 0;
+        for (int i = 0, w = WAVELENGTH_MIN; i < CIE_XBAR.length; i++, w += WAVELENGTH_STEP) {
+            float s = sample(w);
+            X += s * CIE_XBAR[i];
+            Y += s * CIE_YBAR[i];
+            Z += s * CIE_ZBAR[i];
+        }
+        return new XYZColor(X, Y, Z).mul(WAVELENGTH_STEP);
+    }
+}

data/src/main/java/org/sunflow/image/XYZColor.java

@@ -0,0 +1,50 @@
+package org.sunflow.image;
+
+public final class XYZColor {
+
+    private float X, Y, Z;
+
+    public XYZColor() {
+    }
+
+    public XYZColor(float X, float Y, float Z) {
+        this.X = X;
+        this.Y = Y;
+        this.Z = Z;
+    }
+
+    public final float getX() {
+        return X;
+    }
+
+    public final float getY() {
+        return Y;
+    }
+
+    public final float getZ() {
+        return Z;
+    }
+
+    public final XYZColor mul(float s) {
+        X *= s;
+        Y *= s;
+        Z *= s;
+        return this;
+    }
+
+    public final void normalize() {
+        float XYZ = X + Y + Z;
+        if (XYZ < 1e-6f) {
+            return;
+        }
+        float s = 1 / XYZ;
+        X *= s;
+        Y *= s;
+        Z *= s;
+    }
+
+    @Override
+    public final String toString() {
+        return String.format("(%.3f, %.3f, %.3f)", X, Y, Z);
+    }
+}

data/src/main/java/org/sunflow/image/formats/BitmapBlack.java

@@ -0,0 +1,27 @@
+package org.sunflow.image.formats;
+
+import org.sunflow.image.Bitmap;
+import org.sunflow.image.Color;
+
+public class BitmapBlack extends Bitmap {
+
+    @Override
+    public int getWidth() {
+        return 1;
+    }
+
+    @Override
+    public int getHeight() {
+        return 1;
+    }
+
+    @Override
+    public Color readColor(int x, int y) {
+        return Color.BLACK;
+    }
+
+    @Override
+    public float readAlpha(int x, int y) {
+        return 0;
+    }
+}

data/src/main/java/org/sunflow/image/formats/BitmapG8.java

@@ -0,0 +1,36 @@
+package org.sunflow.image.formats;
+
+import org.sunflow.image.Bitmap;
+import org.sunflow.image.Color;
+
+public class BitmapG8 extends Bitmap {
+
+    private int w, h;
+    private byte[] data;
+
+    public BitmapG8(int w, int h, byte[] data) {
+        this.w = w;
+        this.h = h;
+        this.data = data;
+    }
+
+    @Override
+    public int getWidth() {
+        return w;
+    }
+
+    @Override
+    public int getHeight() {
+        return h;
+    }
+
+    @Override
+    public Color readColor(int x, int y) {
+        return new Color((data[x + y * w] & 0xFF) * INV255);
+    }
+
+    @Override
+    public float readAlpha(int x, int y) {
+        return 1;
+    }
+}

data/src/main/java/org/sunflow/image/formats/BitmapGA8.java

@@ -0,0 +1,30 @@
+package org.sunflow.image.formats;
+
+import org.sunflow.image.Bitmap;
+import org.sunflow.image.Color;
+
+public class BitmapGA8 extends Bitmap {
+
+    private int w, h;
+    private byte[] data;
+
+    @Override
+    public int getWidth() {
+        return w;
+    }
+
+    @Override
+    public int getHeight() {
+        return h;
+    }
+
+    @Override
+    public Color readColor(int x, int y) {
+        return new Color((data[2 * (x + y * w) + 0] & 0xFF) * INV255);
+    }
+
+    @Override
+    public float readAlpha(int x, int y) {
+        return (data[2 * (x + y * w) + 1] & 0xFF) * INV255;
+    }
+}

data/src/main/java/org/sunflow/image/formats/BitmapRGB8.java

@@ -0,0 +1,40 @@
+package org.sunflow.image.formats;
+
+import org.sunflow.image.Bitmap;
+import org.sunflow.image.Color;
+
+public class BitmapRGB8 extends Bitmap {
+
+    private int w, h;
+    private byte[] data;
+
+    public BitmapRGB8(int w, int h, byte[] data) {
+        this.w = w;
+        this.h = h;
+        this.data = data;
+    }
+
+    @Override
+    public int getWidth() {
+        return w;
+    }
+
+    @Override
+    public int getHeight() {
+        return h;
+    }
+
+    @Override
+    public Color readColor(int x, int y) {
+        int index = 3 * (x + y * w);
+        float r = (data[index + 0] & 0xFF) * INV255;
+        float g = (data[index + 1] & 0xFF) * INV255;
+        float b = (data[index + 2] & 0xFF) * INV255;
+        return new Color(r, g, b);
+    }
+
+    @Override
+    public float readAlpha(int x, int y) {
+        return 1;
+    }
+}

data/src/main/java/org/sunflow/image/formats/BitmapRGBA8.java

@@ -0,0 +1,40 @@
+package org.sunflow.image.formats;
+
+import org.sunflow.image.Bitmap;
+import org.sunflow.image.Color;
+
+public class BitmapRGBA8 extends Bitmap {
+
+    private int w, h;
+    private byte[] data;
+
+    public BitmapRGBA8(int w, int h, byte[] data) {
+        this.w = w;
+        this.h = h;
+        this.data = data;
+    }
+
+    @Override
+    public int getWidth() {
+        return w;
+    }
+
+    @Override
+    public int getHeight() {
+        return h;
+    }
+
+    @Override
+    public Color readColor(int x, int y) {
+        int index = 4 * (x + y * w);
+        float r = (data[index + 0] & 0xFF) * INV255;
+        float g = (data[index + 1] & 0xFF) * INV255;
+        float b = (data[index + 2] & 0xFF) * INV255;
+        return new Color(r, g, b);
+    }
+
+    @Override
+    public float readAlpha(int x, int y) {
+        return (data[4 * (x + y * w) + 3] & 0xFF) * INV255;
+    }
+}

data/src/main/java/org/sunflow/image/formats/BitmapRGBE.java

@@ -0,0 +1,60 @@
+package org.sunflow.image.formats;
+
+import org.sunflow.image.Bitmap;
+import org.sunflow.image.Color;
+
+public class BitmapRGBE extends Bitmap {
+
+    private int w, h;
+    private int[] data;
+    private static final float[] EXPONENT = new float[256];
+
+    static {
+        EXPONENT[0] = 0;
+        for (int i = 1; i < 256; i++) {
+            float f = 1.0f;
+            int e = i - (128 + 8);
+            if (e > 0) {
+                for (int j = 0; j < e; j++) {
+                    f *= 2.0f;
+                }
+            } else {
+                for (int j = 0; j < -e; j++) {
+                    f *= 0.5f;
+                }
+            }
+            EXPONENT[i] = f;
+        }
+    }
+
+    public BitmapRGBE(int w, int h, int[] data) {
+        this.w = w;
+        this.h = h;
+        this.data = data;
+    }
+
+    @Override
+    public int getWidth() {
+        return w;
+    }
+
+    @Override
+    public int getHeight() {
+        return h;
+    }
+
+    @Override
+    public Color readColor(int x, int y) {
+        int rgbe = data[x + y * w];
+        float f = EXPONENT[rgbe & 0xFF];
+        float r = f * ((rgbe >>> 24) + 0.5f);
+        float g = f * (((rgbe >> 16) & 0xFF) + 0.5f);
+        float b = f * (((rgbe >> 8) & 0xFF) + 0.5f);
+        return new Color(r, g, b);
+    }
+
+    @Override
+    public float readAlpha(int x, int y) {
+        return 1;
+    }
+}

data/src/main/java/org/sunflow/image/formats/BitmapXYZ.java

@@ -0,0 +1,38 @@
+package org.sunflow.image.formats;
+
+import org.sunflow.image.Bitmap;
+import org.sunflow.image.Color;
+import org.sunflow.image.XYZColor;
+
+public class BitmapXYZ extends Bitmap {
+
+    private int w, h;
+    private float[] data;
+
+    public BitmapXYZ(int w, int h, float[] data) {
+        this.w = w;
+        this.h = h;
+        this.data = data;
+    }
+
+    @Override
+    public int getWidth() {
+        return w;
+    }
+
+    @Override
+    public int getHeight() {
+        return h;
+    }
+
+    @Override
+    public Color readColor(int x, int y) {
+        int index = 3 * (x + y * w);
+        return Color.NATIVE_SPACE.convertXYZtoRGB(new XYZColor(data[index], data[index + 1], data[index + 2])).mul(0.1f);
+    }
+
+    @Override
+    public float readAlpha(int x, int y) {
+        return 1;
+    }
+}