/// <summary>
/// Die Gestalt wird nach hinten abgedunkelt.
/// </summary>
protected void DarkenPlane()
{
/* Testweise grau */
float mainDeph = _maxY - _minY;// borderMaxY - borderMinY;
for (int i = 0; i < pData.Width; i++)
{
for (int j = 0; j < pData.Height; j++)
{
FloatVec3 col = _rgbPlane[i, j];
float yd = _heightMap[i, j];
if (yd != float.MinValue)
{
float ydNormalized = (yd - _minY) / mainDeph;
float greyYd = 1.0f - ydNormalized;
greyYd = greyYd * greyYd;
col.X = col.X * ydNormalized + (_backColorRed * greyYd);
col.Y = col.Y * ydNormalized + (_backColorGreen * greyYd);
col.Z = ydNormalized * col.Z + (_backColorBlue * greyYd);
}
else
{
col.X = _backColorRed;
col.Y = _backColorGreen;
col.Z = _backColorBlue;
}
}
}
}
/// <summary>
/// Surface color according to _shininessFactor and _shininess.
/// </summary>
protected virtual FloatVec3 GetLightF(FloatVec3 normal)
{
FloatVec3 retVal = new FloatVec3((float)_backColorRed, (float)_backColorGreen, (float)_backColorBlue);
if (!_useLight)
{
return(new FloatVec3((float)0.5, (float)0.5, (float)0.5));
}
if (normal == null)
{
return(retVal);
}
float weight_shini = (float)_shininessFactor;
float weight_diffuse = 1 - weight_shini;
float norm = (float)Math.Sqrt(normal.X * normal.X + normal.Y * normal.Y + normal.Z * normal.Z);
// Scalar product with (0,-1,0) light ray
float angle = 0;
if (norm == 0)
{
return(retVal);
}
FloatVec3 lightVec = new FloatVec3((float)_lightRay.X, (float)_lightRay.Y, (float)_lightRay.Z);
lightVec.Normalize();
float norm2 = lightVec.Norm;
angle = (float)(Math.Acos((normal.X * lightVec.X + normal.Y * lightVec.Y + normal.Z * lightVec.Z) / (norm * norm2)) / (Math.PI / 2.0));
angle = 1 - angle;
if (angle < 0)
{
angle = 0;
}
if (angle > 1)
{
angle = 1;
}
float light = (float)(weight_diffuse * angle + weight_shini * Math.Pow(angle, _shininess));
if (light < 0)
{
light = 0;
}
if (light > 1)
{
light = 1;
}
retVal.X = light;
retVal.Y = light;
retVal.Z = light;
return(retVal);
}
/// <summary>
/// Create rgb image.
/// </summary>
protected void DrawPlane()
{
_rgbPlane = new FloatVec3[pData.Width, pData.Height];
for (int i = 0; i < pData.Width; i++)
{
for (int j = 0; j < pData.Height; j++)
{
Vec3 rgb = GetRgb(i, j);
if (double.IsNaN(rgb.X) || double.IsNaN(rgb.Y) || double.IsNaN(rgb.Z))
{
}
_rgbPlane[i, j] = new FloatVec3((float)rgb.X, (float)rgb.Y, (float)rgb.Z);
}
}
}
/// <summary>
/// Bildpunkte, die auf Grund fehlender Informationen nicht geladen werden konnten, werden
/// aus den Umgebungsinformationen gemittelt.
/// </summary>
protected void SmoothEmptyPixel()
{
for (int i = 0; i < pData.Width; i++)
{
for (int j = 0; j < pData.Height; j++)
{
PixelInfo pInfo = pData.Points[i, j];
if (pInfo == null)
{
// Dieser Wert ist zu setzen
// Aber nur, wenn es sich nicht um den Hintergrund handelt.
FloatVec3 col = _rgbPlane[i, j];
col.X = _backColorRed;
col.Y = _backColorGreen;
col.Z = _backColorBlue;
float pixelCount = 0;
for (int k = i - 1; k <= i + 1; k++)
{
for (int l = j - 1; l <= j + 1; l++)
{
if (k >= 0 && k < pData.Width && l >= 0 && l < pData.Height && k != i && l != j)
{
PixelInfo pInfo2 = pData.Points[k, l];
if (pInfo2 != null)
{
pixelCount++;
FloatVec3 otherColor = _rgbPlane[k, l];
col.X += otherColor.X;
col.Y += otherColor.Y;
col.Z += otherColor.Z;
}
}
}
}
//pixelCount++; // Etwas dunkler sollte es schon werden
if (pixelCount > 1)
{
col.X /= pixelCount;
col.Y /= pixelCount;
col.Z /= pixelCount;
_isBackground[i, j] = false;
}
}
}
}
}
/// <summary>
/// Creates boundingbox infos.
/// </summary>
protected void CreateStatisticInfo()
{
_minPoint.X = float.MaxValue;
_minPoint.Y = float.MaxValue;
_minPoint.Z = float.MaxValue;
_maxPoint.X = float.MinValue;
_maxPoint.Y = float.MinValue;
_maxPoint.Z = float.MinValue;
for (int i = 0; i < _pictureData.Width; i++)
{
for (int j = 0; j < _pictureData.Height; j++)
{
FloatPixelInfo pInfo = _pictureData.Points[i, j];
if (pInfo != null)
{
FloatVec3 coord = pInfo.Coord;
if (coord.X < _minPoint.X)
{
_minPoint.X = coord.X;
}
if (coord.Y < _minPoint.Y)
{
_minPoint.Y = coord.Y;
}
if (coord.Z < _minPoint.Z)
{
_minPoint.Z = coord.Z;
}
if (coord.X > _maxPoint.X)
{
_maxPoint.X = coord.X;
}
if (coord.Y > _maxPoint.Y)
{
_maxPoint.Y = coord.Y;
}
if (coord.Z > _maxPoint.Z)
{
_maxPoint.Z = coord.Z;
}
}
}
}
}
/// <summary>
/// Add point other.
/// </summary>
public void Add(FloatVec3 other)
{
X += other.X;
Y += other.Y;
Z += other.Z;
}
public FloatVec3(FloatVec3 other)
{
this.X = other.X;
this.Y = other.Y;
this.Z = other.Z;
}
/// <summary>
/// + Operator
/// </summary>
public FloatVec3 Sum(FloatVec3 other)
{
return new FloatVec3(X + other.X, Y + other.Y, Z + other.Z);
}
/// <summary>
/// Subtract componentwise.
/// </summary>
public void Subtract(FloatVec3 other)
{
X -= other.X;
Y -= other.Y;
Z -= other.Z;
}
/// <summary>
/// Compute field of view.
/// </summary>
protected void SmoothPlane()
{
_rgbSmoothPlane1 = new FloatVec3[pData.Width, pData.Height];
_rgbSmoothPlane2 = new FloatVec3[pData.Width, pData.Height];
int intRange = 3;
for (int i = 0; i < pData.Width; i++)
{
for (int j = 0; j < pData.Height; j++)
{
_rgbSmoothPlane2[i, j] = _rgbPlane[i, j];
}
}
if (_ambientIntensity == 0)
{
// No field of view defined:
return;
}
// starts with rgbSmoothPlane2
FloatVec3[,] currentPlane = _rgbSmoothPlane2;
FloatVec3[,] nextPlane = _rgbSmoothPlane1;
// contain the result colors
FloatVec3[,] resultPlane = _rgbSmoothPlane1;
for (int m = 0; m < _ambientIntensity; m++)
{
if (_stopRequest)
{
return;
}
for (int i = 0; i < pData.Width; i++)
{
for (int j = 0; j < pData.Height; j++)
{
double neighborsFound = 0;
FloatPixelInfo pInfo = _pictureData.Points[i, j];
FloatVec3 nColor = new FloatVec3();
float ydNormalized = (float)GetAmbientValue(_heightMap[i, j]);
ydNormalized = (float)Math.Sqrt(ydNormalized);
intRange = 1;
if (intRange == 0)
{
nColor = currentPlane[i, j];
neighborsFound = 1;
}
float sumColor = 0;
// if(pData.Points[i, j]!=null) {
//if (true)
{
for (int k = -intRange; k <= intRange; k++)
{
for (int l = -intRange; l <= intRange; l++)
{
// Center Pixel is also allowed
// if (k != 0 || l != 0) {
int posX = i + k;
int posY = j + l;
if (posX >= 0 && posX < pData.Width && posY >= 0 && posY < pData.Height)
{
FloatVec3 nColor1 = new FloatVec3();
if (true)
// if ( (ylocalDiff > 0) ||(i==posX && j==posY))
// if ((ylocalDiff < 0) || (i == posX && j == posY))
// if(false)
{
//double range = (k * k + l * l) / (intRange * intRange);
int range = (k * k + l * l);
float mult = 1;
if (range == 0)
{
// mult = 0.6;
mult = ydNormalized * 0.3f;
//mult = 0.2;
}
if (range == 1)
{
//mult = 0.25;
mult = (1.0f - ydNormalized) * 0.4f;
//mult = 0.45;
}
if (range == 2)
{
//mult=0.15;
mult = (1.0f - ydNormalized) * 0.3f;
//mult = 0.35;
}
// mult += 0.00001;
FloatPixelInfo pInfo2 = _pictureData.Points[posX, posY];
float amount = 1;
if (pInfo != null && pInfo2 != null)
{
float dist = pInfo.Coord.Dist(pInfo2.Coord);
float dGlobal = _maxPoint.Dist(_minPoint);
dGlobal /= 1500;
if (dist < dGlobal)
{
amount = 1.0f;
}
else
{
// else if (dist > dGlobal && dist < 10.0 * dGlobal)
amount = 1.0f - (dGlobal / dist) / 10.0f;
}
// else
// amount = 0.0;
}
mult *= amount;
// mult *= 1.0/ydNormalized;
sumColor += mult;
FloatVec3 currentColor = currentPlane[posX, posY];
nColor1.X = currentColor.X;
nColor1.Y = currentColor.Y;
nColor1.Z = currentColor.Z;
nColor1 = nColor1.Mult(mult); // Scaling;
nColor.Add(nColor1);
neighborsFound++;
}
}
}
}
}
if (neighborsFound > 1)
{
nColor = nColor.Mult(1 / sumColor);
}
else
{
nColor = currentPlane[i, j];
}
nextPlane[i, j] = nColor;
}
}
resultPlane = nextPlane;
nextPlane = currentPlane;
currentPlane = resultPlane;
}
for (int i = 0; i < pData.Width; i++)
{
for (int j = 0; j < pData.Height; j++)
{
_rgbPlane[i, j] = resultPlane[i, j];
}
}
_rgbSmoothPlane1 = null;
_rgbSmoothPlane2 = null;
return;
}
/// <summary>
/// Return path to other point as vector.
/// </summary>
public FloatVec3 Diff(FloatVec3 other)
{
return(new FloatVec3(this.X - other.X, Y - other.Y, Z - other.Z));
}
/// <summary>
/// Subtract componentwise.
/// </summary>
public void Subtract(FloatVec3 other)
{
X -= other.X;
Y -= other.Y;
Z -= other.Z;
}
/// <summary>
/// Surface color according to _shininessFactor and _shininess.
/// </summary>
protected virtual FloatVec3 GetLightF(FloatVec3 normal)
{
FloatVec3 retVal = new FloatVec3((float)_backColorRed, (float)_backColorGreen, (float)_backColorBlue);
if (!_useLight)
{
return new FloatVec3((float)0.5, (float)0.5, (float)0.5);
}
if (normal == null)
return retVal;
float weight_shini = (float)_shininessFactor;
float weight_diffuse = 1 - weight_shini;
float norm = (float)Math.Sqrt(normal.X * normal.X + normal.Y * normal.Y + normal.Z * normal.Z);
// Scalar product with (0,-1,0) light ray
float angle = 0;
if (norm == 0)
return retVal;
FloatVec3 lightVec = new FloatVec3((float)_lightRay.X, (float)_lightRay.Y, (float)_lightRay.Z);
lightVec.Normalize();
float norm2 = lightVec.Norm;
angle = (float)(Math.Acos((normal.X * lightVec.X + normal.Y * lightVec.Y + normal.Z * lightVec.Z) / (norm * norm2)) / (Math.PI / 2.0));
angle = 1 - angle;
if (angle < 0)
angle = 0;
if (angle > 1)
angle = 1;
float light = (float)(weight_diffuse * angle + weight_shini * Math.Pow(angle, _shininess));
if (light < 0)
light = 0;
if (light > 1)
light = 1;
retVal.X = light;
retVal.Y = light;
retVal.Z = light;
return retVal;
}
/// <summary>
/// Compute field of view.
/// </summary>
protected void SmoothPlane()
{
_rgbSmoothPlane1 = new FloatVec3[pData.Width, pData.Height];
_rgbSmoothPlane2 = new FloatVec3[pData.Width, pData.Height];
int intRange = 3;
for (int i = 0; i < pData.Width; i++)
{
for (int j = 0; j < pData.Height; j++)
{
_rgbSmoothPlane2[i, j] = _rgbPlane[i, j];
}
}
if (_ambientIntensity == 0)
{
// No field of view defined:
return;
}
// starts with rgbSmoothPlane2
FloatVec3[,] currentPlane = _rgbSmoothPlane2;
FloatVec3[,] nextPlane = _rgbSmoothPlane1;
// contain the result colors
FloatVec3[,] resultPlane = _rgbSmoothPlane1;
for (int m = 0; m < _ambientIntensity; m++)
{
if (_stopRequest)
return;
for (int i = 0; i < pData.Width; i++)
{
for (int j = 0; j < pData.Height; j++)
{
double neighborsFound = 0;
FloatPixelInfo pInfo = _pictureData.Points[i, j];
FloatVec3 nColor = new FloatVec3();
float ydNormalized = (float)GetAmbientValue(_heightMap[i, j]);
ydNormalized = (float)Math.Sqrt(ydNormalized);
intRange = 1;
if (intRange == 0)
{
nColor = currentPlane[i, j];
neighborsFound = 1;
}
float sumColor = 0;
// if(pData.Points[i, j]!=null) {
//if (true)
{
for (int k = -intRange; k <= intRange; k++)
{
for (int l = -intRange; l <= intRange; l++)
{
// Center Pixel is also allowed
// if (k != 0 || l != 0) {
int posX = i + k;
int posY = j + l;
if (posX >= 0 && posX < pData.Width && posY >= 0 && posY < pData.Height)
{
FloatVec3 nColor1 = new FloatVec3();
if (true)
// if ( (ylocalDiff > 0) ||(i==posX && j==posY))
// if ((ylocalDiff < 0) || (i == posX && j == posY))
// if(false)
{
//double range = (k * k + l * l) / (intRange * intRange);
int range = (k * k + l * l);
float mult = 1;
if (range == 0)
{
// mult = 0.6;
mult = ydNormalized * 0.3f;
//mult = 0.2;
}
if (range == 1)
{
//mult = 0.25;
mult = (1.0f - ydNormalized) * 0.4f;
//mult = 0.45;
}
if (range == 2)
{
//mult=0.15;
mult = (1.0f - ydNormalized) * 0.3f;
//mult = 0.35;
}
// mult += 0.00001;
FloatPixelInfo pInfo2 = _pictureData.Points[posX, posY];
float amount = 1;
if (pInfo != null && pInfo2 != null)
{
float dist = pInfo.Coord.Dist(pInfo2.Coord);
float dGlobal = _maxPoint.Dist(_minPoint);
dGlobal /= 1500;
if (dist < dGlobal)
amount = 1.0f;
else
// else if (dist > dGlobal && dist < 10.0 * dGlobal)
amount = 1.0f - (dGlobal / dist) / 10.0f;
// else
// amount = 0.0;
}
mult *= amount;
// mult *= 1.0/ydNormalized;
sumColor += mult;
FloatVec3 currentColor = currentPlane[posX, posY];
nColor1.X = currentColor.X;
nColor1.Y = currentColor.Y;
nColor1.Z = currentColor.Z;
nColor1 = nColor1.Mult(mult); // Scaling;
nColor.Add(nColor1);
neighborsFound++;
}
}
}
}
}
if (neighborsFound > 1)
{
nColor = nColor.Mult(1 / sumColor);
}
else
{
nColor = currentPlane[i, j];
}
nextPlane[i, j] = nColor;
}
}
resultPlane = nextPlane;
nextPlane = currentPlane;
currentPlane = resultPlane;
}
for (int i = 0; i < pData.Width; i++)
{
for (int j = 0; j < pData.Height; j++)
{
_rgbPlane[i, j] = resultPlane[i, j];
}
}
_rgbSmoothPlane1 = null;
_rgbSmoothPlane2 = null;
return;
}
protected override Vec3 GetRgbAt(int x, int y)
{
FloatVec3 rgb = _rgbPlane[x, y];
return(new Vec3(rgb.X, rgb.Y, rgb.Z));
}
/// <summary>
/// Weißabgleich und Helligkeitskorrektur.
/// </summary>
protected void NormalizePlane()
{
float maxRed = 0;
float maxGreen = 0;
float maxBlue = 0;
for (int i = 1; i < pData.Width - 1; i++)
{
for (int j = 1; j < pData.Height - 1; j++)
{
if ((_picInfo[i, j] == 0) && (_picInfo[i + 1, j] == 0) && (_picInfo[i - 1, j] == 0) && (_picInfo[i, j - 1] == 0) &&
(_picInfo[i, j + 1] == 0) && (_picInfo[i - 1, j - 1] == 0) && (_picInfo[i - 1, j + 1] == 0) && (_picInfo[i + 1, j - 1] == 0) &&
(_picInfo[i - 1, j + 1] == 0))
{
FloatVec3 col = _rgbPlane[i, j];
if (col.X > maxRed)
{
maxRed = col.X;
}
if (col.Y > maxGreen)
{
maxGreen = col.Y;
}
if (col.Z > maxBlue)
{
maxBlue = col.Z;
}
}
}
}
for (int i = 0; i < pData.Width; i++)
{
for (int j = 0; j < pData.Height; j++)
{
FloatVec3 col = _rgbPlane[i, j];
if (_picInfo[i, j] == 0)
{
if (maxRed > 0)
{
col.X /= maxRed;
}
if (maxGreen > 0)
{
col.Y /= maxGreen;
}
if (maxBlue > 0)
{
col.Z /= maxBlue;
}
}
if (col.X < 0)
{
col.X = 0;
}
if (col.X > 1)
{
col.X = 1;
}
if (col.Y < 0)
{
col.Y = 0;
}
if (col.Y > 1)
{
col.Y = 1;
}
if (col.Z < 0)
{
col.Z = 0;
}
if (col.Z > 1)
{
col.Z = 1;
}
}
}
}
/// <summary>
/// + Operator
/// </summary>
public FloatVec3 Sum(FloatVec3 other)
{
return(new FloatVec3(X + other.X, Y + other.Y, Z + other.Z));
}
/// <summary>
/// Return path to other point as vector.
/// </summary>
public FloatVec3 Diff(FloatVec3 other)
{
return (new FloatVec3(this.X - other.X, Y - other.Y, Z - other.Z));
}
/// <summary>
/// Return distance to other point.
/// </summary>
public float Dist(FloatVec3 other)
{
return(Diff(other).Norm);
}
/// <summary>
/// Return distance to other point.
/// </summary>
public float Dist(FloatVec3 other)
{
return Diff(other).Norm;
}
/// <summary>
/// Add point other.
/// </summary>
public void Add(FloatVec3 other)
{
X += other.X;
Y += other.Y;
Z += other.Z;
}
/// <summary>
/// Create rgb image.
/// </summary>
protected void DrawPlane()
{
_rgbPlane = new FloatVec3[pData.Width, pData.Height];
for (int i = 0; i < pData.Width; i++)
{
for (int j = 0; j < pData.Height; j++)
{
Vec3 rgb = GetRgb(i, j);
if(double.IsNaN(rgb.X)|| double.IsNaN(rgb.Y)|| double.IsNaN(rgb.Z))
{
}
_rgbPlane[i, j] = new FloatVec3((float)rgb.X, (float)rgb.Y, (float)rgb.Z);
}
}
}
/// <summary>
/// Get the color information of the bitmap at (x,y)
/// </summary>
protected Vec3 GetRgb(int x, int y)
{
FloatVec3 retVal = new FloatVec3(0, 0, 1); // blau
FloatPixelInfo pInfo = _pictureData.Points[x, y];
if (pInfo == null)
{
return(new Vec3(_backColorRed, _backColorGreen, _backColorBlue));
}
FloatVec3 light = new FloatVec3(0, 0, 0);
FloatVec3 normal = null;
normal = pInfo.Normal;
if (normal == null)
{
return(new Vec3(0, 0, 0));
}
normal.Normalize();
if (pInfo.Normal != null)
{
light = GetLightF(normal);
}
retVal.X = light.X;
retVal.Y = light.Y;
retVal.Z = light.Z;
retVal.X = (float)(_lightIntensity * retVal.X + (1 - _lightIntensity) * (1 - _shadowPlane[x, y]));
retVal.Y = (float)(_lightIntensity * retVal.Y + (1 - _lightIntensity) * (1 - _shadowPlane[x, y]));
retVal.Z = (float)(_lightIntensity * retVal.Z + (1 - _lightIntensity) * (1 - _shadowPlane[x, y]));
if (retVal.X < 0)
{
retVal.X = 0;
}
if (retVal.Y < 0)
{
retVal.Y = 0;
}
if (retVal.Z < 0)
{
retVal.Z = 0;
}
if (retVal.X > 1)
{
retVal.X = 1;
}
if (retVal.Y > 1)
{
retVal.Y = 1;
}
if (retVal.Z > 1)
{
retVal.Z = 1;
}
float brightLightLevel = (float)_parameters.GetDouble("Renderer.BrightLightLevel");
if (brightLightLevel > 0)
{
retVal.X = (1 - brightLightLevel) * retVal.X + brightLightLevel * light.X * (1 - _shadowPlane[x, y]);
retVal.Y = (1 - brightLightLevel) * retVal.Y + brightLightLevel * light.Y * (1 - _shadowPlane[x, y]);
retVal.Z = (1 - brightLightLevel) * retVal.Z + brightLightLevel * light.Z * (1 - _shadowPlane[x, y]);
}
if (retVal.X < 0)
{
retVal.X = 0;
}
if (retVal.Y < 0)
{
retVal.Y = 0;
}
if (retVal.Z < 0)
{
retVal.Z = 0;
}
if (retVal.X > 1)
{
retVal.X = 1;
}
if (retVal.Y > 1)
{
retVal.Y = 1;
}
if (retVal.Z > 1)
{
retVal.Z = 1;
}
// Add surface color
bool useAdditionalColorinfo = true;
if (_colorIntensity <= 0)
{
useAdditionalColorinfo = false;
}
if (useAdditionalColorinfo && ((pInfo.IsInside && _colorInside) || (!pInfo.IsInside && _colorOutside)))
{
if (pInfo != null && pInfo.AdditionalInfo != null)
{
if (double.IsNaN(pInfo.AdditionalInfo.red))
{
pInfo.AdditionalInfo.red = 0;
}
if (double.IsNaN(pInfo.AdditionalInfo.green))
{
pInfo.AdditionalInfo.green = 0;
}
if (double.IsNaN(pInfo.AdditionalInfo.blue))
{
pInfo.AdditionalInfo.blue = 0;
}
if (pInfo.AdditionalInfo.red < 0)
{
pInfo.AdditionalInfo.green -= pInfo.AdditionalInfo.red;
pInfo.AdditionalInfo.blue -= pInfo.AdditionalInfo.red;
pInfo.AdditionalInfo.red = 0;
}
if (pInfo.AdditionalInfo.green < 0)
{
pInfo.AdditionalInfo.red -= pInfo.AdditionalInfo.green;
pInfo.AdditionalInfo.blue -= pInfo.AdditionalInfo.green;
pInfo.AdditionalInfo.green = 0;
}
if (pInfo.AdditionalInfo.blue < 0)
{
pInfo.AdditionalInfo.red -= pInfo.AdditionalInfo.blue;
pInfo.AdditionalInfo.blue -= pInfo.AdditionalInfo.blue;
pInfo.AdditionalInfo.blue = 0;
}
// Normalise;
float r1 = (float)(_colorFactorRed * Math.Pow(pInfo.AdditionalInfo.red, _colorIntensity));
float g1 = (float)(_colorFactorGreen * Math.Pow(pInfo.AdditionalInfo.green, _colorIntensity));
float b1 = (float)(_colorFactorBlue * Math.Pow(pInfo.AdditionalInfo.blue, _colorIntensity));
if (r1 < 0)
{
r1 = -r1;
}
if (g1 < 0)
{
g1 = -g1;
}
if (b1 < 0)
{
b1 = -b1;
}
// Normalize:
float norm = (float)(Math.Sqrt(r1 * r1 + g1 * g1 + b1 * b1) / Math.Sqrt(2.5));
if (norm > 0)
{
r1 = r1 / norm;
g1 = g1 / norm;
b1 = b1 / norm;
if (_colorThreshold > 0)
{
double thresholdIndicator = Math.Max(Math.Abs(r1 - b1), Math.Max(Math.Abs(r1 - g1), Math.Abs(g1 - b1)));
float lowerThresholdFactor = (float)0.7;
if (thresholdIndicator < _colorThreshold * lowerThresholdFactor)
{
r1 = (float)1;
g1 = (float)1;
b1 = (float)1;
}
else if (thresholdIndicator < _colorThreshold)
{
r1 = (float)(_colorThreshold - thresholdIndicator) * (float)0.5 / lowerThresholdFactor + r1;
g1 = (float)(_colorThreshold - thresholdIndicator) * (float)0.5 / lowerThresholdFactor + g1;
b1 = (float)(_colorThreshold - thresholdIndicator) * (float)0.5 / lowerThresholdFactor + b1;
}
}
}
for (int i = 0; i < 5; i++)
{
if (r1 > 1)
{
b1 += (r1 - 1) / (float)2.0;
g1 += (r1 - 1) / (float)2.0;
r1 = 1;
}
if (b1 > 1)
{
r1 += (b1 - 1) / (float)2.0;
g1 += (b1 - 1) / (float)2.0;
b1 = 1;
}
if (g1 > 1)
{
r1 += (g1 - 1) / (float)2.0;
b1 += (g1 - 1) / (float)2.0;
g1 = 1;
}
}
if (r1 > 1)
{
r1 = 1;
}
if (b1 > 1)
{
b1 = 1;
}
if (g1 > 1)
{
g1 = 1;
}
if (_colorGreyness > 0)
{
r1 = (float)(_colorGreyness + (1 - _colorGreyness) * r1);
g1 = (float)(_colorGreyness + (1 - _colorGreyness) * g1);
b1 = (float)(_colorGreyness + (1 - _colorGreyness) * b1);
}
if (r1 > 1)
{
r1 = 1;
}
if (b1 > 1)
{
b1 = 1;
}
if (g1 > 1)
{
g1 = 1;
}
if (norm != 0)
{
switch (_rgbType)
{
case 1:
retVal.X *= r1;
retVal.Y *= g1;
retVal.Z *= b1;
break;
case 2:
retVal.X *= r1;
retVal.Y *= b1;
retVal.Z *= g1;
break;
case 3:
retVal.X *= g1;
retVal.Y *= r1;
retVal.Z *= b1;
break;
case 4:
retVal.X *= g1;
retVal.Y *= b1;
retVal.Z *= r1;
break;
case 5:
retVal.X *= b1;
retVal.Y *= r1;
retVal.Z *= g1;
break;
case 6:
retVal.X *= b1;
retVal.Y *= g1;
retVal.Z *= r1;
break;
default:
retVal.X *= r1;
retVal.Y *= g1;
retVal.Z *= b1;
break;
}
}
}
}
if (_contrast != 1)
{
retVal.X = (float)Math.Pow(retVal.X, _contrast);
retVal.Y = (float)Math.Pow(retVal.Y, _contrast);
retVal.Z = (float)Math.Pow(retVal.Z, _contrast);
}
if (_brightness > 1)
{
retVal.X *= _brightness;
retVal.Y *= _brightness;
retVal.Z *= _brightness;
}
if (retVal.X < 0)
{
retVal.X = 0;
}
if (retVal.X > 1)
{
retVal.X = 1;
}
if (retVal.Y < 0)
{
retVal.Y = 0;
}
if (retVal.Z < 0)
{
retVal.Z = 0;
}
if (retVal.Y > 1)
{
retVal.Y = 1;
}
if (retVal.Z > 1)
{
retVal.Z = 1;
}
if (pInfo != null && pInfo.AdditionalInfo != null)
{
pInfo.AdditionalInfo.red2 = retVal.X;
pInfo.AdditionalInfo.green2 = retVal.Y;
pInfo.AdditionalInfo.blue2 = retVal.Z;
}
return(new Vec3(retVal.X, retVal.Y, retVal.Z));
}
/// <summary>
/// Get the color information of the bitmap at (x,y)
/// </summary>
protected Vec3 GetRgb(int x, int y)
{
FloatVec3 retVal = new FloatVec3(0, 0, 1); // blau
FloatPixelInfo pInfo = _pictureData.Points[x, y];
if (pInfo == null)
{
return new Vec3(_backColorRed, _backColorGreen, _backColorBlue);
}
FloatVec3 light = new FloatVec3(0, 0, 0);
FloatVec3 normal = null;
normal = pInfo.Normal;
if (normal == null)
return new Vec3(0, 0, 0);
normal.Normalize();
if (pInfo.Normal != null)
{
light = GetLightF(normal);
}
retVal.X = light.X;
retVal.Y = light.Y;
retVal.Z = light.Z;
retVal.X = (float)(_lightIntensity * retVal.X + (1 - _lightIntensity) * (1 - _shadowPlane[x, y]));
retVal.Y = (float)(_lightIntensity * retVal.Y + (1 - _lightIntensity) * (1 - _shadowPlane[x, y]));
retVal.Z = (float)(_lightIntensity * retVal.Z + (1 - _lightIntensity) * (1 - _shadowPlane[x, y]));
if (retVal.X < 0)
retVal.X = 0;
if (retVal.Y < 0)
retVal.Y = 0;
if (retVal.Z < 0)
retVal.Z = 0;
if (retVal.X > 1)
retVal.X = 1;
if (retVal.Y > 1)
retVal.Y = 1;
if (retVal.Z > 1)
retVal.Z = 1;
float brightLightLevel = (float)ParameterDict.Current.GetDouble("Renderer.BrightLightLevel");
if (brightLightLevel > 0)
{
retVal.X = (1 - brightLightLevel) * retVal.X + brightLightLevel * light.X * (1 - _shadowPlane[x, y]);
retVal.Y = (1 - brightLightLevel) * retVal.Y + brightLightLevel * light.Y * (1 - _shadowPlane[x, y]);
retVal.Z = (1 - brightLightLevel) * retVal.Z + brightLightLevel * light.Z * (1 - _shadowPlane[x, y]);
}
if (retVal.X < 0)
retVal.X = 0;
if (retVal.Y < 0)
retVal.Y = 0;
if (retVal.Z < 0)
retVal.Z = 0;
if (retVal.X > 1)
retVal.X = 1;
if (retVal.Y > 1)
retVal.Y = 1;
if (retVal.Z > 1)
retVal.Z = 1;
// Add surface color
bool useAdditionalColorinfo = true;
if (_colorIntensity <= 0)
useAdditionalColorinfo = false;
if (useAdditionalColorinfo && ((pInfo.IsInside && _colorInside) || (!pInfo.IsInside && _colorOutside)))
{
if (pInfo != null && pInfo.AdditionalInfo != null)
{
if (double.IsNaN(pInfo.AdditionalInfo.red))
pInfo.AdditionalInfo.red = 0;
if (double.IsNaN(pInfo.AdditionalInfo.green))
pInfo.AdditionalInfo.green = 0;
if (double.IsNaN(pInfo.AdditionalInfo.blue))
pInfo.AdditionalInfo.blue = 0;
if (pInfo.AdditionalInfo.red < 0)
{
pInfo.AdditionalInfo.green -= pInfo.AdditionalInfo.red;
pInfo.AdditionalInfo.blue -= pInfo.AdditionalInfo.red;
pInfo.AdditionalInfo.red = 0;
}
if (pInfo.AdditionalInfo.green < 0)
{
pInfo.AdditionalInfo.red -= pInfo.AdditionalInfo.green;
pInfo.AdditionalInfo.blue -= pInfo.AdditionalInfo.green;
pInfo.AdditionalInfo.green = 0;
}
if (pInfo.AdditionalInfo.blue < 0)
{
pInfo.AdditionalInfo.red -= pInfo.AdditionalInfo.blue;
pInfo.AdditionalInfo.blue -= pInfo.AdditionalInfo.blue;
pInfo.AdditionalInfo.blue = 0;
}
// Normalise;
float r1 = (float)(_colorFactorRed * Math.Pow(pInfo.AdditionalInfo.red, _colorIntensity));
float g1 = (float)(_colorFactorGreen * Math.Pow(pInfo.AdditionalInfo.green, _colorIntensity));
float b1 = (float)(_colorFactorBlue * Math.Pow(pInfo.AdditionalInfo.blue, _colorIntensity));
if (r1 < 0)
r1 = -r1;
if (g1 < 0)
g1 = -g1;
if (b1 < 0)
b1 = -b1;
// Normalize:
float norm = (float)(Math.Sqrt(r1 * r1 + g1 * g1 + b1 * b1) / Math.Sqrt(2.5));
if (norm > 0)
{
r1 = r1 / norm;
g1 = g1 / norm;
b1 = b1 / norm;
if (_colorThreshold > 0)
{
double thresholdIndicator = Math.Max(Math.Abs(r1 - b1), Math.Max(Math.Abs(r1 - g1), Math.Abs(g1 - b1)));
float lowerThresholdFactor = (float)0.7;
if (thresholdIndicator < _colorThreshold * lowerThresholdFactor)
{
r1 = (float)1;
g1 = (float)1;
b1 = (float)1;
}
else if (thresholdIndicator < _colorThreshold)
{
r1 = (float)(_colorThreshold - thresholdIndicator) * (float)0.5 / lowerThresholdFactor + r1;
g1 = (float)(_colorThreshold - thresholdIndicator) * (float)0.5 / lowerThresholdFactor + g1;
b1 = (float)(_colorThreshold - thresholdIndicator) * (float)0.5 / lowerThresholdFactor + b1;
}
}
}
for (int i = 0; i < 5; i++)
{
if (r1 > 1)
{
b1 += (r1 - 1) / (float)2.0;
g1 += (r1 - 1) / (float)2.0;
r1 = 1;
}
if (b1 > 1)
{
r1 += (b1 - 1) / (float)2.0;
g1 += (b1 - 1) / (float)2.0;
b1 = 1;
}
if (g1 > 1)
{
r1 += (g1 - 1) / (float)2.0;
b1 += (g1 - 1) / (float)2.0;
g1 = 1;
}
}
if (r1 > 1)
r1 = 1;
if (b1 > 1)
b1 = 1;
if (g1 > 1)
g1 = 1;
if (_colorGreyness > 0)
{
r1 = (float)(_colorGreyness + (1 - _colorGreyness) * r1);
g1 = (float)(_colorGreyness + (1 - _colorGreyness) * g1);
b1 = (float)(_colorGreyness + (1 - _colorGreyness) * b1);
}
if (r1 > 1)
r1 = 1;
if (b1 > 1)
b1 = 1;
if (g1 > 1)
g1 = 1;
if (norm != 0)
{
switch (_rgbType)
{
case 1:
retVal.X *= r1;
retVal.Y *= g1;
retVal.Z *= b1;
break;
case 2:
retVal.X *= r1;
retVal.Y *= b1;
retVal.Z *= g1;
break;
case 3:
retVal.X *= g1;
retVal.Y *= r1;
retVal.Z *= b1;
break;
case 4:
retVal.X *= g1;
retVal.Y *= b1;
retVal.Z *= r1;
break;
case 5:
retVal.X *= b1;
retVal.Y *= r1;
retVal.Z *= g1;
break;
case 6:
retVal.X *= b1;
retVal.Y *= g1;
retVal.Z *= r1;
break;
default:
retVal.X *= r1;
retVal.Y *= g1;
retVal.Z *= b1;
break;
}
}
}
}
if (_contrast != 1)
{
retVal.X = (float)Math.Pow(retVal.X, _contrast);
retVal.Y = (float)Math.Pow(retVal.Y, _contrast);
retVal.Z = (float)Math.Pow(retVal.Z, _contrast);
}
if (_brightness > 1)
{
retVal.X *= _brightness;
retVal.Y *= _brightness;
retVal.Z *= _brightness;
}
if (retVal.X < 0)
retVal.X = 0;
if (retVal.X > 1)
retVal.X = 1;
if (retVal.Y < 0)
retVal.Y = 0;
if (retVal.Z < 0)
retVal.Z = 0;
if (retVal.Y > 1)
retVal.Y = 1;
if (retVal.Z > 1)
retVal.Z = 1;
if (pInfo != null && pInfo.AdditionalInfo != null)
{
pInfo.AdditionalInfo.red2 = retVal.X;
pInfo.AdditionalInfo.green2 = retVal.Y;
pInfo.AdditionalInfo.blue2 = retVal.Z;
}
return new Vec3(retVal.X, retVal.Y, retVal.Z);
}
public FloatVec3(FloatVec3 other)
{
this.X = other.X;
this.Y = other.Y;
this.Z = other.Z;
}