private double Eval(TicoFeatures query, TicoFeatures template, List <MinutiaPair> matchingPair)
{
if (matchingPair == null)
{
return(0);
}
var qCenterMtia = matchingPair[matchingPair.Count - 1].QueryMtia;
var tCenterMtia = matchingPair[matchingPair.Count - 1].TemplateMtia;
// Computing query bounding region
var qPolygon = GetBounds(query);
qPolygon =
new PolygonMapper(qCenterMtia, tCenterMtia).Map(qPolygon);
var qBoundingRegion = new FingerprintRegion(qPolygon);
// Computing template bounding region
var tPolygon = GetBounds(template);
var tBoundingRegion = new FingerprintRegion(tPolygon);
// Mapping query minutiae
var mtiaMapper = new MtiaMapper(qCenterMtia, tCenterMtia);
var qMtiae = new List <Minutia>(query.Minutiae.Count);
for (var i = 0; i < query.Minutiae.Count; i++)
{
var mtia = query.Minutiae[i];
qMtiae.Add(mtiaMapper.Map(mtia));
}
var qCount = 0;
for (var i = 0; i < qMtiae.Count; i++)
{
var mtia = qMtiae[i];
if (tBoundingRegion.Contains(mtia))
{
qCount++;
}
}
var tCount = 0;
foreach (var mtia in template.Minutiae)
{
if (qBoundingRegion.Contains(mtia))
{
tCount++;
}
}
double sum = 0;
foreach (var mtiaPair in matchingPair)
{
sum += mtiaPair.MatchingValue;
}
return(qCount >= MtiaCountThr && tCount >= MtiaCountThr ? 1.0 * sum * sum / (qCount * tCount) : 0);
}
private List <Feature> MapMultiple(VectorPrescription prescription)
{
List <Feature> features = new List <Feature>();
double outOfFieldRate = -1.0; // something save to compare with
if (prescription.OutOfFieldRate != null)
{
outOfFieldRate = prescription.OutOfFieldRate.Value.Value;
}
foreach (var shaperate in prescription.RxShapeLookups)
{
int index = 0;
if (shaperate.Shape != null && shaperate.Shape.Polygons.Count > 0)
{
foreach (var adaptPolygon in shaperate.Shape.Polygons)
{
Dictionary <string, object> properties = new Dictionary <string, object>();
RxProductLookup product = prescription.RxProductLookups.Where(r => r.Id.ReferenceId == shaperate.Rates[index].RxProductLookupId).FirstOrDefault();
if (product != null)
{
properties.Add("productId", product.ProductId);
Product adaptProduct = _dataModel.Catalog.Products.Where(p => p.Id.ReferenceId == product.ProductId).FirstOrDefault();
if (adaptProduct != null)
{
properties.Add("productDescription", adaptProduct.Description);
properties.Add("productType", adaptProduct.ProductType.ToString());
}
properties.Add("productCode", product.Representation.Code);
properties.Add("productUom", product.UnitOfMeasure.Code);
}
else
{
properties.Add("productId", shaperate.Rates[index].RxProductLookupId);
}
properties.Add("rate", shaperate.Rates[index++].Rate);
features.Add(new Feature(PolygonMapper.MapPolygon(adaptPolygon, _properties.AffineTransformation), properties));
}
}
}
return(features);
}
private Feature Map(VectorPrescription prescription)
{
Feature feature = null;
Dictionary <string, object> properties = new Dictionary <string, object>();
// SpatialPrescription: not sure these rates mean anything for the geojson output
double outOfFieldRate = -1.0; // something save to compare with
if (prescription.OutOfFieldRate != null)
{
properties.Add("OutOfFieldRate", prescription.OutOfFieldRate.Value.Value);
outOfFieldRate = prescription.OutOfFieldRate.Value.Value;
}
if (prescription.LossOfGpsRate != null)
{
properties.Add("LossOfGpsRate", prescription.LossOfGpsRate.Value.Value);
}
double minRate = double.MaxValue;
double maxRate = -1;
foreach (var shaperate in prescription.RxShapeLookups)
{
if (shaperate.Shape != null && shaperate.Shape.Polygons.Count > 0)
{
foreach (var rate in shaperate.Rates)
{
if (rate.Rate != outOfFieldRate)
{
minRate = Math.Min(minRate, rate.Rate);
maxRate = Math.Max(maxRate, rate.Rate);
}
}
}
else // outoffield
{
foreach (var rate in shaperate.Rates)
{
RxProductLookup product = prescription.RxProductLookups.Where(r => r.Id.ReferenceId == rate.RxProductLookupId).FirstOrDefault();
if (product == null)
{
Console.WriteLine("outoffield rate? " + rate.Rate + " " + rate.RxProductLookupId);
}
else
{
Console.WriteLine("outoffield rate? " + rate.Rate + " " + rate.RxProductLookupId + " " + product.Representation.Code);
}
}
}
}
properties.Add("MinRate", minRate);
properties.Add("MaxRate", maxRate);
List <Dictionary <string, object> > products = new List <Dictionary <string, object> > {
};
foreach (var rxproduct in prescription.RxProductLookups)
{
if (products.Where(p => p.ContainsKey("productId") && (int)p["productId"] == rxproduct.ProductId).FirstOrDefault() == null)
{
Dictionary <string, object> product = new Dictionary <string, object>();
product.Add("productId", rxproduct.ProductId);
product.Add("productCode", rxproduct.Representation.Code);
product.Add("productUom", rxproduct.UnitOfMeasure.Code);
Product adaptProduct = _dataModel.Catalog.Products.Where(p => p.Id.ReferenceId == rxproduct.ProductId).FirstOrDefault();
if (adaptProduct != null)
{
product.Add("productDescription", adaptProduct.Description);
product.Add("productType", adaptProduct.ProductType.ToString()); // or via GetName?
}
products.Add(product);
}
}
properties.Add("Products", products.ToArray());
if (prescription.BoundingBox != null)
{
GeoJSON.Net.Geometry.Polygon polygon = PolygonMapper.MapBoundingBox(prescription.BoundingBox, _properties.AffineTransformation); // MinX, MinY, MaxX, MaxY
feature = new Feature(polygon, properties);
}
else
{
// @ToDo dissolve multipolygon or use Envelope?
// NetTopologySuite.Operation.Union.CascadedPolygonUnion
// NetTopologySuite.Algorithm.MinimumDiameter.GetMinimumRectangle
var polygons = new List <GeoJSON.Net.Geometry.Polygon>();
foreach (var shaperate in prescription.RxShapeLookups)
{
if (shaperate.Shape != null && shaperate.Shape.Polygons.Count > 0)
{
foreach (var adaptPolygon in shaperate.Shape.Polygons)
{
GeoJSON.Net.Geometry.Polygon polygon = PolygonMapper.MapPolygon(adaptPolygon, _properties.AffineTransformation);
if (polygon != null)
{
polygons.Add(polygon);
}
}
}
}
feature = new Feature(new GeoJSON.Net.Geometry.MultiPolygon(polygons), properties);
}
return(feature);
}
private List <Feature> MapMultiple(RasterGridPrescription prescription)
{
List <Feature> features = new List <Feature>();
// Based on Open.Topology.TestRunner.Functions/CreateShapeFunctions.Grid
var grid = new List <BoundingBox>();
int nCellsOnSideX = (int)prescription.ColumnCount;
int nCellsOnSideY = (int)prescription.RowCount;
double cellSizeX = prescription.CellWidth.Value.Value;
double cellSizeY = prescription.CellHeight.Value.Value;
double dMinX, dMinY;
if (prescription.BoundingBox != null)
{
dMinX = prescription.BoundingBox.MinX.Value.Value;
dMinY = prescription.BoundingBox.MinY.Value.Value;
}
else
{
dMinX = prescription.Origin.X;
dMinY = prescription.Origin.Y;
}
for (int j = 0; j < nCellsOnSideY; j++)
{
for (int i = 0; i < nCellsOnSideX; i++)
{
double x1 = dMinX + i * cellSizeX;
double y1 = dMinY + j * cellSizeY;
double x2 = dMinX + (i + 1) * cellSizeX;
double y2 = dMinY + (j + 1) * cellSizeY;
var bbox = new BoundingBox();
bbox.MinY = new NumericRepresentationValue(RepresentationInstanceList.vrLatitude.ToModelRepresentation(), new NumericValue(UnitSystemManager.GetUnitOfMeasure("arcdeg"), y1));
bbox.MinX = new NumericRepresentationValue(RepresentationInstanceList.vrLongitude.ToModelRepresentation(), new NumericValue(UnitSystemManager.GetUnitOfMeasure("arcdeg"), x1));
bbox.MaxY = new NumericRepresentationValue(RepresentationInstanceList.vrLatitude.ToModelRepresentation(), new NumericValue(UnitSystemManager.GetUnitOfMeasure("arcdeg"), y2));
bbox.MaxX = new NumericRepresentationValue(RepresentationInstanceList.vrLongitude.ToModelRepresentation(), new NumericValue(UnitSystemManager.GetUnitOfMeasure("arcdeg"), x2));
grid.Add(bbox);
}
}
int index = 0;
double outOfFieldRate = -1.0; // something save to compare with
if (prescription.OutOfFieldRate != null)
{
outOfFieldRate = prescription.OutOfFieldRate.Value.Value;
}
Console.WriteLine($"PrescriptionMapper outOfFieldRate: {outOfFieldRate}, grids {grid.Count}");
// @ToDo merge adjacent grid boxes with same property values?
foreach (var bbox in grid)
{
// skip outOfField grids
if (prescription.Rates[index].RxRates[0].Rate != outOfFieldRate)
{
GeoJSON.Net.Geometry.IGeometryObject geometry = PolygonMapper.MapBoundingBox(bbox, _properties.AffineTransformation);
Dictionary <string, object> properties = new Dictionary <string, object>();
RxProductLookup product = prescription.RxProductLookups.Where(r => r.Id.ReferenceId == prescription.Rates[index].RxRates[0].RxProductLookupId).FirstOrDefault();
if (product != null)
{
properties.Add("productId", product.ProductId);
Product adaptProduct = _dataModel.Catalog.Products.Where(p => p.Id.ReferenceId == product.ProductId).FirstOrDefault();
if (adaptProduct != null)
{
properties.Add("productDescription", adaptProduct.Description);
properties.Add("productType", adaptProduct.ProductType.ToString()); // or via GetName?
}
properties.Add("productCode", product.Representation.Code);
properties.Add("productUom", product.UnitOfMeasure.Code);
}
else
{
properties.Add("productId", prescription.Rates[index].RxRates[0].RxProductLookupId);
}
properties.Add("rate", prescription.Rates[index].RxRates[0].Rate);
features.Add(new Feature(geometry, properties));
}
index++;
}
return(features);
}
private Feature Map(RasterGridPrescription prescription)
{
GeoJSON.Net.Geometry.IGeometryObject geometry = null;
Dictionary <string, object> properties = new Dictionary <string, object>();
if (prescription.BoundingBox != null)
{
geometry = PolygonMapper.MapBoundingBox(prescription.BoundingBox, _properties.AffineTransformation);
}
else
{
geometry = PointMapper.MapPoint2Point(prescription.Origin, _properties.AffineTransformation);
}
properties.Add("RowCount", prescription.RowCount);
properties.Add("ColumnCount", prescription.ColumnCount);
properties.Add("CellWidth", prescription.CellWidth.Value.Value);
properties.Add("CellHeight", prescription.CellHeight.Value.Value);
//properties.Add("RxProductLookups", prescription.RxProductLookups); // Id, ProductId, Representation, UnitOfMeasure
// SpatialPrescription: not sure these rates mean anything for the geojson output
double outOfFieldRate = -1.0; // something save to compare with
if (prescription.OutOfFieldRate != null)
{
properties.Add("OutOfFieldRate", prescription.OutOfFieldRate.Value.Value);
outOfFieldRate = prescription.OutOfFieldRate.Value.Value;
}
if (prescription.LossOfGpsRate != null)
{
properties.Add("LossOfGpsRate", prescription.LossOfGpsRate.Value.Value);
}
double minRate = double.MaxValue;
double maxRate = -1.0;
foreach (var rate in prescription.Rates)
{
if (rate.RxRates[0].Rate != outOfFieldRate)
{
minRate = Math.Min(minRate, rate.RxRates[0].Rate);
maxRate = Math.Max(maxRate, rate.RxRates[0].Rate);
}
}
properties.Add("MinRate", minRate);
properties.Add("MaxRate", maxRate);
List <Dictionary <string, object> > products = new List <Dictionary <string, object> > {
};
foreach (var rxproduct in prescription.RxProductLookups)
{
if (products.Where(p => p.ContainsKey("productId") && (int)p["productId"] == rxproduct.ProductId).FirstOrDefault() == null)
{
Dictionary <string, object> product = new Dictionary <string, object>();
product.Add("productId", rxproduct.ProductId);
product.Add("productCode", rxproduct.Representation.Code);
product.Add("productUom", rxproduct.UnitOfMeasure.Code);
Product adaptProduct = _dataModel.Catalog.Products.Where(p => p.Id.ReferenceId == rxproduct.ProductId).FirstOrDefault();
if (adaptProduct != null)
{
product.Add("productDescription", adaptProduct.Description);
product.Add("productType", adaptProduct.ProductType.ToString()); // or via GetName?
}
products.Add(product);
}
}
//Array Values = products.Select(x => (Object)x).ToArray(); ;
properties.Add("Products", products.ToArray());
return(new Feature(geometry, properties));
}
private static double Eval(QiFeatures query, QiFeatures template, List <MinutiaPair> matchingPair,
double[,] simArr)
{
if (matchingPair == null)
{
return(0);
}
var qIndex = new Dictionary <Minutia, int>(query.Minutiae.Count);
for (var i = 0; i < query.Minutiae.Count; i++)
{
qIndex.Add(query.Minutiae[i], i);
}
var tIndex = new Dictionary <Minutia, int>(template.Minutiae.Count);
for (var i = 0; i < template.Minutiae.Count; i++)
{
tIndex.Add(template.Minutiae[i], i);
}
var qCenterMtia = matchingPair[matchingPair.Count - 1].QueryMtia;
var tCenterMtia = matchingPair[matchingPair.Count - 1].TemplateMtia;
// Computing query bounding region
var qPolygon = GetBounds(query);
qPolygon =
new PolygonMapper(qCenterMtia, tCenterMtia).Map(qPolygon);
var qBoundingRegion = new FingerprintRegion(qPolygon);
// Computing template bounding region
var tPolygon = GetBounds(template);
var tBoundingRegion = new FingerprintRegion(tPolygon);
// Mapping query minutiae
var mtiaMapper = new MtiaMapper(qCenterMtia, tCenterMtia);
var qMtiae = new List <Minutia>(query.Minutiae.Count);
for (var i = 0; i < query.Minutiae.Count; i++)
{
var mtia = query.Minutiae[i];
qMtiae.Add(mtiaMapper.Map(mtia));
}
var qCount = 0;
for (var i = 0; i < qMtiae.Count; i++)
{
var mtia = qMtiae[i];
if (tBoundingRegion.Contains(mtia))
{
qCount++;
}
}
var tCount = 0;
foreach (var mtia in template.Minutiae)
{
if (qBoundingRegion.Contains(mtia))
{
tCount++;
}
}
double sum = 0;
foreach (var mtiaPair in matchingPair)
{
var i = qIndex[mtiaPair.QueryMtia];
var j = tIndex[mtiaPair.TemplateMtia];
sum += simArr[i, j];
}
return(qCount >= 6 && tCount >= 6 ? 1.0 * sum / Math.Max(qCount, tCount) : 0);
}
