Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1761  997 (N15W47 L=282)  2677    460        3 (1708/ 806-1789/1259) -58.0 W36-W54 N04-N28
      2  223 1203 (S15E59 L=176)   833    136        1 ( 183/1153- 266/1246)  50.2 E62-E55 S18-S12
      3 1529 1181 (N17W25 L=261)   233     70        0 (1509/1138-1542/1217) -53.1 W22-W26 N14-N18
      4 1407 1074 (N07W20 L=256)   526    152        0 (1344/1017-1475/1142) -15.3 W14-W25 N05-N09
      5  373 1048 (S21E42 L=194)   478    136        2 ( 318/1024- 419/1074)   6.3 E46-E37 S22-S19
      6  753 1013 (S14E15 L=221)   511    130        0 ( 692/ 995- 828/1033)  20.6 E18-E10 S16-S12
      7  234  778 (S41E52 L=184)   937    291        1 ( 170/ 743- 303/ 813)  21.8 E59-E44 S44-S37
      8  694  743 (S31E12 L=224)   208     75        0 ( 658/ 717- 732/ 775)  -8.1 E15-E08 S31-S29
      9 1177  690 (S20W19 L=255)   472    122        1 (1113/ 667-1241/ 722)   6.3 W14-W23 S21-S18
     10  640  527 (S47E11 L=225)   579    118        0 ( 594/ 498- 687/ 568)  40.6 E14-E07 S49-S42
  

    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