Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1671 1625 (N30W63 L=314)  1034    104        3 (1639/1601-1720/1666)   3.8 W57-W72 N27-N32
      2 1333 1557 (N29W28 L=279)   902    194        0 (1263/1479-1388/1616) -68.2 W24-W30 N23-N33
      3  942 1455 (N27E00 L=250)   393     87        0 ( 907/1430- 991/1486) -36.9 E02-W02 N24-N28
      4 1675 1446 (N18W52 L=302)   324     32        0 (1649/1426-1697/1462) -26.6 W49-W53 N16-N18
      5  493  744 (S11E38 L=213)   249    124        0 ( 429/ 732- 545/ 766)   6.3 E42-E33 S12-S10
      6 1770  621 (S34W57 L=307)   416     34        0 (1753/ 599-1789/ 647) -26.6 W53-W60 S35-S31
      7 1525  578 (S34W31 L=282)   652    125        1 (1497/ 518-1542/ 637) -77.5 W28-W33 S38-S29
      8  293  508 (S24E68 L=183)   678     92        2 ( 251/ 485- 327/ 536)   6.0 E79-E60 S24-S22
      9  947  475 (S34E13 L=238)   431    122        1 ( 920/ 425- 983/ 529) -90.0 E13-E10 S37-S30
     10  970  255 (S52E19 L=232)   405     74        0 ( 936/ 239-1012/ 269)   7.1 E22-E14 S53-S50
     11  802  174 (S57E45 L=205)   761     65        2 ( 775/ 141- 827/ 209)  26.6 E55-E39 S60-S53
  

    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