Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  607 1625 (N57E16 L=180)   311     40        1 ( 589/1601- 630/1649) -63.4 E18-E13 N54-N58
      2  317 1234 (N35E47 L=148)  1895    423        5 ( 183/1148- 493/1296)  -3.3 E65-E30 N31-N38
      3 1598 1164 (S00W38 L=234)   391    117        0 (1550/1143-1656/1177)  -7.1 W34-W42 S02-N01
      4 1756  913 (S20W45 L=241)   712    166        1 (1700/ 848-1814/ 966)   8.7 W38-W51 S21-S17
      5 1220  902 (S05W08 L=204)   363     69        1 (1196/ 860-1243/ 937) -78.7 W06-W09 S07-S02
      6 1769  780 (S29W45 L=241)   829    180        0 (1732/ 712-1790/ 867)  53.1 W41-W47 S32-S23
      7 1606  771 (S23W30 L=226)   260     63        0 (1577/ 753-1634/ 789) -11.3 W27-W32 S24-S21
      8 1478  658 (S26W19 L=215)   436    101        1 (1439/ 628-1508/ 700)  18.4 W15-W21 S27-S24
      9 1115  556 (S22E07 L=189)   753    273        0 (1005/ 493-1192/ 653)  18.4 E15-W00 S24-S18
     10 1403  400 (S43W08 L=204)  1575    417        2 (1231/ 290-1528/ 515)  15.9 E08-W20 S46-S38
     11  795  341 (S29E36 L=160)   414    124        0 ( 739/ 321- 844/ 365)  -6.3 E39-E30 S30-S26
  

    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