Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1715 1553 (N23W65 L=202)   985    233        1 (1670/1529-1767/1570) -14.0 W59-W71 N20-N24
      2 1724 1401 (N14W58 L=195)  2117    330        4 (1614/1269-1813/1521) -69.8 W51-W62 N04-N23
      3  542 1466 (N42E29 L=108)   281     60        1 ( 512/1436- 577/1495)  18.4 E31-E25 N40-N43
      4  386 1449 (N43E44 L= 94)   634    116        0 ( 350/1405- 444/1483)  29.7 E47-E38 N40-N45
      5 1698 1084 (S03W45 L=183)   543    139        0 (1665/1030-1715/1159)  56.3 W41-W47 S05-N01
      6 1384  956 (S04W20 L=158)   268     83        0 (1341/ 932-1429/ 977)   9.5 W17-W23 S04-S02
      7 1478  669 (S24W23 L=161)   526    130        0 (1418/ 642-1529/ 696) -20.6 W19-W27 S25-S20
  

    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