Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1407 1493 (N36W17 L=173)  3673    790       10 (1028/1340-1701/1655)  -6.1 E14-W42 N30-N39
      2  741 1553 (N25E30 L=127)   501    155        1 ( 682/1503- 785/1607)   5.6 E36-E26 N22-N27
      3 1171 1250 (N15W04 L=161)   235     66        0 (1160/1214-1183/1287) -63.4 W03-W05 N12-N17
      4 1211 1137 (N09W09 L=166)   244    107        1 (1174/1114-1271/1170)   9.5 W06-W12 N07-N10
      5  941 1059 (S01E06 L=151)   337     95        0 ( 897/1031- 974/1090)  45.0 E07-E03 S03-N00
      6  352  765 (S30E42 L=115)   422     79        0 ( 330/ 728- 376/ 794)  51.3 E44-E40 S32-S27
      7 1497  722 (S10W36 L=193)  1237    273        4 (1377/ 691-1605/ 773)  14.0 W27-W43 S12-S05
      8 1745  673 (S08W59 L=215)   646    109        2 (1703/ 644-1769/ 712)  45.0 W55-W61 S09-S05
      9  558  422 (S51E25 L=132)   387     81        0 ( 518/ 386- 602/ 461) -24.0 E29-E20 S53-S48
  

    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