Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  792 1769 (N52E12 L= 69)   875    154        0 ( 725/1726- 873/1814)  19.7 E18-E04 N48-N55
      2  724 1633 (N41E17 L= 64)   932    224        3 ( 633/1600- 824/1666) -16.4 E25-E08 N37-N44
      3  502 1459 (N30E33 L= 48)   949    199        2 ( 442/1400- 555/1503) -38.7 E38-E28 N25-N33
      4  970 1360 (N19E00 L= 80)   189     50        0 ( 950/1331- 992/1388)  56.3 E01-W01 N17-N20
      5 1101 1328 (N16W07 L= 88)   286     80        1 (1060/1300-1142/1351) -26.6 W05-W09 N14-N17
      6  989 1295 (N15W00 L= 81)   245     69        1 ( 952/1272-1028/1322) -45.0 E02-W02 N13-N16
      7  241 1279 (N21E54 L= 27)   431    122        0 ( 210/1258- 277/1307) -29.7 E57-E50 N19-N23
      8  121 1015 (N06E69 L= 11)  1038    102        0 ( 103/ 969- 145/1058)  39.8 E72-E65 N03-N09
      9  748  795 (S13E19 L= 62)   468     86        1 ( 728/ 746- 772/ 842) -90.0 E20-E17 S15-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