[gdal-dev] Using GRASSDataset

[Frank Warmerdam]

 On Wed, Jan 23, 2002 at 02:23:30PM -0500, Aleksey wrote:
> Frank,
> 
> thanks for the clarification. 
> 
> > You should install libgrass 1.0.0, and then re-configure GDAL and check
> > the configure output to ensure it is finding and using the GRASS library.
> 
> This worked. Had to do a clean build of GDAL, but it worked.
> 
> > The path to the file in the cellhd directory normally. I have written up
> > a format page for GRASS support in GDAL and updated the format list. I
> > should have done that a long time ago! See:
> >
> > http://www.remotesensing.org/gdal/frmt_grass.html
> 
> Thanks for the docs. "gdalinfo" is able to open GRASS dataset with both the 
> full path and just the file name in "cellhd". However, it can't report 
> min/max:
> 
> [aleksey at botik PERMANENT]$ gdalinfo -mm 
> /home/Geo/GRASS/data/WNY/PERMANENT/cellhd/t10
> Driver: GRASS/GRASS Database Rasters
> Size is 1066, 1419
> Coordinate System is:
> PROJCS["UTM Zone 17, Northern Hemisphere",
> GEOGCS["NAD27",
> DATUM["North_American_Datum_1927",
> SPHEROID["Clarke 1866",6378206.4,294.978698213898,
> AUTHORITY["EPSG",7008]],
> TOWGS84[-3,142,183,0,0,0,0],
> AUTHORITY["EPSG",6267]],
> PRIMEM["Greenwich",0,
> AUTHORITY["EPSG",8901]],
> UNIT["degree",0.0174532925199433,
> AUTHORITY["EPSG",9108]],
> AXIS["Lat","NORTH"],
> AXIS["Long","EAST"],
> AUTHORITY["EPSG",4267]],
> PROJECTION["Transverse_Mercator"],
> PARAMETER["latitude_of_origin",0],
> PARAMETER["central_meridian",-81],
> PARAMETER["scale_factor",0.9996],
> PARAMETER["false_easting",500000],
> PARAMETER["false_northing",0]]
> Origin = (694765.000000,4721995.000000)
> Pixel Size = (10.000000,-10.000000)
> Corner Coordinates:
> Upper Left ( 694765.000, 4721995.000) ( 78d37'29.65"W, 42d37'39.78"N)
> Lower Left ( 694765.000, 4707805.000) ( 78d37'47.08"W, 42d30'0.12"N)
> Upper Right ( 705425.000, 4721995.000) ( 78d29'42.09"W, 42d37'29.82"N)
> Lower Right ( 705425.000, 4707805.000) ( 78d30'0.47"W, 42d29'50.20"N)
> Center ( 700095.000, 4714900.000) ( 78d33'44.83"W, 42d33'45.05"N)
> Band 1 Block=1066x1 Type=Float32, ColorInterp=Gray
> Min=-4294967295.000/0, Max=4294967295.000/0, Computed Min/Max=nan,nan
> [aleksey at botik PERMANENT]$
> 
> GRASS gives me ("r.info"):
> ...
> Range of data: min = 344.000000 max = 614.599976
> ...
> 
> If this helps, the "t10" dataset was imported into GRASS with "r.in.gdal" 
> from a USGS ASCII DEM format that I bugged you about (and you added, thanks!) 
> not so long ago.
> 
> Aleksey
> 
> 
> 
> 
> Your use of Yahoo! Groups is subject to http://docs.yahoo.com/info/terms/ 
> 
> 
> 

-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up | Frank Warmerdam, warmerdam at p...
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush | Geospatial Programmer for Rent