search

yeesian / ArchGDAL.jl

A high level API for GDAL - Geospatial Data Abstraction Library
https://yeesian.github.io/ArchGDAL.jl/stable/
Other
141 stars 27 forks source link

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

Open igaszczesniak opened 2 years ago

maxfreu commented 2 years ago

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

igaszczesniak commented 2 years ago

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
maxfreu commented 2 years ago

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.

yeesian commented 2 years ago

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?

visr commented 2 years ago

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/

rafaqz commented 2 years ago

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.