Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1293 1832 (N66W12 L= 74)  4422    387        5 (1136/1752-1449/1881)   1.3 E10-W38 N59-N71
      2 1354 1638 (N46W17 L= 79)  2267    391        4 (1176/1567-1479/1676)  -4.1 W00-W28 N42-N50
      3  792 1629 (N35E29 L= 34)   318     70        2 ( 759/1611- 829/1649)  36.9 E30-E26 N33-N37
      4  371 1625 (N28E71 L=352)   726     15        0 ( 350/1607- 392/1645)   8.1 E75-E66 N26-N29
      5  762 1546 (N29E27 L= 35)   390     98        0 ( 736/1495- 790/1587) -35.5 E30-E23 N26-N31
      6  593 1454 (N20E36 L= 26)   636    162        1 ( 556/1403- 659/1520) -28.6 E41-E30 N17-N24
      7  641 1283 (N10E28 L= 34)   510    218        1 ( 590/1185- 688/1345) -41.6 E33-E23 N05-N13
      8 1789 1198 (N22W55 L=117)  1259    203        1 (1718/1121-1839/1261) -28.3 W47-W60 N17-N25
      9  302  852 (S21E50 L= 13)   372     66        0 ( 284/ 818- 317/ 889) -53.1 E51-E48 S22-S18
     10  423  761 (S25E38 L= 24)   441    114        0 ( 406/ 715- 438/ 824)  74.1 E40-E37 S27-S20
     11  815  766 (S18E09 L= 53)   292    100        1 ( 775/ 729- 852/ 800) -26.6 E12-E06 S20-S16
     12  522  672 (S29E30 L= 32)   474    162        0 ( 461/ 634- 576/ 726)  35.0 E35-E25 S32-S24
     13 1141  662 (S20W13 L= 75)   225     80        0 (1103/ 637-1175/ 678) -21.8 W10-W15 S21-S19
     14 1586  643 (S15W45 L=107)   350     46        0 (1560/ 622-1612/ 659)   0.0 W42-W46 S16-S13
     15  917  613 (S27E01 L= 61)   636    193        0 ( 823/ 601-1006/ 628)   9.5 E07-W05 S28-S24
     16  641  578 (S34E21 L= 42)   565    152        0 ( 576/ 551- 711/ 607)  24.4 E26-E15 S36-S30
     17 1421  526 (S25W36 L= 98)   739    259        0 (1335/ 462-1531/ 587)  47.7 W31-W42 S30-S19
     18  855  171 (S69W05 L= 67)  1154     94        1 ( 826/ 129- 881/ 213) -40.9 W00-W15 S77-S63
     19  934  156 (S69W18 L= 81)   528     62        0 ( 914/ 139- 961/ 181) -18.4 W12-W24 S71-S64
  

    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