Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  863 1744 (N59W05 L= 12)   306     38        1 ( 842/1723- 887/1763) -45.0 W02-W07 N56-N60
      2  751 1712 (N59E09 L=359)  1523    223        3 ( 681/1641- 831/1753) -16.4 E17-W00 N52-N61
      3 1221 1709 (N46W35 L= 43)  1106    179        0 (1134/1678-1306/1731) -10.0 W26-W44 N43-N48
      4 1452 1682 (N38W53 L= 61)   805    110        1 (1404/1657-1514/1698) -14.0 W47-W59 N35-N39
      5 1539 1512 (N25W48 L= 55)   318    136        0 (1491/1477-1570/1535) -31.0 W44-W50 N23-N27
      6  620 1338 (N32E22 L=345)   321     69        0 ( 590/1311- 659/1357)  11.3 E24-E19 N30-N33
      7 1344 1258 (N15W24 L= 32)  1248    287        2 (1270/1141-1382/1356)  81.9 W20-W26 N07-N21
      8  773 1241 (N24E12 L=355)   267     97        1 ( 755/1190- 793/1273)  59.0 E13-E10 N20-N25
      9  602 1082 (N17E26 L=342)   509    155        1 ( 534/1049- 665/1114)  12.5 E30-E21 N15-N17
     10 1769  716 (S28W50 L= 58)  2643    212        5 (1721/ 611-1834/ 799) -61.9 W44-W58 S36-S21
     11 1027  261 (S46E18 L=349)   329     40        0 (1005/ 242-1046/ 281)  21.8 E21-E16 S47-S43
     12  942  221 (S48E28 L=339)   529     28        0 ( 919/ 202- 970/ 244) -56.3 E30-E25 S49-S44
  

    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