Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  460 1587 (N28E46 L=273)   534    184        0 ( 393/1568- 518/1609)  -5.2 E51-E40 N26-N29
      2 1582 1463 (N27W35 L=354)   480    150        1 (1512/1444-1641/1478)   5.2 W29-W40 N26-N28
      3  128  886 (S16E68 L=251)  2944    455        6 (  84/ 721- 180/ 989)  86.4 E75-E61 S26-S10
      4 1662  908 (S06W42 L=  1)   307    119        0 (1641/ 839-1678/ 961) -66.8 W40-W43 S10-S03
      5  291  776 (S23E51 L=268)   993    275        3 ( 244/ 677- 324/ 889)  72.9 E56-E47 S29-S16
      6  255  808 (S21E54 L=265)   393    104        3 ( 227/ 782- 279/ 846)  36.9 E56-E51 S22-S18
      7  702  786 (S21E19 L=300)   403     70        2 ( 661/ 759- 737/ 825)  31.0 E21-E16 S22-S18
      8 1861  686 (S16W68 L= 26)   526     85        2 (1838/ 657-1877/ 710)  56.3 W64-W69 S17-S14
      9 1637  598 (S25W47 L=  6)   488    106        0 (1590/ 554-1673/ 647)  66.8 W44-W49 S28-S21
     10 1004  538 (S35W03 L=322)   239    129        0 ( 944/ 484-1064/ 571)  45.0 E01-W06 S39-S32
     11  370  503 (S41E50 L=269)   361     82        0 ( 345/ 483- 398/ 529) -11.3 E52-E47 S42-S39
     12 1629  441 (S35W53 L= 12)  2895    528        7 (1460/ 378-1808/ 507)  28.3 W41-W68 S41-S27
  

    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