Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  903 1788 (N61W16 L=287)   647    117        1 ( 852/1763- 956/1811)   0.0 W08-W24 N58-N63
      2 1389 1717 (N40W52 L=323)   538     88        2 (1351/1687-1426/1759)   8.7 W46-W59 N38-N41
      3  450 1627 (N59E47 L=225)   409     47        1 ( 429/1606- 478/1643)  21.8 E50-E43 N56-N59
      4 1120 1591 (N39W22 L=293)  1052    334        1 (1037/1537-1282/1632) -11.3 W13-W33 N34-N43
      5 1734 1522 (N18W69 L=340)  1174    105        0 (1695/1487-1787/1544)  -2.9 W62-W82 N16-N19
      6  344  840 (N07E49 L=222)   411    103        1 ( 302/ 803- 384/ 868)  15.9 E53-E45 N04-N08
      7 1626  690 (S27W34 L=306)   571     86        2 (1608/ 642-1639/ 742) -90.0 W32-W35 S30-S23
      8 1439  687 (S22W19 L=291)   267     69        0 (1418/ 654-1464/ 726)  45.0 W17-W21 S24-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