--- a
+++ b/Thoracic Organs Segmentation code/DrawCircle.m
@@ -0,0 +1,26 @@
+function [Img] = DrawCircle(R,c1,c2,x,y,lineWidth);
+Img = zeros(x,y);
+theta = [0 : 0.01 : 2 * pi];
+stack1 = 0;
+if nargin == 5
+    
+    lineWidth = 1;
+end
+for i = 1 : length( theta )
+
+        temp_x = 0; temp_y = 0;   
+             
+        temp_x =ceil( R* cos (theta(i)) + c1);
+        
+        temp_y = ceil (R* sin (theta(i)) + c2);
+        stack1 = stack1 + 1;
+        I_new(stack1) = temp_x ;
+        J_new(stack1) = temp_y  ;
+        Img( temp_x , temp_y ) = 1;
+       
+        %plot(I_new,J_new,'.',c1,c2,'r*')
+end
+
+se1 = strel('line',lineWidth,0);
+se2 = strel('line',lineWidth,90);
+Img = imdilate(Img,[se1 se2],'same');