APPENDIX I – The concerning section of the executed IDL code to define the maximum radius of interest and implement the coordinate transformation from Cartesian to cylindrical one for a GDL slice at a particular height (in z-direction): Radius=-1 Angle_step=-1 read_tiff(Slice_image) image = float(image) image = size(image) Image_size = xsize = Image_size(1) ysize = Image_size(2) if EQ then (Radius -1) radius= endif if EQ then begin (Angle_step -1) Angle_step=360/(Radius*2*!PI) endif Angle_range=360/Angle_step Array=fltarr(radius+1,Angle_range) 0 = image > for do begin r=0,radius ;printf, lun, 'Radius with respect to Reference Coordinate:', r r for do begin w=0,(360/Angle_step)-1 ;print, w, r theta=w*Angle_step X=Reference_Point_X_coordinate+r*cos(theta*!PI/180) Y=Reference_Point_Y_coordinate+r*sin(theta*!PI/180) Array(r,w)=interpolate(image,X,Y,/GRID ) endfor endfor
APPENDIX – I The concerning section of the executed IDL code to define the maximum radius of interest and implement the coordinate transformation from Cartesian to cylindrical one for a GDL slice at a particular height (in z-direction): Radius=-1 Angle_step=-1 image = read_tiff(Slice_image) image = float(image) Image_size = size(image) xsize = Image_size(1) ysize = Image_size(2) if (Radius EQ -1) then radius= endif if (Angle_step EQ -1) then begin Angle_step=360/(Radius*2*!PI) endif Angle_range=360/Angle_step Array=fltarr(radius+1,Angle_range) = image > 0 for r=0,radius do begin ;printf, lun, 'Radius with respect to Reference Coordinate:', r r for w=0,(360/Angle_step)-1 do begin ;print, w, r theta=w*Angle_step X=Reference_Point_X_coordinate+r*cos(theta*!PI/180) Y=Reference_Point_Y_coordinate+r*sin(theta*!PI/180) Array(r,w)=interpolate(image,X,Y,/GRID ) endfor endfor
1-s2.0-S0360319917343641-mmc1.docx