I think geoTiffs do not handle categorical variables (can anyone confirm?), neither the writeRAST6() function from spgrass6 (tested yesterday), and that's why I looked for a full GRASS option (avoiding using R to go from ArcGIS to GRASS). So it is true that I can import categorical ArcGIS binary rasters into R, but then I'm stuck into R and cannot move to GRASS...

Manuel

On 02/13/2014 11:58 AM, Nuno Sá wrote:
Out of the box:

Import it in R and save it as a Tiff - use GRASS to open this tiff. Maybe it works ;P


On 13 February 2014 10:50, Benjamin Ducke <bendu...@fastmail.fm <mailto:bendu...@fastmail.fm>> wrote:

    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@lists.osgeo.org <mailto:grass-user@lists.osgeo.org>
    >>>>> http://lists.osgeo.org/mailman/listinfo/grass-user
    >>>>>
    >>>>
    >>>
    >>
    >>
    >
    >



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

    bendu...@fastmail.fm <mailto:bendu...@fastmail.fm>
    _______________________________________________
    grass-user mailing list
    grass-user@lists.osgeo.org <mailto:grass-user@lists.osgeo.org>
    http://lists.osgeo.org/mailman/listinfo/grass-user




--

Nuno César de Sá
+351 91 961 90 37



_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


--

----------------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@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Reply via email to