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

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > GDI/图象编程 > 查看源码
C_QGISAlgorithmLib.cpp
资源名称：QGIS.rar　[点击查看]
上传用户：oybseng
上传日期：2015-04-27
资源大小：7831k
文件大小：10k
源码类别：

GDI/图象编程

开发平台：

Visual C++

C_QGISAlgorithmLib.cpp：源码内容
  1. // QGISAlgorithmLib->cpp: implementation of the CQGISAlgorithmLib class.
  2. //
  3. //////////////////////////////////////////////////////////////////////
  5. #include "stdafx.h"
  7. #include "....QObjectsincludeResource.h"
  8. #include "..includeResource.h"
  9. #include "..includeC_QGISAlgorithmLib.h"
  10. #include "....QObjectsincludeQLineObj.h"
  12. #ifdef _DEBUG
  13. #undef THIS_FILE
  14. static char THIS_FILE[]=__FILE__;
  15. #define new DEBUG_NEW
  16. #endif
  18. //////////////////////////////////////////////////////////////////////
  19. // Construction/Destruction
  20. //////////////////////////////////////////////////////////////////////
  22. CQGISAlgorithmLib::CQGISAlgorithmLib()
  23. {
  25. }
  27. CQGISAlgorithmLib::~CQGISAlgorithmLib()
  28. {
  30. }
  32. //////////////////////////////////////////
  33. ///***************CQGIS****************///
  34. ///函数名称:IsLineCross
  35. ///返回类型:BOOL
  36. ///入口参数:CQPoint p1,CQPoint p2,CQPoint q1,CQPoint q2
  37. ///出口参数:无
  38. ///思路说明:边界矩形，矢量跨立实验
  39. ///***************CQGIS****************///
  40. //////////////////////////////////////////
  41. BOOL CQGISAlgorithmLib::IsLineCross(CQPoint p1,CQPoint p2,CQPoint q1,CQPoint q2)
  42. {
  43. //创建两条直线对象 用于进行边界矩形的判断
  44. if(p1 == q1 || p1 == q2 || p2 == q1 || p2 == q2)
  45. return TRUE; //直接得出直线相交
  47. CQLineObj * pLine1 = new CQLineObj;
  48. pLine1->m_PtList.AddPoint(&p1);
  49. pLine1->m_PtList.AddPoint(&p2);
  51. CQLineObj * pLine2 = new CQLineObj;
  52. pLine2->m_PtList.AddPoint(&q1);
  53. pLine2->m_PtList.AddPoint(&q2);
  55. CBoundaryRect * pRect1 = new CBoundaryRect;
  56. pLine1->GetBoundingRect(pRect1);
  57. CBoundaryRect * pRect2 = new CBoundaryRect;
  58. pLine2->GetBoundingRect(pRect2);
  60. if(!pRect1->IsCross(pRect2))
  61. {
  62. //假如两个边界矩形不相交
  63. delete pRect2;
  64. delete pRect1;
  65. delete pLine2;
  66. delete pLine1;
  67. return FALSE; //不相交
  68. }
  70. //进行跨立实验
  71. double p1q1x = p1.GetX() - q1.GetX();
  72. double p1q1y = p1.GetY() - q1.GetY();
  73. double q2q1x = q2.GetX() - q1.GetX();
  74. double q2q1y = q2.GetY() - q1.GetY();
  75. double p2q1x = p2.GetX() - q1.GetX();
  76. double p2q1y = p2.GetY() - q1.GetY();
  78. double fResult;  //运算结果
  79. double fpq = p1q1x*q2q1y - q2q1x*p1q1y;
  80. double fqp = q2q1x*p2q1y - p2q1x*q2q1y;
  81. fResult = fpq*fqp;
  82. if(fResult>=0)
  83. {
  84. delete pRect2;
  85. delete pRect1;
  86. delete pLine2;
  87. delete pLine1;
  88. return TRUE;
  89. }
  91. delete pRect2;
  92. delete pRect1;
  93. delete pLine2;
  94. delete pLine1;
  95. return FALSE;
  96. }
  98. //////////////////////////////////////////
  99. ///***************CQGIS****************///
  100. ///函数名称:CalLineCrossPoint
  101. ///返回类型:无
  102. ///入口参数:CQPoint p1,CQPoint p2,CQPoint q1,CQPoint q2
  103. ///出口参数:CQPointArray * pResult
  104. ///思路说明:判断两条直线是否相交，若相交则返回其交点
  105. ///***************CQGIS****************///
  106. //////////////////////////////////////////
  107. void CQGISAlgorithmLib::CalLineCrossPoint(CQPoint pp1,CQPoint pp2,CQPoint qq1,CQPoint qq2,CQPointArray * pResult)
  108. {
  109. if(!IsLineCross(pp1,pp2,qq1,qq2))
  110. {
  111. pResult = 0;
  112. return;
  113. }
  115. CQPoint * p1 = new CQPoint(pp1);
  116. CQPoint * p2 = new CQPoint(pp2);
  117. CQPoint * q1 = new CQPoint(qq1);
  118. CQPoint * q2 = new CQPoint(qq2);
  120. //计算两条直线的斜率
  121. //double k1 = (p2.GetY()-p1.GetY())/(p2.GetX()-p1.GetX());
  122. //double k2 = (q2.GetY()-q1.GetY())/(q2.GetX()-q1.GetX());
  124. //L0(p1,p2),L1(q1,q2);
  125. //情况1 L0平行于Y轴
  126. //分两种情况 当L1也与Y轴平行时，则两条直线共线 
  127. //若L1与Y轴不平行，则直接代入方程求解
  128. if(p1->GetX() == p2->GetX())
  129. {
  130. if(q1->GetX() == q2->GetX()) //L1也与Y轴平行
  131. {
  132. CalSameLineCrossPoint(*p1,*p2,*q1,*q2,pResult);
  133. return;
  134. }
  135. else      //L1与Y轴不平行
  136. {
  137. double k2 = (q2->GetY()-q1->GetY())/(q2->GetX()-q1->GetX());
  138. double x2 = p1->GetX();
  139. //y2 = k2x2-k2x1+y1;x1 = q1->GetX(),y1 = q1->getY();
  140. double y2 = k2*x2-k2*q1->GetX()+q1->GetY();
  141. pResult->AddPoint(x2,y2);
  142. delete q2;
  143. }
  144. }
  145. else if(q1->GetX() == q2->GetX())
  146. {
  147. double k1 = (p2->GetY()-p1->GetY())/(p2->GetX()-p1->GetX());
  148. double x1 = q1->GetX(),y1;
  149. y1 = k1*x1-k1*p1->GetX()+p1->GetY();
  150. pResult->AddPoint(x1,y1);
  151. }
  152. else if(p1->GetY() == p2->GetY()) //平行与X 轴
  153. {
  154. if(q1->GetY() == q2->GetY())  //L1也平行与X轴
  155. {
  156. CalSameLineCrossPoint(*p1,*p2,*q1,*q2,pResult);
  157. }
  158. else
  159. {
  160. double k2 = (q2->GetY()-q1->GetY())/(q2->GetX()-q1->GetX());
  161. double y2 = p1->GetY();
  162. double x2 = (y2+k2*q1->GetX()-q1->GetY())/k2;
  163. pResult->AddPoint(x2,y2);
  164. }
  165. }
  166. else if(q1->GetY() == q2->GetY())
  167. {
  168. double k1 = (p2->GetY()-p1->GetY())/(p2->GetX()-p1->GetX());
  169. double y1 = q1->GetY();
  170. double x1 = (y1+k1*p1->GetX()-p1->GetY())/k1;
  171. pResult->AddPoint(x1,y1);
  172. }
  173. else
  174. {
  175. double k1 = (p2->GetY()-p1->GetY())/(p2->GetX()-p1->GetX());
  176. double k2 = (q2->GetY()-q1->GetY())/(q2->GetX()-q1->GetX());
  177. if(k1-k2 == SMALL_NUMBER)
  178. {
  179. CalSameLineCrossPoint(*p1,*p2,*q1,*q1,pResult);
  180. }
  181. else
  182. {
  183. double xResult = k1*p1->GetX()-k2*q1->GetX()+q1->GetY()-p1->GetY();
  184. xResult = xResult/(k1-k2);
  185. double yResult = k1*xResult-k1*p1->GetX()+p1->GetY();
  186. pResult->AddPoint(xResult,yResult);
  187. }
  188. }
  190. delete q2;
  191. delete q1;
  192. delete p2;
  193. delete p1;
  194. }
  196. //////////////////////////////////////////
  197. ///***************CQGIS****************///
  198. ///函数名称:CalSameLineCrossPoint
  199. ///返回类型:无
  200. ///入口参数:CQPoint p1,CQPoint p2,CQPoint q1,CQPoint q2
  201. ///出口参数:CQPointArray * pResult
  202. ///思路说明:
  203. ///***************CQGIS****************///
  204. //////////////////////////////////////////
  205. void CQGISAlgorithmLib::CalSameLineCrossPoint(CQPoint pp1,CQPoint pp2,CQPoint qq1,CQPoint qq2,CQPointArray * pResult)
  206. {
  207. //直接计算 不做条件判断 在调用之前一定要确保两直线是相交的
  208. //共线计算只考虑垂直于X，Y轴的情形
  210. double fy1,fy2;
  211. double fy3,fy4;
  212. double fx1,fx2;
  213. double fx3,fx4;
  215. CQPoint * p1 = new CQPoint(pp1);
  216. CQPoint * p2 = new CQPoint(pp2);
  217. CQPoint * q1 = new CQPoint(qq1);
  218. CQPoint * q2 = new CQPoint(qq2);
  220. //先找出较长的直线
  221. double dd1 = p1->Distance(*p2);
  222. double dd2 = q1->Distance(*q2);
  223. //始终用较长的线段为基准进行判断
  224. if(dd1<dd2)
  225. {
  226. CQPoint tmPoint = *p1;
  227. *p1 = *q1;
  228. *q1 = tmPoint;
  229. tmPoint = *p2;
  230. *p2 = *q2;
  231. *q2 = tmPoint;
  232. }
  234. if(p1->GetX() == p2->GetX() && q1->GetX() == q2->GetX())  //两条直线都与Y轴平行 
  235. {
  236. fy1 = min(p1->GetY(),p2->GetY());
  237. fy2 = max(p1->GetY(),p2->GetY());
  238. fy3 = min(q1->GetY(),q2->GetY());
  239. fy4 = max(q1->GetY(),q2->GetY());
  240. p1->SetPoint(p1->GetX(),fy1);  //Y较小的点
  241. p2->SetPoint(p2->GetX(),fy2);  //Y较大的点
  242. q1->SetPoint(q1->GetX(),fy3);
  243. q2->SetPoint(q2->GetX(),fy4);
  245. if(p1->GetY() > q2->GetY() || p2->GetY() < q1->GetY())
  246. {
  247. pResult = 0;
  248. }
  249. else
  250. {
  251. if(p2->GetY() == q1->GetY())
  252. {
  253. pResult->AddPoint(p2->GetX(),p2->GetY());
  254. }
  255. else if(q1->GetY()<p2->GetY() && q1->GetY() >p1->GetY())
  256. {
  257. if(p2->GetY()>q1->GetY() && p2->GetY()<q2->GetY())
  258. {
  259. pResult->AddPoint(q1->GetX(),q1->GetY());
  260. pResult->AddPoint(p2->GetX(),p2->GetY());
  261. }
  262. else if(p2->GetY()>=q2->GetY())
  263. {
  264. pResult->AddPoint(q1->GetX(),q1->GetY());
  265. pResult->AddPoint(q2->GetY(),q2->GetY());
  266. }
  268. }
  269. else if(*p1 == *q1 && *p2 == *q2)  //两直线重合
  270. {
  271. pResult->AddPoint(p1->GetX(),p1->GetY());
  272. pResult->AddPoint(p2->GetX(),p2->GetY());
  273. }
  274. }
  275. }
  276. else if(p1->GetY() == p2->GetY() && q1->GetY() == q2->GetY())  //两条直线都与X轴平行 
  277. {
  278. fx1 = min(p1->GetX(),p2->GetX());
  279. fx2 = max(p1->GetX(),p2->GetX());
  280. fx3 = min(q1->GetX(),q2->GetX());
  281. fx4 = max(q1->GetX(),q2->GetX());
  282. p1->SetPoint(fx1,p1->GetY());
  283. p2->SetPoint(fx2,p2->GetY());
  284. q1->SetPoint(fx3,q1->GetY());
  285. q2->SetPoint(fx4,q2->GetY());
  286. if(p2->GetX() < q1->GetX() || p1->GetX() > q2->GetX())
  287. {
  288. pResult = 0;
  289. }
  290. else if(q1->GetX()>p1->GetX() && q1->GetX()<p2->GetX())
  291. {
  292. if(p2->GetX()>q1->GetX() && p2->GetX()<q2->GetX())
  293. {
  294. pResult->AddPoint(q1->GetX(),q1->GetY());
  295. pResult->AddPoint(p2->GetX(),p2->GetY());
  296. }
  297. else if(p2->GetX()>=q2->GetX())
  298. {
  299. pResult->AddPoint(q1->GetX(),q1->GetY());
  300. pResult->AddPoint(q2->GetX(),q2->GetY());
  301. }
  302. }
  303. else if(p2->GetX() == q1->GetX())
  304. {
  305. pResult->AddPoint(p2->GetX(),p2->GetY());
  306. }
  307. else if(*p1 == *q1 && *p2 == *q2)
  308. {
  309. pResult->AddPoint(p1->GetX(),p2->GetY());
  310. pResult->AddPoint(p2->GetX(),p2->GetY());
  311. }
  312. }
  313. else   //两条直线只是斜率相等
  314. {
  315. if(p1->GetX()>p2->GetX())
  316. {
  317. CQPoint tmPoint;
  318. tmPoint = *p1;
  319. *p1 = *p2;
  320. *p2 = tmPoint;
  321. }
  322. if(q1->GetX()>q2->GetX())
  323. {
  324. CQPoint tmPoint;
  325. tmPoint = *q1;
  326. *q1 = *q2;
  327. *q2 = tmPoint;
  328. }
  329. if(p2->GetX()<q1->GetX())
  330. {
  331. CQPoint tmPoint = *q1;
  332. *p1 = *q1;
  333. *q1 = tmPoint;
  334. tmPoint = *q2;
  335. *p2 = *q2;
  336. *q2 = tmPoint;
  337. }
  339. if(p2->GetX() < q1->GetX() || p1->GetX() > q2->GetX())
  340. {
  341. pResult = 0;
  342. return;
  343. }
  344. else if(q1->GetX()>p1->GetX() && q1->GetX()<p2->GetX())
  345. {
  346. if(p2->GetX()>q1->GetX() && p2->GetX()<q2->GetX())
  347. {
  348. pResult->AddPoint(q1->GetX(),q1->GetY());
  349. pResult->AddPoint(p2->GetX(),p2->GetY());
  350. }
  351. else if(p2->GetX()>=q2->GetX())
  352. {
  353. pResult->AddPoint(q1->GetX(),q1->GetY());
  354. pResult->AddPoint(q2->GetX(),q2->GetY());
  355. }
  356. }
  357. else if(p2->GetX() == q1->GetX())
  358. {
  359. pResult->AddPoint(p2->GetX(),p2->GetY());
  360. }
  361. else if(*p1 == *q1 && *p2 == *q2)
  362. {
  363. pResult->AddPoint(p1->GetX(),p1->GetY());
  364. pResult->AddPoint(p2->GetX(),p2->GetY());
  365. }
  366. }
  368. delete q2;
  369. delete q1;
  370. delete p2;
  371. delete p1;
  372. }
  374. //弧度转换为度单位
  375. double CQGISAlgorithmLib::RadianToDegree(double fRadian)
  376. {
  377. return fRadian*PIunder180;
  378. }
  380. double CQGISAlgorithmLib::DegreeToRadian(double fDegree)
  381. {
  382. return fDegree*PIover180;
  383. }