Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1043 1494 (N28W09 L=116)  1016    273        0 ( 951/1389-1139/1583) -65.0 W05-W13 N20-N35
      2  700 1437 (N31E14 L= 93)  1259    279        2 ( 647/1320- 791/1540)  47.5 E18-E06 N23-N36
      3 1433 1485 (N21W35 L=141)  1782    430        3 (1232/1437-1567/1528) -14.6 W22-W45 N16-N26
      4  412 1346 (N30E37 L= 70)  1387    277        5 ( 300/1298- 558/1382) -19.2 E48-E25 N24-N33
      5 1039 1226 (N12W04 L=111)   788    302        1 ( 921/1151-1103/1328)  30.3 E03-W09 N08-N17
      6 1786 1064 (S09W51 L=158)  2247    406        5 (1706/ 910-1870/1197) -83.7 W44-W62 S18-N00
      7 1883 1016 (S14W61 L=168)  1372    148        3 (1849/ 954-1919/1089) -90.0 W57-W66 S17-S08
      8  657  961 (N03E23 L= 84)   503    156        0 ( 578/ 946- 727/ 977) -12.5 E27-E18 N01-N03
      9 1671  926 (S15W40 L=147)   233     37        0 (1653/ 905-1686/ 947)  45.0 W38-W41 S16-S14
     10  701  887 (S02E21 L= 86)   227     54        0 ( 680/ 857- 728/ 915)  45.0 E22-E18 S03-N00
     11  374  766 (S04E46 L= 61)  1870    419        1 ( 244/ 611- 468/ 862) -49.4 E56-E42 S14-N01
     12 1497  791 (S20W26 L=132)  1148    305        0 (1379/ 719-1627/ 855)  12.5 W17-W35 S23-S18
     13  620  760 (S08E29 L= 78)   562    184        1 ( 554/ 726- 672/ 805) -23.2 E32-E24 S10-S05
     14  307  629 (S11E57 L= 50)  1204    270        3 ( 272/ 530- 345/ 747)  74.1 E60-E53 S17-S03
     15 1641  721 (S27W37 L=144)   452    118        0 (1591/ 705-1710/ 741)  -5.2 W32-W43 S28-S25
     16  489  435 (S26E50 L= 56)  1005    207        0 ( 425/ 377- 568/ 505)  15.9 E61-E40 S28-S22
     17 1396  300 (S55W18 L=124)  1053    142        0 (1326/ 274-1471/ 333)   3.8 W09-W25 S56-S52
  

    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