Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1298 1832 (N66W13 L= 74)  4070    380        5 (1145/1751-1464/1877)   2.4 E09-W41 N59-N70
      2 1364 1638 (N47W18 L= 79)  2331    399        3 (1177/1569-1490/1675)  -3.9 W00-W29 N43-N50
      3  776 1542 (N29E26 L= 35)   391    116        0 ( 749/1485- 810/1587) -40.6 E29-E21 N25-N31
      4  607 1452 (N21E35 L= 26)   578    161        0 ( 571/1405- 661/1517) -28.6 E40-E29 N18-N24
      5  663 1265 (N10E26 L= 35)   562    195        1 ( 607/1182- 704/1336) -41.6 E31-E22 N05-N13
      6 1798 1197 (N22W56 L=117)  1248    200        0 (1726/1118-1844/1256) -30.3 W48-W61 N17-N24
      7 1857  986 (N09W63 L=124)   418     41        2 (1831/ 970-1881/1005)   0.0 W59-W67 N08-N10
      8  315  846 (S21E48 L= 13)   368     69        1 ( 295/ 815- 337/ 885) -63.4 E50-E46 S22-S18
      9  823  770 (S18E09 L= 52)   327    120        1 ( 782/ 733- 866/ 821) -38.7 E11-E05 S19-S15
     10  435  757 (S25E37 L= 24)   460    114        0 ( 416/ 707- 449/ 816)  82.9 E39-E36 S28-S20
     11  538  669 (S29E29 L= 32)   477    167        0 ( 472/ 628- 591/ 726)  38.7 E34-E24 S32-S24
     12 1157  658 (S20W14 L= 75)   240     77        0 (1120/ 633-1192/ 678) -21.8 W11-W16 S21-S19
     13 1601  639 (S15W46 L=107)   371     51        1 (1571/ 621-1627/ 655) -11.3 W43-W48 S16-S14
     14  942  611 (S27W01 L= 62)   683    238        0 ( 835/ 598-1037/ 651)  22.6 E06-W06 S28-S22
     15  654  575 (S34E20 L= 42)   589    168        0 ( 577/ 547- 725/ 603)  18.4 E26-E14 S36-S30
     16 1487  539 (S23W41 L=102)   967    325        0 (1370/ 486-1571/ 585)  34.7 W33-W50 S28-S19
     17  859  172 (S69W05 L= 66)  1119    100        1 ( 828/ 130- 886/ 212) -36.9 W00-W16 S75-S63
  

    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