spherelight.java

Sunflow是一个照片级的渲染系统

package org.sunflow.core.light;

import org.sunflow.SunflowAPI;
import org.sunflow.core.LightSample;
import org.sunflow.core.LightSource;
import org.sunflow.core.ParameterList;
import org.sunflow.core.Ray;
import org.sunflow.core.Shader;
import org.sunflow.core.ShadingState;
import org.sunflow.core.primitive.Sphere;
import org.sunflow.image.Color;
import org.sunflow.math.Matrix4;
import org.sunflow.math.OrthoNormalBasis;
import org.sunflow.math.Point3;
import org.sunflow.math.Solvers;
import org.sunflow.math.Vector3;

public class SphereLight implements LightSource, Shader {
    private Color radiance;
    private int numSamples;
    private Point3 center;
    private float radius;
    private float r2;

    public SphereLight() {
        radiance = Color.WHITE;
        numSamples = 4;
        center = new Point3();
        radius = r2 = 1;
    }

    public boolean update(ParameterList pl, SunflowAPI api) {
        radiance = pl.getColor("radiance", radiance);
        numSamples = pl.getInt("samples", numSamples);
        radius = pl.getFloat("radius", radius);
        r2 = radius * radius;
        center = pl.getPoint("center", center);
        return true;
    }

    public void init(String name, SunflowAPI api) {
        api.light(name, this);
        api.geometry(name + ".geo", new Sphere());
        api.shader(name + ".shader", this);
        api.parameter("shaders", name + ".shader");
        api.parameter("transform", Matrix4.translation(center.x, center.y, center.z).multiply(Matrix4.scale(radius)));
        api.instance(name + ".instance", name + ".geo");
    }

    public int getNumSamples() {
        return numSamples;
    }

    public int getLowSamples() {
        return 1;
    }

    public boolean isVisible(ShadingState state) {
        return state.getPoint().distanceToSquared(center) > r2;
    }

    public void getSamples(ShadingState state) {
        if (getNumSamples() <= 0)
            return;
        Vector3 wc = Point3.sub(center, state.getPoint(), new Vector3());
        float l2 = wc.lengthSquared();
        if (l2 <= r2)
            return; // inside the sphere?
        // top of the sphere as viewed from the current shading point
        float topX = wc.x + state.getNormal().x * radius;
        float topY = wc.y + state.getNormal().y * radius;
        float topZ = wc.z + state.getNormal().z * radius;
        if (state.getNormal().dot(topX, topY, topZ) <= 0)
            return; // top of the sphere is below the horizon
        float cosThetaMax = (float) Math.sqrt(Math.max(0, 1 - r2 / Vector3.dot(wc, wc)));
        OrthoNormalBasis basis = OrthoNormalBasis.makeFromW(wc);
        int samples = state.getDiffuseDepth() > 0 ? 1 : getNumSamples();
        float scale = (float) (2 * Math.PI * (1 - cosThetaMax));
        Color c = Color.mul(scale / samples, radiance);
        for (int i = 0; i < samples; i++) {
            // random offset on unit square
            double randX = state.getRandom(i, 0, samples);
            double randY = state.getRandom(i, 1, samples);

            // cone sampling
            double cosTheta = (1 - randX) * cosThetaMax + randX;
            double sinTheta = Math.sqrt(1 - cosTheta * cosTheta);
            double phi = randY * 2 * Math.PI;
            Vector3 dir = new Vector3((float) (Math.cos(phi) * sinTheta), (float) (Math.sin(phi) * sinTheta), (float) cosTheta);
            basis.transform(dir);

            // check that the direction of the sample is the same as the
            // normal
            float cosNx = Vector3.dot(dir, state.getNormal());
            if (cosNx <= 0)
                continue;

            float ocx = state.getPoint().x - center.x;
            float ocy = state.getPoint().y - center.y;
            float ocz = state.getPoint().z - center.z;
            float qa = Vector3.dot(dir, dir);
            float qb = 2 * ((dir.x * ocx) + (dir.y * ocy) + (dir.z * ocz));
            float qc = ((ocx * ocx) + (ocy * ocy) + (ocz * ocz)) - r2;
            double[] t = Solvers.solveQuadric(qa, qb, qc);
            if (t == null)
                continue;
            LightSample dest = new LightSample();
            // compute shadow ray to the sampled point
            dest.setShadowRay(new Ray(state.getPoint(), dir));
            // FIXME: arbitrary bias, should handle as in other places
            dest.getShadowRay().setMax((float) t[0] - 1e-3f);
            // prepare sample
            dest.setRadiance(c, c);
            dest.traceShadow(state);
            state.addSample(dest);
        }
    }

    public void getPhoton(double randX1, double randY1, double randX2, double randY2, Point3 p, Vector3 dir, Color power) {
        float z = (float) (1 - 2 * randX2);
        float r = (float) Math.sqrt(Math.max(0, 1 - z * z));
        float phi = (float) (2 * Math.PI * randY2);
        float x = r * (float) Math.cos(phi);
        float y = r * (float) Math.sin(phi);
        p.x = center.x + x * radius;
        p.y = center.y + y * radius;
        p.z = center.z + z * radius;
        OrthoNormalBasis basis = OrthoNormalBasis.makeFromW(new Vector3(x, y, z));
        phi = (float) (2 * Math.PI * randX1);
        float cosPhi = (float) Math.cos(phi);
        float sinPhi = (float) Math.sin(phi);
        float sinTheta = (float) Math.sqrt(randY1);
        float cosTheta = (float) Math.sqrt(1 - randY1);
        dir.x = cosPhi * sinTheta;
        dir.y = sinPhi * sinTheta;
        dir.z = cosTheta;
        basis.transform(dir);
        power.set(radiance);
        power.mul((float) (Math.PI * Math.PI * 4 * r2));
    }

    public float getPower() {
        return radiance.copy().mul((float) (Math.PI * Math.PI * 4 * r2)).getLuminance();
    }

    public Color getRadiance(ShadingState state) {
        if (!state.includeLights())
            return Color.BLACK;
        state.faceforward();
        // emit constant radiance
        return state.isBehind() ? Color.BLACK : radiance;
    }

    public void scatterPhoton(ShadingState state, Color power) {
        // do not scatter photons
    }
}