Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1087 1820 (N48E27 L=359)   383     62        1 (1056/1805-1116/1837)  45.0 E28-E24 N45-N49
      2 1195 1624 (N36E08 L= 18)  1327    290        4 (1124/1500-1260/1737) -39.1 E17-W00 N28-N42
      3  402 1687 (N19E68 L=317)  2618    113        1 ( 359/1639- 452/1730)  69.4 E71-E64 N14-N22
      4  601 1565 (N15E42 L=343)   877    354        2 ( 556/1447- 666/1683) -39.1 E49-E33 N09-N22
      5  417 1502 (N07E54 L=332)   927    194        3 ( 353/1422- 461/1576)  90.0 E57-E50 N01-N12
      6  351 1528 (N07E62 L=324)   635    111        2 ( 312/1485- 385/1556)  51.3 E65-E59 N04-N09
      7 1788 1411 (N43W46 L= 71)   495    105        1 (1766/1386-1815/1429)  39.8 W42-W49 N40-N45
      8 1652 1260 (N27W33 L= 59)   442    129        0 (1623/1205-1671/1325) -60.3 W29-W34 N23-N30
      9 1178 1195 (N09W04 L= 30)   272     49        0 (1155/1165-1196/1230)  90.0 W03-W05 N06-N10
     10  799 1064 (S09E14 L= 12)   186     27        0 ( 784/1042- 815/1091) -45.0 E15-E12 S09-S07
     11 1823  901 (N12W56 L= 82)  2079    342        6 (1774/ 784-1869/1021) -30.5 W49-W66 N06-N17
     12 1440  985 (N04W25 L= 50)   306     66        1 (1407/ 966-1472/1009)  63.4 W22-W26 N02-N06
     13 1489  957 (N04W28 L= 54)   241     19        0 (1470/ 935-1508/ 978)  71.6 W27-W29 N02-N05
     14 1886  926 (N16W63 L= 88)   590    122        0 (1868/ 887-1912/ 951)  90.0 W60-W66 N12-N17
     15 1231  802 (S12W18 L= 43)   534    177        0 (1159/ 740-1321/ 837)  52.1 W15-W22 S16-S07
     16  986  675 (S26W08 L= 34)   675    164        1 ( 895/ 656-1073/ 698)  31.0 W02-W12 S28-S21
     17 1155  524 (S29W24 L= 50)   321     52        0 (1126/ 508-1192/ 548)  11.3 W21-W26 S30-S28
     18  928  477 (S40W12 L= 38)   583    140        2 ( 865/ 452-1009/ 499)  21.8 W07-W17 S42-S36
  

    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