Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1601 1666 (N31W61 L=293)  1008    104        2 (1544/1642-1657/1689)  -3.4 W54-W71 N28-N32
      2 1162 1590 (N33W20 L=252)  1140    301        4 (1032/1548-1303/1624) -25.2 W11-W28 N28-N36
      3 1532 1558 (N25W46 L=278)   494    119        0 (1490/1503-1580/1597) -66.8 W43-W47 N20-N27
      4 1734 1500 (N18W62 L=294)   986    160        1 (1676/1459-1800/1551) -37.9 W58-W69 N14-N21
      5  624 1491 (N36E21 L=211)   681    145        2 ( 553/1446- 693/1517) -24.4 E26-E15 N31-N38
      6  437 1489 (N39E38 L=194)   572     90        1 ( 409/1453- 477/1511)  31.0 E40-E34 N36-N40
      7  494 1149 (N16E32 L=200)   642    119        1 ( 451/1098- 539/1197) -53.1 E34-E28 N11-N19
      8 1787 1061 (S09W52 L=284)   585    102        0 (1745/1020-1817/1098) -45.0 W48-W55 S11-S06
      9  184 1001 (N11E61 L=171)   386     47        0 ( 160/ 989- 212/1022) -18.4 E64-E57 N10-N12
     10 1846  787 (S27W60 L=292)   555     67        1 (1815/ 763-1874/ 814)   8.1 W55-W64 S27-S24
     11 1122  305 (S49E07 L=225)   717    220        0 (1047/ 267-1197/ 387)  23.2 E14-W00 S51-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