[gdal-dev] Using GRASSDataset

[Aleksey]

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