[GRASS-user] r.in.gdal for importing categorical ArcGIS binary rasters
Markus Metz
markus.metz.giswork at gmail.com
Tue Feb 18 12:31:41 PST 2014
On Tue, Feb 18, 2014 at 2:09 PM, manuel.martin
<manuel.martin at orleans.inra.fr> wrote:
> On 02/18/2014 01:56 PM, Markus Metz wrote:
>>
>> 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
>
> Sorry about that, I just sent the landforms/test_lf4.aux.xml to you, but
> unfortunately it does not seem to contain the raster attribute table...:
Right. Unfortunately I can not reproduce the gdalinfo output you sent
previously, and as long as gdal does not see an attribute table, GRASS
can not import the labels in the attribute table. I use gdal 1.10.1,
if that helps.
Markus M
>
>
> GRASS 7.0.svn (lambert93):~/tmp > more
> /home/manuel/Desktop/sftp/Projets/Projets/Zones-humides/base_sig/dsm_data/landforms/test_lf4.aux.xml
> <PAMDataset>
> <SRS>PROJCS["RGF_1993_Lambert_93",GEOGCS["GCS_RGF_1993",DATUM["Reseau_Geodesique_Francais_1993",SPHEROID["GRS_1980",63
> 78137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic_2SP"],PARAM
> ETER["False_Easting",700000.0],PARAMETER["False_Northing",6600000.0],PARAMETER["Central_Meridian",3.0],PARAMETER["Standard_Pa
> rallel_1",44.0],PARAMETER["Standard_Parallel_2",49.0],PARAMETER["Latitude_Of_Origin",46.5],UNIT["Meter",1.0]]</SRS>
> <PAMRasterBand band="1">
> <Description>test_lf4</Description>
> <Histograms>
> <HistItem>
> <HistMin>1</HistMin>
> <HistMax>6</HistMax>
> <BucketCount>256</BucketCount>
> <IncludeOutOfRange>1</IncludeOutOfRange>
> <Approximate>0</Approximate>
> <HistCounts>1313166|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|79933|0|0|0|0|0|0|0|0|0|0|0|0|0|
> 0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|110782|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0
> |0|0|0|0|0|0|0|0|0|0|42435|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|69408|0|0|0|0|0|0|0|0|0|0|0|0|0|0
> |0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|1341719</HistCounts>
> </HistItem>
> </Histograms>
> </PAMRasterBand>
> </PAMDataset>
>
>
>>
>>>
>>> /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------------------
>>>
>
>
> --
>
> ----------------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