Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  952 1559 (N32W01 L=270)   445    108        0 ( 903/1526-1004/1582) -31.0 E02-W04 N29-N34
      2 1736 1493 (N21W59 L=328)  1104    160        3 (1691/1430-1782/1544) -69.4 W55-W62 N16-N24
      3 1306 1473 (N23W23 L=292)   209     55        0 (1271/1453-1335/1494)  11.3 W20-W25 N22-N24
      4 1742 1415 (N16W55 L=325)   467     84        2 (1709/1395-1787/1442) -20.6 W52-W60 N14-N17
      5  251 1321 (N24E54 L=215)   337     46        1 ( 238/1294- 267/1350) -45.0 E56-E52 N22-N26
      6 1201 1173 (N06W12 L=281)   415     84        0 (1157/1139-1238/1199) -45.0 W09-W13 N03-N07
      7 1747 1085 (S04W48 L=317)   319     35        1 (1722/1062-1770/1111)  33.7 W46-W50 S04-S02
      8 1807  893 (S16W54 L=323)  1090    163        1 (1765/ 813-1840/ 964)  77.5 W49-W57 S20-S11
      9  140  805 (S05E70 L=199)  1502    100        2 (  97/ 776- 189/ 838) -15.9 E79-E63 S07-S03
     10 1676  792 (S21W42 L=311)   272     33        0 (1659/ 768-1697/ 818)  56.3 W40-W43 S22-S19
     11 1503  717 (S24W29 L=298)   456    122        0 (1463/ 672-1547/ 776)  50.2 W25-W31 S26-S20
     12 1645  662 (S29W41 L=310)   308     72        0 (1611/ 644-1677/ 680)   9.5 W37-W43 S29-S27
     13  309  599 (S19E57 L=212)  2483    335        7 ( 243/ 524- 395/ 673)  14.7 E67-E47 S23-S14
     14  933  494 (S33E12 L=257)   407     71        1 ( 891/ 458- 963/ 522) -51.3 E14-E10 S35-S30
     15  736  468 (S32E27 L=242)   917    251        1 ( 644/ 408- 859/ 510)  11.3 E36-E16 S35-S30
     16 1179  418 (S41W05 L=274)   578    180        2 (1091/ 379-1259/ 444)   8.7 E02-W11 S43-S39
     17 1417  326 (S52W28 L=297)  1393    196        1 (1337/ 300-1499/ 354)  -3.2 W19-W38 S54-S48
     18 1199  281 (S53W06 L=275)   609     67        1 (1160/ 255-1235/ 318)   7.1 W01-W09 S55-S49
  

    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