Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1320 1242 (N17W14 L=282)  3370    746        4 (1037/ 955-1547/1523) -30.8 E09-W33 N04-N29
      2  309 1343 (N05E56 L=212)   374     53        1 ( 285/1311- 333/1365) -45.0 E58-E54 N02-N06
      3  207 1050 (S14E59 L=209)   470     57        0 ( 182/1030- 226/1074)  36.9 E62-E57 S15-S12
      4  231  956 (S19E55 L=213)   887    230        0 ( 178/ 878- 264/1045)  63.4 E60-E52 S25-S13
      5  127  991 (S19E69 L=200)   976    113        3 ( 110/ 956- 141/1022) -53.1 E71-E66 S21-S17
      6 1576  863 (S00W38 L=306)   269     51        2 (1549/ 841-1605/ 883)   0.0 W35-W40 S01-N00
      7 1320  766 (S11W23 L=291)   689    196        1 (1225/ 724-1415/ 797)  32.5 W17-W28 S14-S07
      8  562  771 (S25E25 L=243)   484    134        2 ( 501/ 740- 624/ 794)   5.7 E30-E20 S25-S23
      9 1055  699 (S20W08 L=277)   743    180        0 ( 973/ 677-1153/ 735)   9.5 W02-W14 S21-S18
     10 1336  308 (S40W43 L=312)  3296    506        3 (1101/ 279-1496/ 364)  12.0 W24-W57 S44-S33
  

    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