Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  972 1914 (N64W45 L= 37)   478     47        1 ( 950/1885- 987/1936)  19.7 W38-W52 N60-N66
      2  810 1836 (N64W16 L=  8)  2142    351        3 ( 685/1740- 928/1913)   7.8 E03-W41 N58-N67
      3  906 1612 (N42W10 L=  2)  1419    326        0 ( 758/1587-1064/1646) -21.8 E01-W19 N35-N45
      4 1677 1540 (N15W60 L= 52)   549     69        0 (1639/1519-1710/1564)   0.0 W55-W64 N14-N16
      5 1190 1499 (N26W23 L= 15)   497     82        1 (1138/1467-1221/1534)   8.1 W19-W26 N24-N28
      6 1818 1429 (N05W68 L= 60)  1925    146        7 (1758/1369-1860/1475) -53.1 W61-W75 N01-N09
      7  370 1414 (N43E40 L=312)   478     67        1 ( 343/1390- 407/1433)   0.0 E42-E35 N41-N44
      8 1496 1303 (N08W35 L= 27)   676    143        0 (1456/1236-1532/1375) -90.0 W33-W36 N03-N12
      9 1780 1226 (S04W54 L= 46)   301     21        1 (1762/1206-1798/1245) -90.0 W52-W55 S05-S02
     10 1161 1084 (N04W09 L=  1)   476    136        0 (1121/1028-1205/1138) -81.9 W07-W10 N00-N07
     11  146 1030 (N23E66 L=286)  5987    245        8 (  80/ 945- 216/1136) -63.4 E80-E56 N17-N29
     12  262 1057 (N23E51 L=301)   324     39        2 ( 240/1032- 291/1078) -21.8 E52-E47 N21-N23
     13  377  917 (N12E43 L=309)   737    174        1 ( 317/ 867- 436/ 986) -56.3 E46-E38 N08-N17
     14 1170  820 (S11W03 L=355)   677    246        0 (1135/ 721-1204/ 951)  59.7 W00-W07 S15-S03
     15  152  706 (N04E75 L=277)  3599    187        5 ( 114/ 643- 193/ 795) -60.9 E87-E69 N00-N09
     16 1741  586 (S42W43 L= 35)   686     99        1 (1691/ 568-1788/ 611) -19.7 W36-W50 S45-S38
     17  548  518 (S14E43 L=309)   746    158        2 ( 518/ 450- 583/ 594)  32.5 E48-E37 S17-S09
     18 1568  508 (S42W24 L= 16)   491     86        0 (1531/ 477-1611/ 540)   7.1 W20-W28 S43-S40
  

    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