#3048: r.out.gdal does not write EPSG code of the projection --------------------------+------------------------- Reporter: sbl | Owner: grass-dev@… Type: defect | Status: new Priority: normal | Milestone: 7.0.5 Component: Raster | Version: svn-trunk Resolution: | Keywords: r.out.gdal CPU: Unspecified | Platform: All --------------------------+-------------------------
Comment (by neteler): I just tried with the NC dataset on my system, no such issue: {{{ GRASS 7.0.5svn (nc_spm_08_grass7):~ > r.out.gdal elevation output=elevation.tif Checking GDAL data type and nodata value... 100% Using GDAL data type <Float32> Exporting raster data to GTiff format... ERROR 6: SetColorTable() only supported for Byte or UInt16 bands in TIFF format. 100% r.out.gdal complete. File <elevation.tif> created. GRASS 7.0.5svn (nc_spm_08_grass7):~ > gdalinfo elevation.tif Driver: GTiff/GeoTIFF Files: elevation.tif Size is 1500, 1350 Coordinate System is: PROJCS["Lambert Conformal Conic", GEOGCS["NAD83", DATUM["North_American_Datum_1983", SPHEROID["GRS 1980",6378137,298.2572221010042, AUTHORITY["EPSG","7019"]], AUTHORITY["EPSG","6269"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4269"]], PROJECTION["Lambert_Conformal_Conic_2SP"], PARAMETER["standard_parallel_1",36.16666666666666], PARAMETER["standard_parallel_2",34.33333333333334], PARAMETER["latitude_of_origin",33.75], PARAMETER["central_meridian",-79], PARAMETER["false_easting",609601.22], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]]] Origin = (630000.000000000000000,228500.000000000000000) Pixel Size = (10.000000000000000,-10.000000000000000) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 630000.000, 228500.000) ( 78d46'27.41"W, 35d48'34.59"N) Lower Left ( 630000.000, 215000.000) ( 78d46'28.63"W, 35d41'16.54"N) Upper Right ( 645000.000, 228500.000) ( 78d36'29.89"W, 35d48'33.08"N) Lower Right ( 645000.000, 215000.000) ( 78d36'32.01"W, 35d41'15.03"N) Center ( 637500.000, 221750.000) ( 78d41'29.49"W, 35d44'54.91"N) Band 1 Block=1500x1 Type=Float32, ColorInterp=Gray Metadata: COLOR_TABLE_RULES_COUNT=5 COLOR_TABLE_RULE_RGB_0=5.557880e+01 7.572900e+01 0 191 191 0 255 0 COLOR_TABLE_RULE_RGB_1=7.572900e+01 9.587920e+01 0 255 0 255 255 0 COLOR_TABLE_RULE_RGB_2=9.587920e+01 1.160290e+02 255 255 0 255 127 0 COLOR_TABLE_RULE_RGB_3=1.160290e+02 1.361800e+02 255 127 0 191 127 63 COLOR_TABLE_RULE_RGB_4=1.361800e+02 1.563300e+02 191 127 63 200 200 200 Generated_with=GRASS GIS 7.0.5svn GRASS 7.0.5svn (nc_spm_08_grass7):~ > gdal-config --version 2.0.2 }}} But doing the same in a UTM32N location, the AUTHORITY is missing: {{{ GRASS 7.0.5svn (utm32n):~ > g.proj -w PROJCS["UTM Zone 32, Northern Hemisphere", GEOGCS["wgs84", DATUM["WGS_1984", SPHEROID["WGS_1984",6378137,298.257223563]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",9], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",0], UNIT["Meter",1]] # pure GDAL: epsg_tr.py 32632 PROJCS["WGS 84 / UTM zone 32N", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.0174532925199433, AUTHORITY["EPSG","9122"]], AUTHORITY["EPSG","4326"]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",9], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]], AXIS["Easting",EAST], AXIS["Northing",NORTH], AUTHORITY["EPSG","32632"]] }}} I see differences in the PROJCS, perhaps the text "warping" partially fails? -- Ticket URL: <https://trac.osgeo.org/grass/ticket/3048#comment:1> GRASS GIS <https://grass.osgeo.org> _______________________________________________ grass-dev mailing list grass-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-dev