Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  572 1765 (N30E56 L=312)  1361    142        2 ( 516/1707- 609/1813) -16.9 E71-E48 N25-N33
      2  738 1579 (N20E31 L=338)   378    120        0 ( 687/1561- 787/1606)   0.0 E35-E27 N18-N21
      3  886 1582 (N23E22 L=347)   294     42        0 ( 862/1567- 917/1604)  14.0 E23-E19 N21-N24
      4 1300 1531 (N30W05 L= 14)   360    178        0 (1231/1476-1387/1565)  -4.8 W00-W12 N28-N30
      5 1571 1327 (N24W27 L= 36)  1674    392        0 (1469/1180-1698/1467) -29.7 W18-W39 N17-N30
      6  348  909 (S27E44 L=325)   794    199        0 ( 280/ 860- 398/ 995)  45.0 E50-E39 S30-S20
      7  455  746 (S35E34 L=335)  1059    292        3 ( 349/ 661- 570/ 810)  35.0 E45-E24 S42-S28
      8 1709  650 (S10W55 L= 64)   337     35        0 (1684/ 632-1729/ 672)  33.7 W53-W57 S11-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