Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  853 1871 (N73W27 L= 73)   589     69        1 ( 833/1837- 871/1904)  19.2 W20-W43 N67-N76
      2 1439 1633 (N32W49 L= 94)   365    120        0 (1418/1583-1476/1665)  15.3 W43-W54 N29-N33
      3 1629 1423 (N15W51 L= 97)   308    102        0 (1584/1399-1667/1444)   0.0 W46-W55 N13-N15
      4  913 1330 (N27W00 L= 45)   661    255        1 ( 862/1236- 971/1422)  39.3 E05-W06 N22-N32
      5  530 1237 (N30E29 L= 17)   533    206        0 ( 494/1139- 550/1321) -76.0 E31-E27 N23-N35
      6 1704 1203 (N01W48 L= 94)   546     54        1 (1670/1178-1730/1234)   0.0 W45-W50 N00-N02
      7 1749 1134 (S04W50 L= 96)   422     71        1 (1728/1106-1778/1162) -53.1 W48-W52 S06-S02
      8 1857 1078 (S11W61 L=106)   496     75        1 (1837/1043-1877/1116)  26.6 W57-W64 S12-S09
      9 1608 1090 (S03W38 L= 83)   350     84        0 (1571/1060-1642/1114) -51.3 W35-W39 S05-N00
     10  659 1062 (N17E22 L= 24)   235     72        0 ( 633/1028- 689/1094)  26.6 E23-E19 N14-N18
     11  413 1011 (N18E40 L=  5)   277     76        1 ( 381/ 981- 439/1044) -59.0 E42-E38 N15-N20
     12 1603  996 (S08W35 L= 81)   251     21        0 (1582/ 979-1621/1017)  18.4 W33-W36 S09-S06
     13 1417  779 (S16W18 L= 64)   378    149        1 (1363/ 750-1487/ 822)   6.3 W14-W23 S16-S14
     14 1688  657 (S32W39 L= 84)   379     75        0 (1653/ 643-1731/ 676)  -7.1 W35-W43 S33-S30
     15 1324  585 (S26W09 L= 54)   426     97        1 (1298/ 552-1345/ 647) -90.0 W07-W10 S28-S21
     16  517  443 (S19E50 L=356)   361     64        1 ( 500/ 410- 533/ 475)  59.0 E52-E48 S21-S16
     17  707  377 (S27E39 L=  7)  1325    260        2 ( 607/ 328- 830/ 434)   7.4 E49-E26 S29-S24
  

    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