Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  645 1769 (N37E50 L=174)  1261    236        3 ( 574/1737- 711/1802)   0.0 E58-E43 N34-N38
      2  853 1645 (N32E28 L=197)  1610    255        3 ( 717/1585- 996/1699)   5.2 E38-E16 N29-N34
      3 1208 1331 (N20W05 L=229)   214     62        1 (1182/1302-1233/1370) -36.9 W02-W06 N18-N21
      4 1869 1064 (N19W60 L=285)   556     65        1 (1854/1031-1891/1097) -63.4 W57-W62 N17-N21
      5 1600  998 (N09W36 L=261)   816    206        1 (1563/ 886-1619/1096) -81.3 W34-W38 N01-N14
      6  500 1068 (S10E33 L=191)   373    111        0 ( 445/1053- 547/1085)   7.1 E37-E29 S11-S08
      7 1322  805 (S08W22 L=247)   289    141        0 (1303/ 743-1346/ 882) -54.5 W19-W25 S11-S04
      8  549  729 (S30E25 L=200)   478    176        2 ( 469/ 714- 629/ 749)   8.7 E31-E18 S31-S27
      9  767  728 (S25E10 L=215)   253     70        1 ( 728/ 710- 808/ 749)   0.0 E12-E06 S25-S23
     10  948  450 (S39W10 L=234)   304     66        0 ( 916/ 428- 967/ 478)  78.7 W07-W11 S41-S36
  

    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