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