On 02/18/2014 01:56 PM, Markus Metz wrote:
On Tue, Feb 18, 2014 at 1:14 PM, manuel.martin
<manuel.mar...@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.mar...@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...:
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.mar...@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.mar...@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@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------------------
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
