Le mercredi 20 mai 2015 01:06:27, EVANS, JAMES R GS-13 USAF ACC 84 RADES/SCZE a écrit : > Claire, > About a year ago I had a similar issue reading NAIP imagery. I was able to > use Global Mapper to convert the images. After that MapServer using GDAL > could serve the images correctly. I don't know if that will help you, but > it worked for me. Converting 4.3 TB of imagery took a while....
This has been fixed in GDAL 2.0dev. Try beta1 : http://download.osgeo.org/gdal/2.0.0beta1/ Otherwise a workaround to avoid creating geotiffs is to gdal_translate to a VRT, and remove the Alpha color interpretation from the 4th band. > > -----Original Message----- > From: [email protected] > [mailto:[email protected]] On Behalf Of Claire Porter Sent: > Tuesday, May 19, 2015 4:54 PM > To: [email protected] > Subject: [gdal-dev] JP2OpenJPEG Driver Reads 4th Band of NITF as Alpha > > Hello, > > I use the OpenJPEG 2.0 driver to decode the JP2 images within NITF files > (1B level NITF images from DigitalGlobe). For 4 band images, the driver > interprets the 4th band as Alpha. When I warp the imagery using gdalwarp, > the 4th band in the output is all zeros. If I convert the image using > gdal_translate, the 11 bit values in the 4th band are truncated to 255. > Using GTiff creation options of PHOTOMETRIC=RGB or MINISBLACK and ALPHA=NO > do not change the results. > > I found this only occurs with the OpenJPEG driver. The Jasper, Kakadu, and > ECW drivers interpret the 4th band as Undefined and the 4th band behaves > normally in gdalwarp and gdal_translate. I've tried gdal versions from > 1.9.0 to 1.11.1. I thought perhaps the issue would be solved in OpenJPEG > 2.1, and so I tried compiling GDAL 1.11.1 with OpenJPEG 2.1.0. But, that > still resulted in the 4th band being read as Alpha. > > I'd really like to continue using OpenJPEG since it will be so well > supported in GDAL 2.0. Anyone know a work around for this? Is this a > GDAL issue or an OpenJPEG issue? > > Thanks, > Claire > > > gdalinfo using the OpenJPEG driver and GDAL 1.11.1: > > > OPENJPEG: info: Start to read j2k main header (0). > OPENJPEG: info: Main header has been correctly decoded. > GDAL: > GDALOpen(/vsisubfile/3935_198497332,output/WV02_12OCT082110317-M1BS-103001 > 001B68B600.ntf, this=0x948760) succeeds as JP2OpenJPEG. GDAL: > NITFDataset::Open() wasn't able to derive a first order geotransform. It > will be returned as GCPs. GDAL: GDALDefaultOverviews::OverviewScan() > GDAL: GDALOpen(output/WV02_12OCT082110317-M1BS-103001001B68B600.ntf, > this=0x8bfd10) succeeds as NITF. Driver: NITF/National Imagery > Transmission Format > Files: output/WV02_12OCT082110317-M1BS-103001001B68B600.ntf > output/WV02_12OCT082110317-M1BS-103001001B68B600.xml > Size is 9216, 8192 > Coordinate System is `' > GCP Projection = > GEOGCS["WGS 84", > DATUM["WGS_1984", > SPHEROID["WGS 84",6378137,298.257223563, ........... > SAMP_OFF=4409 > SAMP_SCALE=4410 > Corner Coordinates: > Upper Left ( 0.0, 0.0) > Lower Left ( 0.0, 8192.0) > Upper Right ( 9216.0, 0.0) > Lower Right ( 9216.0, 8192.0) > Center ( 4608.0, 4096.0) > Band 1 Block=1024x1024 Type=UInt16, ColorInterp=Blue > Overviews: 4608x4096, 2304x2048, 1152x1024, 576x512, 288x256 > Overviews: arbitrary > Image Structure Metadata: > NBITS=11 > Band 2 Block=1024x1024 Type=UInt16, ColorInterp=Green > Overviews: 4608x4096, 2304x2048, 1152x1024, 576x512, 288x256 > Overviews: arbitrary > Image Structure Metadata: > NBITS=11 > Band 3 Block=1024x1024 Type=UInt16, ColorInterp=Red > Overviews: 4608x4096, 2304x2048, 1152x1024, 576x512, 288x256 > Overviews: arbitrary > Image Structure Metadata: > NBITS=11 > Band 4 Block=1024x1024 Type=UInt16, ColorInterp=Alpha > Overviews: 4608x4096, 2304x2048, 1152x1024, 576x512, 288x256 > Overviews: arbitrary > Image Structure Metadata: > NBITS=11 > GDAL: > GDALClose(/vsisubfile/3935_198497332,output/WV02_12OCT082110317-M1BS-10300 > 1001B68B600.ntf, this=0x94d3f0) GDAL: > GDALClose(/vsisubfile/3935_198497332,output/WV02_12OCT082110317-M1BS-10300 > 1001B68B600.ntf, this=0x948760) GDAL: > GDALClose(output/WV02_12OCT082110317-M1BS-103001001B68B600.ntf, > this=0x8bfd10) > > > gdalinfo using ECW and 1.11.1 (Jasper and Kakadu are similar): > > ECW: > NCScbmOpenFileView(/vsisubfile/3935_198497332,Z:\claire\debug_ortho\output > 2\WV02_12OCT082110317-M1BS-103001001B68B600.ntf): eErr=2, will try VSIL > stream. ECW: Got mutex. > ECW: Using filename 'C:\Users\porte254\AppData\Local\Temp\99544_0.ntf' for > temporary directory determination purposes. ECW: Using filename > 'C:\Users\porte254\AppData\Local\Temp\99544_1.ntf' for temporary directory > determination purposes. ECW: Using filename > 'C:\Users\porte254\AppData\Local\Temp\99544_2.ntf' for temporary directory > determination purposes. ECW: Using filename > 'C:\Users\porte254\AppData\Local\Temp\99544_3.ntf' for temporary directory > determination purposes. ECW: Using filename > 'C:\Users\porte254\AppData\Local\Temp\99544_4.ntf' for temporary directory > determination purposes. ECW: Using filename > 'C:\Users\porte254\AppData\Local\Temp\99544_5.ntf' for temporary directory > determination purposes. ECW: Using filename > 'C:\Users\porte254\AppData\Local\Temp\99544_6.ntf' for temporary directory > determination purposes. ECW: Using filename > 'C:\Users\porte254\AppData\Local\Temp\99544_7.ntf' for temporary directory > determination purposes. ECW: Using filename > 'C:\Users\porte254\AppData\Local\Temp\99544_8.ntf' for temporary directory > determination purposes. ECW: FileInfo: SizeXY=9216,8192 Bands=4 > OriginXY=0,0 CellIncrementXY=1,1 > ColorSpace=3, eCellType=1 > > > GDAL: > GDALOpen(/vsisubfile/3935_198497332,Z:\claire\debug_ortho\output2\WV02_12O > CT082110317-M1BS-103001001B68B600.ntf, this=000000000032FBA0) succeeds as > JP2ECW. GDAL: NITFDataset::Open() wasn't able to derive a first order > geotransform. It will be returned as GCPs. > GDAL: GDALDefaultOverviews::OverviewScan() > GDAL: > GDALOpen(Z:\claire\debug_ortho\output2\WV02_12OCT082110317-M1BS-103001001B > 68B600.ntf, this=000000000031EC40) succeeds as NITF. Driver: NITF/National > Imagery Transmission Format > Files: output2\WV02_12OCT082110317-M1BS-103001001B68B600.ntf > output2\WV02_12OCT082110317-M1BS-103001001B68B600.XML > Size is 9216, 8192 > Coordinate System is `' > GCP Projection = > GEOGCS["WGS 84", > DATUM["WGS_1984", > SPHEROID["WGS 84",6378137,298.257223563, > ..................................................... > SAMP_OFF=4409 > SAMP_SCALE=4410 > Corner Coordinates: > Upper Left ( 0.0, 0.0) > Lower Left ( 0.0, 8192.0) > Upper Right ( 9216.0, 0.0) > Lower Right ( 9216.0, 8192.0) > Center ( 4608.0, 4096.0) > Band 1 Block=256x256 Type=UInt16, ColorInterp=Blue > Overviews: 4608x4096, 2304x2048, 1152x1024, 576x512, 288x256 > Image Structure Metadata: > NBITS=11 > Band 2 Block=256x256 Type=UInt16, ColorInterp=Green > Overviews: 4608x4096, 2304x2048, 1152x1024, 576x512, 288x256 > Image Structure Metadata: > NBITS=11 > Band 3 Block=256x256 Type=UInt16, ColorInterp=Red > Overviews: 4608x4096, 2304x2048, 1152x1024, 576x512, 288x256 > Image Structure Metadata: > NBITS=11 > Band 4 Block=256x256 Type=UInt16, ColorInterp=Undefined > Overviews: 4608x4096, 2304x2048, 1152x1024, 576x512, 288x256 > Image Structure Metadata: > NBITS=11 > > _______________________________________________ > gdal-dev mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev -- Spatialys - Geospatial professional services http://www.spatialys.com _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
