[GRASS-user] r.in.gdal for importing categorical ArcGIS binary rasters

[Benjamin Ducke]

Ouch. I am beginning to suspect that the ESRI "field"
concept for raster layers does not translate well
to GRASS/GDAL. Sorry for not being able to help.
Maybe someone else on this list has experience with
this type of data structure?

At least now we know that r.in.gdal starts counting
bands with "1"!

Ben

On 13/02/14 11:14, manuel.martin wrote:
> I had already tried this (with "2") ;-). No luck either. Reading your
> email I tried with "0" (maybe would gdal detect only two bands over
> three). It comes out that only "band=1" works....
> 
> 
> GRASS 7.0.svn (lambert93):~ > r.in.gdal --o
> input=/home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/clc/clc06_250m
> band=0 output=testclc
> WARNING: Raster map <testclc> already exists and will be overwritten
> WARNING: Datum <Not_specified_based_on_GRS_1980_ellipsoid> not recognised
>          by GRASS and no parameters found
> Projection of input dataset and current location appear to match
> ERROR 5:
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/clc/clc06_250m:
> GDALDataset::GetRasterBand(0) - Illegal band #
> 
> ERROR: Selected band (0) does not exist
> 
> 
> 
> 
> On 02/13/2014 10:59 AM, Benjamin Ducke wrote:
>> Hmm, the manual page does not say so explicitely,
>> but maybe band numbering starts with "0"?
>> In that case your band 3 would have index number "2".
>>
>> Ben
>>
>> On 13/02/14 10:48, manuel.martin wrote:
>>> Hi Ben,
>>> Thanks for the reply. I tried with the band option but it looks like
>>> gdal does not detect the multiple fields as multiple bands :
>>>
>>> GRASS 7.0.svn (lambert93):~ > r.in.gdal --o \
>>> input=/home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/clc/clc06_250m
>>>
>>> \
>>>> band=3 output=testclc
>>> WARNING: Raster map <testclc> already exists and will be overwritten
>>> WARNING: Datum <Not_specified_based_on_GRS_1980_ellipsoid> not
>>> recognised
>>>           by GRASS and no parameters found
>>> Projection of input dataset and current location appear to match
>>> ERROR 5:
>>> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/clc/clc06_250m:
>>>
>>> GDALDataset::GetRasterBand(3) - Illegal band #
>>>
>>> ERROR: Selected band (3) does not exist
>>>
>>>
>>> Cheers, Manuel
>>>
>>> On 02/12/2014 08:25 PM, Benjamin Ducke wrote:
>>>> Hi Manuel,
>>>>
>>>> GRASS does not support multi-band rasters, so you'll have
>>>> to import each band as a separate raster map.
>>>>
>>>> r.in.gdal has the "band=" option to specify a band number
>>>> to import.
>>>>
>>>> Best,
>>>>
>>>> Ben
>>>>
>>>> On 12/02/14 17:05, manuel.martin wrote:
>>>>> Hi all,
>>>>>
>>>>> I tried to import an ArcGIS binary raster (corine land cover for the
>>>>> French territory) with three fields on each pixels (VALUE, COUNT and
>>>>> CODE_06) using the *r.in.gdal* command. The import works just fine,
>>>>> except that in the resulting raster, in GRASS, I only get one field,
>>>>> which seems to be the VALUE field (and not the CODE_06 field, which
>>>>> I am
>>>>> interested in, and which is a categorical field, although coded with
>>>>> integers).
>>>>>
>>>>> Is there a way to produce a resulting raster with all the fields of
>>>>> the
>>>>> initial ArcGIS layer, or alternatively to choose one unique field?
>>>>> Also,
>>>>> is there a limitation on categorical fields. For instance, what if
>>>>> instead of integer corine land cover codes I have literal labels, i.e.
>>>>> levels of my categorical field coded with strings?
>>>>>
>>>>> Cheers, Manuel
>>>>>
>>>>> -- 
>>>>>
>>>>> ----------------0----0------
>>>>> INRA - InfoSol
>>>>> Centre de recherche d'Orléans
>>>>> 2163 Avenue de la Pomme de Pin
>>>>> CS 40001 ARDON
>>>>> 45075 ORLEANS Cedex 2
>>>>> tel : (33) (0)2 38 41 48 21
>>>>> fax : (33) (0)2 38 41 78 69
>>>>> http://www.gissol.fr
>>>>> ----0----0------------------
>>>>>
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> grass-user mailing list
>>>>> grass-user at lists.osgeo.org
>>>>> http://lists.osgeo.org/mailman/listinfo/grass-user
>>>>>
>>>>
>>>
>>
>>
> 
> 



-- 
Dr. Benjamin Ducke, M.A.
{*} Geospatial Consultant
{*} GIS Developer

  benducke at fastmail.fm