Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  648 1696 (N62E09 L=245)   376     89        0 ( 629/1655- 679/1729)  35.5 E11-E03 N58-N64
      2 1261 1627 (N34W37 L=290)  1251    236        1 (1151/1575-1361/1671) -53.1 W31-W40 N28-N40
      3  434 1613 (N62E42 L=212)   440     53        1 ( 407/1594- 456/1638)  23.2 E46-E37 N59-N63
      4  304 1153 (N31E48 L=205)   783    154        1 ( 235/1129- 386/1182) -16.4 E57-E40 N28-N33
      5 1774  884 (S23W46 L=300)  1820    282        2 (1695/ 788-1860/ 988)  14.0 W37-W57 S29-S19
      6 1683  714 (S30W36 L=289)   608    116        1 (1633/ 692-1738/ 740) -11.3 W31-W41 S32-S28
      7  678  442 (S19E39 L=215)   330     55        2 ( 654/ 412- 700/ 470)  36.9 E40-E36 S20-S17
      8  873  425 (S25E27 L=227)   615    229        0 ( 815/ 390- 927/ 470) -45.0 E30-E22 S27-S21
      9 1532  419 (S48W20 L=273)   466     96        0 (1502/ 376-1571/ 455)  29.7 W16-W23 S50-S45
     10  524  341 (S23E59 L=195)  3385    248        4 ( 444/ 283- 645/ 412) -47.7 E68-E49 S27-S16
  

    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