source: branches/2817-BinPackingSpeedup/HeuristicLab.Problems.BinPacking.Views/3.3/Container3DView.xaml.cs @ 15307

Last change on this file since 15307 was 15307, checked in by abeham, 4 years ago

#2817:

  • Added checkbox to control showing extreme points in visualization
    • Automatically determine size of extreme point cubes
  • Fixed some bugs in extreme point generation
File size: 20.2 KB
Line 
1/* HeuristicLab
2 * Copyright (C) 2002-2016 Joseph Helm and Heuristic and Evolutionary Algorithms Laboratory (HEAL)
3 *
4 * This file is part of HeuristicLab.
5 *
6 * HeuristicLab is free software: you can redistribute it and/or modify
7 * it under the terms of the GNU General Public License as published by
8 * the Free Software Foundation, either version 3 of the License, or
9 * (at your option) any later version.
10 *
11 * HeuristicLab is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14 * GNU General Public License for more details.
15 *
16 * You should have received a copy of the GNU General Public License
17 * along with HeuristicLab. If not, see<http://www.gnu.org/licenses/> .
18 */
19
20using System;
21using System.Collections.Generic;
22using System.Linq;
23using System.Windows;
24using System.Windows.Controls;
25using System.Windows.Input;
26using System.Windows.Media;
27using System.Windows.Media.Media3D;
28using HeuristicLab.Problems.BinPacking3D;
29
30namespace HeuristicLab.Problems.BinPacking.Views {
31  public partial class Container3DView : UserControl {
32    private static readonly Color[] colors = new[] {
33      Color.FromRgb(0x40, 0x6A, 0xB7),
34      Color.FromRgb(0xB1, 0x6D, 0x01),
35      Color.FromRgb(0x4E, 0x8A, 0x06),
36      Color.FromRgb(0x75, 0x50, 0x7B),
37      Color.FromRgb(0x72, 0x9F, 0xCF),
38      Color.FromRgb(0xA4, 0x00, 0x00),
39      Color.FromRgb(0xAD, 0x7F, 0xA8),
40      Color.FromRgb(0x29, 0x50, 0xCF),
41      Color.FromRgb(0x90, 0xB0, 0x60),
42      Color.FromRgb(0xF5, 0x89, 0x30),
43      Color.FromRgb(0x55, 0x57, 0x53),
44      Color.FromRgb(0xEF, 0x59, 0x59),
45      Color.FromRgb(0xED, 0xD4, 0x30),
46      Color.FromRgb(0x63, 0xC2, 0x16),
47    };
48
49    private static readonly Color hiddenColor = Color.FromArgb(0x1A, 0xAA, 0xAA, 0xAA);
50    private static readonly Color containerColor = Color.FromArgb(0x7F, 0xAA, 0xAA, 0xAA);
51
52    private Point startPos;
53    private bool mouseDown = false;
54    private double startAngleX;
55    private double startAngleY;
56    private int selectedItemKey;
57    private bool showExtremePoints;
58
59    private BinPacking<BinPacking3D.PackingPosition, PackingShape, PackingItem> packing;
60    public BinPacking<BinPacking3D.PackingPosition, PackingShape, PackingItem> Packing {
61      get { return packing; }
62      set {
63        if (packing != value) {
64          this.packing = value;
65          ClearSelection(); // also updates visualization
66        }
67      }
68    }
69
70    private Dictionary<int, DiffuseMaterial> materials;
71
72    public Container3DView() {
73      InitializeComponent();
74      camMain.Position = new Point3D(0.5, 3, 3); // for design time we use a different camera position
75      materials = new Dictionary<int, DiffuseMaterial>();
76      Clear();
77    }
78
79
80    protected override void OnRenderSizeChanged(SizeChangedInfo sizeInfo) {
81      base.OnRenderSizeChanged(sizeInfo);
82      var s = Math.Min(sizeInfo.NewSize.Height, sizeInfo.NewSize.Width);
83      var mySize = new Size(s, s);
84      viewport3D1.RenderSize = mySize;
85    }
86
87    public void SelectItem(int itemKey) {
88      // selection of an item should make all other items semi-transparent
89      selectedItemKey = itemKey;
90      UpdateVisualization();
91    }
92    public void ClearSelection() {
93      // remove all transparency
94      selectedItemKey = -1;
95      UpdateVisualization();
96    }
97
98    private void UpdateVisualization() {
99      Clear();
100      if (packing == null) return; // nothing to display
101
102      var modelGroup = (Model3DGroup)MyModel.Content;
103      var hiddenMaterial = new DiffuseMaterial(new SolidColorBrush(hiddenColor));
104
105      if (selectedItemKey >= 0) {
106        var selectedItem = packing.Items.Single(x => selectedItemKey == x.Key);
107        var selectedPos = packing.Positions[selectedItem.Key];
108
109        var colorIdx = selectedItem.Value.Material;
110        while (colorIdx < 0) colorIdx += colors.Length;
111        colorIdx = colorIdx % colors.Length;
112        var color = colors[colorIdx];
113        var material = new DiffuseMaterial { Brush = new SolidColorBrush(color) };
114        materials[selectedItem.Value.Material] = material;
115
116        var selectedModel = new GeometryModel3D { Geometry = new MeshGeometry3D(), Material = material };
117        AddSolidCube((MeshGeometry3D)selectedModel.Geometry, selectedPos.X, selectedPos.Y, selectedPos.Z,
118          selectedPos.Rotated ? selectedItem.Value.Depth : selectedItem.Value.Width,
119          selectedItem.Value.Height,
120          selectedPos.Rotated ? selectedItem.Value.Width : selectedItem.Value.Depth);
121        modelGroup.Children.Add(selectedModel);
122
123        foreach (var item in packing.Items.Where(x => selectedItemKey != x.Key)) {
124          var position = packing.Positions[item.Key];
125
126          var w = position.Rotated ? item.Value.Depth : item.Value.Width;
127          var h = item.Value.Height;
128          var d = position.Rotated ? item.Value.Width : item.Value.Depth;
129
130          var model = new GeometryModel3D { Geometry = new MeshGeometry3D(), Material = hiddenMaterial };
131          AddWireframeCube((MeshGeometry3D)model.Geometry, position.X, position.Y, position.Z, w, h, d, 1);
132          modelGroup.Children.Add(model);
133        }
134      } else {
135        foreach (var item in packing.Items) {
136          var position = packing.Positions[item.Key];
137
138          var w = position.Rotated ? item.Value.Depth : item.Value.Width;
139          var h = item.Value.Height;
140          var d = position.Rotated ? item.Value.Width : item.Value.Depth;
141
142          var model = new GeometryModel3D { Geometry = new MeshGeometry3D() };
143          DiffuseMaterial material;
144          if (!materials.TryGetValue(item.Value.Material, out material)) {
145            var colorIdx = item.Value.Material;
146            while (colorIdx < 0) colorIdx += colors.Length;
147            colorIdx = colorIdx % colors.Length;
148            var color = colors[colorIdx];
149            material = new DiffuseMaterial { Brush = new SolidColorBrush(color) };
150            materials[item.Value.Material] = material;
151          }
152          var selectedModel = new GeometryModel3D { Geometry = new MeshGeometry3D(), Material = material };
153          AddSolidCube((MeshGeometry3D)selectedModel.Geometry, position.X, position.Y, position.Z, w, h, d);
154          modelGroup.Children.Add(selectedModel);
155        }
156      }
157
158      if (showExtremePoints) {
159        // draw extreme-points
160        var maxMag = (int)Math.Log10(Math.Max(packing.BinShape.Depth, Math.Max(packing.BinShape.Height, packing.BinShape.Width)));
161        var cubeSize = (int)Math.Max(Math.Pow(10, maxMag - 2), 1);
162        foreach (var ep in packing.ExtremePoints) {
163          var epModel = new GeometryModel3D { Geometry = new MeshGeometry3D(), Material = new DiffuseMaterial() { Brush = new SolidColorBrush(Colors.Red) } };
164          AddSolidCube((MeshGeometry3D)epModel.Geometry, ep.X, ep.Y, ep.Z, cubeSize, cubeSize, cubeSize);
165          modelGroup.Children.Add(epModel);
166        }
167      }
168
169      var container = packing.BinShape;
170      var containerModel = new GeometryModel3D(new MeshGeometry3D(), new DiffuseMaterial(new SolidColorBrush(containerColor)));
171      modelGroup.Children.Add(containerModel);
172      AddWireframeCube((MeshGeometry3D)containerModel.Geometry, container.Origin.X - .5, container.Origin.Y - .5, container.Origin.Z - .5, container.Width + 1, container.Height + 1, container.Depth + 1);
173
174      var ratio = Math.Max(container.Width, Math.Max(container.Height, container.Depth));
175      scale.ScaleX = 1.0 / ratio;
176      scale.ScaleY = 1.0 / ratio;
177      scale.ScaleZ = 1.0 / ratio;
178
179      scaleZoom.CenterX = rotateX.CenterX = rotateY.CenterX = container.Width / (2.0 * ratio);
180      scaleZoom.CenterY = rotateX.CenterY = rotateY.CenterY = container.Height / (2.0 * ratio);
181      scaleZoom.CenterZ = rotateX.CenterZ = rotateY.CenterZ = container.Depth / (2.0 * ratio);
182
183      camMain.Position = new Point3D(
184        scaleZoom.CenterX,
185        3,
186        3);
187      camMain.LookDirection = new Vector3D(
188        0,
189        scaleZoom.CenterY - camMain.Position.Y,
190        scaleZoom.CenterZ - camMain.Position.Z);
191    }
192
193
194    private void Clear() {
195      ((Model3DGroup)MyModel.Content).Children.Clear();
196      materials.Clear();
197
198      mouseDown = false;
199      startAngleX = 0;
200      startAngleY = 0;
201    }
202
203    private void Container3DView_MouseMove(object sender, MouseEventArgs e) {
204      if (!mouseDown) return;
205      var pos = e.GetPosition((IInputElement)this);
206
207      ((AxisAngleRotation3D)rotateX.Rotation).Angle = startAngleX + (pos.X - startPos.X) / 4;
208      ((AxisAngleRotation3D)rotateY.Rotation).Angle = startAngleY + (pos.Y - startPos.Y) / 4;
209    }
210
211    private void Container3DView_MouseDown(object sender, MouseButtonEventArgs e) {
212      startAngleX = ((AxisAngleRotation3D)rotateX.Rotation).Angle;
213      startAngleY = ((AxisAngleRotation3D)rotateY.Rotation).Angle;
214      this.startPos = e.GetPosition((IInputElement)this);
215      this.mouseDown = true;
216    }
217
218    private void Container3DView_MouseUp(object sender, MouseButtonEventArgs e) {
219      mouseDown = false;
220    }
221
222    private void Container3DView_OnMouseWheel(object sender, MouseWheelEventArgs e) {
223      if (e.Delta > 0) {
224        scaleZoom.ScaleX *= 1.1;
225        scaleZoom.ScaleY *= 1.1;
226        scaleZoom.ScaleZ *= 1.1;
227      } else if (e.Delta < 0) {
228        scaleZoom.ScaleX /= 1.1;
229        scaleZoom.ScaleY /= 1.1;
230        scaleZoom.ScaleZ /= 1.1;
231      }
232    }
233
234    private void Container3DView_OnMouseEnter(object sender, MouseEventArgs e) {
235      Focus(); // for mouse wheel events
236    }
237    private void showExtremePointsCheckBoxOnChecked(object sender, RoutedEventArgs e) {
238      showExtremePoints = true;
239      UpdateVisualization();
240    }
241    private void showExtremePointsCheckBoxOnUnchecked(object sender, RoutedEventArgs e) {
242      showExtremePoints = false;
243      UpdateVisualization();
244    }
245
246    #region helper for cubes
247    /// <summary>
248    /// Creates a solid cube by adding the respective points and triangles.
249    /// </summary>
250    /// <param name="mesh">The mesh to which points and triangles are added.</param>
251    /// <param name="x">The leftmost point</param>
252    /// <param name="y">The frontmost point</param>
253    /// <param name="z">The lowest point</param>
254    /// <param name="width">The extension to the right</param>
255    /// <param name="height">The extension to the back</param>
256    /// <param name="depth">The extension to the top</param>
257    private void AddSolidCube(MeshGeometry3D mesh, int x, int y, int z, int width, int height, int depth) {
258      // ground
259      mesh.Positions.Add(new Point3D(x, y, z));
260      mesh.Positions.Add(new Point3D(x + width, y, z));
261      mesh.Positions.Add(new Point3D(x + width, y + height, z));
262      mesh.Positions.Add(new Point3D(x, y + height, z));
263      // top
264      mesh.Positions.Add(new Point3D(x, y, z + depth));
265      mesh.Positions.Add(new Point3D(x + width, y, z + depth));
266      mesh.Positions.Add(new Point3D(x + width, y + height, z + depth));
267      mesh.Positions.Add(new Point3D(x, y + height, z + depth));
268
269      // front
270      AddPlane(mesh, 0, 1, 5, 4);
271      // right side
272      AddPlane(mesh, 1, 2, 6, 5);
273      // back
274      AddPlane(mesh, 3, 7, 6, 2);
275      // left side
276      AddPlane(mesh, 0, 4, 7, 3);
277      // top
278      AddPlane(mesh, 4, 5, 6, 7);
279      // bottom
280      AddPlane(mesh, 0, 3, 2, 1);
281    }
282
283    /// <summary>
284    /// Creates a wireframe cube by adding the respective points and triangles.
285    /// </summary>
286    /// <param name="mesh">The mesh to which points and triangles are added.</param>
287    /// <param name="x">The leftmost point</param>
288    /// <param name="y">The frontmost point</param>
289    /// <param name="z">The lowest point</param>
290    /// <param name="width">The extension to the right</param>
291    /// <param name="height">The extension to the back</param>
292    /// <param name="depth">The extension to the top</param>
293    /// <param name="thickness">The thickness of the frame</param>
294    private void AddWireframeCube(MeshGeometry3D mesh, double x, double y, double z, double width, double height, double depth, double thickness = double.NaN) {
295      // default thickness of the wireframe is 5% of smallest dimension
296      if (double.IsNaN(thickness))
297        thickness = Math.Min(width, Math.Min(height, depth)) * 0.05;
298
299      // The cube contains of 8 corner, each corner has 4 points:
300      // 1. The corner point
301      // 2. A point on the edge to the right of the corner
302      // 3. A point on the edge atop or below the corner
303      // 4. A point on the edge to the left of the corner
304
305      // Point 0, Front Left Bottom
306      mesh.Positions.Add(new Point3D(x, y, z));
307      mesh.Positions.Add(new Point3D(x + thickness, y, z));
308      mesh.Positions.Add(new Point3D(x, y, z + thickness));
309      mesh.Positions.Add(new Point3D(x, y + thickness, z));
310
311      // Point 1, Front Right Bottom
312      mesh.Positions.Add(new Point3D(x + width, y, z));
313      mesh.Positions.Add(new Point3D(x + width, y + thickness, z));
314      mesh.Positions.Add(new Point3D(x + width, y, z + thickness));
315      mesh.Positions.Add(new Point3D(x + width - thickness, y, z));
316
317      // Point 2, Back Right Bottom
318      mesh.Positions.Add(new Point3D(x + width, y + height, z));
319      mesh.Positions.Add(new Point3D(x + width - thickness, y + height, z));
320      mesh.Positions.Add(new Point3D(x + width, y + height, z + thickness));
321      mesh.Positions.Add(new Point3D(x + width, y + height - thickness, z));
322
323      // Point 3, Back Left Bottom
324      mesh.Positions.Add(new Point3D(x, y + height, z));
325      mesh.Positions.Add(new Point3D(x, y + height - thickness, z));
326      mesh.Positions.Add(new Point3D(x, y + height, z + thickness));
327      mesh.Positions.Add(new Point3D(x + thickness, y + height, z));
328
329      // Point 4, Front Left Top
330      mesh.Positions.Add(new Point3D(x, y, z + depth));
331      mesh.Positions.Add(new Point3D(x + thickness, y, z + depth));
332      mesh.Positions.Add(new Point3D(x, y, z + depth - thickness));
333      mesh.Positions.Add(new Point3D(x, y + thickness, z + depth));
334
335      // Point 5, Front Right Top
336      mesh.Positions.Add(new Point3D(x + width, y, z + depth));
337      mesh.Positions.Add(new Point3D(x + width, y + thickness, z + depth));
338      mesh.Positions.Add(new Point3D(x + width, y, z + depth - thickness));
339      mesh.Positions.Add(new Point3D(x + width - thickness, y, z + depth));
340
341      // Point 6, Back Right Top
342      mesh.Positions.Add(new Point3D(x + width, y + height, z + depth));
343      mesh.Positions.Add(new Point3D(x + width - thickness, y + height, z + depth));
344      mesh.Positions.Add(new Point3D(x + width, y + height, z + depth - thickness));
345      mesh.Positions.Add(new Point3D(x + width, y + height - thickness, z + depth));
346
347      // Point 7, Back Left Top
348      mesh.Positions.Add(new Point3D(x, y + height, z + depth));
349      mesh.Positions.Add(new Point3D(x, y + height - thickness, z + depth));
350      mesh.Positions.Add(new Point3D(x, y + height, z + depth - thickness));
351      mesh.Positions.Add(new Point3D(x + thickness, y + height, z + depth));
352
353      // Point 0, non-edge
354      mesh.Positions.Add(new Point3D(x + thickness, y, z + thickness));
355      mesh.Positions.Add(new Point3D(x, y + thickness, z + thickness));
356      mesh.Positions.Add(new Point3D(x + thickness, y + thickness, z));
357
358      // Point 1, non-edge
359      mesh.Positions.Add(new Point3D(x + width, y + thickness, z + thickness));
360      mesh.Positions.Add(new Point3D(x + width - thickness, y, z + thickness));
361      mesh.Positions.Add(new Point3D(x + width - thickness, y + thickness, z));
362
363      // Point 2, non-edge
364      mesh.Positions.Add(new Point3D(x + width - thickness, y + height, z + thickness));
365      mesh.Positions.Add(new Point3D(x + width, y + height - thickness, z + thickness));
366      mesh.Positions.Add(new Point3D(x + width - thickness, y + height - thickness, z));
367
368      // Point 3, non-edge
369      mesh.Positions.Add(new Point3D(x, y + height - thickness, z + thickness));
370      mesh.Positions.Add(new Point3D(x + thickness, y + height, z + thickness));
371      mesh.Positions.Add(new Point3D(x + thickness, y + height - thickness, z));
372
373      // Point 4, non-edge
374      mesh.Positions.Add(new Point3D(x + thickness, y, z + depth - thickness));
375      mesh.Positions.Add(new Point3D(x, y + thickness, z + depth - thickness));
376      mesh.Positions.Add(new Point3D(x + thickness, y + thickness, z + depth));
377
378      // Point 5, non-edge
379      mesh.Positions.Add(new Point3D(x + width, y + thickness, z + depth - thickness));
380      mesh.Positions.Add(new Point3D(x + width - thickness, y, z + depth - thickness));
381      mesh.Positions.Add(new Point3D(x + width - thickness, y + thickness, z + depth));
382
383      // Point 6, non-edge
384      mesh.Positions.Add(new Point3D(x + width - thickness, y + height, z + depth - thickness));
385      mesh.Positions.Add(new Point3D(x + width, y + height - thickness, z + depth - thickness));
386      mesh.Positions.Add(new Point3D(x + width - thickness, y + height - thickness, z + depth));
387
388      // Point 7, non-edge
389      mesh.Positions.Add(new Point3D(x, y + height - thickness, z + depth - thickness));
390      mesh.Positions.Add(new Point3D(x + thickness, y + height, z + depth - thickness));
391      mesh.Positions.Add(new Point3D(x + thickness, y + height - thickness, z + depth));
392
393      // Draw the 24 corner plates
394      AddPlane(mesh, 0, 1, 32, 2);
395      AddPlane(mesh, 0, 2, 33, 3);
396      AddPlane(mesh, 0, 3, 34, 1);
397
398      AddPlane(mesh, 4, 6, 36, 7);
399      AddPlane(mesh, 4, 5, 35, 6);
400      AddPlane(mesh, 4, 7, 37, 5);
401
402      AddPlane(mesh, 8, 10, 39, 11);
403      AddPlane(mesh, 8, 9, 38, 10);
404      AddPlane(mesh, 8, 11, 40, 9);
405
406      AddPlane(mesh, 12, 13, 41, 14);
407      AddPlane(mesh, 12, 14, 42, 15);
408      AddPlane(mesh, 12, 15, 43, 13);
409
410      AddPlane(mesh, 16, 18, 44, 17);
411      AddPlane(mesh, 16, 19, 45, 18);
412      AddPlane(mesh, 16, 17, 46, 19);
413
414      AddPlane(mesh, 20, 23, 48, 22);
415      AddPlane(mesh, 20, 22, 47, 21);
416      AddPlane(mesh, 20, 21, 49, 23);
417
418      AddPlane(mesh, 24, 27, 51, 26);
419      AddPlane(mesh, 24, 26, 50, 25);
420      AddPlane(mesh, 24, 25, 52, 27);
421
422      AddPlane(mesh, 28, 31, 54, 30);
423      AddPlane(mesh, 28, 30, 53, 29);
424      AddPlane(mesh, 28, 29, 55, 31);
425
426      // Draw the connecting plates
427      // on the bottom
428      AddPlane(mesh, 1, 7, 36, 32);
429      AddPlane(mesh, 1, 34, 37, 7);
430
431      AddPlane(mesh, 5, 11, 39, 35);
432      AddPlane(mesh, 5, 37, 40, 11);
433
434      AddPlane(mesh, 9, 15, 42, 38);
435      AddPlane(mesh, 9, 40, 43, 15);
436
437      AddPlane(mesh, 13, 3, 33, 41);
438      AddPlane(mesh, 13, 43, 34, 3);
439
440      // between bottom and top
441      AddPlane(mesh, 2, 32, 44, 18);
442      AddPlane(mesh, 2, 18, 45, 33);
443
444      AddPlane(mesh, 6, 22, 48, 36);
445      AddPlane(mesh, 6, 35, 47, 22);
446
447      AddPlane(mesh, 10, 26, 51, 39);
448      AddPlane(mesh, 10, 38, 50, 26);
449
450      AddPlane(mesh, 14, 30, 54, 42);
451      AddPlane(mesh, 14, 41, 53, 30);
452
453      // on the top
454      AddPlane(mesh, 17, 44, 48, 23);
455      AddPlane(mesh, 17, 23, 49, 46);
456
457      AddPlane(mesh, 21, 47, 51, 27);
458      AddPlane(mesh, 21, 27, 52, 49);
459
460      AddPlane(mesh, 25, 50, 54, 31);
461      AddPlane(mesh, 25, 31, 55, 52);
462
463      AddPlane(mesh, 29, 19, 46, 55);
464      AddPlane(mesh, 29, 53, 45, 19);
465    }
466
467    /// <summary>
468    /// Adds a plane by two triangles. The indices of the points have to be given
469    /// in counter-clockwise sequence.
470    /// </summary>
471    /// <param name="mesh">The mesh to add the triangles to</param>
472    /// <param name="a">The index of the first point</param>
473    /// <param name="b">The index of the second point</param>
474    /// <param name="c">The index of the third point</param>
475    /// <param name="d">The index of the fourth point</param>
476    private void AddPlane(MeshGeometry3D mesh, int a, int b, int c, int d) {
477      // two triangles form a plane
478      mesh.TriangleIndices.Add(a);
479      mesh.TriangleIndices.Add(b);
480      mesh.TriangleIndices.Add(d);
481      mesh.TriangleIndices.Add(c);
482      mesh.TriangleIndices.Add(d);
483      mesh.TriangleIndices.Add(b);
484    }
485    #endregion
486  }
487}
Note: See TracBrowser for help on using the repository browser.