[GRASS-user] r.in.gdal for importing categorical ArcGIS binary rasters
manuel.martin
manuel.martin at orleans.inra.fr
Thu Feb 13 04:32:55 PST 2014
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 <benducke at fastmail.fm
> <mailto:benducke at 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 at lists.osgeo.org <mailto: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 <mailto:benducke at fastmail.fm>
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org <mailto:grass-user at 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 at 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------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20140213/aa6c6d7d/attachment-0001.html>
More information about the grass-user
mailing list