Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1010 1876 (N60W33 L=225)  1746    236        2 ( 935/1810-1116/1937)   7.3 W19-W58 N55-N62
      2 1225 1741 (N42W36 L=228)   315     48        0 (1201/1714-1251/1762)  18.4 W33-W39 N39-N42
      3 1327 1732 (N38W44 L=235)   380     52        1 (1295/1717-1360/1751)  -9.5 W40-W46 N37-N39
      4  252 1161 (N25E52 L=139)   745    158        1 ( 238/1090- 273/1247) -68.2 E54-E49 N20-N30
      5 1675 1017 (S12W41 L=232)   280     61        0 (1655/ 991-1701/1048) -63.4 W39-W42 S13-S09
      6  831  722 (S11E17 L=174)   359    148        0 ( 792/ 658- 859/ 792) -90.0 E18-E16 S15-S06
      7 1884  713 (S36W67 L=259)   711     52        2 (1852/ 692-1915/ 739)  -6.7 W60-W77 S37-S33
      8 1134  670 (S21E01 L=191)   593    266        0 (1034/ 596-1240/ 738)  12.1 E08-W06 S22-S17
      9  248  665 (S04E63 L=128)   387     26        0 ( 231/ 643- 265/ 690)  21.8 E65-E60 S05-S02
     10  886  355 (S36E27 L=165)  1203    319        2 ( 744/ 336-1047/ 385) -24.2 E35-E14 S39-S30
  

    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