Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1150 1880 (N54E30 L= 42)   518     40        1 (1130/1858-1178/1900)  90.0 E33-E27 N51-N57
      2  637 1791 (N29E57 L= 15)   421     59        1 ( 614/1768- 662/1812) -18.4 E61-E52 N27-N30
      3  585 1654 (N18E49 L= 23)   913    166        2 ( 547/1579- 652/1710) -20.6 E55-E39 N14-N21
      4 1221 1231 (N10W05 L= 77)   254     94        1 (1203/1184-1244/1271) -78.7 W04-W06 N07-N12
      5  323 1178 (S15E50 L= 22)   408    128        2 ( 304/1114- 337/1229) -66.8 E52-E47 S19-S12
      6  172 1111 (S22E66 L=  6)   656     65        0 ( 152/1076- 191/1142)  63.4 E68-E63 S23-S19
      7  201 1041 (S25E60 L= 11)  1822    189        4 ( 159/ 966- 238/1127)  55.0 E66-E56 S30-S20
      8  373 1012 (S24E42 L= 30)   983    344        0 ( 247/ 924- 508/1087)  35.5 E53-E32 S31-S16
      9 1074  996 (S07W03 L= 75)   305    114        1 (1044/ 957-1120/1024)  45.0 W01-W06 S09-S05
  

    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