开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 中间件编程 > 查看源码
Smooth.cpp
资源名称：08.zip　[点击查看]
上传用户：ynjin1970
上传日期：2014-10-13
资源大小：6438k
文件大小：10k
源码类别：

中间件编程

开发平台：

Visual C++

Smooth.cpp：源码内容
  1. // Smooth.cpp: implementation of the CSmooth class.
  2. //
  3. //////////////////////////////////////////////////////////////////////
  5. #include "stdafx.h"
  6. #include "Onecontour.h"
  7. #include <math.h>
  8. #include "Smooth.h"
  10. #ifdef _DEBUG
  11. #undef THIS_FILE
  12. static char THIS_FILE[]=__FILE__;
  13. #define new DEBUG_NEW
  14. #endif
  16. //////////////////////////////////////////////////////////////////////
  17. // Construction/Destruction
  18. //////////////////////////////////////////////////////////////////////
  20. CSmooth::CSmooth(XYZ* points, int count)
  21. {
  22. // 初始哈原始点和点数量
  23. points_ = points;
  24. count_ = count;
  26. ptArrayCount = 0;
  28. // other initial
  29. ssum = 0;
  30. dts=20;
  32. }
  34. CSmooth::~CSmooth()
  35. {
  37. }
  39. /*
  40. 获得平滑过后的点
  41. */
  42. void CSmooth::GetSmoothedPoints(XYZ **points, int& count)
  43. {
  44. count = ptArrayCount;
  45. *points = ptArray;
  46. }
  49. //平滑由多个点组成的折线, 并把这些平滑的线画在DC上
  50. void CSmooth::Smooth()
  51. {
  52. double x, y, s(0);
  53. int i;
  54. int MAXSIZE = 512; // ptArray的最大大小，动态增加
  55. XYZ *tmpXYZ;
  57. ptArray = new XYZ[MAXSIZE]; // 不释放，将在哪儿释放呢？
  59. // bakup points and count for remove repeat points later.
  60. int count = count_;
  61. XYZ* points = new XYZ[count];
  63. for(i=0; i<count; i++)
  64. points[i] = points_[i];
  65. RemoveRepeatPoints(points, count);//规范化处理，剔除重复的点
  67. ssum = 0;
  68. ptArray[ptArrayCount++] = points[0]; //往数组中添加一个元素
  70. // 点的个数不够,退出平滑处理
  71. if(count < 2)
  72. {
  73. delete []points; return;
  74. }
  76. // 两个点。
  77. if(count == 2) 
  78. {
  79. ptArray[ptArrayCount++] = points[1];
  80. delete []points; return;
  81. } 
  83. // 只有三个点
  84. if(count==3) 
  85. {
  86. CalcPara(0, true, points);
  87. ax0 = ax; bx0 = bx; cx0 = cx;
  88. ay0 = ay; by0 = by; cy0 = cy;
  89. si00 = si0; si10 = si1; si20 = si2; si30 = si3;
  90. do 
  91. {
  92. s += dts;
  93. if (s>=si20) s = si20;
  94. x = x1(s); y = y1(s);
  95. ptArray[ptArrayCount++] = XYZ(x, y);
  96. if ( ptArrayCount >= MAXSIZE )
  97. {
  98. tmpXYZ = ptArray;
  99. ptArray = new XYZ[MAXSIZE + 512];
  100. memcpy(ptArray, tmpXYZ, sizeof(XYZ) * MAXSIZE);
  101. delete []tmpXYZ;
  102. MAXSIZE += 512;
  103. }
  104. } while(s<si20);
  105. delete []points; 
  106. return;
  107. }
  109. //有四个或更多的点
  110. //绘制第一个和第二个点之间的曲线
  111. CalcPara(0, true, points);
  112. ax0 = ax; bx0 = bx; cx0 = cx;
  113. ay0 = ay; by0 = by; cy0 = cy;
  114. si00 = si0; si10 = si1; si20 = si2;
  115. do 
  116. {
  117. s += dts;
  118. if (s>=si10) s=si10;
  119. x = x1(s); y = y1(s);
  120. ptArray[ptArrayCount++] = XYZ(x, y);
  121. if ( ptArrayCount >= MAXSIZE )
  122. {
  123. tmpXYZ = ptArray;
  124. ptArray = new XYZ[MAXSIZE + 512];
  125. memcpy(ptArray, tmpXYZ, sizeof(XYZ) * MAXSIZE);
  126. delete []tmpXYZ;
  127. MAXSIZE += 512;
  128. }
  129. } while(s<si10);
  130. //绘制中间点的曲线
  131. for (i=0; i<count-3; i++)
  132. {
  133. CalcPara(i, true, points);
  134. ax0 = ax; bx0 = bx; cx0 = cx;
  135. ay0 = ay; by0 = by; cy0 = cy;
  136. si00 = si0; si10 = si1; si20 = si2; si30 = si3;
  137. ssum = si10;
  138. CalcPara(i+1, false, points);
  140. s = si10;
  141. do
  142. {
  143. s += dts;
  144. if ( s >= si20) s = si20;
  145. x = a(s)*x1(s)+(1-a(s))*x2(s);
  146. y = a(s)*y1(s)+(1-a(s))*y2(s);
  147. ptArray[ptArrayCount++] = XYZ(x, y);
  148. ASSERT(ptArrayCount <= 10000);
  149. if ( ptArrayCount >= MAXSIZE )
  150. {
  151. tmpXYZ = ptArray;
  152. ptArray = new XYZ[MAXSIZE + 512];
  153. memcpy(ptArray, tmpXYZ, sizeof(XYZ) * MAXSIZE);
  154. delete []tmpXYZ;
  155. MAXSIZE += 512;
  156. }
  157. }
  158. while(s<si20);
  159. }
  160. //绘制最后的点与倒数第二个点之间的曲线
  161. s = si20;
  162. do 
  163. {
  164. s += dts;
  165. if (s >= si2) s = si2;
  166. x = x2(s); y = y2(s);
  167. ptArray[ptArrayCount++] = XYZ(x, y);
  168. if ( ptArrayCount >= MAXSIZE )
  169. {
  170. tmpXYZ = ptArray;
  171. ptArray = new XYZ[MAXSIZE + 512];
  172. memcpy(ptArray, tmpXYZ, sizeof(XYZ) * MAXSIZE);
  173. delete []tmpXYZ;
  174. MAXSIZE += 512;
  175. }
  176. } while(s<si2);
  177. delete []points;
  178. }
  180. double CSmooth::x1(const double s)
  181. {
  182. return (ax0+bx0+cx0)*s*s-(ax0*(si10+si20)+bx0*(si00+si20)+cx0*(si00+si10))*s+(ax0*si10*si20+bx0*si20*si00+cx0*si00*si10);
  183. }
  185. double CSmooth::y1(const double s)
  186. {
  187. return (ay0+by0+cy0)*s*s-(ay0*(si10+si20)+by0*(si00+si20)+cy0*(si00+si10))*s+(ay0*si10*si20+by0*si20*si00+cy0*si00*si10);
  188. }
  190. double CSmooth::x2(const double s)
  191. {
  192. return (ax+bx+cx)*s*s-(ax*(si1+si2)+bx*(si0+si2)+cx*(si0+si1))*s+(ax*si1*si2+bx*si2*si0+cx*si0*si1);
  193. }
  195. double CSmooth::y2(const double s)
  196. {
  197. return (ay+by+cy)*s*s-(ay*(si1+si2)+by*(si0+si2)+cy*(si0+si1))*s+(ay*si1*si2+by*si2*si0+cy*si0*si1);
  198. }
  200. double CSmooth::a(const double s)
  201. {
  202. return (si20-s)/(si20-si10);
  203. }
  205. //算法中的处理
  206. void CSmooth::CalcPara(const int i, bool bCalcS, const XYZ* points)
  207. {
  208. si0 = ssum;
  209. si1=si0+sqrt((points[i+1].x-points[i].x)*(points[i+1].x-points[i].x)+(points[i+1].y-points[i].y)*(points[i+1].y-points[i].y));
  210. si2=si1+sqrt((points[i+2].x-points[i+1].x)*(points[i+2].x-points[i+1].x)+(points[i+2].y-points[i+1].y)*(points[i+2].y-points[i+1].y));
  211. // si3=si2+sqrt((points[i+2].x-points[i+3].x)*(points[i+2].x-points[i+3].x)+(points[i+2].y-points[i+3].y)*(points[i+2].y-points[i+3].y));
  212. ax = points[i].x/((si0-si1)*(si0-si2));
  213. bx = points[i+1].x/((si1-si0)*(si1-si2));
  214. cx = points[i+2].x/((si2-si0)*(si2-si1));
  215. ay = points[i].y/((si0-si1)*(si0-si2));
  216. by = points[i+1].y/((si1-si0)*(si1-si2));
  217. cy = points[i+2].y/((si2-si0)*(si2-si1));
  218. }
  220. //消除重复的点
  221. void CSmooth::RemoveRepeatPoints(XYZ* points, int& count)
  222. {
  223. for (int i= 0 ; i<count-1; i++)
  224. {
  225. if (points[i] == points[i+1])
  226. {
  227. for ( int j = i; j<count-1; j++)
  228. points[j] = points[j+1];
  229. count--;
  230. i--;
  231. }
  232. }
  233. }
  236. //画平滑线
  237. void CSmooth::DrawSmoothLines(CDC *dc)
  238. {
  239. /* CPoint point;
  240. point = ptArray.ElementAt(0);
  241. dc->MoveTo(point);
  242. for (int i=1; i<ptArray.GetSize(); i++)
  243. {
  244. point = ptArray.ElementAt(i);
  245. dc->LineTo(point);
  246. }
  247. */}
  249. //标注
  250. void CSmooth::DrawLable(CDC *dc, double zvalue)
  251. {
  252. // dc->Ellipse(ptArray.GetAt(0).x-50,ptArray.GetAt(0).y-50,ptArray.GetAt(0).x+50,ptArray.GetAt(0).y+50);
  253. //if (zvalue!=1750) return;
  254. /* int i;
  255. int kkkmmm=0;
  256. double bzjd_m;
  257. double dX1, dY1, dX2, dY2, dX3, dY3;
  258. const double pie = 3.1416;
  259. CString TextString;
  260. LOGFONT myFont;
  261.     double r1;
  262. bool bSearched;
  264. myFont.lfHeight = 780;
  265.     myFont.lfWidth = 460;
  266.     myFont.lfEscapement = 0;
  267.     myFont.lfWeight = 400;
  268.     myFont.lfItalic = 0;
  269.     myFont.lfUnderline = 0;
  270.     myFont.lfStrikeOut = 0;
  271.     myFont.lfOutPrecision = 0;
  272.     myFont.lfClipPrecision = 0;
  273.     myFont.lfQuality = 0;
  274.     myFont.lfPitchAndFamily = 0;
  275.     myFont.lfCharSet = 0;
  276. strcpy(myFont.lfFaceName, "Arial");            
  278. //寻找标注点的位置及字体的标注角度
  279.     r1 = 0;
  280.     int nnk = 0;
  281.     double RS = 0;
  282.     i = ptArray.GetSize() - 1;
  283.     TextString.Format("%7.2f", zvalue);
  285. int iWidth = TextString.GetLength()*myFont.lfWidth;
  286. int iHeight = myFont.lfHeight;
  287.     int iNeedDots = (int)floor(iWidth/dts)+1;
  289. if (i >= 10)
  290. {
  291. for (int t = 5; t<=i-10; t++)
  292. {
  293. bSearched = false;//假设没有找到。。。
  294. int iik=t;
  295. dX1 = ptArray.GetAt(iik).x;
  296. dY1 = ptArray.GetAt(iik).y;
  297. dX2 = ptArray.GetAt(iik+1).x;
  298. dY2 = ptArray.GetAt(iik+1).y;
  299. dX3 = ptArray.GetAt(iik+2).x;
  300. dY3 = ptArray.GetAt(iik+2).y;
  301. double dist=sqrt((dX1-dX2)*(dX1-dX2)+(dY1-dY2)*(dY1-dY2))+sqrt((dX3-dX2)*(dX3-dX2)+(dY3-dY2)*(dY3-dY2));
  302. double mmk, kkm, xc, yc, ro;
  303. do 
  304. {
  305.   double pf1 = dX1*dX1 + dY1*dY1 - dX2*dX2 - dY2*dY2;
  306.   double pf2 = dX2*dX2 + dY2*dY2 - dX3*dX3 - dY3*dY3;
  307.   if (dY1 ==dY2) goto lbl100;
  308.   mmk = pf2 - ((dY2 - dY3) / (dY1 - dY2)) * pf1;
  309.   kkm = 2 * (dX2 - dX3 - (dY2 - dY3) * (dX1 - dX2) / (dY1 - dY2));
  310.   if (fabs(kkm) < 0.0001) goto lbl100;
  311.   xc = mmk / kkm;//（xc,yc）为园心，ro为半径
  312.   yc = (pf2 - 2 * (dX1 - dX2) * xc) / (2 * (dY1 - dY2));
  313.   ro = sqrt((yc - dX2)*(yc - dX2) + (yc - dY2)*(yc - dY2));
  314.   if (ro<myFont.lfHeight*5) {bSearched=false;break;}
  316. lbl100:   //如果满足要求的话，继续。。。
  317.   bSearched=true;
  318.   if (dist>iWidth) break;
  319.   iik++;
  320.   if (iik+2>i-10) {bSearched=false; break;}//点的个数不够，停止追踪
  321.   dX1=dX2;dY1=dY2; dX2=dX3; dY2=dY3;
  322.   dX3 = ptArray.GetAt(iik+2).x;
  323.   dY3 = ptArray.GetAt(iik+2).y;
  324.   dist+=sqrt((dX3-dX2)*(dX3-dX2)+(dY3-dY2)*(dY3-dY2));
  325. // dc->Ellipse(ptArray.GetAt(iik+2).x-10,ptArray.GetAt(iik+2).y-10,ptArray.GetAt(iik+2).x+10,ptArray.GetAt(iik+2).y+10);
  326. }while(1);
  327. if (bSearched == true) {kkkmmm = t; break;}
  328. }
  329.         if (bSearched == false || kkkmmm==0) 
  330. {
  331. return;
  332. }
  333. dc->Ellipse(ptArray.GetAt(kkkmmm).x-40,ptArray.GetAt(kkkmmm).y-40,ptArray.GetAt(kkkmmm).x+40,ptArray.GetAt(kkkmmm).y+40);
  335. double xzwz = ptArray.GetAt(kkkmmm).x;
  336. double yzwz = ptArray.GetAt(kkkmmm).y;
  338. double dltx = ptArray.GetAt(kkkmmm).x - ptArray.GetAt(kkkmmm+iNeedDots).x;
  339. double dlty = ptArray.GetAt(kkkmmm).y - ptArray.GetAt(kkkmmm+iNeedDots).y;
  340. double fzqz;
  341.    
  342.         if (dltx != 0)
  343. {
  344. fzqz = -atan(dlty/dltx);
  345. if (dlty == 0) fzqz = 0;
  346. }
  347.         else
  348.              double fzqz = 3 * pie / 2;
  349.         
  350.         bzjd_m = -fzqz;
  351.         if ((fzqz<pie/2 && fzqz>0) || (fzqz<-3*pie/2 && fzqz>-2*pie) || (fzqz>pie && fzqz<3*pie/2) || (fzqz>-pie && fzqz<-pie/2))
  352. {
  353. // xzwz += (sqrt(iHeight*iHeight+iWidth*iWidth)/2 * cos(fzqz - atan(iHeight/iWidth)));
  354. // yzwz += (sqrt(iHeight*iHeight+iWidth*iWidth)/2 * sin(fzqz - atan(iHeight/iWidth)));
  355. //xzwz = xzwz - abs(myFont.lfWidth * cos(fzqz)) - abs(myFont.lfHeight * sin(fzqz) / 2);
  356. // yzwz = yzwz - abs(myFont.lfHeight * cos(fzqz) / 2) + abs(myFont.lfWidth * sin(fzqz));
  357. }
  358.         else
  359. {
  360. // xzwz -= (sqrt(iHeight*iHeight+iWidth*iWidth)/2 * cos(fzqz - atan(iHeight/iWidth)));
  361. // yzwz += (sqrt(iHeight*iHeight+iWidth*iWidth)/2 * sin(fzqz - atan(iHeight/iWidth)));
  362. // xzwz = xzwz - abs(myFont.lfWidth * cos(fzqz)) + abs(myFont.lfHeight * sin(fzqz) / 2);
  363. // yzwz = yzwz - abs(myFont.lfHeight * cos(fzqz) / 2) - abs(myFont.lfWidth * sin(fzqz));
  364. }
  366.         myFont.lfEscapement =(long)( fzqz * 1800.0 / pie);
  367.        //标注          
  368.         CFont NewFont;
  369. NewFont.CreateFontIndirect(&myFont);
  370.         CFont *OrgFont = dc->SelectObject(&NewFont);
  371.         dc->TextOut((int)xzwz, (int)yzwz, TextString, 5); //TextString.GetLength());
  372.         dc->SelectObject(OrgFont);
  373. NewFont.DeleteObject();
  374. OrgFont->DeleteObject();
  375. } 
  376. */}