ArchGDAL.jl icon indicating copy to clipboard operation
ArchGDAL.jl copied to clipboard

Unknown GCI: ArchGDAL.GDALColorInterp[ArchGDAL.GCI_Undefined, ArchGDAL.GCI_GrayIndex]. while using ArchGDAL.imread(dataset)

Open igaszczesniak opened this issue 3 years ago • 6 comments
trafficstars

igaszczesniak avatar Aug 19 '22 09:08 igaszczesniak

Please provide a description of the error, some minimum failing code and if possible sample data.

maxfreu avatar Aug 19 '22 11:08 maxfreu

So, I would like to display the multiple channel (13) image using the ArchGDAL package. I can read the file but later ArchGDAL.imread(dataset) command does not work;

image

igaszczesniak avatar Aug 19 '22 11:08 igaszczesniak

I have no solution to your plotting problem using ArchGDAL, but I once wrote a plotting function for the Rasters package, which builds on top of ArchGDAL. You find it here https://github.com/rafaqz/Rasters.jl/issues/182#issuecomment-1027817429.

maxfreu avatar Aug 19 '22 12:08 maxfreu

Hi, thanks for filing the issue! I think we encouraged it at https://github.com/yeesian/ArchGDAL.jl/blob/d6edcde5ff88d9eac837e5b41e2f77bb2303f3b5/src/raster/images.jl#L177-L179 and appreciate you taking the time.

In the case of your dataset, it is unclear how you want to display the multi-channel (13) raster, so the package was unable to do so. If you wish to visualize them, you can possibly try reading and displaying them a channel at a time?

yeesian avatar Aug 19 '22 13:08 yeesian

It might be nice to add something like ArchGDAL.imread(path; red=3, blue=5, green=6)

https://www.esri.com/arcgis-blog/products/product/imagery/band-combinations-for-landsat-8/

visr avatar Aug 19 '22 14:08 visr

It might be nice to add something like ArchGDAL.imread(path; red=3, blue=5, green=6)

https://www.esri.com/arcgis-blog/products/product/imagery/band-combinations-for-landsat-8/

Was thinking the same thing recently.

rafaqz avatar Aug 19 '22 14:08 rafaqz