Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1524 1635 (N43W45 L=131)  2068    319        4 (1461/1569-1568/1701)  33.7 W36-W52 N38-N48
      2 1082 1506 (N35W04 L= 90)   499    166        1 (1060/1435-1101/1572) -84.3 W02-W05 N29-N39
      3  741 1511 (N35E21 L= 65)   443     77        2 ( 716/1478- 765/1546) -51.3 E23-E19 N32-N37
      4 1688 1321 (N21W49 L=135)   382    102        0 (1651/1279-1722/1350) -45.0 W45-W51 N18-N22
      5  731 1307 (N21E20 L= 66)   244     74        0 ( 716/1266- 745/1342) -68.2 E20-E18 N18-N23
      6 1063 1190 (N14W02 L= 88)   285    155        0 (1050/1118-1086/1269)  77.5 W01-W04 N09-N18
      7 1680 1234 (N15W46 L=132)   757    186        1 (1624/1195-1743/1264)  10.3 W41-W52 N13-N17
      8 1280 1197 (N14W16 L=102)   543    121        1 (1235/1136-1312/1259) -76.0 W13-W18 N10-N18
      9 1719 1107 (N07W48 L=134)   806    263        4 (1626/1077-1800/1159) -17.4 W40-W56 N05-N10
     10  773 1127 (N10E16 L= 70)   288    105        1 ( 721/1111- 830/1147)  -8.1 E19-E12 N08-N11
     11  662  982 (N01E23 L= 63)   396    137        0 ( 607/ 952- 719/1030)  29.7 E26-E19 S01-N03
     12  291  827 (S10E53 L= 33)   456    120        1 ( 265/ 798- 317/ 869)  38.7 E55-E49 S11-S07
     13  221  786 (S13E61 L= 25)   535    155        0 ( 182/ 767- 257/ 814) -11.3 E66-E56 S14-S11
     14 1721  678 (S20W51 L=137)  1389    264        2 (1646/ 580-1768/ 812)  68.2 W46-W54 S26-S11
     15  633  586 (S25E27 L= 59)   479    120        1 ( 573/ 567- 685/ 607)   7.1 E31-E23 S26-S23
     16 1622  561 (S27W45 L=131)   315     38        1 (1599/ 546-1646/ 580)   0.0 W42-W47 S28-S25
     17 1199  369 (S41W14 L=100)  1707    391        6 (1124/ 303-1310/ 475)  51.8 W08-W21 S47-S32
     18  859  362 (S42E14 L= 72)   503     92        1 ( 820/ 339- 890/ 398)   9.5 E17-E11 S43-S38
  

    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