Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1364 1374 (N15W29 L=205)   308     97        0 (1321/1359-1415/1388) -18.4 W26-W32 N13-N16
      2  819 1308 (N25E05 L=171)  1091    224        1 ( 710/1248- 927/1352) -42.3 E11-W00 N19-N29
      3  606 1139 (N20E22 L=154)   204     54        0 ( 577/1122- 636/1157) -33.7 E23-E20 N18-N21
      4 1393  833 (S16W17 L=193)  1131    293        2 (1313/ 731-1441/ 946) -86.2 W13-W19 S22-S07
      5 1307  634 (S25W07 L=183)   269     34        1 (1284/ 609-1333/ 654) -56.3 W05-W08 S27-S23
      6 1053  441 (S31E14 L=162)   292     93        0 (1005/ 427-1106/ 462) -15.9 E17-E10 S32-S29
      7 1518  355 (S54W19 L=195)   673     98        1 (1488/ 330-1551/ 387) -53.1 W16-W24 S58-S50
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y