Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1444 1821 (N38W73 L=150)   654     67        0 (1427/1784-1457/1852)   5.7 W65-W85 N36-N38
      2 1462 1730 (N33W61 L=137)  1118    104        2 (1415/1693-1505/1764) -90.0 W57-W65 N29-N35
      3  447 1457 (N48E33 L= 43)   358     74        1 ( 413/1434- 475/1481)  15.9 E36-E29 N46-N49
      4 1831 1115 (S12W57 L=133)   511     69        1 (1806/1089-1857/1151) -45.0 W54-W59 S13-S09
      5 1745  974 (S17W44 L=121)  1329    215        2 (1720/ 873-1764/1093)  65.2 W40-W47 S22-S09
      6  599 1018 (N17E25 L= 51)   266     66        1 ( 573/ 994- 638/1041) -45.0 E27-E23 N15-N18
      7  340  860 (N13E49 L= 28)   602     88        2 ( 311/ 812- 362/ 894) -68.2 E50-E46 N10-N15
      8  176  734 (N08E72 L=  4)  1351    132        1 ( 125/ 693- 215/ 797) -26.6 E81-E66 N05-N11
      9  396  746 (N05E47 L= 29)   289     50        2 ( 375/ 721- 418/ 766)  45.0 E48-E45 N04-N06
     10 1412  306 (S53W05 L= 81)   391     66        0 (1382/ 290-1445/ 318)  -8.1 W01-W08 S54-S50
     11  892  270 (S36E35 L= 42)   574    109        0 ( 848/ 235- 935/ 314)  18.4 E40-E28 S37-S33
  

    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