Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  751 1760 (N53E32 L=159)   434     57        1 ( 713/1739- 781/1776)   8.1 E36-E29 N51-N54
      2 1585 1678 (N48W57 L=248)   714     60        1 (1551/1654-1620/1697)   7.1 W50-W66 N45-N49
      3 1203 1522 (N35W12 L=203)   672    178        0 (1131/1469-1255/1566) -24.0 W06-W15 N31-N38
      4  402 1491 (N30E52 L=139)  1951    329        4 ( 269/1442- 536/1542)  -4.4 E76-E37 N26-N32
      5 1249 1328 (N22W14 L=205)   477    108        1 (1225/1269-1273/1379)  74.1 W12-W15 N18-N25
      6  274 1164 (N08E54 L=137)   828    177        1 ( 252/1097- 304/1256) -55.0 E58-E50 N04-N14
      7 1661 1052 (N05W43 L=234)  1527    327        2 (1578/ 927-1705/1181) -61.9 W36-W46 S02-N13
      8  887 1096 (N07E09 L=183)   220    122        0 ( 840/1061- 914/1133) -38.7 E11-E06 N04-N09
      9  604 1088 (N06E27 L=164)   265    115        0 ( 564/1054- 665/1115)  29.7 E29-E22 N03-N07
     10 1874 1061 (N06W65 L=256)   506     68        2 (1855/1028-1890/1093) -63.4 W62-W67 N03-N07
     11 1841  834 (S08W62 L=253)   546     45        0 (1823/ 805-1857/ 866) -63.4 W59-W64 S10-S06
     12 1562  761 (S13W37 L=228)   346     46        1 (1541/ 736-1586/ 792) -76.0 W35-W38 S14-S10
     13  254  699 (S21E59 L=132)   890    144        1 ( 203/ 652- 296/ 742)  22.4 E70-E53 S24-S17
     14 1760  606 (S23W60 L=251)  2942    281        5 (1670/ 525-1842/ 734)  54.5 W52-W66 S28-S14
     15  577  644 (S22E30 L=161)   371    100        0 ( 561/ 587- 596/ 688)  74.1 E31-E27 S26-S19
     16  661  597 (S25E24 L=167)   549    152        0 ( 588/ 566- 711/ 642)  29.1 E29-E20 S27-S21
     17  852  583 (S25E11 L=181)   927     80        0 ( 787/ 544- 932/ 636) -24.0 E14-E05 S28-S21
     18 1274  218 (S55W32 L=223)  2743    233        7 (1171/ 186-1393/ 257)  11.8 W19-W48 S59-S50
  

    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