Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  729 1782 (N32E49 L=230)   535    125        1 ( 679/1724- 765/1829) -63.4 E53-E47 N26-N36
      2 1065 1705 (N37E21 L=258)  1949    415        1 ( 928/1558-1147/1772) -22.2 E31-E09 N28-N41
      3  533 1698 (N21E58 L=221)  1167    168        3 ( 474/1655- 589/1765)  90.0 E63-E54 N17-N26
      4  592 1676 (N21E51 L=229)   373     37        1 ( 568/1656- 613/1692)  33.7 E53-E48 N19-N22
      5  733 1618 (N21E37 L=242)   516     95        1 ( 715/1559- 750/1665) -54.5 E39-E34 N16-N23
      6  320 1270 (S09E53 L=226)   812    167        4 ( 252/1236- 354/1331)  40.6 E59-E49 S12-S05
      7  168 1139 (S20E67 L=212)  4402    415        7 (  97/ 961- 218/1295)  57.5 E83-E61 S31-S09
      8 1595 1239 (N24W29 L=308)   266     75        0 (1570/1216-1625/1277) -21.8 W26-W31 N22-N25
      9  273 1178 (S16E55 L=224)   579     69        2 ( 250/1143- 301/1214) -45.0 E57-E51 S17-S13
     10 1634  941 (N08W39 L=318)   499    204        0 (1611/ 864-1657/1048) -68.2 W36-W40 N03-N13
     11  313  766 (S40E46 L=233)   878    130        4 ( 261/ 724- 371/ 795)  23.2 E53-E39 S43-S36
     12 1331  738 (S12W26 L=305)   314    116        0 (1285/ 701-1379/ 765)  56.3 W24-W28 S15-S09
     13 1866  700 (N04W73 L=352)   679     49        2 (1849/ 678-1884/ 732) -71.6 W69-W79 N02-N05
  

    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