Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  594 1759 (N58E31 L=160)   988    172        1 ( 524/1721- 675/1797)   6.0 E40-E20 N55-N60
      2 1104 1694 (N42W18 L=209)  1187    219        3 (1005/1669-1202/1725) -14.9 W10-W26 N38-N43
      3 1287 1619 (N33W29 L=220)   767    137        0 (1239/1579-1358/1643)   5.7 W24-W34 N30-N35
      4 1653 1533 (N22W55 L=246)  1032    141        1 (1585/1497-1722/1568) -26.6 W49-W60 N18-N24
      5 1639  793 (S22W38 L=229)   493    146        0 (1610/ 736-1663/ 868)  63.4 W35-W40 S25-S17
      6 1019  680 (S20E05 L=186)   747    245        1 ( 973/ 584-1045/ 784)  85.6 E08-E03 S26-S13
      7 1431  605 (S31W22 L=213)   429     98        1 (1380/ 583-1472/ 636)  15.9 W17-W24 S32-S29
      8  403  537 (S21E53 L=138)   445    103        0 ( 386/ 489- 424/ 586)  40.6 E56-E48 S23-S17
      9 1475  384 (S49W28 L=219)  1121    160        1 (1403/ 353-1548/ 419)   0.0 W20-W35 S52-S45
  

    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