IEnumerator AddLines(List <Reading> values, List <DateTime> dates, int track, List <int> allIndices, List <Color> allColours, List <Vector3> allVerts)
{
tracksToUpdate -= 1;
double max = values.Max(val => val.value);
double min = values.Min(val => val.value);
DateTime start = dates.Min();
DateTime finsih = dates.Max();
float theta;
System.Drawing.Color col;
List <Vector3> tempVerts;
for (int i = 0; i < values.Count; i++)
{
theta = getTheta(dates[i], start, finsih - start);
col = cMap.GetColor((double)theta);
getLineOnTrack(values[i].value, max, min, theta * 2 * Mathf.PI, track, Vector3.zero, out tempVerts);
for (int j = 0; j < tempVerts.Count; j++)
{
allIndices.Add(allVerts.Count + j);
allColours.Add(new Color(col.R * 1.0f / 255, col.G * 1.0f / 255, col.B * 1.0f / 255));
}
allVerts.AddRange(tempVerts);
yield return(null);
}
tracksToUpdate += 1;
UnityEngine.Debug.Log("Added mesh line");
}
private void ShowSurface2D(Func <double, double, double> function)
{
var bmp = new Bitmap(SurfaceRect.Width, SurfaceRect.Height);
for (var x = SurfaceRect.Left; x < SurfaceRect.Right; x++)
{
for (var y = SurfaceRect.Top; y < SurfaceRect.Bottom; y++)
{
var z = function(x, y);
bmp.SetPixel(x - SurfaceRect.Left, y - SurfaceRect.Top, _cmap.GetColor(z));
}
}
_surfacePanel.BackgroundImage = bmp;
}
/// <inheritdoc />
public void AttachLayerData(float[] values, ColorMap colormap)
{
int length = values.Length;
double min = double.MaxValue;
double max = double.MinValue;
// Find min and max of given values.
for (int i = 0; i < length; i++)
{
if (values[i] < min)
{
min = values[i];
}
if (values[i] > max)
{
max = values[i];
}
}
var colorData = new Color[length];
for (int i = 0; i < length; i++)
{
colorData[i] = colormap.GetColor(values[i], min, max);
}
colors = new ColorBuffer(colorData, 1);
}
public void GetColor_should_interpolate_between_black_and_the_entries()
{
var map = new ColorMap(new ColorMapEntry(0.5, Colors.White), new ColorMapEntry(0.75, Colors.Red));
map.GetColor(0).Should().Be(Colors.Black);
map.GetColor(0.1).Should().Be(new Color(0.2, 0.2, 0.2));
map.GetColor(0.3).Should().Be(new Color(0.6, 0.6, 0.6));
map.GetColor(0.5).Should().Be(Colors.White);
map.GetColor(0.65).Should().Be(new Color(1, 0.4, 0.4));
map.GetColor(0.75).Should().Be(Colors.Red);
map.GetColor(0.85).Should().Be(new Color(0.6, 0, 0));
map.GetColor(1).Should().Be(Colors.Black);
}
public static void GenerateStretchRep(int resolution, float thickness, float offset, List <float> lengths, List <float> energies, Mesh mesh)
{
//Get colors
List <Color> energyColors = new List <Color>();
float minEnergy = energies[0];
float maxEnergy = energies[0];
float energy;
for (int energyNum = 1; energyNum < energies.Count; energyNum++)
{
energy = energies[energyNum];
if (energy < minEnergy)
{
minEnergy = energy;
}
if (energy > maxEnergy)
{
maxEnergy = energy;
}
}
for (int energyNum = 0; energyNum < energies.Count; energyNum++)
{
energy = energies[energyNum];
energyColors.Add(ColorMap.GetColor((energy - minEnergy) / (maxEnergy - minEnergy)));
}
List <Vector3> vertices = new List <Vector3>();
List <Face> faces = new List <Face>();
List <Color> colors = new List <Color>();
int numSegments = lengths.Count - 1;
for (int segmentNum = 0; segmentNum < numSegments; segmentNum++)
{
AddSegment(
resolution,
lengths[segmentNum],
lengths[segmentNum + 1],
thickness,
offset,
energyColors[segmentNum],
energyColors[segmentNum + 1],
vertices,
faces,
colors
);
}
mesh.vertices = vertices.ToArray();
mesh.triangles = FacesToTriangles(faces);
mesh.colors = colors.ToArray();
}
public void HandlesSimpleCase(int val, byte red)
{
// Arrange
map.AddSample(0);
map.AddSample(128);
// Act
var color = map.GetColor(val);
// Assert
color.R.Should().Be(red);
color.G.Should().Be(0);
color.B.Should().Be(0);
}
private void DrawModulationSpectraHerz(List <float[]> spectra)
{
var minValue = spectra.SelectMany(s => s).Min();
var maxValue = spectra.SelectMany(s => s).Max();
var cmap = new ColorMap("blues", minValue, maxValue);
var g = modulationSpectrumPanel.CreateGraphics();
g.Clear(Color.White);
var spectrumBitmap = new Bitmap(spectra.Count, spectra[0].Length);
for (var i = 0; i < spectra.Count; i++)
{
for (var j = 0; j < spectra[i].Length; j++)
{
spectrumBitmap.SetPixel(i, spectra[i].Length - 1 - j, cmap.GetColor(spectra[i][j]));
}
}
g.DrawImage(spectrumBitmap, 25, 25, modulationSpectrumPanel.Width - 25, modulationSpectrumPanel.Height - 25);
}
/// <summary>
/// フィールド値を描画する
/// </summary>
/// <param name="g"></param>
/// <param name="ofs"></param>
/// <param name="delta"></param>
/// <param name="regionSize"></param>
/// <param name="colorMap"></param>
/// <param name="valueDv"></param>
public override void DrawField(Graphics g, Size ofs, Size delta, Size regionSize, FemElement.FieldDV fieldDv, FemElement.ValueDV valueDv, ColorMap colorMap)
{
//base.DrawField(g, ofs, delta, regionSize, colorMap);
if (_Nodes == null || _FValues == null || _RotXFValues == null || _RotYFValues == null || _PoyntingXFValues == null || _PoyntingYFValues == null)
{
return;
}
Complex[] tagtValues = null;
if (fieldDv == FemElement.FieldDV.Field)
{
tagtValues = _FValues;
}
else if (fieldDv == FemElement.FieldDV.RotX)
{
tagtValues = _RotXFValues;
}
else if (fieldDv == FemElement.FieldDV.RotY)
{
tagtValues = _RotYFValues;
}
else
{
return;
}
const int ndim = Constants.CoordDim2D; //2; // 座標の次元数
const int vertexCnt = Constants.QuadVertexCnt; //3; // 四角形形の頂点の数(2次要素でも同じ)
//const int nodeCnt = Constants.QuadNodeCnt_SecondOrder_Type2; //8; // 四角形2次要素
int nodeCnt = NodeNumbers.Length;
if (nodeCnt != Constants.QuadNodeCnt_SecondOrder_Type2 && nodeCnt != Constants.QuadNodeCnt_FirstOrder)
{
return;
}
// 四角形節点座標を取得
double[][] pp = new double[nodeCnt][];
for (int ino = 0; ino < pp.GetLength(0); ino++)
{
FemNode node = _Nodes[ino];
System.Diagnostics.Debug.Assert(node.Coord.Length == ndim);
pp[ino] = new double[ndim];
pp[ino][0] = node.Coord[0] * delta.Width + ofs.Width;
pp[ino][1] = regionSize.Height - node.Coord[1] * delta.Height + ofs.Height;
}
// 四角形内部を四角形で分割
// 要素節点座標( 局所r,s成分 )
// s
// |
// 3+ 6 +2
// | | |
// ---7---+---5-->r
// | | |
// 0+ 4 +1
// |
//
double[][] n_pts =
{
// r, s
new double[] {-1.0, -1.0}, //0
new double[] { 1.0, -1.0}, //1
new double[] { 1.0, 1.0}, //2
new double[] {-1.0, 1.0}, //3
new double[] { 0, -1.0}, //4
new double[] { 1.0, 0}, //5
new double[] { 0, 1.0}, //6
new double[] {-1.0, 0}, //7
};
int ndiv = this.IsCoarseFieldMesh ? (Constants.TriDrawFieldMshDivCnt / 2) : Constants.TriDrawFieldMshDivCnt;
double defdr = 2.0 / (double)ndiv;
double defds = defdr;
for (int i1 = 0; i1 < ndiv; i1++)
{
double r = - 1.0 + i1 * defdr;
double rNext = r + defdr;
for (int i2 = 0; i2 < ndiv; i2++)
{
double s = -1.0 + i2 * defds;
double sNext = s + defds;
// 四角形の頂点
const int rectVCnt = 4;
double[][] rect_local_p = new double[rectVCnt][]
{
new double[]{r , s },
new double[]{rNext, s },
new double[]{rNext, sNext},
new double[]{r , sNext}
};
double[][] rectpp = new double[rectVCnt][];
for (int ino = 0; ino < rectVCnt; ino++)
{
double work_r = rect_local_p[ino][0];
double work_s = rect_local_p[ino][1];
double xx = 0.0;
double yy = 0.0;
for (int k = 0; k < vertexCnt; k++)
{
double ri = n_pts[k][0];
double si = n_pts[k][1];
xx += pp[k][0] * 0.25 * (1 + ri * work_r) * (1 + si * work_s);
yy += pp[k][1] * 0.25 * (1 + ri * work_r) * (1 + si * work_s);
}
rectpp[ino] = new double[] { xx, yy };
}
// 表示する位置
double[] disp_p = new double[] { (rect_local_p[0][0] + rect_local_p[1][0]) * 0.5, (rect_local_p[0][1] + rect_local_p[3][1]) * 0.5 };
// 表示する値
Complex cvalue = new Complex(0.0, 0.0);
// 表示する位置の形状関数値
double[] workN = new double[nodeCnt];
if (nodeCnt == Constants.QuadNodeCnt_FirstOrder)
{
double work_r = disp_p[0];
double work_s = disp_p[1];
for (int i = 0; i < 4; i++)
{
// 節点の局所座標
double ri = n_pts[i][0];
double si = n_pts[i][1];
workN[i] = 0.25 * (1.0 + ri * work_r) * (1.0 + si * work_s);
}
}
else
{
double work_r = disp_p[0];
double work_s = disp_p[1];
// 節点0～3 : 四角形の頂点
for (int i = 0; i < 4; i++)
{
// 節点の局所座標
double ri = n_pts[i][0];
double si = n_pts[i][1];
// 形状関数N
workN[i] = 0.25 * (1.0 + ri * work_r) * (1.0 + si * work_s) * (ri * work_r + si * work_s - 1.0);
}
// 節点4,6 : r方向辺上中点
foreach (int i in new int[] { 4, 6 })
{
// 節点の局所座標
double ri = n_pts[i][0];
double si = n_pts[i][1];
// 形状関数N
workN[i] = 0.5 * (1.0 - work_r * work_r) * (1.0 + si * work_s);
}
// 節点5,7 : s方向辺上中点
foreach (int i in new int[] { 5, 7 })
{
// 節点の局所座標
double ri = n_pts[i][0];
double si = n_pts[i][1];
// 形状関数N
workN[i] = 0.5 * (1.0 + ri * work_r) * (1.0 - work_s * work_s);
}
}
for (int k = 0; k < nodeCnt; k++)
{
cvalue += tagtValues[k] * workN[k];
}
// 四角形の頂点（描画用）
Point[] rectp = new Point[rectVCnt];
for (int ino = 0; ino < rectVCnt; ino++)
{
rectp[ino] = new Point((int)rectpp[ino][0], (int)rectpp[ino][1]);
}
try
{
// 表示する値
double showValue = 0.0;
if (valueDv == ValueDV.Real)
{
showValue = cvalue.Real;
}
else if (valueDv == ValueDV.Imaginary)
{
showValue = cvalue.Imaginary;
}
else
{
// 既定値は絶対値
showValue = Complex.Abs(cvalue);
}
// 塗りつぶし色の取得
Color fillColor = colorMap.GetColor(showValue);
// 塗りつぶし
using (Brush brush = new SolidBrush(fillColor))
{
g.FillPolygon(brush, rectp);
}
}
catch (Exception exception)
{
System.Diagnostics.Debug.WriteLine(exception.Message + " " + exception.StackTrace);
}
}
}
}
public void GetColor_should_return_black_for_an_empty_map()
{
var map = new ColorMap();
map.GetColor(0).Should().Be(Colors.Black);
}
/// <summary>
/// フィールド値を描画する
/// </summary>
/// <param name="g"></param>
/// <param name="ofs"></param>
/// <param name="delta"></param>
/// <param name="regionSize"></param>
/// <param name="fieldDv"></param>
/// <param name="valueDv"></param>
/// <param name="colorMap"></param>
/*
public override void DrawField(Graphics g, Size ofs, Size delta, Size regionSize, FemElement.FieldDV fieldDv, FemElement.ValueDV valueDv, ColorMap colorMap)
{
//base.DrawField(g, ofs, delta, regionSize, colorMap);
if (_Nodes == null || _FValues == null || _RotXFValues == null || _RotYFValues == null || _PoyntingXFValues == null || _PoyntingYFValues == null)
{
return;
}
Complex[] tagtValues = null;
if (fieldDv == FemElement.FieldDV.Field)
{
tagtValues = _FValues;
}
else if (fieldDv == FemElement.FieldDV.RotX)
{
tagtValues = _RotXFValues;
}
else if (fieldDv == FemElement.FieldDV.RotY)
{
tagtValues = _RotYFValues;
}
else
{
return;
}
const int ndim = Constants.CoordDim2D; //2; // 座標の次元数
const int vertexCnt = Constants.TriVertexCnt; //3; // 三角形の頂点の数(2次要素でも同じ)
//const int nodeCnt = Constants.TriNodeCnt_SecondOrder; //6; // 三角形2次要素
int nodeCnt = NodeNumbers.Length;
if (nodeCnt != Constants.TriNodeCnt_SecondOrder && nodeCnt != Constants.TriNodeCnt_FirstOrder)
{
return;
}
// 三角形の節点座標を取得
double[][] pp = new double[nodeCnt][];
for (int ino = 0; ino < pp.GetLength(0); ino++)
{
FemNode node = _Nodes[ino];
System.Diagnostics.Debug.Assert(node.Coord.Length == ndim);
pp[ino] = new double[ndim];
pp[ino][0] = node.Coord[0] * delta.Width + ofs.Width;
pp[ino][1] = regionSize.Height - node.Coord[1] * delta.Height + ofs.Height;
}
// 下記分割ロジックの原点となる頂点
// 頂点0固定で計算していたが、原点の内角が直角のとき長方形メッシュになるので原点を2（頂点を0,1,2としたとき）にする
int orginVertexNo = 2;
// 内角が最大の頂点を取得し、その頂点を原点とする(後のロジックは原点が頂点を0,1,2としたとき、2になっている
{
double minCosth = double.MaxValue;
int minCosthVertexNo = 0;
for (int ino = 0; ino < vertexCnt; ino++)
{
const int vecCnt = 2;
double[][] vec = new double[vecCnt][] { new double[ndim]{0, 0}, new double[ndim]{0, 0} };
double[] len = new double[vecCnt];
double costh;
{
int n1 = ino;
int n2 = (ino + 1) % 3;
int n3 = (ino + 2) % 3;
vec[0][0] = pp[n2][0] - pp[n1][0];
vec[0][1] = pp[n2][1] - pp[n1][1];
vec[1][0] = pp[n3][0] - pp[n1][0];
vec[1][1] = pp[n3][1] - pp[n1][1];
len[0] = FemMeshLogic.GetDistance(pp[n1], pp[n2]);
len[1] = FemMeshLogic.GetDistance(pp[n1], pp[n3]);
costh = (vec[0][0] * vec[1][0] + vec[0][1] * vec[1][1]) / (len[0] * len[1]);
if (costh < minCosth)
{
minCosth = costh;
minCosthVertexNo = ino;
}
}
}
orginVertexNo = (minCosthVertexNo + 2) % 3;
}
// 三角形内部を四角形で分割
// 面積座標L1方向分割数
//int ndiv = 4;
int ndiv = Constants.TriDrawFieldMshDivCnt;
double defdL1 = 1.0 / (double)ndiv;
double defdL2 = defdL1;
for (int i1 = 0; i1 < ndiv; i1++)
{
double vL1 = i1 * defdL1;
double vL1Next = (i1 + 1) * defdL1;
if (i1 == ndiv - 1)
{
vL1Next = 1.0;
}
double vL2max = 1.0 - vL1;
if (vL2max < 0.0)
{
// ERROR
System.Diagnostics.Debug.WriteLine("logic error vL2max = {0}", vL2max);
continue;
}
double fdiv2 = (double)ndiv * vL2max;
int ndiv2 = (int)fdiv2;
if (fdiv2 - (double)ndiv2 > Constants.PrecisionLowerLimit)
{
ndiv2++;
}
for (int i2 = 0; i2 < ndiv2; i2++)
{
double vL2 = i2 * defdL2;
double vL2Next = (i2 + 1) * defdL2;
if (i2 == ndiv2 - 1)
{
vL2Next = vL2max;
}
double vL3 = 1.0 - vL1 - vL2;
if (vL3 < 0.0)
{
// ERROR
System.Diagnostics.Debug.WriteLine("logic error vL3 = {0}", vL3);
continue;
}
// 四角形の頂点
const int rectVCnt = 4;
double[][] rectLi = new double[rectVCnt][]
{
new double[]{vL1 , vL2 , 0},
new double[]{vL1Next, vL2 , 0},
new double[]{vL1Next, vL2Next, 0},
new double[]{vL1 , vL2Next, 0}
};
if ((i1 == ndiv - 1) || (i2 == ndiv2 - 1))
{
for (int k = 0; k < 3; k++)
{
rectLi[2][k] = rectLi[3][k];
}
}
double[][] rectpp = new double[rectVCnt][];
for (int ino = 0; ino < rectVCnt; ino++)
{
if (rectLi[ino][0] < 0.0)
{
rectLi[ino][0] = 0.0;
System.Diagnostics.Debug.WriteLine("logical error rectLi[{0}][0] = {1}", ino, rectLi[ino][0]);
}
if (rectLi[ino][0] > 1.0)
{
rectLi[ino][0] = 1.0;
System.Diagnostics.Debug.WriteLine("logical error rectLi[{0}][0] = {1}", ino, rectLi[ino][0]);
}
if (rectLi[ino][1] < 0.0)
{
rectLi[ino][1] = 0.0;
System.Diagnostics.Debug.WriteLine("logical error rectLi[{0}][1] = {1}", ino, rectLi[ino][1]);
}
if (rectLi[ino][1] > (1.0 - rectLi[ino][0])) // L2最大値(1 - L1)チェック
{
rectLi[ino][1] = 1.0 - rectLi[ino][0];
}
rectLi[ino][2] = 1.0 - rectLi[ino][0] - rectLi[ino][1];
if (rectLi[ino][2] < 0.0)
{
System.Diagnostics.Debug.WriteLine("logical error rectLi[{0}][2] = {1}", ino, rectLi[ino][2]);
}
}
for (int ino = 0; ino < rectVCnt; ino++)
{
double[] vLpp = rectLi[ino];
double xx = 0.0;
double yy = 0.0;
for (int k = 0; k < vertexCnt; k++)
{
xx += pp[k][0] * vLpp[(k + orginVertexNo) % vertexCnt];
yy += pp[k][1] * vLpp[(k + orginVertexNo) % vertexCnt];
}
rectpp[ino] = new double[] { xx, yy };
}
// 表示する位置
double[] vLi = new double[] { (rectLi[0][0] + rectLi[1][0]) * 0.5, (rectLi[0][1] + rectLi[3][1]) * 0.5, 0 };
if (vLi[0] < 0.0)
{
vLi[0] = 0.0;
}
if (vLi[0] > 1.0)
{
vLi[0] = 1.0;
}
if (vLi[1] < 0.0)
{
vLi[1] = 0.0;
}
if (vLi[1] > (1.0 - vLi[0]))
{
vLi[1] = (1.0 - vLi[0]);
}
vLi[2] = 1.0 - vLi[0] - vLi[1];
if (vLi[2] < 0.0)
{
System.Diagnostics.Debug.WriteLine("logic error vLi[2] = {0}", vLi[2]);
}
// 表示する値
Complex cvalue = new Complex(0.0, 0.0);
// 表示する位置の形状関数値
double[] vNi = null;
double[] shiftedLi = new double[vertexCnt];
for (int i = 0; i < vertexCnt; i++)
{
shiftedLi[i] = vLi[(i + orginVertexNo) % vertexCnt];
}
if (nodeCnt == Constants.TriNodeCnt_FirstOrder)
{
vNi = new double[]
{
shiftedLi[0],
shiftedLi[1],
shiftedLi[2]
};
}
else
{
vNi = new double[]
{
shiftedLi[0] * (2.0 * shiftedLi[0] - 1.0),
shiftedLi[1] * (2.0 * shiftedLi[1] - 1.0),
shiftedLi[2] * (2.0 * shiftedLi[2] - 1.0),
4.0 * shiftedLi[0] * shiftedLi[1],
4.0 * shiftedLi[1] * shiftedLi[2],
4.0 * shiftedLi[2] * shiftedLi[0],
};
}
for (int k = 0; k < nodeCnt; k++)
{
cvalue += tagtValues[k] * vNi[k];
}
// 四角形の頂点（描画用）
Point[] rectp = new Point[rectVCnt];
for (int ino = 0; ino < rectVCnt; ino++)
{
rectp[ino] = new Point((int)rectpp[ino][0], (int)rectpp[ino][1]);
}
try
{
// 表示する値
double showValue = 0.0;
if (valueDv == ValueDV.Real)
{
showValue = cvalue.Real;
}
else if (valueDv == ValueDV.Imaginary)
{
showValue = cvalue.Imaginary;
}
else
{
// 既定値は絶対値
showValue = Complex.Abs(cvalue);
}
// 塗りつぶし色の取得
Color fillColor = colorMap.GetColor(showValue);
// 塗りつぶし
using (Brush brush = new SolidBrush(fillColor))
{
g.FillPolygon(brush, rectp);
}
}
catch (Exception exception)
{
System.Diagnostics.Debug.WriteLine(exception.Message + " " + exception.StackTrace);
}
}
}
}
*/
public override void DrawField(Graphics g, Size ofs, Size delta, Size regionSize, FemElement.FieldDV fieldDv, FemElement.ValueDV valueDv, ColorMap colorMap)
{
//base.DrawField(g, ofs, delta, regionSize, colorMap);
if (_Nodes == null || _FValues == null || _RotXFValues == null || _RotYFValues == null || _PoyntingXFValues == null || _PoyntingYFValues == null)
{
return;
}
Complex[] tagtValues = null;
if (fieldDv == FemElement.FieldDV.Field)
{
tagtValues = _FValues;
}
else if (fieldDv == FemElement.FieldDV.RotX)
{
tagtValues = _RotXFValues;
}
else if (fieldDv == FemElement.FieldDV.RotY)
{
tagtValues = _RotYFValues;
}
else
{
return;
}
const int ndim = Constants.CoordDim2D; //2; // 座標の次元数
const int vertexCnt = Constants.TriVertexCnt; //3; // 三角形の頂点の数(2次要素でも同じ)
//const int nodeCnt = Constants.TriNodeCnt_SecondOrder; //6; // 三角形2次要素
int nodeCnt = NodeNumbers.Length;
if (nodeCnt != Constants.TriNodeCnt_SecondOrder && nodeCnt != Constants.TriNodeCnt_FirstOrder)
{
return;
}
// 三角形の節点座標を取得
double[][] pp = new double[nodeCnt][];
for (int ino = 0; ino < pp.GetLength(0); ino++)
{
FemNode node = _Nodes[ino];
System.Diagnostics.Debug.Assert(node.Coord.Length == ndim);
pp[ino] = new double[ndim];
pp[ino][0] = node.Coord[0] * delta.Width + ofs.Width;
pp[ino][1] = regionSize.Height - node.Coord[1] * delta.Height + ofs.Height;
}
// 長方形描画領域のリスト
IList<double[][]> rectLiList = _RectLiList;
// 描画ロジック上の原点となる頂点
int orginVertexNo = _OrginVertexNo;
// 四角形の頂点
const int rectVCnt = 4;
foreach (double[][] rectLi in rectLiList)
{
double[][] rectpp = new double[rectVCnt][];
for (int ino = 0; ino < rectVCnt; ino++)
{
double[] vLpp = rectLi[ino];
double xx = 0.0;
double yy = 0.0;
for (int k = 0; k < vertexCnt; k++)
{
xx += pp[k][0] * vLpp[(k + orginVertexNo) % vertexCnt];
yy += pp[k][1] * vLpp[(k + orginVertexNo) % vertexCnt];
}
rectpp[ino] = new double[] { xx, yy };
}
// 表示する位置
double[] vLi = new double[] { (rectLi[0][0] + rectLi[1][0]) * 0.5, (rectLi[0][1] + rectLi[3][1]) * 0.5, 0 };
if (vLi[0] < 0.0)
{
vLi[0] = 0.0;
}
if (vLi[0] > 1.0)
{
vLi[0] = 1.0;
}
if (vLi[1] < 0.0)
{
vLi[1] = 0.0;
}
if (vLi[1] > (1.0 - vLi[0]))
{
vLi[1] = (1.0 - vLi[0]);
}
vLi[2] = 1.0 - vLi[0] - vLi[1];
if (vLi[2] < 0.0)
{
System.Diagnostics.Debug.WriteLine("logic error vLi[2] = {0}", vLi[2]);
}
// 表示する値
Complex cvalue = new Complex(0.0, 0.0);
// 表示する位置の形状関数値
double[] vNi = null;
double[] shiftedLi = new double[vertexCnt];
for (int i = 0; i < vertexCnt; i++)
{
shiftedLi[i] = vLi[(i + orginVertexNo) % vertexCnt];
}
if (nodeCnt == Constants.TriNodeCnt_FirstOrder)
{
vNi = new double[]
{
shiftedLi[0],
shiftedLi[1],
shiftedLi[2]
};
}
else
{
vNi = new double[]
{
shiftedLi[0] * (2.0 * shiftedLi[0] - 1.0),
shiftedLi[1] * (2.0 * shiftedLi[1] - 1.0),
shiftedLi[2] * (2.0 * shiftedLi[2] - 1.0),
4.0 * shiftedLi[0] * shiftedLi[1],
4.0 * shiftedLi[1] * shiftedLi[2],
4.0 * shiftedLi[2] * shiftedLi[0],
};
}
for (int k = 0; k < nodeCnt; k++)
{
cvalue += tagtValues[k] * vNi[k];
}
// 四角形の頂点（描画用）
Point[] rectp = new Point[rectVCnt];
for (int ino = 0; ino < rectVCnt; ino++)
{
rectp[ino] = new Point((int)rectpp[ino][0], (int)rectpp[ino][1]);
}
try
{
// 表示する値
double showValue = 0.0;
if (valueDv == ValueDV.Real)
{
showValue = cvalue.Real;
}
else if (valueDv == ValueDV.Imaginary)
{
showValue = cvalue.Imaginary;
}
else
{
// 既定値は絶対値
showValue = Complex.Abs(cvalue);
}
// 塗りつぶし色の取得
Color fillColor = colorMap.GetColor(showValue);
// 塗りつぶし
using (Brush brush = new SolidBrush(fillColor))
{
g.FillPolygon(brush, rectp);
}
}
catch (Exception exception)
{
System.Diagnostics.Debug.WriteLine(exception.Message + " " + exception.StackTrace);
}
}
}
protected override void OnPaint(PaintEventArgs e)
{
base.OnPaint(e);
if (_spectrogram == null)
{
return;
}
var g = e.Graphics;
g.Clear(Color.White);
var width = Math.Min(Width, _spectrogram.Count);
var spectrogramBitmap = new Bitmap(width, _spectrogram[0].Length);
var realPos = 0;
var startPos = -AutoScrollPosition.X;
for (var i = startPos; i < startPos + spectrogramBitmap.Width; i++, realPos++)
{
for (var j = 0; j < _spectrogram[i].Length; j++)
{
spectrogramBitmap.SetPixel(realPos, _spectrogram[i].Length - 1 - j, _cmap.GetColor(_spectrogram[i][j]));
}
}
g.DrawImage(spectrogramBitmap, 0, 0);
if (_markline != null)
{
var pen = new Pen(Color.DeepPink, _marklineThickness);
realPos = 1;
for (var i = startPos + 1; i < startPos + spectrogramBitmap.Width; i++, realPos++)
{
g.DrawLine(pen, realPos - 1, _spectrogram[i].Length - 1 - (int)(_markline[i - 1]),
realPos, _spectrogram[i].Length - 1 - (int)(_markline[i]));
}
pen.Dispose();
}
}
/// <summary>
/// Get color.
/// </summary>
/// <param name="value">Value to return color from map.</param>
/// <returns>Color from color map.</returns>
public Color GetColor(float value)
{
return(ColorMap.GetColor(value));
}
public void Run(TimeMapOptions options)
{
// Set up the world, which will be used to find the files in the world.
var world = options.GetWorld();
// To properly scale the histogram, all chunks must be scanned. A 2D sparse matrix is used
// to store the inhabited times for each chunk. While that is being populated, the histogram
// is also created.
var matrix = new Sparse2DMatrix <int, int, long>();
var colorMap = new ColorMap();
foreach (var regionPath in world.ListRegionPaths())
{
var region = new Region(regionPath);
// TODO - region should expose a Coordinate2D, not separate X and Z values
var regionPt = new Coordinate2D(region.X, region.Z, CoordinateType2D.Region);
Console.WriteLine("Region {0}, X={1}, Z={2}:", regionPath, region.X, region.Z);
var found = 0;
for (var x = 0; x < 32; x++)
{
for (var z = 0; z < 32; z++)
{
var pt = new Coordinate2D(x, z, CoordinateType2D.ChunkWithinRegion).ToChunk(regionPt);
if (region.HasChunk(pt))
{
found += 1;
var chunk = region.GetChunk(pt);
var time = chunk.InhabitedTime;
matrix[pt.X, pt.Z] = time;
colorMap.AddSample(time);
}
}
}
Console.WriteLine(" {0} chunks", found);
}
// Set up the transformation from chunk coords to pixel coords
var dx = matrix.MaxX - matrix.MinX;
var dy = matrix.MaxY - matrix.MinY;
var transform = new Transform2D();
transform.InputRange(matrix.MinX, matrix.MaxX, matrix.MinY, matrix.MaxY);
transform.OutputRange(0, dx - 1, 0, dy - 1);
// Use the data to create the image
using (var image = new Image <Rgba32>(dx * options.ChunkSize, dy * options.ChunkSize))
{
for (var x = matrix.MinX; x < matrix.MaxX; x++)
{
for (var y = matrix.MinY; y < matrix.MaxY; y++)
{
var val = matrix[x, y];
var black = new Rgba32(0, 0, 0, 255);
Rgba32 color;
if (val == 0)
{
color = black;
}
else
{
color = colorMap.GetColor(val);
}
var px = options.ChunkSize * (int)transform.TransformX(x);
var py = options.ChunkSize * (int)transform.TransformY(y);
for (var xx = 0; xx < options.ChunkSize; xx++)
{
for (var yy = 0; yy < options.ChunkSize; yy++)
{
image[px + xx, py + yy] = color;
}
}
}
}
image.SaveAsPng(options.OutputPath);
Console.WriteLine("Wrote image to {0}.", options.OutputPath);
}
}
public void PlotLines(List <Wave> waves)
{
m_Model.Series.Clear();
m_Model.InvalidatePlot(true);
if (waves.Count != 0)
{
var minW = waves.Min(x => x.lambda) * 0.95;
var maxW = waves.Max(x => x.lambda) * 1.05;
m_Model.Axes[0].Minimum = minW;
m_Model.Axes[0].Maximum = maxW;
}
if (m_Limits != null)
{
SetWaveLimits(m_Limits);
}
var index = 0;
foreach (var groupedWave in waves.GroupBy(x => x.Order))
{
var waveSeries = new LinearBarSeries()
{
BarWidth = 0.0001,
StrokeThickness = 2,
StrokeColor = OxyColor.Parse(ColorMap.GetColor(index)).ToColor()
};
++index;
waveSeries.ItemsSource = groupedWave.Select(x =>
{
var dh = (x.dlambda * 1e-3) / 2.0;
return(new DataPoint(x.lambda, x.intensity));
}).OrderBy(x => x.X).ToList();
m_Model.Series.Add(waveSeries);
m_Model.InvalidatePlot(true);
}
var lineSeries = new LineSeries()
{
LineStyle = LineStyle.None,
MarkerFill = OxyColors.Blue.ToColor(),
MarkerSize = 5,
MarkerType = MarkerType.Circle
};
lineSeries.YAxisKey = "WaveReflectivity";
lineSeries.XAxisKey = "WaveKey";
List <DataPoint> pointList = new List <DataPoint>();
foreach (var wave in waves)
{
pointList.Add(new DataPoint(wave.lambda, wave.Efficiency));
}
lineSeries.ItemsSource = pointList;
lineSeries.Items.Refresh();
m_Model.Series.Add(lineSeries);
m_Model.InvalidatePlot(true);
}