Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1373 1749 (N48W48 L=201)   362     55        0 (1345/1722-1396/1778)  15.3 W43-W54 N46-N49
      2 1501 1755 (N45W64 L=217)   560     59        0 (1468/1736-1527/1775)   3.8 W57-W72 N44-N46
      3  652 1696 (N56E28 L=125)   476     68        0 ( 627/1664- 680/1722)  53.1 E30-E25 N52-N58
      4  423 1671 (N55E63 L= 90)   954     98        1 ( 372/1649- 474/1690)  -3.6 E82-E50 N53-N56
      5  761 1380 (N30E14 L=139)   218     64        0 ( 737/1356- 778/1411)  45.0 E15-E12 N28-N32
      6  599 1357 (N30E26 L=126)  1041    212        2 ( 500/1326- 701/1376)  -3.6 E34-E18 N28-N31
      7  509 1112 (N16E33 L=120)   564    196        0 ( 470/1036- 544/1190)  71.6 E35-E30 N11-N20
      8  941  993 (N05E05 L=147)   947    379        0 ( 859/ 834- 981/1123) -74.5 E09-E02 S05-N13
      9 1636  960 (S06W39 L=192)   588    103        3 (1604/ 920-1674/1006) -56.3 W36-W41 S08-S02
     10  788  766 (S07E17 L=135)   558    245        0 ( 692/ 737- 903/ 783)  -4.4 E23-E10 S09-S06
     11 1828  724 (S25W60 L=212)   964     94        0 (1807/ 677-1851/ 769)  90.0 W56-W62 S27-S21
     12 1765  671 (S27W53 L=205)   632     91        1 (1742/ 624-1789/ 715)  63.4 W50-W54 S30-S24
     13 1294  621 (S22W13 L=165)   953    178        3 (1239/ 558-1346/ 694)  41.2 W08-W16 S26-S18
     14  733  624 (S16E23 L=129)   353     52        1 ( 708/ 599- 763/ 656)  36.9 E25-E21 S17-S13
     15 1444  483 (S35W24 L=177)  2994    496        3 (1290/ 402-1592/ 630)  14.6 W12-W35 S40-S26
     16  831  470 (S27E20 L=132)   458    109        0 ( 794/ 447- 880/ 500) -31.0 E22-E16 S29-S24
     17  627  398 (S31E38 L=114)   615    108        1 ( 589/ 371- 669/ 432)  18.4 E42-E33 S33-S28
     18 1472  281 (S55W35 L=187)   835     88        2 (1416/ 264-1515/ 303)  -7.6 W26-W42 S57-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