[gdal-dev] How do I select 3 of 5 bands for output?

[Stephen Woodbridge]

Hi again,

I have a script that was processing MrSID files and need to convert it 
to support jp2 files also. It seems pretty straight forward except the 
jp2 files have bands R,G,B,IR and I need to chop these into GeoTif tiles 
with only R,G,B bands.

The script was iterating over the src image and chopping it into tiles 
using the following command:

system("gdalwarp -srcnodata 0 -dstnodata 0 -s_srs +init=epsg:$epsg 
-t_srs EPSG:4326 -te $xm $ym $xx $yx -rb -wm 250 --config 
GDAL_ONE_BIG_READ ON -co "TILED=YES" $infile $work/$outd/$outf\n");

1) How does this behave with R,G,B,IR as input and how can I get R,G,B 
only into the GeoTif?

2) Is the option -rb still valid? or should I change this to "-r bilinear"?

Thanks for the help,
   -Steve