public List <MyLib.Point> findObjectinarea(TMBR mbr, int node_id, List <MyLib.Point> array, int idx) { double eps = Dbscanwithrtree.eps; MyLib.Point f = new MyLib.Point(mbr.Left.X + eps, mbr.Left.Y - eps, -1); for (int i = 0; i < FNodeArr[node_id].FObject.Length; i++) { if (FNodeArr[node_id].FObject[i].idx != idx && f.Distance(FNodeArr[node_id].FObject[i]) <= eps) { array.Add(FNodeArr[node_id].FObject[i]); } } return(array); }
/// <summary> /// метод для вставки обьекта в дерево /// </summary> /// <param name="obj"></param> public void insertObject(MyLib.Point obj) { int node_id = FRoot; Console.WriteLine(ct++); node_id = chooseSubtree(obj, node_id); if (FNodeArr[node_id].FObject.Length < MAX_M) //если количество объектов в узле меньше максимального { FNodeArr[node_id].setObject(FNodeArr[node_id].FObject.Length, obj); updateMBR(node_id); } else//если количество объектов достигло максимально допустимого { splitNodeRStar(node_id, obj); //делим узел } }
private int chooseSubtree(MyLib.Point obj, int node_id) { //Console.WriteLine("chooseSubtree work"); int i, id_child; int[] idChild_area = new int[0]; int id_zero; double dspace; if (FNodeArr[node_id].isLeaf) { return(node_id); } Array.Resize(ref idChild_area, FNodeArr[node_id].FChildren.Length); dspace = double.MaxValue; id_zero = 0; for (i = 0; i < FNodeArr[node_id].FChildren.Length; i++) { idChild_area[i] = FNodeArr[node_id].FChildren[i]; } dspace = double.MaxValue; for (i = 0; i < idChild_area.Length; i++) { id_child = idChild_area[i]; if (FNodeArr[id_child].Area() < dspace) { id_zero = idChild_area[i]; dspace = FNodeArr[id_child].Area(); } } node_id = id_zero; chooseSubtree(obj, node_id); return(node_id); }
static bool Check(TRtree tree, MyLib.Point point) { ForRstartree.Point lefttop = new ForRstartree.Point(point.x - eps, point.y + eps); ForRstartree.Point rightdown = new ForRstartree.Point(point.x + eps, point.y - eps); ForRstartree.TMBR findzone = new ForRstartree.TMBR(lefttop, rightdown);//Получили квадрат описывающий круг с центром в текущей точке List <MyLib.Point> list = new List <MyLib.Point>(0); List <MyLib.Point> find = new List <MyLib.Point>(0); find = tree.findobjectinarea(findzone, tree.FRoot, list, point.idx);//получили лист точек, которые находятся в eps-окрестности текущей(ее не берем) for (int i = 0; i < tree.FNodeArr.Length; i++) { tree.FNodeArr[i].IsVisited = false; } point.findpoint = find; if (find.Count >= minPTS) { return(true); } else { return(false); } }
/// <summary> /// TAXis /// </summary> /// <param name="obj"></param> /// <param name="node_id"></param> /// <returns></returns> #region private TAxis chooseSplitAxis(MyLib.Point obj, int node_id) { //Console.WriteLine("chooseSplitAxis1 work"); MyLib.Point[] arr_obj; int i, j, k, idx; TRNode node_1, node_2; double perimetr_min, perimetr; TAxis result = new TAxis(); arr_obj = new MyLib.Point[MAX_M + 1]; if (!FNodeArr[node_id].isLeaf) { return(0); } for (i = 0; i < FNodeArr[node_id].FObject.Length; i++) { arr_obj[i] = FNodeArr[node_id].FObject[i]; } arr_obj[arr_obj.Length - 1] = obj; node_1 = new TRNode(); node_2 = new TRNode(); perimetr_min = double.MaxValue; for (i = 0; i <= 1; i++) { perimetr = 0; for (j = 0; j <= 1; j++) { node_1.clearObjects(); node_2.clearObjects(); QuickSort(arr_obj, 0, arr_obj.Length - 1, (TAxis)i); for (k = 1; k <= MAX_M - MIN_M * 2 + 2; k++) // высчитваем периметры { idx = 0; while (idx < ((MIN_M - 1) + k)) { node_1.setObject(idx, arr_obj[idx]); idx++; } for (; idx < arr_obj.Length; idx++) { node_2.setObject(idx - ((MIN_M - 1) + k), arr_obj[idx]); } updateMBR(node_1); updateMBR(node_2); perimetr = perimetr + ((node_1.mbr.Right.X - node_1.mbr.Left.X) * 2 + (node_2.mbr.Left.Y - node_2.mbr.Right.Y) * 2); } } if (perimetr <= perimetr_min) { result = (TAxis)i; perimetr_min = perimetr; } } return(result); }
/// <summary> /// разделяет узел на 2 в соответствии с алгоритмами R*-tree (node_id = ссылка на узел для разделения obj = обьект для вставки) /// </summary> #region private void splitNodeRStar(int node_id, MyLib.Point obj) { //Console.WriteLine("splitNodeRStar work1 " + FHeight); TAxis axe; int parent_id, new_child_id; TRNode node_1, node_2, node_1_min, node_2_min; int i, j, k; MyLib.Point[] arr_obj; double area_min, Area; if (!FNodeArr[node_id].isLeaf) { return; } if (isRoot(node_id)) { //Console.WriteLine("it works"); parent_id = newNode(); FNodeArr[FRoot].Parent = parent_id; FNodeArr[parent_id].setChild(0, FRoot); FNodeArr[parent_id].Level = FNodeArr[FNodeArr[parent_id].getChild(0)].Level + 1; FRoot = parent_id; FHeight++; } else { parent_id = newNode(); FNodeArr[parent_id].Parent = FNodeArr[node_id].Parent; FNodeArr[FNodeArr[parent_id].Parent].setChild(FNodeArr[FNodeArr[parent_id].Parent].FChildren.Length, parent_id); FNodeArr[node_id].Parent = parent_id; FNodeArr[parent_id].setChild(0, node_id); /* * parent_id = FNodeArr[node_id].Parent; * Console.WriteLine("f**k); " + parent_id); * FNodeArr[FNodeArr[node_id].Parent].FisLeaf = false; */ } arr_obj = new MyLib.Point[MAX_M + 1]; for (i = 0; i < arr_obj.Length - 1; i++) { arr_obj[i] = FNodeArr[node_id].getObject(i); } arr_obj[arr_obj.Length - 1] = obj; node_1_min = new TRNode(); node_2_min = new TRNode(); node_1 = new TRNode(); node_2 = new TRNode(); axe = chooseSplitAxis(obj, node_id); area_min = double.MaxValue; for (i = 0; i <= 1; i++) { QuickSort(arr_obj, 0, arr_obj.Length - 1, axe); for (k = MIN_M - 1; k <= MAX_M - MIN_M; k++) { node_1.clearObjects(); node_2.clearObjects(); j = 0; while (j <= k) { node_1.setObject(j, arr_obj[j]); j++; } for (j = k; j < arr_obj.Length - 1; j++) { node_2.setObject(j - k, arr_obj[j + 1]); } updateMBR(node_1); updateMBR(node_2); Area = node_1.Area() + node_2.Area(); if (Area < area_min) { node_1_min.copy(node_1); node_2_min.copy(node_2); area_min = Area; } } } node_1_min.Level = 0; node_2_min.Level = 0; FNodeArr[node_id].copy(node_1_min); FNodeArr[node_id].Parent = parent_id; updateMBR(node_id); new_child_id = newNode(); FNodeArr[new_child_id].copy(node_2_min); FNodeArr[new_child_id].Parent = parent_id; updateMBR(new_child_id); if (FNodeArr[parent_id].FChildren.Length < MAX_M) { FNodeArr[parent_id].setChild(FNodeArr[parent_id].FChildren.Length, new_child_id); updateMBR(parent_id); FNodeArr[parent_id].FisLeaf = false; } else { splitNodeRStar(parent_id, new_child_id); } }