Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  215 1279 (N25E57 L=291)   427     53        0 ( 177/1267- 244/1296)  -5.7 E63-E53 N23-N25
      2 1111 1245 (N11W08 L=356)   176     73        1 (1078/1226-1147/1265) -18.4 W06-W09 N10-N12
      3  515 1191 (N15E30 L=318)   500    124        2 ( 464/1155- 573/1235)  14.0 E33-E25 N13-N17
      4 1667 1097 (S04W42 L= 30)  1648    364        5 (1557/ 943-1736/1202) -58.0 W34-W46 S13-N03
      5 1935 1067 (S08W70 L= 58)   578     64        1 (1901/1040-1955/1097)  12.5 W64-W74 S09-S06
      6  134  942 (N05E67 L=281)  1453    162        3 (  63/ 891- 196/ 985) -18.4 E83-E59 N01-N08
      7 1407  944 (S10W22 L= 10)   231     59        2 (1383/ 916-1430/ 970)  33.7 W20-W23 S11-S08
      8 1895  899 (S18W64 L= 52)   561    122        0 (1874/ 839-1919/ 951)  66.8 W61-W67 S21-S14
      9 1741  818 (S21W47 L= 35)  1505    257        1 (1703/ 710-1772/ 947)  70.3 W44-W49 S27-S13
     10  622  671 (S16E30 L=318)   697    150        1 ( 604/ 599- 644/ 759)  66.0 E32-E27 S20-S11
     11  975  588 (S26E09 L=339)  2697    671        6 ( 677/ 476-1258/ 721)   8.3 E31-W10 S30-S21
     12 1637  575 (S36W41 L= 29)  2808    440        8 (1447/ 466-1761/ 691)  18.4 W24-W51 S41-S28
     13  407  608 (S17E48 L=300)   677    166        1 ( 376/ 523- 436/ 668)  78.7 E49-E45 S23-S13
     14  373  483 (S25E58 L=290)  1469    162        2 ( 295/ 452- 436/ 534) -26.6 E63-E51 S26-S20
     15  487  368 (S34E55 L=293)  1134    212        0 ( 449/ 290- 528/ 444)  20.6 E69-E45 S38-S29
     16 1238  411 (S42W07 L=355)   275     80        1 (1204/ 393-1279/ 431) -18.4 W04-W10 S43-S39
     17  702  204 (S50E50 L=298)  1200    115        2 ( 654/ 178- 755/ 236)   0.0 E57-E42 S52-S47
  

    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