Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1448 1746 (N35W52 L=179)  1533    229        4 (1380/1704-1520/1778)  -4.4 W44-W58 N32-N37
      2  331 1376 (N38E43 L= 84)   529    133        1 ( 279/1354- 386/1397)   0.0 E47-E37 N37-N39
      3 1696 1152 (S05W45 L=172)  1802    398        3 (1635/ 958-1762/1303)  38.3 W36-W55 S14-N01
      4  117  949 (N15E69 L= 58)  1263    211        6 (  82/ 915- 153/ 984) -32.0 E73-E64 N12-N17
      5 1686  823 (S24W38 L=165)   453    119        2 (1630/ 801-1727/ 849)   0.0 W33-W42 S25-S22
      6 1266  515 (S33W04 L=130)   952    175        2 (1206/ 449-1298/ 596)  45.0 E01-W07 S36-S28
      7  827  445 (S28E27 L=100)   748    178        0 ( 798/ 377- 852/ 534)  45.0 E31-E22 S31-S21
      8  382  451 (S18E60 L= 67)   458     83        1 ( 367/ 423- 407/ 473)  36.9 E62-E57 S19-S16
      9  386  370 (S24E65 L= 62)   648    152        1 ( 368/ 345- 407/ 391)  23.2 E68-E61 S25-S22
  

    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