Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  321 1538 (N07E66 L=246)  1592    160        1 ( 271/1471- 366/1618)  90.0 E69-E63 N02-N12
      2  773 1565 (N19E31 L=281)   249     78        1 ( 740/1542- 807/1596)  -9.5 E34-E28 N17-N19
      3  945 1425 (N15E16 L=296)   331    105        0 ( 917/1379- 966/1476) -71.6 E17-E14 N11-N17
      4 1524 1080 (N12W27 L=339)   344     92        1 (1498/1034-1558/1128) -26.6 W24-W30 N09-N13
      5 1641 1097 (N16W35 L=347)   323     89        0 (1615/1067-1679/1122)  51.3 W33-W37 N14-N19
      6 1696  881 (N06W45 L=357)  1575    199        2 (1652/ 787-1727/ 981) -45.0 W40-W49 N01-N11
      7 1621  951 (N08W37 L=349)   432     83        3 (1587/ 915-1659/ 975)  45.0 W34-W39 N05-N09
      8 1838  837 (N09W60 L= 12)   544    128        0 (1813/ 793-1858/ 882) -20.6 W56-W64 N07-N10
      9  545  708 (S37E22 L=290)   518    137        1 ( 494/ 675- 611/ 732)  36.9 E26-E17 S39-S33
     10 1675  440 (S17W66 L= 18)  1097    246        1 (1648/ 390-1692/ 508) -24.4 W60-W72 S18-S13
  

    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