public static PointPolar PtdDecartToPolar(PointD decartPointD, PointD center, PointD zeroAnglePoint) { double diffX = zeroAnglePoint.X - center.X; double diffY = zeroAnglePoint.Y - center.Y; double zeroAnglePointR = Math.Sqrt(diffX * diffX + diffY * diffY); double zeroAngle = 0.0d; if (zeroAnglePointR != 0.0d) { zeroAngle = Math.Acos(diffX / zeroAnglePointR); if (diffY >= 0) { zeroAngle = 2.0d * Math.PI - zeroAngle; } } PointPolar polarPointD = new PointPolar(); diffX = decartPointD.X - center.X; diffY = decartPointD.Y - center.Y; polarPointD.R = Math.Sqrt(diffX * diffX + diffY * diffY); polarPointD.Phi = Math.Acos(diffX / polarPointD.R); if (diffY >= 0) { polarPointD.Phi = 2.0d * Math.PI - polarPointD.Phi; } polarPointD.Phi -= zeroAngle; polarPointD.CropAngle(true); return(polarPointD); }
public static PointD PtdPolarToCartesian(PointPolar polarPointD, PointD centerCartesian, double zeroAngleInCartesian) { PointD cartesianPointD = new PointD(); cartesianPointD.X = polarPointD.R * Math.Cos(polarPointD.Phi + zeroAngleInCartesian) + centerCartesian.X; cartesianPointD.Y = -polarPointD.R * Math.Sin(polarPointD.Phi + zeroAngleInCartesian) + centerCartesian.Y; return(cartesianPointD); }
public static PointD PtdPolarToCartesian(PointPolar polarPointD, PointD centerCartesian, PointD zeroAnglePointInCartesian) { double diffX = zeroAnglePointInCartesian.X - centerCartesian.X; double diffY = zeroAnglePointInCartesian.Y - centerCartesian.Y; double zeroAnglePointR = Math.Sqrt(diffX * diffX + diffY * diffY); double zeroAngle = 0.0d; if (zeroAnglePointR != 0.0d) { zeroAngle = Math.Acos(diffX / zeroAnglePointR); if (diffY >= 0) { zeroAngle += Math.PI; } } PointD decartPointD = new PointD(); decartPointD.X = polarPointD.R * Math.Cos(polarPointD.Phi + zeroAngle) + centerCartesian.X; decartPointD.Y = -polarPointD.R * Math.Sin(polarPointD.Phi + zeroAngle) + centerCartesian.Y; return(decartPointD); }
public static List <PointPolar> ListDecartToPolar(List <PointD> decartPointdList, PointD center, PointD zeroAnglePoint, out double zeroAngleValue) { double diffX = zeroAnglePoint.X - center.X; double diffY = zeroAnglePoint.Y - center.Y; double zeroAnglePointR = Math.Sqrt(diffX * diffX + diffY * diffY); double zeroAngle = 0.0d; if (zeroAnglePointR != 0.0d) { zeroAngle = Math.Acos(diffX / zeroAnglePointR); if (diffY >= 0) { zeroAngle = 2.0d * Math.PI - zeroAngle; } } zeroAngleValue = zeroAngle; List <PointPolar> retList = new List <PointPolar>(); foreach (PointD curPointDecart in decartPointdList) { PointPolar polarPointD = new PointPolar(); diffX = curPointDecart.X - center.X; diffY = curPointDecart.Y - center.Y; polarPointD.R = Math.Sqrt(diffX * diffX + diffY * diffY); polarPointD.Phi = Math.Acos(diffX / polarPointD.R); if (diffY >= 0) { polarPointD.Phi = 2.0d * Math.PI - polarPointD.Phi; } polarPointD.Phi -= zeroAngle; polarPointD.CropAngle(true); retList.Add(polarPointD); } return(retList); }
/// <summary> /// Gets the local minimums distribution. /// </summary> /// <param name="dmFieldData">The dm field data.</param> /// <param name="dimensionNumber">The dimension number: /// 1 - rows (angle) /// 2 - columns (distance) /// </param> /// <returns>DenseMatrix.</returns> //public static DenseMatrix GetLocalMinimumsDistribution(DenseMatrix dmFieldData, PointD sunCenterPoint, PointD imageCenterPoint, double imageRadius, int imageHeight, double imageCircleCropFactor = 0.9d, int dimensionNumber = 1) public static List <Point3D> GetLocalMinimumsDistribution(DenseMatrix dmFieldData, RoundData sunDiskData, RoundData imageRoundData, int imageHeight, double imageCircleCropFactor = 0.9d) { // DenseMatrix dmFieldminimumsData = DenseMatrix.Create(dmFieldData.RowCount, dmFieldData.ColumnCount, 0.0d); List <Point3D> lRetPoints = new List <Point3D>(); double imageRadius = imageRoundData.DRadius; PointD imageCenterPoint = imageRoundData.pointDCircleCenter(); PointPolar imageCenterPointRelatedToSunCenter = new PointPolar(imageCenterPoint - sunDiskData.pointDCircleCenter(), true); double distanceSunCenterToImageCenter = PointD.Distance(imageCenterPoint, sunDiskData.pointDCircleCenter()); #region // obsolete //if (dimensionNumber == 1) //{ #endregion // obsolete for (int i = 0; i < dmFieldData.RowCount; i++) { bool itsTheCropCase = false; //если направлени на кроп кадра - то не берем в расмотрение double currentAngle = ((double)i / (double)(dmFieldData.RowCount - 1)) * 2.0d * Math.PI; LineDescription2D line, lineMargin; if (currentAngle < Math.PI) { //верхняя половина, смотрим направление на y=0.0d line = new LineDescription2D(sunDiskData.pointDCircleCenter(), new Vector2D(Math.Cos(currentAngle), -Math.Sin(currentAngle))); lineMargin = new LineDescription2D(new PointD(0.0d, 0.0d), new Vector2D(1.0d, 0.0d)); } else { line = new LineDescription2D(sunDiskData.pointDCircleCenter(), new Vector2D(Math.Cos(currentAngle), Math.Sin(currentAngle))); lineMargin = new LineDescription2D(new PointD(0.0d, imageHeight), new Vector2D(1.0d, 0.0d)); } PointD crossPointD = LineDescription2D.CrossPoint(line, lineMargin); if (crossPointD.Distance(imageCenterPoint) < imageRadius) { itsTheCropCase = true; } #region // obsolete //double yMargin = 0.0d; //double xMargin = sunCenterPoint.X + (yMargin - sunCenterPoint.Y) / Math.Tan(currentAngle); //double dx = xMargin - imageCenterPoint.X; //double dy = yMargin - imageCenterPoint.Y; //if (Math.Sqrt(dx * dx + dy * dy) < imageRadius) itsTheCropCase = true; #endregion // obsolete #region //obsolete //else //{ // //нижняя половина, смотрим направление на y=imageHeight // double yMargin = (double)imageHeight; // double xMargin = sunCenterPoint.X + (yMargin - sunCenterPoint.Y) / Math.Tan(currentAngle); // double dx = xMargin - imageCenterPoint.X; // double dy = yMargin - imageCenterPoint.Y; // if (Math.Sqrt(dx * dx + dy * dy) < imageRadius) itsTheCropCase = true; //} #endregion //obsolete //Если слишком близко к краю изображения - тоже исключаем. Минимум должен лежать не ближе, например, 1/15 //DenseMatrix dmSlicedDataMatrix = (DenseMatrix)dmFieldData.SubMatrix(i, 1, 0, dmFieldData.ColumnCount); DenseVector dvRowDataVector = (DenseVector)dmFieldData.EnumerateRows().ElementAt(i); #region // debug plotting //dvRowDataVector.SaveVectorDataAsImagePlot( // "D:\\_gulevlab\\SkyImagesAnalysis_appData\\patent-samples\\result.2015-03-24\\img-2014-09-20T16-03-58devID1\\dvRowDataVector-plot-image-" + // i.ToString("D03") + "-step1.png"); #endregion // debug plotting dvRowDataVector.MapIndexedInplace((idx, x) => ((x == 0.0d) || (idx < sunDiskData.DRadius * 1.5d)) ? (1.0d) : (x)); #region // debug plotting //dvRowDataVector.SaveVectorDataAsImagePlot( // "D:\\_gulevlab\\SkyImagesAnalysis_appData\\patent-samples\\result.2015-03-24\\img-2014-09-20T16-03-58devID1\\dvRowDataVector-plot-image-" + // i.ToString("D03") + "-step2.png"); #endregion // debug plotting double phiFromImageCenterToDirection = imageCenterPointRelatedToSunCenter.Phi - currentAngle; double distanceToImageMargin = distanceSunCenterToImageCenter * Math.Cos(phiFromImageCenterToDirection) + Math.Sqrt(imageRadius * imageRadius - distanceSunCenterToImageCenter * distanceSunCenterToImageCenter * Math.Sin(phiFromImageCenterToDirection) * Math.Sin(phiFromImageCenterToDirection)); dvRowDataVector.MapIndexedInplace( (idx, x) => ((double)idx / distanceToImageMargin >= imageCircleCropFactor) ? (1.0d) : (x)); #region // debug plotting //dvRowDataVector.SaveVectorDataAsImagePlot( // "D:\\_gulevlab\\SkyImagesAnalysis_appData\\patent-samples\\result.2015-03-24\\img-2014-09-20T16-03-58devID1\\dvRowDataVector-plot-image-" + // i.ToString("D03") + "-step3.png"); #endregion // debug plotting double minValue = dvRowDataVector.Minimum(); int minValueIndex = dvRowDataVector.MinimumIndex(); //if (!itsTheCropCase) dmFieldminimumsData[i, minValueIndex] = minValue; if ((!itsTheCropCase) && ((double)minValueIndex > sunDiskData.DRadius)) { lRetPoints.Add(new Point3D(currentAngle, minValueIndex, minValue)); } else { continue; } } #region // obsolete //} //else if (dimensionNumber == 2) //{ // for (int i = 0; i < dmFieldData.ColumnCount; i++) // { // DenseMatrix dmSlicedDataMatrix = (DenseMatrix)dmFieldData.SubMatrix(0, dmFieldData.RowCount, i, 1); // DenseVector dvSlicedDataVector = DenseVector.OfEnumerable(dmSlicedDataMatrix.Values); // dvSlicedDataVector.MapInplace(new Func<double, double>(x => (x == 0.0d) ? (1.0d) : (x))); // double minValue = dvSlicedDataVector.Minimum(); // int minValueIndex = dvSlicedDataVector.MinimumIndex(); // dmFieldminimumsData[minValueIndex, i] = minValue; // } //} #endregion // obsolete //return dmFieldminimumsData; return(lRetPoints); }