[WIP] Handle TIFFs with more than one image

src/GeoTIFF.jl

@@ -4,7 +4,7 @@
 
 module GeoTIFF
 
-using TiffImages: AbstractTIFF, DenseTaggedImage, WidePixel
+using TiffImages: AbstractTIFF, DenseTaggedImage, StridedTaggedImage, WidePixel
 using ColorTypes: Colorant, Gray
 using MappedArrays: mappedarray
 using StaticArrays: SVector, SMatrix, SA

src/image.jl

@@ -15,7 +15,7 @@ to get the metadata and the image with corrected axes, respectively.
 * The [`GeoTIFF.image`](@ref) function is necessary because 
   the GeoTIFF format swaps the order of the image axes;
 """
-struct GeoTIFFImage{T,N,I<:AbstractTIFF{T,N}} <: AbstractArray{T,N}
+struct GeoTIFFImage{T,I<:AbstractMatrix{T}} <: AbstractMatrix{T}
   tiff::I
   metadata::Metadata
 end
@@ -59,4 +59,11 @@ channel(geotiff::GeoTIFFImage, i) = mappedarray(c -> channel(c, i), image(geotif
 Base.size(geotiff::GeoTIFFImage) = size(geotiff.tiff)
 Base.getindex(geotiff::GeoTIFFImage, i...) = getindex(geotiff.tiff, i...)
 Base.setindex!(geotiff::GeoTIFFImage, v, i...) = setindex!(geotiff.tiff, v, i...)
-Base.IndexStyle(::Type{GeoTIFFImage{T,N,I}}) where {T,N,I} = IndexStyle(I)
+Base.IndexStyle(::Type{GeoTIFFImage{T,I}}) where {T,I} = IndexStyle(I)
+
+struct GeoTIFFImages{I<:GeoTIFFImage} <: AbstractVector{I}
+  geotiffs::Vector{I}
+end
+
+Base.size(geotiff::GeoTIFFImages) = size(geotiff.geotiffs)
+Base.getindex(geotiff::GeoTIFFImages, i) = getindex(geotiff.geotiffs, i)

src/load.jl

@@ -15,17 +15,24 @@ The keyword arguments are forward to the `TiffImages.load` function.
   * [`ModelTransformation`](@ref) with `A` as identity matrix and `b` as vector of zeros.
 """
 function load(fname; kwargs...)
+  geotiff(img, ifd) = GeoTIFFImage(img, _getmetadata(ifd))
   tiff = TiffImages.load(fname; kwargs...)
-  metadata = _getmetadata(tiff)
-  GeoTIFFImage(tiff, metadata)
+  ifds = TiffImages.ifds(tiff)
+  if tiff isa StridedTaggedImage
+    GeoTIFFImages(geotiff.(tiff, ifds))
+  elseif ndims(tiff) == 3
+    imgs = eachslice(tiff, dims=3)
+    GeoTIFFImages(geotiff.(imgs, ifds))
+  else
+    geotiff(tiff, ifds)
+  end
 end
 
 # -----------------
 # HELPER FUNCTIONS
 # -----------------
 
-function _getmetadata(tiff)
-  ifd = TiffImages.ifds(tiff)
+function _getmetadata(ifd)
   geokeydirectory = _gettag(ifd, GeoKeyDirectoryTag, GeoKeyDirectory)
   geodoubleparams = _gettag(ifd, GeoDoubleParamsTag, GeoDoubleParams)
   geoasciiparams = _gettag(ifd, GeoAsciiParamsTag, GeoAsciiParams)