朗伯着色和阴影问题（光线追踪）

Question

我做了一个ray tracer:

我一直在实现 Blinn-Phong 算法，但是在特定情况下渲染时遇到了一些问题(otherwise it works quite well)

• 当球体悬停在空中时，朗伯着色无法正常工作（似乎仅在球体高于光源位置时才会发生）-Screenshot

• 当光源低于球体时 - Screenshot

• 当阴影投射到球体上时（也可以在第一个示例中看到）-Screenshot

• 当场景中有超过 1 个平面时（仅当阴影打开时）-Screenshot

我假设我的阴影和朗伯着色算法有问题，但我不确定是什么问题，有人可以帮忙吗？

代码（没有外部库，直接编译）：

Color finalColor;
Color ambient;
Color diffuse;
FPType lambertian;
FPType phong;
Color specular;
bool shadowed = false;

// Ambient
if(AMBIENT_ON)
{
    ambient =  closestObjectMaterial.GetColor() * AMBIENT_LIGHT * closestObjectMaterial.GetAmbient();
    finalColor += ambient;
}

// Shadows, Diffuse, Specular
for(const auto &lightSource : lightSources)
{
    Vector lightDir = (lightSource->GetPosition() - point); // Calculate the directional vector towards the lightSource
    FPType distance = lightDir.Magnitude();
    lightDir = lightDir.Normalize();
    lambertian = closestObjectNormal.Dot(lightDir);

    // Shadows
    if(SHADOWS_ON && lambertian > 0)
    {
        Ray shadowRay(point, lightDir); // Cast a ray from the first intersection to the light

        std::vector<FPType> secondaryIntersections;
        for(const auto &sceneObject : sceneObjects)
        {
            secondaryIntersections.push_back(sceneObject->GetIntersection(shadowRay));
        }

        for(const auto &secondaryIntersection : secondaryIntersections)
        {
            if(secondaryIntersection > TOLERANCE) // If shadow ray intersects with some object along the way
            {
                shadowed = true;
                finalColor *= closestObjectMaterial.GetDiffuse() * AMBIENT_LIGHT;
                break;
            }
        }
    }

    // Diffuse
    if(DIFFUSE_ON && shadowed == false)
    {
        diffuse = closestObjectMaterial.GetColor().Average(lightSource->GetColor()) * closestObjectMaterial.GetDiffuse() * lightSource->GetIntensity() * std::fmax(lambertian, 0) / distance;
        finalColor += diffuse;
    }
    if(shadowed == false && SPECULAR_ON)
    {
        // Specular
        if(closestObjectMaterial.GetSpecular() > 0 && closestObjectMaterial.GetSpecular() <= 1)
        {
            Vector V = -sceneDirection;
            // Blinn-Phong
            Vector H = (lightDir + V).Normalize();
            FPType NdotH = closestObjectNormal.Dot(H);

            phong = pow(NdotH, 300);
            specular = lightSource->GetColor() * std::fmax(0, phong) * lightSource->GetIntensity() / distance; // closestObjectMaterial.GetSpecular(); add or no?
            finalColor += specular;
        }
    }
}

FPType Sphere::GetIntersection(Ray ray)
{
    Vector length = center - ray.GetOrigin(); // Length of the vector between the center and the ray origin (hypotenuse)
    FPType tca = length.Dot(ray.GetDirection()); // opposide side

    if(tca < 0) // No intersection registered
        return -1;

if(tca > 0) // Intersection registered
{
    FPType a = sqrt(length.Dot(length) - tca*tca); // Adjacent side (a = sqrt(c²-b²))

    if(a > radius || a < 0)
        return -1;

    FPType thc = sqrt(radius*radius - a*a); // the line between 2 intersection points / 2

    FPType primaryIntersection;
    primaryIntersection = tca - thc;
    if(primaryIntersection > 0)
        return primaryIntersection;
    else
    {
        FPType secondaryIntersection = thc + tca;
        return secondaryIntersection;
    }
}
return -1;

}

• Lines 97-132

• To set the scene up

Answer 1

我不确定为什么您的 Sphere-Ray 交集有两个平方根，您想根据 this 测试我的实现吗？

FPType Sphere::GetIntersection(Ray ray)
{
  Vector delta = ray.GetOrigin() - center;
  Vector dir = ray.GetDirection();

  //Quadratic equation describing the distance along ray to intersection
  FPType a = dir.Dot(dir);
  FPType b = dir.Dot(delta); //removed factor of 2 later divide by a, NOT 2a
  FPType c = delta.Dot(delta) - radius*radius;

  FPType descrim = b*b - a*c;
  if (descrim < FPType(0)) {
    return -1;
  }
  //Find solutions to quadratic equation
  descrim = sqrt(descrim) / a;
  b = -b / a;

  FPType intersection0 = b - descrim
  if(intersection0 >= FPType(0)) {
    return intersection0;
  }

  FPType intersection1 = b + descrim
  if(intersection1 >= FPType(0)) {
    return intersection1;
  }

  //Both solutions were negative
  return -1;
}