Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1320 1655 (N32W41 L=102)   552    155        0 (1261/1616-1380/1680) -11.3 W35-W45 N30-N33
      2 1389 1465 (N20W35 L= 96)   233     46        1 (1368/1443-1406/1490)  18.4 W33-W36 N18-N20
      3  778 1161 (N20E10 L= 51)  1624    457        4 ( 598/1095-1004/1233) -39.0 E20-W01 N11-N28
      4 1733  967 (S18W42 L=103)  1770    333        5 (1687/ 783-1758/1175)  57.0 W35-W49 S27-S06
      5 1817 1117 (S12W54 L=115)   661     69        1 (1790/1088-1841/1147) -76.0 W51-W56 S14-S10
      6 1314 1106 (N03W18 L= 79)   296     50        1 (1291/1075-1339/1129)  18.4 W16-W19 N01-N04
      7  273 1025 (N24E50 L= 12)   431     91        2 ( 250/ 980- 294/1072)  68.2 E52-E47 N21-N26
      8 1510  699 (S26W20 L= 81)  1816    451        3 (1386/ 581-1567/ 948)  46.6 W10-W29 S32-S13
      9  828  750 (S04E18 L= 43)   205     52        0 ( 805/ 715- 852/ 777) -78.7 E19-E17 S06-S01
     10  967  697 (S10E12 L= 49)   468    168        2 ( 882/ 665-1047/ 742) -11.3 E16-E06 S12-S07
     11 1054  610 (S18E10 L= 51)   343    148        1 (1027/ 542-1077/ 677)  49.4 E13-E06 S21-S14
  

    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