Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1799 1302 (N33W51 L=199)   471     22        1 (1778/1282-1819/1324)  21.8 W48-W53 N31-N34
      2 1614 1192 (N21W34 L=183)   718    135        1 (1573/1129-1649/1257) -35.5 W30-W37 N17-N24
      3 1834 1207 (N27W54 L=203)   511     73        1 (1813/1185-1854/1236) -26.6 W52-W56 N25-N29
      4 1852 1116 (N22W57 L=206)   479     66        1 (1829/1093-1877/1136)  36.9 W54-W60 N20-N23
      5  377 1004 (S16E42 L=107)   279     45        1 ( 354/ 983- 395/1032)  56.3 E43-E40 S17-S14
      6  461  886 (S22E33 L=115)   372    153        0 ( 400/ 863- 518/ 906)  18.4 E38-E29 S22-S19
      7  346  836 (S27E42 L=106)   581    147        1 ( 315/ 773- 372/ 889)  63.4 E45-E40 S31-S23
      8 1171  810 (S11W13 L=162)   308    104        0 (1116/ 788-1223/ 830)  38.7 W10-W15 S13-S09
      9  940  688 (S24W02 L=151)   221     75        0 ( 902/ 675- 979/ 706)  36.9 W00-W04 S24-S21
     10 1523  479 (S23W48 L=196)  1808    382        2 (1376/ 399-1659/ 603)  57.7 W40-W53 S31-S12
     11  997  469 (S37W12 L=161)   526    143        0 ( 920/ 452-1072/ 483)  21.8 W06-W17 S38-S34
  

    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