Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  341 1543 (N07E64 L=246)  1808    211        3 ( 290/1465- 382/1624) -80.5 E68-E60 N01-N13
      2  781 1569 (N19E31 L=279)   313    121        1 ( 730/1533- 827/1602)  -6.3 E35-E26 N17-N20
      3  974 1410 (N15E14 L=296)   284    107        1 ( 946/1364- 999/1461)  90.0 E15-E12 N11-N18
      4 1552 1069 (N12W29 L=339)   425    137        1 (1524/1023-1583/1109) -31.0 W26-W32 N10-N14
      5 1662 1081 (N16W37 L=347)   339    103        1 (1633/1047-1698/1108)  59.0 W35-W39 N13-N18
      6 1718  859 (N06W47 L=357)  1482    258        5 (1667/ 772-1749/ 949) -45.0 W42-W53 N01-N10
      7 1856  823 (N09W63 L= 13)   637    130        1 (1827/ 780-1882/ 855) -14.0 W58-W68 N07-N10
      8  554  692 (S37E21 L=289)   470    134        1 ( 505/ 659- 621/ 720)  36.9 E24-E16 S40-S34
      9 1687  439 (S16W67 L= 17)  1061    242        3 (1666/ 384-1700/ 504) -31.0 W62-W73 S18-S12
  

    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