Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1148 1876 (N59W04 L=301)   400     41        0 (1119/1862-1175/1887)  23.2 W00-W07 N57-N60
      2  993 1863 (N55E13 L=285)   548     87        2 ( 948/1838-1038/1882)  12.5 E17-E08 N52-N57
      3  801 1855 (N53E32 L=265)   474     49        2 ( 778/1826- 824/1880) -71.6 E35-E29 N49-N55
      4 1133 1642 (N36W03 L=300)   237     59        1 (1110/1606-1147/1671) -59.0 W01-W04 N33-N38
      5  383 1624 (N31E59 L=239)   687     83        1 ( 357/1584- 413/1653)  90.0 E63-E56 N28-N33
      6  740 1564 (N28E24 L=274)   580    180        1 ( 681/1511- 816/1587)  26.6 E27-E18 N23-N29
      7 1411 1448 (N24W23 L=320)   284    110        0 (1361/1422-1459/1470) -15.9 W19-W26 N22-N25
      8  904  824 (S17E06 L=291)   270    100        1 ( 880/ 777- 936/ 864)  71.6 E07-E04 S20-S14
      9  636  744 (S24E24 L=274)   571    157        0 ( 576/ 707- 719/ 770)  21.8 E28-E18 S26-S21
     10 1008  487 (S39W03 L=300)   594    123        4 ( 972/ 440-1061/ 529)  56.3 W00-W06 S42-S35
  

    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