Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  656 1738 (N26E49 L= 31)   543     75        2 ( 640/1705- 674/1787) -33.7 E54-E45 N23-N29
      2 1052 1648 (N30E17 L= 62)   527    167        0 ( 992/1580-1098/1693) -20.0 E23-E12 N27-N32
      3  717 1637 (N21E38 L= 42)   311     80        1 ( 670/1625- 750/1659)  36.9 E40-E35 N18-N22
      4  772 1576 (N18E31 L= 48)   297     72        1 ( 748/1551- 802/1610) -18.4 E34-E28 N16-N19
      5 1657 1057 (N13W37 L=116)   241      0       -1 (1643/1006-1680/1135) -45.0 W34-W39 N10-N16
      6  160 1033 (S26E66 L= 13)   742     83        2 ( 114/1015- 199/1056)   0.0 E75-E60 S26-S24
      7 1729  657 (S05W57 L=137)   636     60        2 (1704/ 624-1756/ 694) -53.1 W54-W60 S07-S03
      8 1414  435 (S27W46 L=125)   908    279        0 (1279/ 395-1509/ 501)  57.3 W39-W50 S34-S20
      9 1640  472 (S18W61 L=141)   449     65        0 (1612/ 449-1676/ 488)  45.0 W59-W64 S19-S15
     10 1620  380 (S23W69 L=148)  2262    156        5 (1576/ 327-1683/ 428) -82.9 W63-W77 S25-S17
     11 1385  303 (S36W56 L=135)   391     68        0 (1357/ 274-1415/ 333) -68.2 W54-W57 S37-S32
  

    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