[GRASS-user] r.in.gdal for importing categorical ArcGIS binary rasters
Markus Metz
markus.metz.giswork at gmail.com
Tue Feb 18 04:56:33 PST 2014
On Tue, Feb 18, 2014 at 1:14 PM, manuel.martin
<manuel.martin at orleans.inra.fr> wrote:
> On 02/18/2014 12:47 PM, Markus Metz wrote:
>>
>> On Tue, Feb 18, 2014 at 11:43 AM, manuel.martin
>> <manuel.martin at orleans.inra.fr> wrote:
>>>
>>> Hi, I just compiled grass r59073 and gave a try.
>>> I could not get the attributes imported along with the raster, but maybe
>>> am
>>> I missing something from the new revision of r.in.gdal (although from the
>>> code it looks that the import of values from the attribute table should
>>> be
>>> automatic?).
>>
>> Yes, r.in.gdal automatically imports whatever it finds as color tables
>> or attribute labels.
>>
>>> Markus, I just sent you the ESRI layer I am trying to import, in case you
>>> would like to try as well.
>>
>> GDAL does not see an attribute table in the data you sent to me, thus
>> GRASS can not import any attributes. If there is really an attribute
>> table, it is not recognized by GDAL.
>
> Then maybe this is a gdal problem : I'm saying this because through the
> command gdalinfo, it looks like gdal detects an attribute table (see below).
>
> Cheers
>
>
> GRASS 7.0.svn (lambert93):~/tmp > gdalinfo
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4
> Driver: AIG/Arc/Info Binary Grid
> Files:
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4
This file is missing in the archive you sent:
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4.aux.xml
I am pretty sure that it contains the raster attribute table.
Markus M
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4/hdr.adf
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4/dblbnd.adf
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4/w001001.adf
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4/w001001x.adf
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4/sta.adf
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4/vat.adf
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4/log
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4/prj.adf
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4/metadata.xml
> Size is 2504, 2089
> Coordinate System is:
> PROJCS["unnamed",
> GEOGCS["Unknown datum based upon the GRS 1980 ellipsoid",
> DATUM["Not_specified_based_on_GRS_1980_ellipsoid",
> SPHEROID["GRS 1980",6378137,298.257222101,
> AUTHORITY["EPSG","7019"]],
> AUTHORITY["EPSG","6019"]],
> PRIMEM["Greenwich",0,
> AUTHORITY["EPSG","8901"]],
> UNIT["degree",0.0174532925199433,
> AUTHORITY["EPSG","9122"]],
> AUTHORITY["EPSG","4019"]],
> PROJECTION["Lambert_Conformal_Conic_2SP"],
> PARAMETER["standard_parallel_1",44],
> PARAMETER["standard_parallel_2",49],
> PARAMETER["latitude_of_origin",46.5],
> PARAMETER["central_meridian",3],
> PARAMETER["false_easting",700000],
> PARAMETER["false_northing",6600000],
> UNIT["METERS",1]]
> Origin = (198837.500000000000000,6808087.500000000000000)
> Pixel Size = (50.000000000000000,-50.000000000000000)
> Corner Coordinates:
> Upper Left ( 198837.500, 6808087.500) ( 3d45' 3.55"W, 48d10'48.77"N)
> Lower Left ( 198837.500, 6703637.500) ( 3d37'59.07"W, 47d14'36.54"N)
> Upper Right ( 324037.500, 6808087.500) ( 2d 4'11.49"W, 48d15'52.06"N)
> Lower Right ( 324037.500, 6703637.500) ( 1d58'52.06"W, 47d19'34.69"N)
> Center ( 261437.500, 6755862.500) ( 2d51'32.26"W, 47d45'23.95"N)
> Band 1 Block=256x16 Type=Byte, ColorInterp=Undefined
> Description = test_lf4
> Min=1.000 Max=6.000
> NoData Value=255
> <GDALRasterAttributeTable>
> <FieldDefn index="0">
> <Name>VALUE</Name>
> <Type>0</Type>
> <Usage>5</Usage>
> </FieldDefn>
> <FieldDefn index="1">
> <Name>COUNT</Name>
> <Type>0</Type>
> <Usage>1</Usage>
> </FieldDefn>
> <FieldDefn index="2">
> <Name>CLASS_NAME</Name>
> <Type>2</Type>
> <Usage>0</Usage>
> </FieldDefn>
> <Row index="0">
> <F>1</F>
> <F>1313166</F>
> <F>Valleys</F>
> </Row>
> <Row index="1">
> <F>2</F>
> <F>79933</F>
> <F>Lower Slopes</F>
> </Row>
> <Row index="2">
> <F>3</F>
> <F>110782</F>
> <F>Gentle Slopes</F>
> </Row>
> <Row index="3">
> <F>4</F>
> <F>42435</F>
> <F>Steep Slopes</F>
> </Row>
> <Row index="4">
> <F>5</F>
> <F>69408</F>
> <F>Upper Slopes</F>
> </Row>
> <Row index="5">
> <F>6</F>
> <F>1341719</F>
> <F>Ridges</F>
> </Row>
> </GDALRasterAttributeTable>
>
>
>>
>> Markus M
>>
>>> Thank you in advance
>>>
>>> Manuel
>>>
>>>
>>>
>>> On 02/17/2014 02:48 PM, Markus Metz wrote:
>>>>
>>>> On Fri, Feb 14, 2014 at 5:41 PM, manuel.martin
>>>> <manuel.martin at orleans.inra.fr> wrote:
>>>>>
>>>>> On 02/14/2014 05:08 PM, Markus Metz wrote:
>>>>>>
>>>>>> On Wed, Feb 12, 2014 at 5:05 PM, manuel.martin
>>>>>> <manuel.martin at orleans.inra.fr> 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).
>>>>>>
>>>>>> r.in.gdal will soon support import of category labels such as those
>>>>>> stored in the CODE_06 field. You can then inspect the labels with
>>>>>> r.category. GDAL supports reading and writing of raster attribute
>>>>>> tables, so far only the color rules stored in a raster attribute table
>>>>>> are imported, but I plan to add import of labels to r.in.gdal. Note
>>>>>> that this applies only to GRASS 7.
>>>>>>
>>>>>> Markus M
>>>>>
>>>>>
>>>>> Good to hear that ;-). I will wait for the next release of r.in.gdal
>>>>
>>>> Done in r59073.
>>>>
>>>> Markus M
>>>>
>>>>
>>>>> (meanwhile I did import clc vector layers with v.in.ogr and later
>>>>> managed
>>>>> the category labels with v.to.rast).
>>>>>
>>>>> All the best, Manuel
>>>>>
>>>>>>> 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
>>>>>
>>>>>
>>>>>
>>>>> --
>>>>>
>>>>> ----------------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------------------
>>>>>
>>>
>>> --
>>>
>>> ----------------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------------------
>>>
>
>
> --
>
> ----------------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------------------
>
More information about the grass-user
mailing list