Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1684 1530 (N15W59 L= 84)   951    349        2 (1570/1501-1784/1568) -25.3 W51-W70 N11-N20
      2 1510 1415 (N13W40 L= 65)   241     77        0 (1468/1382-1543/1450) -59.0 W38-W41 N10-N16
      3 1797 1374 (N03W63 L= 87)   796    205        2 (1751/1326-1843/1418)   3.6 W55-W71 N01-N04
      4 1455 1092 (S03W26 L= 51)   312    107        0 (1422/1045-1488/1156) -81.9 W25-W27 S06-N01
      5 1901 1054 (S17W62 L= 87)   674    117        2 (1890/ 991-1919/1103)  56.3 W59-W64 S20-S14
      6  711  874 (N02E21 L=  4)   220     61        0 ( 691/ 846- 735/ 897) -90.0 E22-E19 N00-N03
      7  640  843 (N02E26 L=358)   227     76        1 ( 607/ 823- 686/ 858) -11.3 E28-E23 N00-N02
      8  460  582 (S09E47 L=338)   553    177        5 ( 397/ 528- 510/ 632) -49.4 E50-E43 S12-S05
      9  404  487 (S14E56 L=329)  1569    280        2 ( 313/ 425- 502/ 555) -20.6 E66-E50 S19-S09
  

    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