opencv 实现的椭圆检测的源代码

1. / readbmp.cpp : Defines the entry point for the console application. 2. / 3. 4. #include “stdafx.h“ 5. #include 6. 7. #include “ipl.h“ 8. #include “cv.h“ 9. /#include “_cv.h“ 10. /#include “vlgrfmts.h“ 11. 12. /- 13. /UTILITY FUNCTIONS 14. / 15. / Read in a BMP and convert it into an IPL 16. / fin The file name of the BMP with path 17. / plannar If set, convert the image to plannar 18. / flip If set, mirror the image around the horizontal axis, set 19. / IPL_ORIGIN_BL) 20. IplImage * readBMP2IPL(char *fin, int plannar = 0, int flip = 0) 21. 22. int filter; 23. int width=0,height=0,color=55,hdr_ret=0,close_ret = 0,rd_ret = 0; 24. IplImage *I; 25. filter = gr_fmt_find_filter(fin ); /Get file type 26. hdr_ret = gr_fmt_read_header( filter, /image p arams 27. I = cvCreateImageHeader( cvSize(width, height), IPL_DEPTH_8U, 3); /creat e image header 28. I- align = IPL_ALIGN_QWORD; /fix the problem that cvCreateImage uses align 4 29. I- depth != IPL_DEPTH_32F ? iplAllocateImage( I, 0, 0 ) :iplAllocateImageFP( I, 0, 0 ); 30. rd_ret = gr_fmt_read_data(filter, I-imageData, I- widthStep, color); /read in the data, origin TL 31. close_ret = gr_fmt_close_filter( filter ); 32. if(plannar) /If plannar order desired, turn it into plannar order 33. 34. IplImage *Ip = iplCreateImageHeader(3, 0, IPL_DEPTH_8U, “RGB“,“BGR“, IPL_DATA_ORDER_PLANE, I- origin, IPL_ALIGN_QWORD, width, height, NULL,NULL,NULL,NULL); 35. iplAllocateImage(Ip, 0,0); 36. iplConvert(I,Ip); 37. cvReleaseImage( 38. I = Ip; 39. 40. if(flip) /User wants IPL_ORIGIN_BL 41. 42. iplMirror(I, I, 0); /flipporuni 43. I-origin = IPL_ORIGIN_BL; 44. 45. return I; 46. 47. 48. / 49. / Write IPL image to disk as a BMP image (not optimized but step by step sim plistic) 50. / fout Pathfilename to write to 51. / I IPL image to be written 52. / 53. / Image fomat may be any of: 54. / dataOrder Plane or Pixel 55. / 1 channel (GRAY) or 3 (RGB) 56. / 8U, 8S, 32F 57. / Origin TL or BL (if BL, image is flipped to conform to TL bmps) 58. / 59. / If image is not RGB or GRAY, RGB is simply imposed on the Image. E.g. if i ts YUV, 60. / the image will be stored as BGR with B=Y, G=U, and R=V 61. void writeIPL2BMP(char *fout, IplImage *I) 62. 63. IplImage *Isav,*Itmp; 64. Isav = iplCloneImage(I); 65. /Convert to pixel order if needed 66. if(I-dataOrder = IPL_DATA_ORDER_PLANE) 67. 68. Itmp = iplCreateImageHeader(I-nChannels, 69. I-alphaChannel, I-depth, I-colorModel, 70. I-channelSeq, IPL_DATA_ORDER_PIXEL, I-origin, 71. I-align, 72. I-width, I-height, NULL, NULL,NULL,NULL); 73. if(I-depth = IPL_DEPTH_32F) 74. iplAllocateImageFP(Itmp, 1,0.0); 75. else 76. iplAllocateImage(Itmp, 1,0); 77. iplConvert(I, Itmp); 78. iplDeallocate (I, IPL_IMAGE_ALL); 79. I = Itmp; 80. 81. /Image is now pixel order. Convert to 8U if needed. 82. if(I-depth = IPL_DEPTH_8S) 83. I-depth = IPL_DEPTH_8U; 84. else if(I-depth = IPL_DEPTH_32F) 85. 86. float min,max; 87. iplMinMaxFP (I, 88. Itmp = iplCreateImageHeader(I-nChannels, 89. I-alphaChannel, IPL_DEPTH_8U, I-colorModel, 90. I-channelSeq, IPL_DATA_ORDER_PIXEL, I-origin, 91. I-align, 92. I-width, I-height, NULL, NULL,NULL,NULL); 93. iplAllocateImage(Itmp, 1,0); 94. if(min = max) max = (float)(min + 1.0); 95. iplScaleFP (I, Itmp, min, max); 96. iplDeallocate (I, IPL_IMAGE_ALL); 97. I = Itmp; 98. 99. /Image is now pixel order, 8U. Convert to RGB if needed 100. if(I-nChannels = 3) 101. 102. strcpy(I-colorModel,“RGB“); 103. strcpy(I-channelSeq,“BGR“); 104. 105. else if(I-nChannels = 1) 106. 107. Itmp = iplCreateImageHeader(3, 108. I-alphaChannel, IPL_DEPTH_8U, “RGB“, 109. “BGR“, IPL_DATA_ORDER_PIXEL, I-origin, I-align, 110. I-width, I-height, NULL, NULL,NULL,NULL); 111. iplAllocateImage(Itmp, 1,0); 112. iplGrayToColor (I, Itmp, 1.0, 1.0, 1.0); 113. iplDeallocate (I, IPL_IMAGE_ALL); 114. I = Itmp; 115. 116. else 117. return; 118. /Image is now pixel order, 8U, RGB, flip it if the origin is BL 119. if(I-origin = IPL_ORIGIN_BL) 120. 121. iplMirror(I, I, 0); 122. 123. /FINALLY, WRITE THE IMAGE: 124. int write_ret = gr_fmt_write_image( I-imageData, I-widthStep, 125. I-width, I-height, I-colorModel0 =R?1:0,fout, “bmp“); 126. /Return image to its original state 127. iplDeallocate(I, IPL_IMAGE_ALL); 128. I = Isav; 129. 130. 131./- 132. 133.void main() 134. 135. IplImage * input_img=NULL; 136. int plannar=0; 137. int flip=0; 138. char fin=“circles.bmp“; 139. char fout=“result.bmp“; 140. 141. 142. /作一些声明 143. CvMemStorage *storage; 144. int header_size, i, count,length, width; 145. CvSeq *contour; 146. CvBox2D * myBox; /用于画圆和椭圆 147. CvPoint *PointArray; 148. CvPoint2D32f *PointArray32f; 149. CvPoint myCenter; 150. 151. / 分配内存 152. myBox= (CvBox2D *) malloc(sizeof(CvBox2D); 153. header_size = sizeof(CvContour); 154. storage = cvCreateMemStorage (1000); / For FindContours. 155. 156. / 读入图象 input_img 并转换为 IPLIMAGE 格式 157. /已经定义 char fin=“circles.bmp“; 158. input_img=readBMP2IPL(fin,plannar, flip); 159. 160. coutnChannelscenter.x; 162. myCenter.y = (int) myBox-center.y; 163. 164. /HEIGHT AND WITH (THE 2 AXES OF THE ELLIPSE) 165. length =(int)myBox-size.height; 166. width = (int)myBox-size.width; 167. 168. /AND THE ANGLE 169. float myAngle= myBox-angle; 170. 171. / 检查这个椭圆是否在图象内 172. if (myCenter.x 0) 176. 177. /画椭圆 178. cvEllipse(input_img, 179. myCenter, 180. cvSize (int)width/2,(int)length/2) , 181. myBox-angle, 182. 0,360,RGB(255,0,255),1); 183. 184. 185. 186. free(PointArray32f ); 187. fre