Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  576 1765 (N30E56 L=312)  1628    146        3 ( 518/1707- 617/1816) -18.4 E71-E47 N25-N34
      2  812 1579 (N21E26 L=342)   825    241        2 ( 690/1558- 924/1607)  10.6 E35-E19 N18-N24
      3 1293 1535 (N30W05 L= 13)   374    157        1 (1233/1492-1371/1565)  -5.7 W00-W10 N28-N30
      4 1593 1302 (N23W29 L= 38)  1891    466        1 (1477/1129-1711/1458) -32.5 W18-W41 N15-N30
      5  243  999 (S23E55 L=313)   330     35        0 ( 222/ 979- 263/1021)  21.8 E57-E52 S24-S21
      6  352  905 (S27E43 L=325)   859    197        0 ( 283/ 857- 403/ 989)  39.3 E50-E39 S30-S21
      7  193  925 (S28E61 L=307)   383     60        0 ( 164/ 912- 226/ 944)   0.0 E65-E56 S28-S26
      8  458  742 (S35E33 L=335)  1179    306        1 ( 353/ 652- 575/ 810)  36.9 E44-E24 S43-S28
      9  752  784 (S27E11 L=357)   192     57        1 ( 731/ 762- 775/ 805)  71.6 E12-E09 S28-S25
  

    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