예제 #1
0
        /// <summary>
        ///追踪处理单元中的射线
        /// </summary>
        private void TracingUnitRay(RayTracingModel unit, Node tx, Terrain ter, ReceiveBall rxBall, City cityBuilding, List <FrequencyBand> txFrequencyBand)
        {
            double rayBeamAngle = unit.AngleOfRayBeam * Math.PI / 180;

            //追踪射线1
            if (unit.FirstRay.WhetherHaveTraced == false)//若该射线之前还没有追踪
            {
                RayInfo     newRay   = new RayInfo(unit.FirstRay.Origin, unit.FirstRay.SVector);
                List <Path> pathtemp = this.GetSingleRayPaths(tx, ter, rxBall, cityBuilding, newRay, rayBeamAngle, txFrequencyBand);
                unit.FirstRay.Flag = this.JudgeIfArriveRx(pathtemp);
                unit.FirstRay.Path = pathtemp;
                unit.FirstRay.WhetherHaveTraced = true;
            }
            //追踪射线2
            if (unit.SecondRay.WhetherHaveTraced == false)//若该射线之前还没有追踪
            {
                RayInfo     newRay   = new RayInfo(unit.SecondRay.Origin, unit.SecondRay.SVector);
                List <Path> pathtemp = this.GetSingleRayPaths(tx, ter, rxBall, cityBuilding, newRay, rayBeamAngle, txFrequencyBand);
                unit.SecondRay.Flag             = this.JudgeIfArriveRx(pathtemp);
                unit.SecondRay.Path             = pathtemp;
                unit.FirstRay.WhetherHaveTraced = true;
            }
            //追踪射线3
            if (unit.ThirdRay.WhetherHaveTraced == false)//若该射线之前还没有追踪
            {
                RayInfo     newRay   = new RayInfo(unit.ThirdRay.Origin, unit.ThirdRay.SVector);
                List <Path> pathtemp = this.GetSingleRayPaths(tx, ter, rxBall, cityBuilding, newRay, rayBeamAngle, txFrequencyBand);
                unit.ThirdRay.Flag = this.JudgeIfArriveRx(pathtemp);
                unit.ThirdRay.Path = pathtemp;
                unit.FirstRay.WhetherHaveTraced = true;
            }
            unit.HaveTraced = true;
            return;
        }
예제 #2
0
        /// <summary>
        ///获取每个三角形细分的顶点
        /// </summary>
        protected List <Point> GetSubDivisionPoints(RayTracingModel unit, int TessellationFrequency)
        {
            List <Point> vertex = new List <Point>();
            Point        param1 = new Point(unit.FirstRay.SVector);
            Point        param2 = new Point(unit.SecondRay.SVector);
            Point        param3 = new Point(unit.ThirdRay.SVector);

            vertex.Add(param1);
            vertex.Add(param2);
            vertex.Add(param3);
            SpectVector n1     = new SpectVector(param1, param2);
            SpectVector n2     = new SpectVector(param1, param3);
            double      length = param1.GetDistance(param2) / TessellationFrequency;

            Point[] p1 = new Point[TessellationFrequency];
            Point[] p2 = new Point[TessellationFrequency];
            for (int i = 0; i < TessellationFrequency; i++)
            {
                p1[i] = new Point(param1.X + (i + 1) * (n1.a / Math.Sqrt(Math.Pow(n1.a, 2) + Math.Pow(n1.b, 2) + Math.Pow(n1.c, 2)) * length), param1.Y + (i + 1) * (n1.b / Math.Sqrt(Math.Pow(n1.a, 2) + Math.Pow(n1.b, 2) + Math.Pow(n1.c, 2)) * length), param1.Z + (i + 1) * (n1.c / Math.Sqrt(Math.Pow(n1.a, 2) + Math.Pow(n1.b, 2) + Math.Pow(n1.c, 2)) * length));
                p2[i] = new Point(param1.X + (i + 1) * (n2.a / Math.Sqrt(Math.Pow(n2.a, 2) + Math.Pow(n2.b, 2) + Math.Pow(n2.c, 2)) * length), param1.Y + (i + 1) * (n2.b / Math.Sqrt(Math.Pow(n2.a, 2) + Math.Pow(n2.b, 2) + Math.Pow(n2.c, 2)) * length), param1.Z + (i + 1) * (n2.c / Math.Sqrt(Math.Pow(n2.a, 2) + Math.Pow(n2.b, 2) + Math.Pow(n2.c, 2)) * length));
                vertex.Add(p1[i]);
                vertex.Add(p2[i]);
            }
            for (int j = 0; j < TessellationFrequency; j++)
            {
                SpectVector m = new SpectVector(p1[j], p2[j]);
                GetTrianglePoint(m, p1[j], p2[j], length, vertex);
            }
            return(vertex);
        }
예제 #3
0
        /// <summary>
        ///将初始三角形上的各条射线分成处理单元
        /// </summary>
        protected List <RayTracingModel> HandleInitialRayToUnit(List <ClassNewRay>[] raysPath)
        {
            try
            {
                List <RayTracingModel> triangleUnits = new List <RayTracingModel>();
                for (int i = 0; i < raysPath.Length - 1; i++)
                {
                    for (int j = 0; j < raysPath[i].Count; j++)
                    {
                        RayTracingModel param = new RayTracingModel(raysPath[i][j], raysPath[i + 1][j], raysPath[i + 1][j + 1]);
                        param.HaveTraced = true;
                        triangleUnits.Add(param);
                    }
                }
                if (raysPath.Length > 2)//若三角形平分大于2层
                {
                    for (int m = raysPath.Length - 1; m > 2; m--)
                    {
                        for (int n = 1; n < raysPath[m].Count - 1; n++)
                        {
                            RayTracingModel param = new RayTracingModel(raysPath[m][n], raysPath[m - 1][n - 1], raysPath[m - 1][n]);
                            param.HaveTraced = true;
                            triangleUnits.Add(param);
                        }
                    }
                }

                return(triangleUnits);
            }
            catch (Exception)
            {
                LogFileManager.ObjLog.error("在分配初始射线成处理单元时出错");
                return(new List <RayTracingModel>());
            }
        }
예제 #4
0
        /// <summary>
        ///细分处理单元
        /// </summary>
        protected List <RayTracingModel> DivideRayUnit(RayTracingModel unit)
        {
            if (unit == null)
            {
                LogFileManager.ObjLog.error("在细分处理单元时,输入的处理单元为null");
                return(new List <RayTracingModel>());
            }
            List <RayTracingModel> divideUnits = new List <RayTracingModel>();

            //3条射线都没有到达接收机
            if (unit.FirstRay.Flag == false && unit.SecondRay.Flag == false && unit.ThirdRay.Flag == false)
            {
                return(divideUnits);
            }
            //只有射线1到达接收机
            else if (unit.FirstRay.Flag == true && unit.SecondRay.Flag == false && unit.ThirdRay.Flag == false)
            {
                divideUnits = GetFirstRayUnits(unit);
                return(divideUnits);
            }
            //只有射线2到达接收机
            else if (unit.FirstRay.Flag == false && unit.SecondRay.Flag == true && unit.ThirdRay.Flag == false)
            {
                divideUnits = GetSecongRayUnits(unit);
                return(divideUnits);
            }
            //只有射线3到达接收机
            else if (unit.FirstRay.Flag == false && unit.SecondRay.Flag == false && unit.ThirdRay.Flag == true)
            {
                divideUnits = GetThirdRayUnits(unit);
                return(divideUnits);
            }
            //射线1和射线2到达接收机
            else if (unit.FirstRay.Flag == true && unit.SecondRay.Flag == true && unit.ThirdRay.Flag == false)
            {
                divideUnits = GetFirstAndSecondRaysUnits(unit);
                return(divideUnits);
            }
            //射线1和射线3到达接收机
            else if (unit.FirstRay.Flag == true && unit.SecondRay.Flag == false && unit.ThirdRay.Flag == true)
            {
                divideUnits = GetFirstAndThirdRaysUnits(unit);
                return(divideUnits);
            }
            //射线2和射线3到达接收机
            else if (unit.FirstRay.Flag == false && unit.SecondRay.Flag == true && unit.ThirdRay.Flag == true)
            {
                divideUnits = GetSecongdAndThirdRaysUnits(unit);
                return(divideUnits);
            }
            else
            {
                divideUnits = GetAllRaysUnits(unit);
                return(divideUnits);
            }
        }
예제 #5
0
        /// <summary>
        ///获取射线3的细分单元
        /// </summary>
        protected List <RayTracingModel> GetThirdRayUnits(RayTracingModel unit)
        {
            ClassNewRay            ray3        = (ClassNewRay)unit.ThirdRay.Clone();
            ClassNewRay            ray4        = new ClassNewRay(unit.ThirdRay.Origin, GetMiddleVectorOfTwoVectors(unit.FirstRay.SVector, unit.SecondRay.SVector), false);
            ClassNewRay            ray5        = new ClassNewRay(unit.ThirdRay.Origin, GetMiddleVectorOfTwoVectors(unit.FirstRay.SVector, unit.ThirdRay.SVector), false);
            ClassNewRay            ray6        = new ClassNewRay(unit.ThirdRay.Origin, GetMiddleVectorOfTwoVectors(unit.SecondRay.SVector, unit.ThirdRay.SVector), false);
            List <RayTracingModel> divideUnits = new List <RayTracingModel>();

            divideUnits.Add(new RayTracingModel(ray3, ray5, ray6));
            divideUnits.Add(new RayTracingModel(ray4, ray5, ray6));
            divideUnits.Add(new RayTracingModel(unit.FirstRay, ray4, ray5, true));
            divideUnits.Add(new RayTracingModel(unit.SecondRay, ray4, ray6, true));
            return(divideUnits);
        }
예제 #6
0
 /// <summary>
 ///提取路径
 /// </summary>
 protected void ExtractPathOfUnit(RayTracingModel unit, List <Path> pathInfo)
 {
     if (unit.FirstRay.Path != null && unit.FirstRay.Path.Count != 0)
     {
         pathInfo.AddRange(unit.FirstRay.Path);
     }
     if (unit.SecondRay.Path != null && unit.SecondRay.Path.Count != 0)
     {
         pathInfo.AddRange(unit.SecondRay.Path);
     }
     if (unit.ThirdRay.Path != null && unit.ThirdRay.Path.Count != 0)
     {
         pathInfo.AddRange(unit.ThirdRay.Path);
     }
 }
예제 #7
0
        /// <summary>
        ///初始阶段在每个三角形上发出一定的射线
        /// </summary>
        /// <param name="tx">发射机</param>
        /// <param name="ter">地形</param>
        /// <param name="rxBall">接收球</param>
        ///  <param name="cityBuilding">建筑物</param>
        ///   <param name="unit">射线处理单元</param>
        ///  <param name="TessellationFrequency">镶嵌次数,等于三角形每条边的细分次数</param>
        ///  <param name="divideAngle">射线束夹角</param>
        ///  <param name="TxFrequencyBand">发射机频段信息</param>
        /// <returns></returns>
        private List <RayTracingModel> GetInitialTrianglePath(Node tx, Terrain ter, ReceiveBall rxBall, City cityBuilding, RayTracingModel unit, int TessellationFrequency, double divideAngle, List <FrequencyBand> txFrequencyBand)
        {
            Point       vertex1  = new Point(unit.FirstRay.SVector);  //三角形的顶点
            Point       vertex2  = new Point(unit.SecondRay.SVector); //三角形的顶点
            Point       vertex3  = new Point(unit.ThirdRay.SVector);  //三角形的顶点
            SpectVector vector12 = new SpectVector(vertex1, vertex2);
            SpectVector vector23 = new SpectVector(vertex2, vertex3);
            double      length   = vertex1.GetDistance(vertex2) / TessellationFrequency; //每段的距离

            Point[] pointOfEdge12 = new Point[TessellationFrequency];                    //存放棱上的点,不包括顶点1
            //    p1[0] = vertex1;
            for (int i = 0; i < TessellationFrequency; i++)
            {
                pointOfEdge12[i] =
                    new Point(vertex1.X + (i + 1) * (vector12.a / Math.Sqrt(Math.Pow(vector12.a, 2) + Math.Pow(vector12.b, 2) + Math.Pow(vector12.c, 2)) * length),
                              vertex1.Y + (i + 1) * (vector12.b / Math.Sqrt(Math.Pow(vector12.a, 2) + Math.Pow(vector12.b, 2) + Math.Pow(vector12.c, 2)) * length),
                              vertex1.Z + (i + 1) * (vector12.c / Math.Sqrt(Math.Pow(vector12.a, 2) + Math.Pow(vector12.b, 2) + Math.Pow(vector12.c, 2)) * length));
            }
            List <Point>[] vertexPoints = new List <Point> [TessellationFrequency + 1];//存放三角形内每条平行边上的点
            vertexPoints[0] = new List <Point> {
                vertex1
            };                                               //把顶点1放到数组第一位
            for (int j = 1; j <= TessellationFrequency; j++) //得到三角形切分后每条边上的点
            {
                vertexPoints[j] = new List <Point>();
                vertexPoints[j].Add(pointOfEdge12[j - 1]);
                this.GetTriangleInternalEdgePoint(vector23, pointOfEdge12[j - 1], length, j, vertexPoints[j]);
            }
            List <ClassNewRay>[]   raysPath      = this.GetEachTriangleRay(tx, ter, rxBall, cityBuilding, TessellationFrequency, divideAngle, txFrequencyBand, vertexPoints); //根据得到的点构造射线
            List <RayTracingModel> triangleUnits = this.HandleInitialRayToUnit(raysPath);                                                                                     //得到初始追踪完的三角形处理单元

            return(triangleUnits);
        }