Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  546 1125 (N23E28 L=318)   622    193        3 ( 460/1094- 638/1153) -26.6 E33-E21 N20-N26
      2 1547  986 (S08W30 L= 16)   458    125        0 (1505/ 939-1603/1011)   7.1 W26-W34 S09-S06
      3  851  868 (N02E14 L=333)   810    239        3 ( 824/ 749- 884/ 968)  71.6 E17-E11 S04-N08
      4 1687  810 (S23W37 L= 23)   888    245        0 (1639/ 693-1713/ 937)  57.3 W32-W41 S29-S15
      5 1388  835 (S12W17 L=  3)   712    173        1 (1329/ 753-1418/ 920)  48.8 W11-W19 S16-S07
      6 1239  712 (S16W05 L=351)   347    110        1 (1191/ 681-1294/ 737)   8.1 W01-W08 S16-S14
      7 1402  623 (S26W13 L=360)   764    215        0 (1309/ 550-1482/ 719)  21.8 W05-W20 S28-S22
      8 1718  662 (S34W39 L= 25)   594    102        1 (1673/ 627-1759/ 699)   6.3 W34-W43 S36-S32
      9 1011  486 (S25E15 L=332)   291     68        1 ( 984/ 459-1046/ 516)  18.4 E17-E11 S25-S23
     10  962  279 (S39E27 L=319)   544    183        1 ( 937/ 223- 994/ 339)  39.3 E32-E21 S43-S34
  

    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