Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1145 1759 (N42E17 L=198)   407     45        1 (1123/1728-1164/1785) -38.7 E19-E14 N39-N44
      2  757 1712 (N27E39 L=176)   345     74        0 ( 725/1673- 777/1740) -20.6 E43-E35 N24-N28
      3 1595 1638 (N51W25 L=240)   437     47        0 (1571/1614-1624/1656)  36.9 W21-W29 N47-N54
      4 1653 1632 (N54W33 L=248)   708     62        0 (1629/1601-1686/1652)  36.9 W28-W40 N49-N58
      5 1786 1117 (N21W47 L=262)   317     63        2 (1763/1085-1810/1149) -21.8 W44-W49 N19-N22
      6  419 1068 (S19E39 L=176)   316     60        0 ( 392/1041- 446/1090)  63.4 E40-E36 S21-S17
      7  371  909 (S30E41 L=175)   762    123        3 ( 351/ 840- 395/ 985)  90.0 E42-E38 S33-S25
      8  249  834 (S37E54 L=161)   791    153        1 ( 193/ 775- 300/ 886) -22.6 E61-E47 S39-S33
      9  315  605 (S51E50 L=165)   493     59        0 ( 291/ 581- 340/ 636) -18.4 E53-E46 S52-S48
     10 1085  251 (S51W40 L=256)   891    115        0 (1041/ 210-1113/ 307) -48.4 W35-W45 S54-S45
     11 1187  248 (S47W48 L=263)   360     50        0 (1155/ 234-1218/ 262)  21.8 W44-W51 S48-S45
  

    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