Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1213 1786 (N47E15 L=107)   358     41        1 (1196/1761-1235/1811) -31.0 E18-E13 N44-N49
      2  643 1721 (N24E49 L= 74)   481     88        1 ( 628/1678- 661/1766) -36.9 E54-E45 N21-N27
      3  546 1677 (N19E53 L= 69)   628    129        1 ( 515/1628- 585/1720) -81.9 E55-E51 N15-N22
      4 1531 1688 (N53W16 L=139)   409     40        1 (1504/1671-1561/1706)  40.6 W13-W20 N50-N56
      5  484 1652 (N16E57 L= 65)   510     56        1 ( 463/1625- 506/1678)  90.0 E59-E54 N14-N17
      6  611 1510 (N10E39 L= 84)   457     76        1 ( 579/1475- 642/1544) -36.9 E41-E35 N08-N12
      7 1324 1523 (N31W03 L=126)   273     46        0 (1300/1497-1350/1540)  45.0 W01-W04 N28-N32
      8 1818 1213 (N29W49 L=172)   428    111        0 (1775/1174-1858/1248)  -5.7 W43-W54 N28-N29
      9 1179  438 (S35W31 L=153)   597    144        0 (1149/ 370-1211/ 506) -26.6 W25-W37 S37-S31
     10 1350  400 (S31W44 L=167)   293     57        0 (1325/ 376-1383/ 423)  90.0 W42-W45 S32-S28
  

    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