Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1255 1742 (N45E09 L= 14)  2621    430        3 (1137/1586-1412/1885) -24.0 E29-W07 N37-N55
      2  618 1693 (N22E48 L=335)  2463    365        8 ( 501/1604- 743/1771)  45.0 E55-E41 N15-N28
      3 1024 1680 (N32E21 L=  2)   532    191        1 ( 963/1599-1049/1745) -33.7 E28-E16 N26-N35
      4  672 1586 (N16E38 L=345)   319     53        3 ( 648/1565- 696/1618) -21.8 E40-E35 N14-N17
      5 1369 1382 (N23W10 L= 33)   220     91        0 (1344/1356-1396/1426) -21.8 W07-W12 N21-N25
      6 1770 1237 (N28W43 L= 67)   443     86        1 (1738/1210-1798/1283) -26.6 W39-W46 N26-N30
      7  197 1017 (S27E61 L=322)   564    123        1 ( 154/ 977- 241/1064)  26.6 E67-E55 S29-S23
      8  664  932 (S21E19 L=  5)   240     76        1 ( 619/ 911- 689/ 951)  36.9 E21-E16 S23-S19
      9 1724  649 (S06W57 L= 80)   432     37        2 (1702/ 628-1747/ 673)   0.0 W54-W60 S06-S04
     10 1283  416 (S33W39 L= 62)   304     68        1 (1262/ 380-1301/ 452) -53.1 W36-W41 S34-S30
     11 1311  190 (S44W67 L= 90)   625     75        1 (1277/ 158-1342/ 219)   3.4 W60-W77 S45-S42
  

    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