Skip to content

rgdal, multiband CCITT compression

6 messages · Oliver Soong, Michael Sumner

Original
Oliver Soong
# 
Hello,

I can't seem to write a multi-band 1-bit CCITT-compressed GeoTIFF
(using any of the CCITT compression methods).  The file writes, but
when reloaded the data are wrong (I can't tell for sure whether this
is a read or write problem).  I can write a 1-band 1-bit
CCITT-compressed GeoTIFF fine, and I can write a multi-band 1-bit
deflate-compressed GeoTIFF fine.  This is on Windows XP x86, R 2.12.0,
rgdal 0.6-28, sp 0.9-72, GDAL 1.7.1, PROJ.4 4.7.1.

On a similar note, only GDAL seems to be able to handle any multi-band
1-bit GeoTIFFs it produces.  ArcGIS 10 explicitly states (somewhere
deep in their documentation) that multi-band 1-bit TIFFs are not
supported.  Even without compression, ENVI throws an error about an
array with a corrupted descriptor, then crashes out.

Is this a bug, or a limitation of TIFF, CCITT, or GDAL that I'm not aware of?

Thoughts are appreciated,
Oliver
Oliver Soong
# 
Here's an example:




require(rgdal)
formals(paste)$sep <- ""

SP27GTIF <- readGDAL(
  system.file("pictures/SP27GTIF.TIF", package = "rgdal")[1],
  output.dim=c(100,100)
)
img_1bit_1band <- SP27GTIF
img_1bit_1band at data$band1 <- as.integer(SP27GTIF at data$band1 > 85)
img_1bit_2band <- SP27GTIF
img_1bit_2band at data <- data.frame(
  band1 = img_1bit_1band at data$band1,
  band2 = as.integer(SP27GTIF at data$band1 > 170)
)

for(img in c("img_1bit_1band", "img_1bit_2band")) {
  cat(paste("*** Correct Stats: ", img, "\n"))
  print(sapply(get(img)@data, table))
  for(compress in c("NONE", "DEFLATE", "LZW", "CCITTRLE",
                    "CCITTFAX3", "CCITTFAX4")) {
    file <- paste(img, "_", compress, ".tif")
    writeGDAL(get(img), file, type = "Byte",
      options = c("NBITS=1", paste("COMPRESS=", compress)))
    cat(paste("Compression: ", compress, " from ", file, "\n"))
    reread <- readGDAL(file, silent = TRUE)
    print(sapply(reread at data, table))
  }
  cat("\n")
}





The expected output is:

*** Correct Stats: img_1bit_1band
  band1
0  3086
1  6914
Compression: NONE from img_1bit_1band_NONE.tif
  band1
0  3086
1  6914
Compression: DEFLATE from img_1bit_1band_DEFLATE.tif
  band1
0  3086
1  6914
Compression: LZW from img_1bit_1band_LZW.tif
  band1
0  3086
1  6914
Compression: CCITTRLE from img_1bit_1band_CCITTRLE.tif
  band1
0  3086
1  6914
Compression: CCITTFAX3 from img_1bit_1band_CCITTFAX3.tif
  band1
0  3086
1  6914
Compression: CCITTFAX4 from img_1bit_1band_CCITTFAX4.tif
  band1
0  3086
1  6914

*** Correct Stats: img_1bit_2band
  band1 band2
0  3086  8353
1  6914  1647
Compression: NONE from img_1bit_2band_NONE.tif
  band1 band2
0  3086  8353
1  6914  1647
Compression: DEFLATE from img_1bit_2band_DEFLATE.tif
  band1 band2
0  3086  8353
1  6914  1647
Compression: LZW from img_1bit_2band_LZW.tif
  band1 band2
0  3086  8353
1  6914  1647
Compression: CCITTRLE from img_1bit_2band_CCITTRLE.tif
  band1 band2
0  6131  8783
1  3869  1217
Compression: CCITTFAX3 from img_1bit_2band_CCITTFAX3.tif
  band1 band2
0  6131  8783
1  3869  1217
Compression: CCITTFAX4 from img_1bit_2band_CCITTFAX4.tif
  band1 band2
0  6131  8783
1  3869  1217



Everything looks fine except CCITT compression with 2 bands.

I wasn't able to dig up Windows GDAL 1.7.1 binaries, and the 1.7.0
binaries from FW_Tools seem to give similar behavior, although it's a
little hard to check because neither ArcGIS nor ENVI seems to like any
of the multi-band files.
Oliver Soong
# 
Sorry, the word "expected" was a bad choice.  The output is what I
get.  I'm highlighting two tables that should be the same, but aren't.
 I appreciate you taking the time to poke around with me.

*** Correct Stats: img_1bit_2band
 band1 band2
0  3086  8353
1  6914  1647
Compression: CCITTRLE from img_1bit_2band_CCITTRLE.tif
 band1 band2
0  6131  8783
1  3869  1217
On Thu, Nov 4, 2010 at 5:14 PM, Michael Sumner <mdsumner at gmail.com> wrote: