Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  867 1687 (N47W03 L=105)   257     29        1 ( 848/1668- 893/1704)   0.0 W01-W05 N45-N48
      2 1354 1629 (N31W37 L=138)   278     55        0 (1333/1600-1377/1657)  76.0 W35-W38 N29-N33
      3 1639 1540 (N20W55 L=157)   416     33        0 (1615/1518-1661/1566) -56.3 W53-W57 N18-N22
      4  560 1388 (N32E25 L= 77)   434    199        0 ( 504/1320- 633/1425)  24.0 E29-E18 N28-N33
      5 1527  997 (S09W30 L=131)   388     91        0 (1505/ 953-1545/1045)  68.2 W28-W31 S11-S06
      6  301  865 (N04E51 L= 51)   980    175        0 ( 230/ 827- 355/ 928)  12.1 E58-E44 N02-N07
      7 1076  842 (S10E00 L=101)   710    197        1 (1025/ 756-1130/ 919)  48.8 E04-W03 S14-S06
      8  483  746 (S06E39 L= 62)   537    113        1 ( 457/ 691- 529/ 782)  29.7 E42-E35 S08-S03
      9 1818  703 (S33W57 L=158)   944     91        1 (1785/ 659-1845/ 747) -38.7 W51-W63 S37-S29
     10 1717  675 (S33W44 L=146)   582    173        0 (1659/ 612-1762/ 740)  41.2 W39-W48 S36-S29
     11  499  616 (S13E42 L= 60)   765    110        2 ( 472/ 579- 538/ 663) -90.0 E44-E38 S16-S10
  

    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