Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1507 1759 (N45W66 L=217)   639     63        1 (1473/1739-1534/1779)   0.0 W58-W74 N43-N46
      2  661 1699 (N56E27 L=125)   501     71        1 ( 626/1668- 694/1727)  23.2 E30-E23 N53-N58
      3  427 1671 (N55E62 L= 90)  1046     98        1 ( 372/1649- 479/1691)  -1.8 E82-E50 N53-N56
      4  779 1385 (N31E13 L=139)   247     70        0 ( 754/1358- 798/1418)  45.0 E14-E11 N28-N32
      5  612 1357 (N30E25 L=126)  1111    216        2 ( 510/1328- 718/1378)  -3.6 E33-E17 N28-N32
      6  517 1102 (N15E32 L=119)   717    250        0 ( 458/1004- 559/1192)  65.6 E37-E28 N09-N20
      7  961  993 (N05E04 L=147)  1080    399        0 ( 877/ 819-1004/1128) -78.1 E07-E01 S06-N13
      8 1653  960 (S06W40 L=192)   644    109        2 (1616/ 910-1694/1004) -50.2 W37-W43 S09-S03
      9  834  769 (S08E14 L=137)   514    231        1 ( 735/ 735- 941/ 787) -14.0 E20-E08 S10-S06
     10 1835  726 (S25W61 L=212)   947     94        1 (1817/ 679-1856/ 775) -60.3 W58-W63 S28-S21
     11 1765  659 (S28W53 L=205)   608     99        2 (1744/ 612-1790/ 706)  80.5 W51-W55 S30-S24
     12 1311  625 (S22W14 L=165)   914    174        3 (1254/ 564-1360/ 696)  53.1 W10-W17 S26-S18
     13  744  627 (S16E22 L=129)   290     47        1 ( 723/ 603- 765/ 658)  36.9 E24-E20 S17-S13
     14 1463  489 (S35W26 L=177)  3139    485        4 (1305/ 406-1608/ 633)  20.0 W13-W37 S40-S26
     15  844  473 (S27E19 L=132)   411    106        0 ( 812/ 451- 891/ 501) -21.8 E21-E15 S29-S25
     16  627  377 (S33E39 L=112)  1095    183        3 ( 575/ 298- 683/ 434)  29.1 E50-E32 S39-S28
     17 1471  287 (S54W34 L=186)  1086    101        2 (1414/ 267-1523/ 312) -12.5 W25-W43 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