On Tue, Feb 18, 2014 at 2:09 PM, manuel.martin <manuel.mar...@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.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...:
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.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