Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1161 1641 (N36W29 L=103)   418     82        0 (1127/1620-1206/1660) -18.4 W26-W32 N33-N38
      2 1683 1599 (N17W65 L=138)   707    174        0 (1633/1582-1737/1620) -26.6 W60-W70 N14-N19
      3 1435 1404 (N14W35 L=108)   668    198        3 (1339/1347-1507/1442) -52.1 W30-W38 N09-N18
      4  617 1104 (N19E22 L= 51)   383     78        3 ( 573/1077- 659/1122) -38.7 E24-E19 N17-N21
      5 1899 1077 (S18W62 L=135)  1400    168        3 (1856/1035-1933/1121) -53.1 W57-W65 S21-S13
      6 1521  819 (S20W24 L= 97)   387    133        0 (1493/ 785-1548/ 858)  26.6 W21-W25 S22-S17
      7 1627  635 (S35W29 L=102)  1956    283        2 (1535/ 517-1682/ 767)  40.6 W20-W34 S41-S27
      8 1092  527 (S25E10 L= 63)  1100    312        4 (1052/ 397-1142/ 641)  70.6 E13-E04 S35-S18
      9  416  587 (S05E51 L= 22)   355     71        2 ( 400/ 546- 435/ 625)  53.1 E52-E48 S07-S03
     10  769  415 (S24E34 L= 40)   962    232        1 ( 710/ 344- 798/ 502)  39.3 E40-E29 S27-S18
     11  497  413 (S17E53 L= 20)  1828    242        6 ( 448/ 335- 549/ 478)  74.7 E56-E48 S23-S12
  

    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