Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1213 1595 (N39W01 L=195)   305     52        1 (1191/1576-1241/1626) -45.0 W00-W03 N37-N41
      2  286 1425 (N10E64 L=130)  1101    157        6 ( 239/1360- 353/1467)  58.0 E72-E57 N05-N13
      3 1538 1360 (N29W29 L=223)   410     69        3 (1505/1324-1572/1386)   0.0 W26-W31 N26-N31
      4  122 1041 (S16E76 L=118)  1616    123        3 ( 107/ 993- 146/1113) -74.1 E79-E69 S18-S11
      5 1567 1031 (N09W34 L=228)   602    143        1 (1511/1010-1643/1068)  16.7 W30-W40 N07-N10
      6  559  849 (S20E27 L=167)  2002    498        4 ( 389/ 723- 739/1018)  47.4 E40-E16 S31-S06
      7 1708  650 (S10W56 L=250)   369     52        1 (1682/ 624-1724/ 670) -18.4 W52-W58 S11-S08
      8 1444  621 (S17W35 L=229)   303     32        0 (1421/ 605-1472/ 643) -14.0 W33-W37 S18-S15
      9 1791  579 (S13W72 L=266)   927     77        0 (1756/ 563-1835/ 596)  10.6 W66-W83 S14-S11
     10 1309  500 (S27W31 L=224)   240     91        0 (1262/ 488-1346/ 513)  18.4 W27-W33 S28-S25
     11 1520  457 (S26W50 L=243)   350    104        0 (1486/ 430-1549/ 505)  90.0 W48-W51 S28-S22
     12 1042  232 (S55W27 L=221)  1750    329        3 ( 896/ 209-1195/ 255)  22.9 W13-W39 S60-S48
  

    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