public static RayImage ArrayLogRange(double[,] imagedata_in, ColorGradient c)
{
int ni = imagedata_in.GetLength(0);
int nj = imagedata_in.GetLength(1);
double[,] imagedata = new double[ni, nj];
for (int i = 0; i < ni; i++)
{
for (int j = 0; j < nj; j++)
{
imagedata[i, j] = Math.Log10(Math.Abs(imagedata_in[i, j]) + 1e-10);
}
}
RayImage output = new RayImage(ni, nj);
double minn = double.PositiveInfinity;
double maxx = double.NegativeInfinity;
for (int i = 0; i < ni; i++)
{
for (int j = 0; j < nj; j++)
{
maxx = (imagedata[i, j] > maxx) ? imagedata[i, j] : maxx;
minn = (imagedata[i, j] < minn) ? imagedata[i, j] : minn;
}
}
for (int i = 0; i < ni; i++)
{
for (int j = 0; j < nj; j++)
{
output.SetPixel(i, j, c.GetColor(imagedata[i, j], minn, maxx));
}
}
return(output);
}
public static RayImage ArrayRangeXY(double[,] imagedata, ColorGradient c)
{
int ni = imagedata.GetLength(0);
int nj = imagedata.GetLength(1);
RayImage output = new RayImage(ni, nj);
double minn = double.PositiveInfinity;
double maxx = double.NegativeInfinity;
for (int i = 0; i < ni; i++)
{
for (int j = 0; j < nj; j++)
{
maxx = (imagedata[i, j] > maxx) ? imagedata[i, j] : maxx;
minn = (imagedata[i, j] < minn) ? imagedata[i, j] : minn;
}
}
for (int i = 0; i < ni; i++)
{
for (int j = 0; j < nj; j++)
{
Triple col = c.GetColor(imagedata[i, j], minn, maxx);
output.SetPixelXY(i, j, col);
}
}
return(output);
}
public static void RdrTimeMaps()
{
int nx = 800;
int ny = 600;
double radius = 31;
double cam_z = 30;
int N = 60;
double dtheta = 2 * Math.PI / N;
for (int i = 0; i < N; i++)
{
Info.WriteLine(i);
double theta = 0.5 * Math.PI + i * dtheta;
Camera c = new Camera(new Triple(-radius * Math.Cos(-theta), radius * Math.Sin(-theta), cam_z), 0, 0, nx, ny, 0.9);
c.AzimuthAngle = theta;
Background basic = new Background();
GlobalLightSource light = new GlobalLightSource(new Triple(1, 2, -4));
Scene main_scene = new Scene(basic, c);
Stl stl_subject = new Stl("stl/cat.stl");
FacetBody subject = stl_subject.ToFacetBody(new Triple(0, 0, 25));
RectangularPrism floor = new RectangularPrism(new Triple(0, 0, 21.3), 190, 190, 0.3);
subject.BodyOpticalProperties.BaseColor = new Triple(0.2, 0.2, 0.9);
subject.BodyOpticalProperties.IsReflective = true;
subject.BodyOpticalProperties.Reflectivity = 0.21;
main_scene.AddBody(floor);
main_scene.AddBody(subject);
main_scene.AddLight(light);
double[,] dist, ts;
int[,] id;
RayImage picture = main_scene.Render(out dist, out id, out ts);
RayImage dist_picture = RayImage.ArrayLogRangeXY(dist, ColorGradient.Jedi());
RayImage t_picture = RayImage.ArrayLogRangeXY(ts, ColorGradient.Jedi());
picture.Save("outputdata/images/" + i.ToString().PadLeft(N.ToString().Length, '0') + ".png");
dist_picture.Save("outputdata/dists/" + i.ToString().PadLeft(N.ToString().Length, '0') + ".png");
t_picture.Save("outputdata/times/" + i.ToString().PadLeft(N.ToString().Length, '0') + ".png");
}
/*double[,] a = new double[nx, ny];
* for (int x = 0; x < nx; x++)
* {
* for (int y = 0; y < ny; y++)
* {
* a[x,y] = (double)(x+y);
* }
* }
* RayImage test = RayImage.ArrayRangeXY(a, ColorGradient.Rgb());
* test.Save("frames/s.png");*/
}
public static void Liberty()
{
int nx = 1950;
int ny = 1000;
double radius = 280;
double cam_z = 26;
double theta = 0.5 * Math.PI;
Camera c = new Camera(new Triple(-radius * Math.Cos(-theta), radius * Math.Sin(-theta), cam_z), 0, 0, nx, ny, 0.9);
c.AzimuthAngle = theta;
Background basic = new Background();
GlobalLightSource light = new GlobalLightSource(new Triple(1, 2, -4));
Scene main_scene = new Scene(basic, c);
main_scene.DoShadows = true;
Info.WriteLine("Importing...");
Stl stl_subject = new Stl("stl/liberty.stl");
Info.WriteLine("Done importing.");
FacetBody subject = stl_subject.ToFacetBody(new Triple(0, 0, 250));
RectangularPrism floor = new RectangularPrism(new Triple(0, 0, subject.ZminGlobal - 0.15), 190, 190, 0.3);
c.Position.Z = subject.ZminGlobal;
c.ElevationAngle = 0.78 * Math.PI / 3;
subject.BodyOpticalProperties.BaseColor = new Triple(0.8, 0.5, 0.8);
subject.BodyOpticalProperties.IsReflective = true;
subject.BodyOpticalProperties.Reflectivity = 0.21;
//main_scene.AddBody(floor);
main_scene.AddBody(subject);
main_scene.AddLight(light);
double[,] dist, ts;
int[,] id;
Info.WriteLine("Rendering...");
RayImage picture = main_scene.Render(out dist, out id, out ts);
RayImage dist_picture = RayImage.ArrayLogRangeXY(dist, ColorGradient.Jedi());
RayImage t_picture = RayImage.ArrayLogRangeXY(ts, ColorGradient.Jedi());
Info.WriteLine("Done.");
picture.Save("frames/lib.png");
dist_picture.Save("frames/libdists.png");
t_picture.Save("frames/libtimes.png");
Utils.WriteCsv("frames/dists.csv", dist);
Utils.WriteCsv("frames/times.csv", ts);
}
public static void DestroyComputer()
{
int nx = 1950;
int ny = 1000;
double radius = 25;
double cam_z = 26;
double theta = 0.5 * Math.PI;
Camera c = new Camera(new Triple(-radius * Math.Cos(-theta), radius * Math.Sin(-theta), cam_z), 0, 0, nx, ny, 0.9);
c.AzimuthAngle = theta;
Background basic = new Background();
basic.HasFloor = false;
basic.SkyColor = new Triple(1, 1, 1);
GlobalLightSource light = new GlobalLightSource(new Triple(1, 2, -4));
Scene main_scene = new Scene(basic, c);
main_scene.DoShadows = true;
Info.WriteLine("Importing...");
Stl stl_subject = new Stl("stl/car-a.stl", true);
Info.WriteLine("Done importing.");
FacetBody subject = stl_subject.ToFacetBody(new Triple(0, 0, 25));
RectangularPrism floor = new RectangularPrism(new Triple(0, 0, subject.ZminGlobal - 0.15), 190, 190, 0.3);
subject.BodyOpticalProperties.BaseColor = new Triple(0.5, 0.5, 0.5);
subject.BodyOpticalProperties.IsReflective = false;
subject.BodyOpticalProperties.Reflectivity = 0.21;
//main_scene.AddBody(floor);
main_scene.AddBody(subject);
main_scene.AddLight(light);
double[,] dist, ts;
int[,] id;
Info.WriteLine("Rendering...");
RayImage picture = main_scene.Render(out dist, out id, out ts);
RayImage dist_picture = RayImage.ArrayLogRangeXY(dist, ColorGradient.Jedi());
RayImage t_picture = RayImage.ArrayLogRangeXY(ts, ColorGradient.Jedi());
Info.WriteLine("Done.");
picture.Save("cfdstuff/car.png");
dist_picture.Save("cfdstuff/cardists.png");
t_picture.Save("cfdstuff/cartimes.png");
}
