<div dir="ltr">Nike,<div><br></div><div>The gdalinfo report on the file is as follows.  It is an RGBA file - that is it has four bands which are red, green, blue and alpha.  It appears you are configured to only work with single band pseudo-colored images (images with a palette).  I suppose you will need to generalize your code for this case. </div>
<div><br></div><div><div>gdalinfo jrc_21.tiff </div><div>Driver: GTiff/GeoTIFF</div><div>Files: jrc_21.tiff</div><div>Size is 1009, 753</div><div>Coordinate System is:</div><div>PROJCS["WGS 84 / UTM zone 32N",</div>
<div>    GEOGCS["WGS 84",</div><div>        DATUM["WGS_1984",</div><div>            SPHEROID["WGS 84",6378137,298.257223563,</div><div>                AUTHORITY["EPSG","7030"]],</div>
<div>            AUTHORITY["EPSG","6326"]],</div><div>        PRIMEM["Greenwich",0],</div><div>        UNIT["degree",0.0174532925199433],</div><div>        AUTHORITY["EPSG","4326"]],</div>
<div>    PROJECTION["Transverse_Mercator"],</div><div>    PARAMETER["latitude_of_origin",0],</div><div>    PARAMETER["central_meridian",9],</div><div>    PARAMETER["scale_factor",0.9996],</div>
<div>    PARAMETER["false_easting",500000],</div><div>    PARAMETER["false_northing",0],</div><div>    UNIT["metre",1,</div><div>        AUTHORITY["EPSG","9001"]],</div><div>
    AUTHORITY["EPSG","32632"]]</div><div>Origin = (470363.734899999981280,5073657.708399999886751)</div><div>Pixel Size = (2.060455896927651,-2.061088977423639)</div><div>Metadata:</div><div>  AREA_OR_POINT=Area</div>
<div>  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)</div><div>  TIFFTAG_XRESOLUTION=96</div><div>  TIFFTAG_YRESOLUTION=96</div><div>Image Structure Metadata:</div><div>  INTERLEAVE=PIXEL</div><div>Corner Coordinates:</div><div>
Upper Left  (  470363.735, 5073657.708) (  8d37' 6.69"E, 45d48'57.04"N)</div><div>Lower Left  (  470363.735, 5072105.708) (  8d37' 7.04"E, 45d48' 6.76"N)</div><div>Upper Right (  472442.735, 5073657.708) (  8d38'43.03"E, 45d48'57.35"N)</div>
<div>Lower Right (  472442.735, 5072105.708) (  8d38'43.35"E, 45d48' 7.07"N)</div><div>Center      (  471403.235, 5072881.708) (  8d37'55.03"E, 45d48'32.06"N)</div><div>Band 1 Block=1009x2 Type=Byte, ColorInterp=Red</div>
<div>  NoData Value=0</div><div>Band 2 Block=1009x2 Type=Byte, ColorInterp=Green</div><div>  NoData Value=0</div><div>Band 3 Block=1009x2 Type=Byte, ColorInterp=Blue</div><div>  NoData Value=0</div><div>Band 4 Block=1009x2 Type=Byte, ColorInterp=Alpha</div>
<div>  NoData Value=0</div></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Jul 15, 2013 at 4:07 PM, Nik Sands <span dir="ltr"><<a href="mailto:nixanz@nixanz.com" target="_blank">nixanz@nixanz.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi!<br>
<br>
I'm using the code below as part of a routine to warp images.  For most images it works fine.  However, for one of my users, their images are being converted from full colour to black and white which is not supposed to happen.  I've narrowed down the issue to the code below where 'hCT' remains NULL.  An example of one such source image can be downloaded at:<br>

<br>
        <a href="https://dl.dropboxusercontent.com/u/12436846/jrc_21.tiff" target="_blank">https://dl.dropboxusercontent.com/u/12436846/jrc_21.tiff</a><br>
<br>
I guess I must be doing something wrong in my code, but can't figure it out.  Could somebody please help me to understand how to do this correctly?<br>
<br>
Ie, how to I get the colour table from the image linked to above, and set that colour table for my destination image?<br>
<br>
Cheers,<br>
Nik.<br>
<br>
<br>
<br>
GDALColorTableH hCT = GDALGetRasterColorTable( GDALGetRasterBand(hSrcDS,1) );<br>
<br>
if ( hCT == NULL )<br>
        NSLog(@"Failed to get colour table.");<br>
else<br>
        GDALSetRasterColorTable( GDALGetRasterBand(hDstDS,1), hCT );<br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>---------------------------------------+--------------------------------------<br>I set the clouds in motion - turn up   | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a><br>
light and sound - activate the windows | <a href="http://pobox.com/~warmerdam" target="_blank">http://pobox.com/~warmerdam</a><br>and watch the world go round - Rush    | Geospatial Software Developer<br>
</div>