Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1174 1677 (N38E02 L=243)   723    221        2 (1074/1623-1270/1718) -10.6 E10-W06 N35-N39
      2 1029 1529 (N24E08 L=237)   437    165        0 ( 956/1476-1096/1569) -16.7 E13-E03 N21-N26
      3  352  898 (S23E44 L=201)   651    193        0 ( 301/ 817- 382/ 978)  63.4 E49-E41 S27-S17
      4 1020  924 (S12W01 L=247)   216    142        0 ( 966/ 883-1081/ 962)  50.2 E01-W04 S15-S09
      5 1092  843 (S16W07 L=252)   204     95        0 (1048/ 822-1122/ 873)  36.9 W04-W08 S17-S13
      6 1024  689 (S26W06 L=251)  2488    911        7 ( 636/ 583-1423/ 835)  28.0 E20-W27 S36-S10
      7 1678  593 (S19W55 L=300)  1299    268        4 (1584/ 536-1749/ 676)  45.0 W48-W59 S23-S12
      8 1055  195 (S62W32 L=277)   932    126        2 (1004/ 154-1094/ 247) -13.4 W21-W42 S65-S57
  

    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