On 02/18/2014 09:31 PM, Markus Metz wrote:
On Tue, Feb 18, 2014 at 2:09 PM, manuel.martin
<[email protected]> wrote:
On 02/18/2014 01:56 PM, Markus Metz wrote:
On Tue, Feb 18, 2014 at 1:14 PM, manuel.martin
<[email protected]> wrote:
On 02/18/2014 12:47 PM, Markus Metz wrote:
On Tue, Feb 18, 2014 at 11:43 AM, manuel.martin
<[email protected]> 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
I am running GDAL 1.10.0, maybe that explains the difference. I will let you know when I will have 1.10.1 running if I still get the same gdalinfo results.

Manuel

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[&quot;RGF_1993_Lambert_93&quot;,GEOGCS[&quot;GCS_RGF_1993&quot;,DATUM[&quot;Reseau_Geodesique_Francais_1993&quot;,SPHEROID[&quot;GRS_1980&quot;,63
78137.0,298.257222101]],PRIMEM[&quot;Greenwich&quot;,0.0],UNIT[&quot;Degree&quot;,0.0174532925199433]],PROJECTION[&quot;Lambert_Conformal_Conic_2SP&quot;],PARAM
ETER[&quot;False_Easting&quot;,700000.0],PARAMETER[&quot;False_Northing&quot;,6600000.0],PARAMETER[&quot;Central_Meridian&quot;,3.0],PARAMETER[&quot;Standard_Pa
rallel_1&quot;,44.0],PARAMETER[&quot;Standard_Parallel_2&quot;,49.0],PARAMETER[&quot;Latitude_Of_Origin&quot;,46.5],UNIT[&quot;Meter&quot;,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
<[email protected]> wrote:
On 02/14/2014 05:08 PM, Markus Metz wrote:
On Wed, Feb 12, 2014 at 5:05 PM, manuel.martin
<[email protected]> 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
[email protected]
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------------------



--

----------------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
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user

Reply via email to