protected override void OnPaint(PaintEventArgs e) { chart3D1.C3Axes.XMin = -3; chart3D1.C3Axes.XMax = 3; chart3D1.C3Axes.YMin = -3; chart3D1.C3Axes.YMax = 3; chart3D1.C3Axes.ZMin = -8; chart3D1.C3Axes.ZMax = 8; chart3D1.C3Axes.XTick = 1; chart3D1.C3Axes.YTick = 1; chart3D1.C3Axes.ZTick = 4; chart3D1.C3DataSeries.XDataMin = chart3D1.C3Axes.XMin; chart3D1.C3DataSeries.YDataMin = chart3D1.C3Axes.YMin; chart3D1.C3DataSeries.XSpacing = 0.1f; chart3D1.C3DataSeries.YSpacing = 0.1f; chart3D1.C3DataSeries.XNumber = Convert.ToInt16((chart3D1.C3Axes.XMax - chart3D1.C3Axes.XMin) / chart3D1.C3DataSeries.XSpacing) + 1; chart3D1.C3DataSeries.YNumber = Convert.ToInt16((chart3D1.C3Axes.YMax - chart3D1.C3Axes.YMin) / chart3D1.C3DataSeries.YSpacing) + 1; Point3[,] pts = new Point3[chart3D1.C3DataSeries.XNumber, chart3D1.C3DataSeries.YNumber]; for (int i = 0; i < chart3D1.C3DataSeries.XNumber; i++) { for (int j = 0; j < chart3D1.C3DataSeries.YNumber; j++) { float x = chart3D1.C3DataSeries.XDataMin + i * chart3D1.C3DataSeries.XSpacing; float y = chart3D1.C3DataSeries.YDataMin + j * chart3D1.C3DataSeries.YSpacing; double zz = 3 * Math.Pow((1 - x), 2) * Math.Exp(-x * x - (y + 1) * (y + 1)) - 10 * (0.2 * x - Math.Pow(x, 3) - Math.Pow(y, 5)) * Math.Exp(-x * x - y * y) - 1 / 3 * Math.Exp(-(x + 1) * (x + 1) - y * y); float z = (float)zz; pts[i, j] = new Point3(x, y, z, 1); } } chart3D1.C3DataSeries.PointArray = pts; }
//метод для построения графика private void AddData() { FunctionGraphik.C3Axes.XMin = (float)minX; //задаем границы для графика FunctionGraphik.C3Axes.XMax = (float)maxX; FunctionGraphik.C3Axes.YMin = (float)minY; FunctionGraphik.C3Axes.YMax = (float)maxY; FunctionGraphik.C3Labels.XLabel = "Х"; //именуем оси графика FunctionGraphik.C3Labels.YLabel = "Y"; FunctionGraphik.C3Labels.ZLabel = "Z"; FunctionGraphik.C3Labels.Title = "График функции"; //заголовок для графика FunctionGraphik.C3Axes.XTick = (float)((FunctionGraphik.C3Axes.XMax - FunctionGraphik.C3Axes.XMin) / 4); //для осей задаем количество меток FunctionGraphik.C3Axes.YTick = (float)((FunctionGraphik.C3Axes.YMax - FunctionGraphik.C3Axes.YMin) / 4); FunctionGraphik.C3DataSeries.XDataMin = FunctionGraphik.C3Axes.XMin; //задаем минимальные границы, чтобы график при повороте не ускакивал FunctionGraphik.C3DataSeries.YDataMin = FunctionGraphik.C3Axes.YMin; FunctionGraphik.C3DataSeries.XSpacing = (float)((FunctionGraphik.C3Axes.XMax - FunctionGraphik.C3Axes.XMin) / 20); //шаг, с которым мы увеличиваем координаты FunctionGraphik.C3DataSeries.YSpacing = (float)((FunctionGraphik.C3Axes.YMax - FunctionGraphik.C3Axes.YMin) / 20); float min = (float)Fun(FunctionGraphik.C3Axes.XMin, FunctionGraphik.C3Axes.YMin); //для поиска минимума/максимума для координаты Z float max = (float)Fun(FunctionGraphik.C3Axes.XMin, FunctionGraphik.C3Axes.YMin); for (int i = 0; i < 20; i++) { //задаем матрицу float x = FunctionGraphik.C3DataSeries.XDataMin + i * FunctionGraphik.C3DataSeries.XSpacing; //изменяем координату по Х for (int j = 0; j < 20; j++) { float y = FunctionGraphik.C3DataSeries.YDataMin + j * FunctionGraphik.C3DataSeries.YSpacing; //изменяем координату по Y float z = (float)Fun(x, y); //задаем координату Z pts[i, j] = new Point3(x, y, z, 1); //добавляем точку в матрицу if (min > z) min = z; //ищем минимум if (max < z) max = z; //ищем максимум } } FunctionGraphik.C3Axes.ZMin = (float)min; //задаем границы оси Z FunctionGraphik.C3Axes.ZMax = (float)max; FunctionGraphik.C3Axes.ZTick = (float)((max - min) / 4); //для осей задаем количество меток FunctionGraphik.C3ViewAngle.Elevation = PovorotElevation.Value; //задаем поворот графика FunctionGraphik.C3ViewAngle.Azimuth = PovorotAzimuth.Value; //задаем поворот графика FunctionGraphik.C3DataSeries.PointArray = pts; //"даем" графику матрицу точек }
private void AddTicks(Graphics g, Axes ax, ViewAngle va, ChartLabels cl) { Matrix3 m = Matrix3.AzimuthElevation(va.Elevation, va.Azimuth); Point3[] pta = new Point3[2]; Point3[] pts = CoordinatesOfChartBox(ax, va); ; Pen aPen = new Pen(ax.AxisStyle.LineColor, ax.AxisStyle.Thickness); aPen.DashStyle = ax.AxisStyle.Pattern; // Add x ticks: float offset = (ax.YMax - ax.YMin) / 30.0f; float ticklength = offset; for (float x = ax.XMin; x <= ax.XMax; x = x + ax.XTick) { if (va.Elevation >= 0) { if (va.Azimuth >= -90 && va.Azimuth < 90) ticklength = -offset; } else if (va.Elevation < 0) { if ((va.Azimuth >= -180 && va.Azimuth < -90) || va.Azimuth >= 90 && va.Azimuth <= 180) ticklength = -(ax.YMax - ax.YMin) / 30; } pta[0] = new Point3(x, pts[1].Y + ticklength, pts[1].Z, pts[1].W); pta[1] = new Point3(x, pts[1].Y, pts[1].Z, pts[1].W); for (int i = 0; i < pta.Length; i++) { pta[i].Transform(m, chart3d, ax, this, cl); } g.DrawLine(aPen, pta[0].X, pta[0].Y, pta[1].X, pta[1].Y); } // Add y ticks: offset = (ax.XMax - ax.XMin) / 30.0f; ticklength = offset; for (float y = ax.YMin; y <= ax.YMax; y = y + ax.YTick) { pts = CoordinatesOfChartBox(ax, va); ; if (va.Elevation >= 0) { if (va.Azimuth >= -180 && va.Azimuth < 0) ticklength = -offset; } else if (va.Elevation < 0) { if (va.Azimuth >= 0 && va.Azimuth < 180) ticklength = -offset; } pta[0] = new Point3(pts[1].X + ticklength, y, pts[1].Z, pts[1].W); pta[1] = new Point3(pts[1].X, y, pts[1].Z, pts[1].W); for (int i = 0; i < pta.Length; i++) { pta[i].Transform(m, chart3d, ax, this, cl); } g.DrawLine(aPen, pta[0].X, pta[0].Y, pta[1].X, pta[1].Y); } float xoffset = (ax.XMax - ax.XMin) / 45.0f; float yoffset = (ax.YMax - ax.YMin) / 20.0f; float xticklength = xoffset; float yticklength = yoffset; for (float z = ax.ZMin; z <= ax.ZMax; z = z + ax.ZTick) { if (va.Elevation >= 0) { if (va.Azimuth >= -180 && va.Azimuth < -90) { xticklength = 0; yticklength = yoffset; } else if (va.Azimuth >= -90 && va.Azimuth < 0) { xticklength = xoffset; yticklength = 0; } else if (va.Azimuth >= 0 && va.Azimuth < 90) { xticklength = 0; yticklength = -yoffset; } else if (va.Azimuth >= 90 && va.Azimuth <= 180) { xticklength = -xoffset; yticklength = 0; } } else if (va.Elevation < 0) { if (va.Azimuth >= -180 && va.Azimuth < -90) { yticklength = 0; xticklength = xoffset; } else if (va.Azimuth >= -90 && va.Azimuth < 0) { yticklength = -yoffset; xticklength = 0; } else if (va.Azimuth >= 0 && va.Azimuth < 90) { yticklength = 0; xticklength = -xoffset; } else if (va.Azimuth >= 90 && va.Azimuth <= 180) { yticklength = yoffset; xticklength = 0; } } pta[0] = new Point3(pts[2].X, pts[2].Y, z, pts[2].W); pta[1] = new Point3(pts[2].X + yticklength, pts[2].Y + xticklength, z, pts[2].W); for (int i = 0; i < pta.Length; i++) { pta[i].Transform(m, chart3d, ax, this, cl); } g.DrawLine(aPen, pta[0].X, pta[0].Y, pta[1].X, pta[1].Y); } aPen.Dispose(); }
private void AddLabels(Graphics g, Axes ax, ViewAngle va, ChartLabels cl) { Matrix3 m = Matrix3.AzimuthElevation(va.Elevation, va.Azimuth); Point3 pt = new Point3(); Point3[] pts = CoordinatesOfChartBox(ax, va); SolidBrush aBrush = new SolidBrush(cl.LabelFontColor); StringFormat sf = new StringFormat(); sf.Alignment = StringAlignment.Center; // Add x tick labels: float offset = (ax.YMax - ax.YMin) / 20; float labelSpace = offset; for (float x = ax.XMin + ax.XTick; x < ax.XMax; x = x + ax.XTick) { if (va.Elevation >= 0) { if (va.Azimuth >= -90 && va.Azimuth < 90) labelSpace = -offset; } else if (va.Elevation < 0) { if ((va.Azimuth >= -180 && va.Azimuth < -90) || va.Azimuth >= 90 && va.Azimuth <= 180) labelSpace = -offset; } pt = new Point3(x, pts[1].Y + labelSpace, pts[1].Z, pts[1].W); pt.Transform(m, chart3d, ax, this, cl); g.DrawString(x.ToString(), cl.TickFont, aBrush, new PointF(pt.X, pt.Y), sf); } // Add y tick labels: offset = (ax.XMax - ax.XMin) / 20; labelSpace = offset; for (float y = ax.YMin + ax.YTick; y < ax.YMax; y = y + ax.YTick) { pts = CoordinatesOfChartBox(ax, va); if (va.Elevation >= 0) { if (va.Azimuth >= -180 && va.Azimuth < 0) labelSpace = -offset; } else if (va.Elevation < 0) { if (va.Azimuth >= 0 && va.Azimuth < 180) labelSpace = -offset; } pt = new Point3(pts[1].X + labelSpace, y, pts[1].Z, pts[1].W); pt.Transform(m, chart3d, ax, this, cl); g.DrawString(y.ToString(), cl.TickFont, aBrush, new PointF(pt.X, pt.Y), sf); } // Add z tick labels: float xoffset = (ax.XMax - ax.XMin) / 30.0f; float yoffset = (ax.YMax - ax.YMin) / 15.0f; float xlabelSpace = xoffset; float ylabelSpace = yoffset; SizeF s = g.MeasureString("A", cl.TickFont); for (float z = ax.ZMin; z <= ax.ZMax; z = z + ax.ZTick) { sf.Alignment = StringAlignment.Far; pts = CoordinatesOfChartBox(ax, va); if (va.Elevation >= 0) { if (va.Azimuth >= -180 && va.Azimuth < -90) { xlabelSpace = 0; ylabelSpace = yoffset; } else if (va.Azimuth >= -90 && va.Azimuth < 0) { xlabelSpace = xoffset; ylabelSpace = 0; } else if (va.Azimuth >= 0 && va.Azimuth < 90) { xlabelSpace = 0; ylabelSpace = -yoffset; } else if (va.Azimuth >= 90 && va.Azimuth <= 180) { xlabelSpace = -xoffset; ylabelSpace = 0; } } else if (va.Elevation < 0) { if (va.Azimuth >= -180 && va.Azimuth < -90) { ylabelSpace = 0; xlabelSpace = xoffset; } else if (va.Azimuth >= -90 && va.Azimuth < 0) { ylabelSpace = -yoffset; xlabelSpace = 0; } else if (va.Azimuth >= 0 && va.Azimuth < 90) { ylabelSpace = 0; xlabelSpace = -xoffset; } else if (va.Azimuth >= 90 && va.Azimuth <= 180) { ylabelSpace = yoffset; xlabelSpace = 0; } } pt = new Point3(pts[2].X + ylabelSpace, pts[2].Y + xlabelSpace, z, pts[2].W); pt.Transform(m, chart3d, ax, this, cl); g.DrawString(z.ToString(), cl.TickFont, aBrush, new PointF(pt.X - labelSpace, pt.Y - s.Height / 2), sf); } // Add Title: sf.Alignment = StringAlignment.Center; aBrush = new SolidBrush(cl.TitleColor); if (cl.Title != "No Title") { g.DrawString(cl.Title, cl.TitleFont, aBrush, new PointF(chart3d.Width / 2, chart3d.Height / 30), sf); } aBrush.Dispose(); // Add x axis label: offset = (ax.YMax - ax.YMin) / 3; labelSpace = offset; sf.Alignment = StringAlignment.Center; aBrush = new SolidBrush(cl.LabelFontColor); float offset1 = (ax.XMax - ax.XMin) / 10; float xc = offset1; if (va.Elevation >= 0) { if (va.Azimuth >= -90 && va.Azimuth < 90) labelSpace = -offset; if (va.Azimuth >= 0 && va.Azimuth <= 180) xc = -offset1; } else if (va.Elevation < 0) { if ((va.Azimuth >= -180 && va.Azimuth < -90) || va.Azimuth >= 90 && va.Azimuth <= 180) labelSpace = -offset; if (va.Azimuth >= -180 && va.Azimuth <= 0) xc = -offset1; } Point3[] pta = new Point3[2]; pta[0] = new Point3(ax.XMin, pts[1].Y + labelSpace, pts[1].Z, pts[1].W); pta[1] = new Point3((ax.XMin + ax.XMax) / 2 - xc, pts[1].Y + labelSpace, pts[1].Z, pts[1].W); pta[0].Transform(m, chart3d, ax, this, cl); pta[1].Transform(m, chart3d, ax, this, cl); float theta = (float)Math.Atan((pta[1].Y - pta[0].Y) / (pta[1].X - pta[0].X)); theta = theta * 180 / (float)Math.PI; GraphicsState gs = g.Save(); g.TranslateTransform(pta[1].X, pta[1].Y); g.RotateTransform(theta); g.DrawString(cl.XLabel, cl.LabelFont, aBrush, new PointF(0, 0), sf); g.Restore(gs); // Add y axis label: offset = (ax.XMax - ax.XMin) / 3; offset1 = (ax.YMax - ax.YMin) / 5; labelSpace = offset; float yc = ax.YTick; if (va.Elevation >= 0) { if (va.Azimuth >= -180 && va.Azimuth < 0) labelSpace = -offset; if (va.Azimuth >= -90 && va.Azimuth <= 90) yc = -offset1; } else if (va.Elevation < 0) { yc = -offset1; if (va.Azimuth >= 0 && va.Azimuth < 180) labelSpace = -offset; if (va.Azimuth >= -90 && va.Azimuth <= 90) yc = offset1; } pta[0] = new Point3(pts[1].X + labelSpace, ax.YMin, pts[1].Z, pts[1].W); pta[1] = new Point3(pts[1].X + labelSpace, (ax.YMin + ax.YMax) / 2 + yc, pts[1].Z, pts[1].W); pta[0].Transform(m, chart3d, ax, this, cl); pta[1].Transform(m, chart3d, ax, this, cl); theta = (float)Math.Atan((pta[1].Y - pta[0].Y) / (pta[1].X - pta[0].X)); theta = theta * 180 / (float)Math.PI; gs = g.Save(); g.TranslateTransform(pta[1].X, pta[1].Y); g.RotateTransform(theta); g.DrawString(cl.YLabel, cl.LabelFont, aBrush, new PointF(0, 0), sf); g.Restore(gs); // Add z axis labels: float zticklength = 10; labelSpace = -1.3f * offset; offset1 = (ax.ZMax - ax.ZMin) / 8; float zc = -offset1; for (float z = ax.ZMin; z < ax.ZMax; z = z + ax.ZTick) { SizeF size = g.MeasureString(z.ToString(), cl.TickFont); if (zticklength < size.Width) zticklength = size.Width; } float zlength = -zticklength; if (va.Elevation >= 0) { if (va.Azimuth >= -180 && va.Azimuth < -90) { zlength = -zticklength; labelSpace = -1.3f * offset; zc = -offset1; } else if (va.Azimuth >= -90 && va.Azimuth < 0) { zlength = zticklength; labelSpace = 2 * offset / 3; zc = offset1; } else if (va.Azimuth >= 0 && va.Azimuth < 90) { zlength = zticklength; labelSpace = 2 * offset / 3; zc = -offset1; } else if (va.Azimuth >= 90 && va.Azimuth <= 180) { zlength = -zticklength; labelSpace = -1.3f * offset; zc = offset1; } } else if (va.Elevation < 0) { if (va.Azimuth >= -180 && va.Azimuth < -90) { zlength = -zticklength; labelSpace = -1.3f * offset; zc = offset1; } else if (va.Azimuth >= -90 && va.Azimuth < 0) { zlength = zticklength; labelSpace = 2 * offset / 3; zc = -offset1; } else if (va.Azimuth >= 0 && va.Azimuth < 90) { zlength = zticklength; labelSpace = 2 * offset / 3; zc = offset1; } else if (va.Azimuth >= 90 && va.Azimuth <= 180) { zlength = -zticklength; labelSpace = -1.3f * offset; zc = -offset1; } } pta[0] = new Point3(pts[2].X - labelSpace, pts[2].Y, (ax.ZMin + ax.ZMax) / 2 + zc, pts[2].W); pta[0].Transform(m, chart3d, ax, this, cl); gs = g.Save(); g.TranslateTransform(pta[0].X - zlength, pta[0].Y); g.RotateTransform(270); g.DrawString(cl.ZLabel, cl.LabelFont, aBrush, new PointF(0, 0), sf); g.Restore(gs); }
private void AddGrids(Graphics g, Axes ax, ViewAngle va, Grid gd, ChartLabels cl) { Matrix3 m = Matrix3.AzimuthElevation(va.Elevation, va.Azimuth); Point3[] pta = new Point3[3]; Point3[] pts = CoordinatesOfChartBox(ax, va); Pen aPen = new Pen(gd.GridStyle.LineColor, gd.GridStyle.Thickness); aPen.DashStyle = gd.GridStyle.Pattern; // Draw x gridlines: if (gd.IsXGrid) { for (float x = ax.XMin; x <= ax.XMax; x = x + ax.XTick) { pts = CoordinatesOfChartBox(ax, va); pta[0] = new Point3(x, pts[1].Y, pts[1].Z, pts[1].W); if (va.Elevation >= 0) { if ((va.Azimuth >= -180 && va.Azimuth < -90) || (va.Azimuth >= 0 && va.Azimuth < 90)) { pta[1] = new Point3(x, pts[0].Y, pts[1].Z, pts[1].W); pta[2] = new Point3(x, pts[0].Y, pts[3].Z, pts[1].W); } else { pta[1] = new Point3(x, pts[2].Y, pts[1].Z, pts[1].W); pta[2] = new Point3(x, pts[2].Y, pts[3].Z, pts[1].W); } } else if (va.Elevation < 0) { if ((va.Azimuth >= -180 && va.Azimuth < -90) || (va.Azimuth >= 0 && va.Azimuth < 90)) { pta[1] = new Point3(x, pts[2].Y, pts[1].Z, pts[1].W); pta[2] = new Point3(x, pts[2].Y, pts[3].Z, pts[1].W); } else { pta[1] = new Point3(x, pts[0].Y, pts[1].Z, pts[1].W); pta[2] = new Point3(x, pts[0].Y, pts[3].Z, pts[1].W); } } for (int i = 0; i < pta.Length; i++) { pta[i].Transform(m, chart3d, ax, this, cl); } g.DrawLine(aPen, pta[0].X, pta[0].Y, pta[1].X, pta[1].Y); g.DrawLine(aPen, pta[1].X, pta[1].Y, pta[2].X, pta[2].Y); } } // Draw y gridlines: if (gd.IsYGrid) { for (float y = ax.YMin; y <= ax.YMax; y = y + ax.YTick) { pts = CoordinatesOfChartBox(ax, va); pta[0] = new Point3(pts[1].X, y, pts[1].Z, pts[1].W); if (va.Elevation >= 0) { if ((va.Azimuth >= -180 && va.Azimuth < -90) || (va.Azimuth >= 0 && va.Azimuth < 90)) { pta[1] = new Point3(pts[2].X, y, pts[1].Z, pts[1].W); pta[2] = new Point3(pts[2].X, y, pts[3].Z, pts[1].W); } else { pta[1] = new Point3(pts[0].X, y, pts[1].Z, pts[1].W); pta[2] = new Point3(pts[0].X, y, pts[3].Z, pts[1].W); } } if (va.Elevation < 0) { if ((va.Azimuth >= -180 && va.Azimuth < -90) || (va.Azimuth >= 0 && va.Azimuth < 90)) { pta[1] = new Point3(pts[0].X, y, pts[1].Z, pts[1].W); pta[2] = new Point3(pts[0].X, y, pts[3].Z, pts[1].W); } else { pta[1] = new Point3(pts[2].X, y, pts[1].Z, pts[1].W); pta[2] = new Point3(pts[2].X, y, pts[3].Z, pts[1].W); } } for (int i = 0; i < pta.Length; i++) { pta[i].Transform(m, chart3d, ax, this, cl); } g.DrawLine(aPen, pta[0].X, pta[0].Y, pta[1].X, pta[1].Y); g.DrawLine(aPen, pta[1].X, pta[1].Y, pta[2].X, pta[2].Y); } } // Draw Z gridlines: if (gd.IsZGrid) { for (float z = ax.ZMin; z <= ax.ZMax; z = z + ax.ZTick) { pts = CoordinatesOfChartBox(ax, va); pta[0] = new Point3(pts[2].X, pts[2].Y, z, pts[2].W); if (va.Elevation >= 0) { if ((va.Azimuth >= -180 && va.Azimuth < -90) || (va.Azimuth >= 0 && va.Azimuth < 90)) { pta[1] = new Point3(pts[2].X, pts[0].Y, z, pts[1].W); pta[2] = new Point3(pts[0].X, pts[0].Y, z, pts[1].W); } else { pta[1] = new Point3(pts[0].X, pts[2].Y, z, pts[1].W); pta[2] = new Point3(pts[0].X, pts[1].Y, z, pts[1].W); } } if (va.Elevation < 0) { if ((va.Azimuth >= -180 && va.Azimuth < -90) || (va.Azimuth >= 0 && va.Azimuth < 90)) { pta[1] = new Point3(pts[0].X, pts[2].Y, z, pts[1].W); pta[2] = new Point3(pts[0].X, pts[0].Y, z, pts[1].W); } else { pta[1] = new Point3(pts[2].X, pts[0].Y, z, pts[1].W); pta[2] = new Point3(pts[0].X, pts[0].Y, z, pts[1].W); } } for (int i = 0; i < pta.Length; i++) { pta[i].Transform(m, chart3d, ax, this, cl); } g.DrawLine(aPen, pta[0].X, pta[0].Y, pta[1].X, pta[1].Y); g.DrawLine(aPen, pta[1].X, pta[1].Y, pta[2].X, pta[2].Y); } } }
private void AddXYColor(Graphics g, DataSeries ds, ChartStyle cs, ChartStyle2D cs2d, Axes ax, ViewAngle va, ChartLabels cl) { Pen aPen = new Pen(ds.LineStyle.LineColor, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; SolidBrush aBrush = new SolidBrush(Color.White); PointF[] pta = new PointF[4]; Point3[,] pts = ds.PointArray; Matrix3 m = new Matrix3(); // Find the minumum and maximum z values: float zmin = ds.ZDataMin(); float zmax = ds.ZDataMax(); // Draw surface on the XY plane: if (!IsInterp) { for (int i = 0; i < pts.GetLength(0) - 1; i++) { for (int j = 0; j < pts.GetLength(1) - 1; j++) { pta[0] = cs2d.Point2D(new PointF(pts[i, j].X, pts[i, j].Y), ax); pta[1] = cs2d.Point2D(new PointF(pts[i, j + 1].X, pts[i, j + 1].Y), ax); pta[2] = cs2d.Point2D(new PointF(pts[i + 1, j + 1].X, pts[i + 1, j + 1].Y), ax); pta[3] = cs2d.Point2D(new PointF(pts[i + 1, j].X, pts[i + 1, j].Y), ax); Color color = AddColor(cs, pts[i, j], zmin, zmax, ax, va, cl); aBrush = new SolidBrush(color); g.FillPolygon(aBrush, pta); if (ds.LineStyle.IsVisible) { g.DrawPolygon(aPen, pta); } } } } // Draw refined surface: else if (IsInterp) { for (int i = 0; i < pts.GetLength(0) - 1; i++) { for (int j = 0; j < pts.GetLength(1) - 1; j++) { Point3[] points = new Point3[4]; points[0] = pts[i, j]; points[1] = pts[i, j + 1]; points[2] = pts[i + 1, j + 1]; points[3] = pts[i + 1, j]; Interp(g, cs, cs2d, m, points, zmin, zmax, 2, ax, va, cl); pta[0] = cs2d.Point2D(new PointF(pts[i, j].X, pts[i, j].Y), ax); pta[1] = cs2d.Point2D(new PointF(pts[i, j + 1].X, pts[i, j + 1].Y), ax); pta[2] = cs2d.Point2D(new PointF(pts[i + 1, j + 1].X, pts[i + 1, j + 1].Y), ax); pta[3] = cs2d.Point2D(new PointF(pts[i + 1, j].X, pts[i + 1, j].Y), ax); if (ds.LineStyle.IsVisible) { g.DrawPolygon(aPen, pta); } } } } }
public Point3 Spherical(float r, float theta, float phi) { Point3 pt = new Point3(); float snt = (float)Math.Sin(theta * Math.PI / 180); float cnt = (float)Math.Cos(theta * Math.PI / 180); float snp = (float)Math.Sin(phi * Math.PI / 180); float cnp = (float)Math.Cos(phi * Math.PI / 180); pt.X = r * snt * cnp; pt.Y = r * cnt; pt.Z = -r * snt * snp; pt.W = 1; return pt; }
public void Peak3D(DataSeries ds, Axes ax) { ax.XMin = -3; ax.XMax = 3; ax.YMin = -3; ax.YMax = 3; ax.ZMin = -8; ax.ZMax = 8; ax.XTick = 1; ax.YTick = 1; ax.ZTick = 4; ds.XDataMin = ax.XMin; ds.YDataMin = ax.YMin; ds.XSpacing = 0.3f; ds.YSpacing = 0.3f; ds.XNumber = Convert.ToInt16((ax.XMax - ax.XMin) / ds.XSpacing) + 1; ds.YNumber = Convert.ToInt16((ax.YMax - ax.YMin) / ds.YSpacing) + 1; Point3[,] pts = new Point3[ds.XNumber, ds.YNumber]; for (int i = 0; i < ds.XNumber; i++) { for (int j = 0; j < ds.YNumber; j++) { float x = ds.XDataMin + i * ds.XSpacing; float y = ds.YDataMin + j * ds.YSpacing; double zz = 3 * Math.Pow((1 - x), 2) * Math.Exp(-x * x - (y + 1) * (y + 1)) - 10 * (0.2 * x - Math.Pow(x, 3) - Math.Pow(y, 5)) * Math.Exp(-x * x - y * y) - 1 / 3 * Math.Exp(-(x + 1) * (x + 1) - y * y); float z = (float)zz; pts[i, j] = new Point3(x, y, z, 1); } } ds.PointArray = pts; }
private void AddLine(Graphics g, DataSeries ds, ChartStyle cs, Axes ax, ViewAngle va, ChartLabels cl) { Pen aPen = new Pen(ds.LineStyle.LineColor, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; Matrix3 m = Matrix3.AzimuthElevation(va.Elevation, va.Azimuth); Point3[] pts = new Point3[ds.PointList.Count]; // Find zmin and zmax values, then perform transformation on points: float zmin = 0; float zmax = 0; for (int i = 0; i < pts.Length; i++) { pts[i] = (Point3)ds.PointList[i]; zmin = Math.Min(zmin, pts[i].Z); zmax = Math.Max(zmax, pts[i].Z); pts[i].Transform(m, chart3d, ax, cs, cl); } // Draw line: if (ds.LineStyle.IsVisible == true) { for (int i = 1; i < pts.Length; i++) { if (IsColorMap) { Color color = AddColor(cs, pts[i], zmin, zmax, ax, va, cl); aPen = new Pen(color, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; } g.DrawLine(aPen, pts[i - 1].X, pts[i - 1].Y, pts[i].X, pts[i].Y); } } aPen.Dispose(); }
private void AddContour3D(Graphics g, DataSeries ds, ChartStyle cs, ChartStyle2D cs2d, Axes ax, ViewAngle va, ChartLabels cl) { Pen aPen = new Pen(ds.LineStyle.LineColor, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; Point3[] pta = new Point3[2]; Point3[,] pts = ds.PointArray; Matrix3 m = Matrix3.AzimuthElevation(va.Elevation, va.Azimuth); // Find the minumum and maximum z values: float zmin = ds.ZDataMin(); float zmax = ds.ZDataMax(); float[] zlevels = new float[numberContours]; for (int i = 0; i < numberContours; i++) { zlevels[i] = zmin + i * (zmax - zmin) / (numberContours - 1); } int i0, i1, i2, j0, j1, j2; float zratio = 1; // Draw contour on the XY plane: for (int i = 0; i < pts.GetLength(0) - 1; i++) { for (int j = 0; j < pts.GetLength(1) - 1; j++) { for (int k = 0; k < numberContours; k++) { if (IsColorMap && ChartType != ChartTypeEnum.FillContour && ChartType != ChartTypeEnum.SurfaceFillContour) { Color color = AddColor(cs, pts[i, j], zmin, zmax, ax, va, cl); aPen = new Pen(color, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; } // Left triangle: i0 = i; j0 = j; i1 = i; j1 = j + 1; i2 = i + 1; j2 = j + 1; if ((zlevels[k] >= pts[i0, j0].Z && zlevels[k] < pts[i1, j1].Z || zlevels[k] < pts[i0, j0].Z && zlevels[k] >= pts[i1, j1].Z) && (zlevels[k] >= pts[i1, j1].Z && zlevels[k] < pts[i2, j2].Z || zlevels[k] < pts[i1, j1].Z && zlevels[k] >= pts[i2, j2].Z)) { zratio = (zlevels[k] - pts[i0, j0].Z) / (pts[i1, j1].Z - pts[i0, j0].Z); pta[0] = new Point3(pts[i0, j0].X, (1 - zratio) * pts[i0, j0].Y + zratio * pts[i1, j1].Y, ax.ZMin, 1); zratio = (zlevels[k] - pts[i1, j1].Z) / (pts[i2, j2].Z - pts[i1, j1].Z); pta[1] = new Point3((1 - zratio) * pts[i1, j1].X + zratio * pts[i2, j2].X, pts[i1, j1].Y, ax.ZMin, 1); pta[0].Transform(m, chart3d, ax, cs, cl); pta[1].Transform(m, chart3d, ax, cs, cl); g.DrawLine(aPen, pta[0].X, pta[0].Y, pta[1].X, pta[1].Y); } else if ((zlevels[k] >= pts[i0, j0].Z && zlevels[k] < pts[i2, j2].Z || zlevels[k] < pts[i0, j0].Z && zlevels[k] >= pts[i2, j2].Z) && (zlevels[k] >= pts[i1, j1].Z && zlevels[k] < pts[i2, j2].Z || zlevels[k] < pts[i1, j1].Z && zlevels[k] >= pts[i2, j2].Z)) { zratio = (zlevels[k] - pts[i0, j0].Z) / (pts[i2, j2].Z - pts[i0, j0].Z); pta[0] = new Point3((1 - zratio) * pts[i0, j0].X + zratio * pts[i2, j2].X, (1 - zratio) * pts[i0, j0].Y + zratio * pts[i2, j2].Y, ax.ZMin, 1); zratio = (zlevels[k] - pts[i1, j1].Z) / (pts[i2, j2].Z - pts[i1, j1].Z); pta[1] = new Point3((1 - zratio) * pts[i1, j1].X + zratio * pts[i2, j2].X, pts[i1, j1].Y, ax.ZMin, 1); pta[0].Transform(m, chart3d, ax, cs, cl); pta[1].Transform(m, chart3d, ax, cs, cl); g.DrawLine(aPen, pta[0].X, pta[0].Y, pta[1].X, pta[1].Y); } else if ((zlevels[k] >= pts[i0, j0].Z && zlevels[k] < pts[i1, j1].Z || zlevels[k] < pts[i0, j0].Z && zlevels[k] >= pts[i1, j1].Z) && (zlevels[k] >= pts[i0, j0].Z && zlevels[k] < pts[i2, j2].Z || zlevels[k] < pts[i0, j0].Z && zlevels[k] >= pts[i2, j2].Z)) { zratio = (zlevels[k] - pts[i0, j0].Z) / (pts[i1, j1].Z - pts[i0, j0].Z); pta[0] = new Point3(pts[i0, j0].X, (1 - zratio) * pts[i0, j0].Y + zratio * pts[i1, j1].Y, ax.ZMin, 1); zratio = (zlevels[k] - pts[i0, j0].Z) / (pts[i2, j2].Z - pts[i0, j0].Z); pta[1] = new Point3(pts[i0, j0].X * (1 - zratio) + pts[i2, j2].X * zratio, pts[i0, j0].Y * (1 - zratio) + pts[i2, j2].Y * zratio, ax.ZMin, 1); pta[0].Transform(m, chart3d, ax, cs, cl); pta[1].Transform(m, chart3d, ax, cs, cl); g.DrawLine(aPen, pta[0].X, pta[0].Y, pta[1].X, pta[1].Y); } // right triangle: i0 = i; j0 = j; i1 = i + 1; j1 = j; i2 = i + 1; j2 = j + 1; if ((zlevels[k] >= pts[i0, j0].Z && zlevels[k] < pts[i1, j1].Z || zlevels[k] < pts[i0, j0].Z && zlevels[k] >= pts[i1, j1].Z) && (zlevels[k] >= pts[i1, j1].Z && zlevels[k] < pts[i2, j2].Z || zlevels[k] < pts[i1, j1].Z && zlevels[k] >= pts[i2, j2].Z)) { zratio = (zlevels[k] - pts[i0, j0].Z) / (pts[i1, j1].Z - pts[i0, j0].Z); pta[0] = new Point3(pts[i0, j0].X * (1 - zratio) + pts[i1, j1].X * zratio, pts[i0, j0].Y, ax.ZMin, 1); zratio = (zlevels[k] - pts[i1, j1].Z) / (pts[i2, j2].Z - pts[i1, j1].Z); pta[1] = new Point3(pts[i1, j1].X, pts[i1, j1].Y * (1 - zratio) + pts[i2, j2].Y * zratio, ax.ZMin, 1); pta[0].Transform(m, chart3d, ax, cs, cl); pta[1].Transform(m, chart3d, ax, cs, cl); g.DrawLine(aPen, pta[0].X, pta[0].Y, pta[1].X, pta[1].Y); } else if ((zlevels[k] >= pts[i0, j0].Z && zlevels[k] < pts[i2, j2].Z || zlevels[k] < pts[i0, j0].Z && zlevels[k] >= pts[i2, j2].Z) && (zlevels[k] >= pts[i1, j1].Z && zlevels[k] < pts[i2, j2].Z || zlevels[k] < pts[i1, j1].Z && zlevels[k] >= pts[i2, j2].Z)) { zratio = (zlevels[k] - pts[i0, j0].Z) / (pts[i2, j2].Z - pts[i0, j0].Z); pta[0] = new Point3(pts[i0, j0].X * (1 - zratio) + pts[i2, j2].X * zratio, pts[i0, j0].Y * (1 - zratio) + pts[i2, j2].Y * zratio, ax.ZMin, 1); zratio = (zlevels[k] - pts[i1, j1].Z) / (pts[i2, j2].Z - pts[i1, j1].Z); pta[1] = new Point3(pts[i1, j1].X, pts[i1, j1].Y * (1 - zratio) + pts[i2, j2].Y * zratio, ax.ZMin, 1); pta[0].Transform(m, chart3d, ax, cs, cl); pta[1].Transform(m, chart3d, ax, cs, cl); g.DrawLine(aPen, pta[0].X, pta[0].Y, pta[1].X, pta[1].Y); } else if ((zlevels[k] >= pts[i0, j0].Z && zlevels[k] < pts[i1, j1].Z || zlevels[k] < pts[i0, j0].Z && zlevels[k] >= pts[i1, j1].Z) && (zlevels[k] >= pts[i0, j0].Z && zlevels[k] < pts[i2, j2].Z || zlevels[k] < pts[i0, j0].Z && zlevels[k] >= pts[i2, j2].Z)) { zratio = (zlevels[k] - pts[i0, j0].Z) / (pts[i1, j1].Z - pts[i0, j0].Z); pta[0] = new Point3(pts[i0, j0].X * (1 - zratio) + pts[i1, j1].X * zratio, pts[i0, j0].Y, ax.ZMin, 1); zratio = (zlevels[k] - pts[i0, j0].Z) / (pts[i2, j2].Z - pts[i0, j0].Z); pta[1] = new Point3(pts[i0, j0].X * (1 - zratio) + pts[i2, j2].X * zratio, pts[i0, j0].Y * (1 - zratio) + pts[i2, j2].Y * zratio, ax.ZMin, 1); pta[0].Transform(m, chart3d, ax, cs, cl); pta[1].Transform(m, chart3d, ax, cs, cl); g.DrawLine(aPen, pta[0].X, pta[0].Y, pta[1].X, pta[1].Y); } } } } aPen.Dispose(); }
private Color AddColor(ChartStyle cs, Point3 pt, float zmin, float zmax, Axes ax, ViewAngle va, ChartLabels cl) { int colorLength = CMap.GetLength(0); int cindex = (int)Math.Round((colorLength * (pt.Z - zmin) + (zmax - pt.Z)) / (zmax - zmin)); if (cindex < 1) cindex = 1; if (cindex > colorLength) cindex = colorLength; Color color = Color.FromArgb(CMap[cindex - 1, 0], CMap[cindex - 1, 1], CMap[cindex - 1, 2], CMap[cindex - 1, 3]); return color; }
public void AddColorBar(Graphics g, DataSeries ds, ChartStyle cs, ChartStyle2D cs2d, Axes ax, ViewAngle va, ChartLabels cl) { if (cs.IsColorBar && IsColorMap) { Pen aPen = new Pen(Color.Black, 1); SolidBrush aBrush = new SolidBrush(cl.TickFontColor); StringFormat sFormat = new StringFormat(); sFormat.Alignment = StringAlignment.Near; SizeF size = g.MeasureString("A", cl.TickFont); int x, y, width, height; Point3[] pts = new Point3[64]; PointF[] pta = new PointF[4]; float zmin, zmax; if (ChartType == ChartTypeEnum.Slice) { zmin = ds.VDataMin(); zmax = ds.VDataMax(); } else { zmin = ds.ZDataMin(); zmax = ds.ZDataMax(); } float dz = (zmax - zmin) / 63; if (ChartType == ChartTypeEnum.Contour || ChartType == ChartTypeEnum.FillContour || ChartType == ChartTypeEnum.XYColor) { x = 5 * cs2d.ChartArea.Width / 6; y = cs2d.PlotArea.Y; width = cs2d.ChartArea.Width / 25; height = cs2d.PlotArea.Height; // Add color bar: for (int i = 0; i < 64; i++) { pts[i] = new Point3(x, y, zmin + i * dz, 1); } for (int i = 0; i < 63; i++) { Color color = AddColor(cs, pts[i], zmin, zmax, ax, va, cl); aBrush = new SolidBrush(color); float y1 = y + height - (pts[i].Z - zmin) * height / (zmax - zmin); float y2 = y + height - (pts[i + 1].Z - zmin) * height / (zmax - zmin); pta[0] = new PointF(x, y2); pta[1] = new PointF(x + width, y2); pta[2] = new PointF(x + width, y1); pta[3] = new PointF(x, y1); g.FillPolygon(aBrush, pta); } g.DrawRectangle(aPen, x, y, width, height); // Add ticks and labels to the color bar: float ticklength = 0.1f * width; for (float z = zmin; z <= zmax; z = z + (zmax - zmin) / 6) { float yy = y + height - (z - zmin) * height / (zmax - zmin); g.DrawLine(aPen, x, yy, x + ticklength, yy); g.DrawLine(aPen, x + width, yy, x + width - ticklength, yy); g.DrawString((Math.Round(z, 2)).ToString(), cl.TickFont, Brushes.Black, new PointF(x + width + 5, yy - size.Height / 2), sFormat); } } else { x = 5 * cs2d.ChartArea.Width / 6; y = cs2d.ChartArea.Height / 10; width = cs2d.ChartArea.Width / 25; height = 8 * cs2d.ChartArea.Height / 10; // Add color bar: for (int i = 0; i < 64; i++) { pts[i] = new Point3(x, y, zmin + i * dz, 1); } for (int i = 0; i < 63; i++) { Color color = AddColor(cs, pts[i], zmin, zmax, ax, va, cl); aBrush = new SolidBrush(color); float y1 = y + height - (pts[i].Z - zmin) * height / (zmax - zmin); float y2 = y + height - (pts[i + 1].Z - zmin) * height / (zmax - zmin); pta[0] = new PointF(x, y2); pta[1] = new PointF(x + width, y2); pta[2] = new PointF(x + width, y1); pta[3] = new PointF(x, y1); g.FillPolygon(aBrush, pta); } g.DrawRectangle(aPen, x, y, width, height); // Add ticks and labels to the color bar: float ticklength = 0.1f * width; for (float z = zmin; z <= zmax; z = z + (zmax - zmin) / 6) { float yy = y + height - (z - zmin) * height / (zmax - zmin); g.DrawLine(aPen, x, yy, x + ticklength, yy); g.DrawLine(aPen, x + width, yy, x + width - ticklength, yy); g.DrawString((Math.Round(z, 2)).ToString(), cl.TickFont, Brushes.Black, new PointF(x + width + 5, yy - size.Height / 2), sFormat); } } } }
private void Interp(Graphics g, ChartStyle cs, ChartStyle2D cs2d, Matrix3 m, Point3[] pta, float zmin, float zmax, int flag, Axes ax, ViewAngle va, ChartLabels cl) { SolidBrush aBrush = new SolidBrush(Color.Black); PointF[] points = new PointF[4]; int npoints = NumberInterp; Point3[,] pts = new Point3[npoints + 1, npoints + 1]; Point3[,] pts1 = new Point3[npoints + 1, npoints + 1]; float x0 = pta[0].X; float y0 = pta[0].Y; float x1 = pta[2].X; float y1 = pta[2].Y; float dx = (x1 - x0) / npoints; float dy = (y1 - y0) / npoints; float C00 = pta[0].Z; float C10 = pta[3].Z; float C11 = pta[2].Z; float C01 = pta[1].Z; float x, y, C; Color color; if (flag == 1) // For Surface chart { for (int i = 0; i <= npoints; i++) { x = x0 + i * dx; for (int j = 0; j <= npoints; j++) { y = y0 + j * dy; C = (y1 - y) * ((x1 - x) * C00 + (x - x0) * C10) / (x1 - x0) / (y1 - y0) + (y - y0) * ((x1 - x) * C01 + (x - x0) * C11) / (x1 - x0) / (y1 - y0); pts[i, j] = new Point3(x, y, C, 1); pts[i, j].Transform(m, chart3d, ax, cs, cl); } } for (int i = 0; i < npoints; i++) { for (int j = 0; j < npoints; j++) { color = AddColor(cs, pts[i, j], zmin, zmax, ax, va, cl); aBrush = new SolidBrush(color); points[0] = new PointF(pts[i, j].X, pts[i, j].Y); points[1] = new PointF(pts[i + 1, j].X, pts[i + 1, j].Y); points[2] = new PointF(pts[i + 1, j + 1].X, pts[i + 1, j + 1].Y); points[3] = new PointF(pts[i, j + 1].X, pts[i, j + 1].Y); g.FillPolygon(aBrush, points); aBrush.Dispose(); } } } else if (flag == 2) // For XYColor chart { for (int i = 0; i <= npoints; i++) { x = x0 + i * dx; for (int j = 0; j <= npoints; j++) { y = y0 + j * dy; C = (y1 - y) * ((x1 - x) * C00 + (x - x0) * C10) / (x1 - x0) / (y1 - y0) + (y - y0) * ((x1 - x) * C01 + (x - x0) * C11) / (x1 - x0) / (y1 - y0); pts[i, j] = new Point3(x, y, C, 1); } } for (int i = 0; i < npoints; i++) { for (int j = 0; j < npoints; j++) { color = AddColor(cs, pts[i, j], zmin, zmax, ax, va, cl); aBrush = new SolidBrush(color); points[0] = cs2d.Point2D(new PointF(pts[i, j].X, pts[i, j].Y), ax); points[1] = cs2d.Point2D(new PointF(pts[i + 1, j].X, pts[i + 1, j].Y), ax); points[2] = cs2d.Point2D(new PointF(pts[i + 1, j + 1].X, pts[i + 1, j + 1].Y), ax); points[3] = cs2d.Point2D(new PointF(pts[i, j + 1].X, pts[i, j + 1].Y), ax); g.FillPolygon(aBrush, points); aBrush.Dispose(); } } } else if (flag == 3) // For XYColor3D chart { for (int i = 0; i <= npoints; i++) { x = x0 + i * dx; for (int j = 0; j <= npoints; j++) { y = y0 + j * dy; C = (y1 - y) * ((x1 - x) * C00 + (x - x0) * C10) / (x1 - x0) / (y1 - y0) + (y - y0) * ((x1 - x) * C01 + (x - x0) * C11) / (x1 - x0) / (y1 - y0); pts1[i, j] = new Point3(x, y, C, 1); pts[i, j] = new Point3(x, y, ax.ZMin, 1); pts[i, j].Transform(m, chart3d, ax, cs, cl); } } for (int i = 0; i < npoints; i++) { for (int j = 0; j < npoints; j++) { color = AddColor(cs, pts1[i, j], zmin, zmax, ax, va, cl); aBrush = new SolidBrush(color); points[0] = new PointF(pts[i, j].X, pts[i, j].Y); points[1] = new PointF(pts[i + 1, j].X, pts[i + 1, j].Y); points[2] = new PointF(pts[i + 1, j + 1].X, pts[i + 1, j + 1].Y); points[3] = new PointF(pts[i, j + 1].X, pts[i, j + 1].Y); g.FillPolygon(aBrush, points); aBrush.Dispose(); } } } }
private void DrawBar(Graphics g, DataSeries ds, ChartStyle cs, Matrix3 m, Point3 pt, float xlength, float ylength, float zorign, float zmax, float zmin, Axes ax, ViewAngle va, ChartLabels cl) { Pen aPen = new Pen(ds.LineStyle.LineColor, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; Color color = AddColor(cs, pt, zmin, zmax, ax, va, cl); SolidBrush aBrush = new SolidBrush(color); Point3[] pts = new Point3[8]; Point3[] pts1 = new Point3[8]; Point3[] pt3 = new Point3[4]; PointF[] pta = new PointF[4]; pts[0] = new Point3(pt.X - xlength, pt.Y - ylength, zorign, 1); pts[1] = new Point3(pt.X - xlength, pt.Y + ylength, zorign, 1); pts[2] = new Point3(pt.X + xlength, pt.Y + ylength, zorign, 1); pts[3] = new Point3(pt.X + xlength, pt.Y - ylength, zorign, 1); pts[4] = new Point3(pt.X + xlength, pt.Y - ylength, pt.Z, 1); pts[5] = new Point3(pt.X + xlength, pt.Y + ylength, pt.Z, 1); pts[6] = new Point3(pt.X - xlength, pt.Y + ylength, pt.Z, 1); pts[7] = new Point3(pt.X - xlength, pt.Y - ylength, pt.Z, 1); for (int i = 0; i < pts.Length; i++) { pts1[i] = new Point3(pts[i].X, pts[i].Y, pts[i].Z, 1); pts[i].Transform(m, chart3d, ax, cs, cl); } int[] nconfigs = new int[8]; if (ds.BarStyle.IsBarSingleColor) { pta[0] = new PointF(pts[4].X, pts[4].Y); pta[1] = new PointF(pts[5].X, pts[5].Y); pta[2] = new PointF(pts[6].X, pts[6].Y); pta[3] = new PointF(pts[7].X, pts[7].Y); g.FillPolygon(aBrush, pta); g.DrawPolygon(aPen, pta); pta[0] = new PointF(pts[0].X, pts[0].Y); pta[1] = new PointF(pts[1].X, pts[1].Y); pta[2] = new PointF(pts[2].X, pts[2].Y); pta[3] = new PointF(pts[3].X, pts[3].Y); g.FillPolygon(aBrush, pta); g.DrawPolygon(aPen, pta); if (va.Azimuth >= -180 && va.Azimuth < -90) { nconfigs = new int[8] { 1, 2, 5, 6, 1, 0, 7, 6 }; } else if (va.Azimuth >= -90 && va.Azimuth < 0) { nconfigs = new int[8] { 1, 0, 7, 6, 0, 3, 4, 7 }; } else if (va.Azimuth >= 0 && va.Azimuth < 90) { nconfigs = new int[8] { 0, 3, 4, 7, 2, 3, 4, 5 }; } else if (va.Azimuth >= 90 && va.Azimuth < 180) { nconfigs = new int[8] { 2, 3, 4, 5, 1, 2, 5, 6 }; } pta[0] = new PointF(pts[nconfigs[0]].X, pts[nconfigs[0]].Y); pta[1] = new PointF(pts[nconfigs[1]].X, pts[nconfigs[1]].Y); pta[2] = new PointF(pts[nconfigs[2]].X, pts[nconfigs[2]].Y); pta[3] = new PointF(pts[nconfigs[3]].X, pts[nconfigs[3]].Y); g.FillPolygon(aBrush, pta); g.DrawPolygon(aPen, pta); pta[0] = new PointF(pts[nconfigs[4]].X, pts[nconfigs[4]].Y); pta[1] = new PointF(pts[nconfigs[5]].X, pts[nconfigs[5]].Y); pta[2] = new PointF(pts[nconfigs[6]].X, pts[nconfigs[6]].Y); pta[3] = new PointF(pts[nconfigs[7]].X, pts[nconfigs[7]].Y); g.FillPolygon(aBrush, pta); g.DrawPolygon(aPen, pta); } else if (!ds.BarStyle.IsBarSingleColor && IsColorMap) { pta[0] = new PointF(pts[4].X, pts[4].Y); pta[1] = new PointF(pts[5].X, pts[5].Y); pta[2] = new PointF(pts[6].X, pts[6].Y); pta[3] = new PointF(pts[7].X, pts[7].Y); g.FillPolygon(aBrush, pta); g.DrawPolygon(aPen, pta); pta[0] = new PointF(pts[0].X, pts[0].Y); pta[1] = new PointF(pts[1].X, pts[1].Y); pta[2] = new PointF(pts[2].X, pts[2].Y); pta[3] = new PointF(pts[3].X, pts[3].Y); color = AddColor(cs, pts[0], zmin, zmax, ax, va, cl); aBrush = new SolidBrush(color); g.FillPolygon(aBrush, pta); g.DrawPolygon(aPen, pta); float dz = (zmax - zmin) / 63; if (pt.Z < zorign) dz = -dz; int nz = (int)((pt.Z - zorign) / dz) + 1; if (nz < 1) nz = 1; float z = zorign; if (va.Azimuth >= -180 && va.Azimuth < -90) { nconfigs = new int[4] { 1, 2, 1, 0 }; } else if (va.Azimuth >= -90 && va.Azimuth < 0) { nconfigs = new int[4] { 1, 0, 0, 3 }; } else if (va.Azimuth >= 0 && va.Azimuth < 90) { nconfigs = new int[4] { 0, 3, 2, 3 }; } else if (va.Azimuth >= 90 && va.Azimuth <= 180) { nconfigs = new int[4] { 2, 3, 1, 2 }; } for (int i = 0; i < nz; i++) { z = zorign + i * dz; pt3[0] = new Point3(pts1[nconfigs[0]].X, pts1[nconfigs[0]].Y, z, 1); pt3[1] = new Point3(pts1[nconfigs[1]].X, pts1[nconfigs[1]].Y, z, 1); pt3[2] = new Point3(pts1[nconfigs[1]].X, pts1[nconfigs[1]].Y, z + dz, 1); pt3[3] = new Point3(pts1[nconfigs[0]].X, pts1[nconfigs[0]].Y, z + dz, 1); for (int j = 0; j < pt3.Length; j++) { pt3[j].Transform(m, chart3d, ax, cs, cl); } pta[0] = new PointF(pt3[0].X, pt3[0].Y); pta[1] = new PointF(pt3[1].X, pt3[1].Y); pta[2] = new PointF(pt3[2].X, pt3[2].Y); pta[3] = new PointF(pt3[3].X, pt3[3].Y); color = AddColor(cs, pt3[0], zmin, zmax, ax, va, cl); aBrush = new SolidBrush(color); g.FillPolygon(aBrush, pta); } pt3[0] = new Point3(pts1[nconfigs[0]].X, pts1[nconfigs[0]].Y, zorign, 1); pt3[1] = new Point3(pts1[nconfigs[1]].X, pts1[nconfigs[1]].Y, zorign, 1); pt3[2] = new Point3(pts1[nconfigs[1]].X, pts1[nconfigs[1]].Y, pt.Z, 1); pt3[3] = new Point3(pts1[nconfigs[0]].X, pts1[nconfigs[0]].Y, pt.Z, 1); for (int j = 0; j < pt3.Length; j++) { pt3[j].Transform(m, chart3d, ax, cs, cl); } pta[0] = new PointF(pt3[0].X, pt3[0].Y); pta[1] = new PointF(pt3[1].X, pt3[1].Y); pta[2] = new PointF(pt3[2].X, pt3[2].Y); pta[3] = new PointF(pt3[3].X, pt3[3].Y); g.DrawPolygon(aPen, pta); for (int i = 0; i < nz; i++) { z = zorign + i * dz; pt3[0] = new Point3(pts1[nconfigs[2]].X, pts1[nconfigs[2]].Y, z, 1); pt3[1] = new Point3(pts1[nconfigs[3]].X, pts1[nconfigs[3]].Y, z, 1); pt3[2] = new Point3(pts1[nconfigs[3]].X, pts1[nconfigs[3]].Y, z + dz, 1); pt3[3] = new Point3(pts1[nconfigs[2]].X, pts1[nconfigs[2]].Y, z + dz, 1); for (int j = 0; j < pt3.Length; j++) { pt3[j].Transform(m, chart3d, ax, cs, cl); } pta[0] = new PointF(pt3[0].X, pt3[0].Y); pta[1] = new PointF(pt3[1].X, pt3[1].Y); pta[2] = new PointF(pt3[2].X, pt3[2].Y); pta[3] = new PointF(pt3[3].X, pt3[3].Y); color = AddColor(cs, pt3[0], zmin, zmax, ax, va, cl); aBrush = new SolidBrush(color); g.FillPolygon(aBrush, pta); } pt3[0] = new Point3(pts1[nconfigs[2]].X, pts1[nconfigs[2]].Y, zorign, 1); pt3[1] = new Point3(pts1[nconfigs[3]].X, pts1[nconfigs[3]].Y, zorign, 1); pt3[2] = new Point3(pts1[nconfigs[3]].X, pts1[nconfigs[3]].Y, pt.Z, 1); pt3[3] = new Point3(pts1[nconfigs[2]].X, pts1[nconfigs[2]].Y, pt.Z, 1); for (int j = 0; j < pt3.Length; j++) { pt3[j].Transform(m, chart3d, ax, cs, cl); } pta[0] = new PointF(pt3[0].X, pt3[0].Y); pta[1] = new PointF(pt3[1].X, pt3[1].Y); pta[2] = new PointF(pt3[2].X, pt3[2].Y); pta[3] = new PointF(pt3[3].X, pt3[3].Y); g.DrawPolygon(aPen, pta); } aPen.Dispose(); aBrush.Dispose(); }
private void AddXYColor3D(Graphics g, DataSeries ds, ChartStyle cs, ChartStyle2D cs2d, Axes ax, ViewAngle va, ChartLabels cl) { Pen aPen = new Pen(ds.LineStyle.LineColor, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; SolidBrush aBrush = new SolidBrush(Color.White); PointF[] pta = new PointF[4]; Point3[,] pts = ds.PointArray; Point3[,] pts1 = new Point3[pts.GetLength(0), pts.GetLength(1)]; Point3[,] pts2 = new Point3[pts.GetLength(0), pts.GetLength(1)]; Matrix3 m = Matrix3.AzimuthElevation(va.Elevation, va.Azimuth); // Find the minumum and maximum z values: float zmin = ds.ZDataMin(); float zmax = ds.ZDataMax(); // Perform transformation on points: for (int i = 0; i < pts.GetLength(0); i++) { for (int j = 0; j < pts.GetLength(1); j++) { // Make a deep copy the points array: pts1[i, j] = new Point3(pts[i, j].X, pts[i, j].Y, ax.ZMin, 1); pts2[i, j] = new Point3(pts[i, j].X, pts[i, j].Y, ax.ZMin, 1); pts1[i, j].Transform(m, chart3d, ax, cs, cl); } } // Draw surface on the XY plane: if (!IsInterp) { for (int i = 0; i < pts.GetLength(0) - 1; i++) { for (int j = 0; j < pts.GetLength(1) - 1; j++) { pta[0] = new PointF(pts1[i, j].X, pts1[i, j].Y); pta[1] = new PointF(pts1[i, j + 1].X, pts1[i, j + 1].Y); pta[2] = new PointF(pts1[i + 1, j + 1].X, pts1[i + 1, j + 1].Y); pta[3] = new PointF(pts1[i + 1, j].X, pts1[i + 1, j].Y); Color color = AddColor(cs, pts[i, j], zmin, zmax, ax, va, cl); aBrush = new SolidBrush(color); g.FillPolygon(aBrush, pta); if (ds.LineStyle.IsVisible) { //g.DrawPolygon(aPen, pta); } } } } // Draw refined surface: else if (IsInterp) { for (int i = 0; i < pts1.GetLength(0) - 1; i++) { for (int j = 0; j < pts1.GetLength(1) - 1; j++) { Point3[] points = new Point3[4]; points[0] = pts[i, j]; points[1] = pts[i, j + 1]; points[2] = pts[i + 1, j + 1]; points[3] = pts[i + 1, j]; Interp(g, cs, cs2d, m, points, zmin, zmax, 3, ax, va, cl); } } } }
private Point3[] CoordinatesOfChartBox(Axes ax, ViewAngle va) { // Create coordinate of the axes: Point3[] pta = new Point3[8]; pta[0] = new Point3(ax.XMax, ax.YMin, ax.ZMin, 1); pta[1] = new Point3(ax.XMin, ax.YMin, ax.ZMin, 1); pta[2] = new Point3(ax.XMin, ax.YMax, ax.ZMin, 1); pta[3] = new Point3(ax.XMin, ax.YMax, ax.ZMax, 1); pta[4] = new Point3(ax.XMin, ax.YMin, ax.ZMax, 1); pta[5] = new Point3(ax.XMax, ax.YMin, ax.ZMax, 1); pta[6] = new Point3(ax.XMax, ax.YMax, ax.ZMax, 1); pta[7] = new Point3(ax.XMax, ax.YMax, ax.ZMin, 1); Point3[] pts = new Point3[4]; int[] npts = new int[4] { 0, 1, 2, 3 }; if (va.Elevation >= 0) { if (va.Azimuth >= -180 && va.Azimuth < -90) { npts = new int[4] { 1, 2, 7, 6 }; } else if (va.Azimuth >= -90 && va.Azimuth < 0) { npts = new int[4] { 0, 1, 2, 3 }; } else if (va.Azimuth >= 0 && va.Azimuth < 90) { npts = new int[4] { 7, 0, 1, 4 }; } else if (va.Azimuth >= 90 && va.Azimuth <= 180) { npts = new int[4] { 2, 7, 0, 5 }; } } else if (va.Elevation < 0) { if (va.Azimuth >= -180 && va.Azimuth < -90) { npts = new int[4] { 1, 0, 7, 6 }; } else if (va.Azimuth >= -90 && va.Azimuth < 0) { npts = new int[4] { 0, 7, 2, 3 }; } else if (va.Azimuth >= 0 && va.Azimuth < 90) { npts = new int[4] { 7, 2, 1, 4 }; } else if (va.Azimuth >= 90 && va.Azimuth <= 180) { npts = new int[4] { 2, 1, 0, 5 }; } } for (int i = 0; i < 4; i++) { pts[i] = pta[npts[i]]; } return pts; }
public void SinROverR3D(DataSeries ds, Axes ax) { ax.XMin = -8; ax.XMax = 8; ax.YMin = -8; ax.YMax = 8; ax.ZMin = -0.5f; ax.ZMax = 1; ax.XTick = 4; ax.YTick = 4; ax.ZTick = 0.5f; ds.XDataMin = ax.XMin; ds.YDataMin = ax.YMin; ds.XSpacing = 0.5f; ds.YSpacing = 0.5f; ds.XNumber = Convert.ToInt16((ax.XMax - ax.XMin) / ds.XSpacing) + 1; ds.YNumber = Convert.ToInt16((ax.YMax - ax.YMin) / ds.YSpacing) + 1; Point3[,] pts = new Point3[ds.XNumber, ds.YNumber]; for (int i = 0; i < ds.XNumber; i++) { for (int j = 0; j < ds.YNumber; j++) { float x = ds.XDataMin + i * ds.XSpacing; float y = ds.YDataMin + j * ds.YSpacing; float r = (float)Math.Sqrt(x * x + y * y) + 0.000001f; float z = (float)Math.Sin(r) / r; pts[i, j] = new Point3(x, y, z, 1); } } ds.PointArray = pts; }
private void AddMeshZ(Graphics g, DataSeries ds, ChartStyle cs, Axes ax, ViewAngle va, ChartLabels cl) { Pen aPen = new Pen(ds.LineStyle.LineColor, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; SolidBrush aBrush = new SolidBrush(Color.White); Matrix3 m = Matrix3.AzimuthElevation(va.Elevation, va.Azimuth); PointF[] pta = new PointF[4]; Point3[,] pts = ds.PointArray; Point3[,] pts1 = new Point3[pts.GetLength(0), pts.GetLength(1)]; Color color; // Find the minumum and maximum z values: float zmin = ds.ZDataMin(); float zmax = ds.ZDataMax(); for (int i = 0; i < pts.GetLength(0); i++) { for (int j = 0; j < pts.GetLength(1); j++) { // Make a deep copy the points array: pts1[i, j] = new Point3(pts[i, j].X, pts[i, j].Y, pts[i, j].Z, 1); // Perform transformations on points: pts[i, j].Transform(m, chart3d, ax, cs, cl); } } for (int i = 0; i < pts.GetLength(0) - 1; i++) { for (int j = 0; j < pts.GetLength(1) - 1; j++) { int ii = i; if (va.Azimuth >= -180 && va.Azimuth < 0) { ii = pts.GetLength(0) - 2 - i; } pta[0] = new PointF(pts[ii, j].X, pts[ii, j].Y); pta[1] = new PointF(pts[ii, j + 1].X, pts[ii, j + 1].Y); pta[2] = new PointF(pts[ii + 1, j + 1].X, pts[ii + 1, j + 1].Y); pta[3] = new PointF(pts[ii + 1, j].X, pts[ii + 1, j].Y); g.FillPolygon(aBrush, pta); if (IsColorMap) { color = AddColor(cs, pts[ii, j], zmin, zmax, ax, va, cl); aPen = new Pen(color, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; } g.DrawPolygon(aPen, pta); } } Point3[] pt3 = new Point3[4]; for (int i = 0; i < pts1.GetLength(0); i++) { int jj = pts1.GetLength(0) - 1; if (va.Elevation >= 0) { if (va.Azimuth >= -90 && va.Azimuth <= 90) jj = 0; } else if (va.Elevation < 0) { jj = 0; if (va.Azimuth >= -90 && va.Azimuth <= 90) jj = pts1.GetLength(0) - 1; } if (i < pts1.GetLength(0) - 1) { pt3[0] = new Point3(pts1[i, jj].X, pts1[i, jj].Y, pts1[i, jj].Z, 1); pt3[1] = new Point3(pts1[i + 1, jj].X, pts1[i + 1, jj].Y, pts1[i + 1, jj].Z, 1); pt3[2] = new Point3(pts1[i + 1, jj].X, pts1[i + 1, jj].Y, ax.ZMin, 1); pt3[3] = new Point3(pts1[i, jj].X, pts1[i, jj].Y, ax.ZMin, 1); for (int k = 0; k < 4; k++) pt3[k].Transform(m, chart3d, ax, cs, cl); pta[0] = new PointF(pt3[0].X, pt3[0].Y); pta[1] = new PointF(pt3[1].X, pt3[1].Y); pta[2] = new PointF(pt3[2].X, pt3[2].Y); pta[3] = new PointF(pt3[3].X, pt3[3].Y); g.FillPolygon(aBrush, pta); if (IsColorMap) { color = AddColor(cs, pt3[0], zmin, zmax, ax, va, cl); aPen = new Pen(color, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; } g.DrawPolygon(aPen, pta); } } for (int j = 0; j < pts1.GetLength(1); j++) { int ii = 0; if (va.Elevation >= 0) { if (va.Azimuth >= 0 && va.Azimuth <= 180) ii = pts1.GetLength(1) - 1; } else if (va.Elevation < 0) { if (va.Azimuth >= -180 && va.Azimuth <= 0) ii = pts1.GetLength(1) - 1; } if (j < pts1.GetLength(1) - 1) { pt3[0] = new Point3(pts1[ii, j].X, pts1[ii, j].Y, pts1[ii, j].Z, 1); pt3[1] = new Point3(pts1[ii, j + 1].X, pts1[ii, j + 1].Y, pts1[ii, j + 1].Z, 1); pt3[2] = new Point3(pts1[ii, j + 1].X, pts1[ii, j + 1].Y, ax.ZMin, 1); pt3[3] = new Point3(pts1[ii, j].X, pts1[ii, j].Y, ax.ZMin, 1); for (int k = 0; k < 4; k++) pt3[k].Transform(m, chart3d, ax, cs, cl); pta[0] = new PointF(pt3[0].X, pt3[0].Y); pta[1] = new PointF(pt3[1].X, pt3[1].Y); pta[2] = new PointF(pt3[2].X, pt3[2].Y); pta[3] = new PointF(pt3[3].X, pt3[3].Y); g.FillPolygon(aBrush, pta); if (IsColorMap) { color = AddColor(cs, pt3[0], zmin, zmax, ax, va, cl); aPen = new Pen(color, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; } g.DrawPolygon(aPen, pta); } } aPen.Dispose(); aBrush.Dispose(); }
public Point3 Cylindrical(float r, float theta, float y) { Point3 pt = new Point3(); float sn = (float)Math.Sin(theta * Math.PI / 180); float cn = (float)Math.Cos(theta * Math.PI / 180); pt.X = r * cn; pt.Y = y; pt.Z = -r * sn; pt.W = 1; return pt; }
public void AddPoint(Point3 pt) { PointList.Add(pt); }
public Point4(Point3 pt3, float v) { point3 = pt3; V = v; }
private void AddWaterfall(Graphics g, DataSeries ds, ChartStyle cs, Axes ax, ViewAngle va, ChartLabels cl) { Pen aPen = new Pen(ds.LineStyle.LineColor, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; SolidBrush aBrush = new SolidBrush(Color.White); Matrix3 m = Matrix3.AzimuthElevation(va.Elevation, va.Azimuth); Point3[,] pts = ds.PointArray; Point3[] pt3 = new Point3[pts.GetLength(0) + 2]; PointF[] pta = new PointF[pts.GetLength(0) + 2]; Color color; // Find the minumum and maximum z values: float zmin = ds.ZDataMin(); float zmax = ds.ZDataMax(); for (int j = 0; j < pts.GetLength(1); j++) { int jj = j; if (va.Elevation >= 0) { if (va.Azimuth >= -90 && va.Azimuth < 90) { jj = pts.GetLength(1) - 1 - j; } } else if (va.Elevation < 0) { jj = pts.GetLength(1) - 1 - j; if (va.Azimuth >= -90 && va.Azimuth < 90) { jj = j; } } for (int i = 0; i < pts.GetLength(0); i++) { pt3[i + 1] = pts[i, jj]; if (i == 0) { pt3[0] = new Point3(pt3[i + 1].X, pt3[i + 1].Y, ax.ZMin, 1); } if (i == pts.GetLength(0) - 1) { pt3[pts.GetLength(0) + 1] = new Point3(pt3[i + 1].X, pt3[i + 1].Y, ax.ZMin, 1); } } for (int i = 0; i < pt3.Length; i++) { pt3[i].Transform(m, chart3d, ax, cs, cl); pta[i] = new PointF(pt3[i].X, pt3[i].Y); } g.FillPolygon(aBrush, pta); for (int i = 1; i < pt3.Length; i++) { if (IsColorMap) { color = AddColor(cs, pt3[i], zmin, zmax, ax, va, cl); aPen = new Pen(color, ds.LineStyle.Thickness); aPen.DashStyle = ds.LineStyle.Pattern; } g.DrawLine(aPen, pta[i - 1], pta[i]); } } aPen.Dispose(); aBrush.Dispose(); }