In addition, this is the output when CPL_DEBUG is ON: [Fri May 19 09:05:44 2006] [error] [client 145.50.90.64] ECW: NCScbmOpenFileView(/data/geopakhuis/aerogrid_nl2005/Aerogrid_NL2005_40cm.ecw): eErr = 0 [Fri May 19 09:05:44 2006] [error] [client 145.50.90.64] ECW: FileInfo: SizeXY=690000,800000 Bands=3 [Fri May 19 09:05:44 2006] [error] [client 145.50.90.64] OriginXY=8000,624000 CellIncrementXY=0.4,-0.4 [Fri May 19 09:05:44 2006] [error] [client 145.50.90.64] [Fri May 19 09:05:44 2006] [error] [client 145.50.90.64] ECW: projection=LOCAL, datum=WGS84 [Fri May 19 09:05:44 2006] [error] [client 145.50.90.64] GDAL: GDALOpen(/data/geopakhuis/aerogrid_nl2005/Aerogrid_NL2005_40cm.ecw) succeeds as ECW. [Fri May 19 09:05:44 2006] [error] [client 145.50.90.64] [Fri May 19 09:05:44 2006] [error] [client 145.50.90.64] ECWRasterBand: RasterIO(nXOff=443174,nYOff=749014,nXSize=1254,nYSize=751 -> 1252x750) [Fri May 19 09:05:45 2006] [error] [client 145.50.90.64] GDAL: GDALClose(/data/geopakhuis/aerogrid_nl2005/Aerogrid_NL2005_40cm.ecw) [Fri May 19 09:05:45 2006] [error] [client 145.50.90.64] [Fri May 19 09:05:45 2006] [error] [client 145.50.90.64] GDAL: GDALDriver::Create(MEM,msSaveImageGDAL_temp,1252,750,4,Byte,(nil)) [Fri May 19 09:05:45 2006] [error] [client 145.50.90.64] ECW: GetGDTProjDat(28992) = NUTM31/RD [Fri May 19 09:05:45 2006] [error] [client 145.50.90.64] ECW: Writing with PROJ=NUTM31, DATUM=RD [Fri May 19 09:06:00 2006] [error] [client 145.50.90.64] ECW: NCScbmOpenFileView(/data/OGC_UMN_services/tmp/ms_tmp/114802234574260.ecw): eErr = 0 [Fri May 19 09:06:00 2006] [error] [client 145.50.90.64] ECW: FileInfo: SizeXY=1252,750 Bands=4 [Fri May 19 09:06:00 2006] [error] [client 145.50.90.64] OriginXY=185270,324394 CellIncrementXY=0.40016,-0.4 [Fri May 19 09:06:00 2006] [error] [client 145.50.90.64] [Fri May 19 09:06:00 2006] [error] [client 145.50.90.64] ECW: projection=NUTM31, datum=RD [Fri May 19 09:06:00 2006] [error] [client 145.50.90.64] GDAL: GDALOpen(/data/OGC_UMN_services/tmp/ms_tmp/114802234574260.ecw) succeeds as ECW. [Fri May 19 09:06:00 2006] [error] [client 145.50.90.64] [Fri May 19 09:06:00 2006] [error] [client 145.50.90.64] GDAL: GDALClose(msSaveImageGDAL_temp) [Fri May 19 09:06:00 2006] [error] [client 145.50.90.64] [Fri May 19 09:06:00 2006] [error] [client 145.50.90.64] GDAL: GDALClose(/data/OGC_UMN_services/tmp/ms_tmp/114802234574260.ecw) [Fri May 19 09:06:00 2006] [error] [client 145.50.90.64] [Fri May 19 09:06:00 2006] [error] [client 145.50.90.64] GDAL: GDALDeregister_GTiff() called.
Best regards, Bart > Hi list, > > using Mapserver 4.8.3, GDAL 1.3.2 final and the latest ECW SDK we are > running into a strange issue. This is on Red Hat Enterprise 3. > > We use Mapserver WCS to output rasters in ECW, Geotiff and MrSID format. > If we use a Geotiff as a source, all output from Mapserver WCS is okay. If > we use an aerial photo in ECW format as input, all output (ECW, MrSID as > well as Geotiff) comes out black and white instead of colour. WMS output > is just fine. Any ideas what could be causing this? > > Btw, a direct gdal_translate works fine and outputs in colour: > > [EMAIL PROTECTED] gdelivery]# gdal_translate -of "Gtiff" -projwin 185269.9 > 324394.4 185771.2 324093.7 > /data/geopakhuis/aerogrid_nl2005/Aerogrid_NL2005_40cm.ecw > /tmp/testbart.tif > Input file size is 690000, 800000 > Computed -srcwin 443174 749014 1253 752 from projected window. > 0...10...20...30...40...50...60...70...80...90...100 - done. > > These are the output formats used: > > OUTPUTFORMAT > NAME "GEOTIFF" > DRIVER "GDAL/GTiff" > MIMETYPE "image/tiff" > IMAGEMODE "RGBA" > EXTENSION "tif" > FORMATOPTION "COMPRESS=PACKBITS" > END > > OUTPUTFORMAT > NAME "ECW" > DRIVER "GDAL/ECW" > MIMETYPE "image/ecw" > IMAGEMODE "RGBA" > EXTENSION "ecw" > END > > OUTPUTFORMAT > NAME "MrSID" > DRIVER "GDAL/MrSID" > MIMETYPE "image/mrsid" > IMAGEMODE "RGBA" > EXTENSION "sid" > END > > This is the gdalinfo output from the source ECW file: > > [EMAIL PROTECTED] gdelivery]# gdalinfo > /data/geopakhuis/aerogrid_nl2005/Aerogrid_NL2005_40cm.ecw > Driver: ECW/ERMapper Compressed Wavelets > Size is 690000, 800000 > Coordinate System is: > LOCAL_CS["LOCAL", > UNIT["Meter",1]] > Origin = (8000.000000,624000.200000) > Pixel Size = (0.40000000,-0.40000000) > Corner Coordinates: > Upper Left ( 8000.000, 624000.200) > Lower Left ( 8000.000, 304000.200) > Upper Right ( 284000.000, 624000.200) > Lower Right ( 284000.000, 304000.200) > Center ( 146000.000, 464000.200) > Band 1 Block=690000x1 Type=Byte, ColorInterp=Red > Overviews: arbitrary > Band 2 Block=690000x1 Type=Byte, ColorInterp=Green > Overviews: arbitrary > Band 3 Block=690000x1 Type=Byte, ColorInterp=Blue > Overviews: arbitrary > > Thanks in advance. > > Best regards, > Bart >
