protected override List <ImagePoint> GetPointsInternal()
{
MyMatrix image = mat_provider.GetMatrix(), dx_mat = dx_provider.GetMatrix(), dy_mat = dy_provider.GetMatrix();
List <ImagePoint <Descriptor> > ans = new List <ImagePoint <Descriptor> >();
foreach (ImagePoint p in points_provider.GetPoints())
{
double
ox = dx_mat.sample(p.x, p.y),
oy = dy_mat.sample(p.x, p.y);
double len = Math.Sqrt(ox * ox + oy * oy);
ox /= len; oy /= len;
double[] vec_flat = new double[8 * 8];
for (int i = 0; i < 40; i++)
{
for (int j = 0; j < 40; j++)
{
// x: oy-ox
// y: -ox-oy
vec_flat[i / 5 * 8 + j / 5] += image.sample(p.x + (j - 19.5) * (oy - ox), p.y + (i - 19.5) * (-ox - oy));
}
}
double mu = vec_flat.Sum() / vec_flat.Length;
double ro = Math.Sqrt(vec_flat.Sum(v => v * v) / vec_flat.Length - mu * mu);
vec_flat = vec_flat.Select(v => (v - mu) / ro).ToArray();
double[,] vec = new double[8, 8];
for (int i = 0; i < 64; i++)
{
vec[i / 8, i % 8] = vec_flat[i];
}
ans.Add(new ImagePoint <Descriptor>(p.x, p.y, p.importance, new Descriptor(vec)));
}
return(ans.Cast <ImagePoint>().ToList());
}