Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  513 1714 (N42E54 L=107)  1640    260        2 ( 445/1679- 580/1739)   4.1 E62-E48 N39-N44
      2 1107 1558 (N37W01 L=161)  1010    213        2 ( 988/1533-1204/1581)  11.3 E07-W08 N35-N39
      3 1131 1270 (N18W05 L=165)   332    139        1 (1081/1220-1190/1291)  39.8 W01-W08 N14-N19
      4  885 1071 (N04E09 L=152)   202    106        0 ( 852/1024- 930/1126)  60.3 E10-E06 N00-N07
      5 1509 1016 (N05W31 L=192)   833    177        1 (1457/ 932-1557/1094)  56.3 W28-W34 N00-N09
      6 1790 1027 (N07W55 L=215)   778    221        3 (1747/ 973-1844/1073) -32.0 W50-W61 N04-N09
      7  771  938 (S05E15 L=145)   366     84        0 ( 725/ 914- 811/ 964)  38.7 E17-E12 S07-S03
      8  209  910 (S13E60 L=101)   516     54        0 ( 184/ 885- 231/ 938) -36.9 E62-E57 S13-S10
      9 1527  736 (S12W36 L=196)   494    145        2 (1465/ 713-1597/ 764)  18.4 W31-W40 S13-S10
  

    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