Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1535 1532 (N25W49 L=199)  1666    261        1 (1432/1443-1614/1606) -70.3 W45-W52 N18-N32
      2  408 1435 (N44E43 L=107)   470     73        1 ( 384/1409- 449/1458)   0.0 E46-E37 N41-N45
      3 1171 1303 (N20W15 L=164)   257    130        1 (1130/1250-1222/1343) -66.8 W13-W16 N16-N23
      4 1816 1205 (S01W60 L=210)   420     26        0 (1799/1180-1835/1224)  11.3 W57-W62 S01-N00
      5 1589 1003 (S06W35 L=184)   463     97        2 (1544/ 975-1637/1037) -39.8 W32-W38 S08-S03
      6 1664  917 (S13W39 L=189)   393    100        1 (1639/ 872-1680/ 966)  51.3 W36-W41 S15-S10
      7  286  703 (S01E59 L= 90)   904     89        2 ( 253/ 657- 315/ 741) -68.2 E61-E55 S04-N01
      8 1699  684 (S29W41 L=191)   967    122        2 (1650/ 636-1739/ 736)  39.8 W36-W45 S32-S26
      9  276  551 (S11E69 L= 80)   791     74        4 ( 250/ 517- 301/ 578) -63.4 E72-E66 S13-S09
     10 1217  523 (S28W02 L=152)   247     57        0 (1198/ 491-1236/ 552)  45.0 W00-W03 S30-S26
     11  464  508 (S16E50 L= 99)   400     70        0 ( 452/ 472- 476/ 546)  45.0 E52-E47 S18-S13
  

    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