On Wed, 2009-07-08 at 10:57 -0500, stephen sefick wrote: > The data I would like to import in GRASS is here: > http://csat.er.usgs.gov/statewide/layers/dem250.html > > I have tried r.in.gdal with no sucess. I will provide anything that > is helpfull. This is my second post, so I don't know what would be > helpful. > thanks for the help, >
Hi. I've downloaded the GRID [1] file. "gdalinfo" works (with some reported errors) on the "dem250/" directory which is stored within "workspl"[2]. So I created a location based on the data itself (using the dem250 directory as a source) and imported in grass with r.in.gdal (again using the dem250 directory as input). --- # enter in _some_ location and create another-one based on the data g.proj -c location=albers_us georef=/geo/geodata/world/us/worksp1/dem250 # exit from current location, enter in newly created one grass64 /geo/grassdb/global/albers_us/PERMANENT/ # import data r.in.gdal in=/geo/geodata/world/us/workspl/dem250 out=dem250 --- There is something strange with the values though (they expand from min = -32687 max = 32645). Don't have the time to dig further, maybe there are details in the meta-data about it (!?). Good Luck, Nikos --- [1] http://csat.er.usgs.gov/download/statewide/dem250grid.zip [2] gdalinfo dem250/ ERROR 3: Attempt to read past EOF in dem250//../info/arc.dir. ERROR 4: Failed to open table .VAT Driver: AIG/Arc/Info Binary Grid Files: dem250/ dem250/dblbnd.adf dem250/hdr.adf dem250/log dem250/prj.adf dem250/sta.adf dem250/vat.adf dem250/w001001.adf dem250/w001001x.adf Size is 7495, 8622 Coordinate System is: PROJCS["unnamed", GEOGCS["NAD83", DATUM["North_American_Datum_1983", SPHEROID["GRS 1980",6378137,298.257222101, AUTHORITY["EPSG","7019"]], TOWGS84[0,0,0,0,0,0,0], AUTHORITY["EPSG","6269"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.0174532925199433, AUTHORITY["EPSG","9108"]], AUTHORITY["EPSG","4269"]], PROJECTION["Albers_Conic_Equal_Area"], PARAMETER["standard_parallel_1",29.5], PARAMETER["standard_parallel_2",45.5], PARAMETER["latitude_of_center",23], PARAMETER["longitude_of_center",-83.5], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["METERS",1]] Origin = (-190574.428637420991436,1327208.631963560357690) Pixel Size = (60.000000000000000,-60.000000000000000) Corner Coordinates: Upper Left ( -190574.429, 1327208.632) ( 85d36'18.54"W, 34d59'5.20"N) Lower Left ( -190574.429, 809888.632) ( 85d29'8.74"W, 30d20'49.23"N) Upper Right ( 259125.571, 1327208.632) ( 80d38'16.82"W, 34d58'7.57"N) Lower Right ( 259125.571, 809888.632) ( 80d48'1.00"W, 30d19'54.45"N) Center ( 34275.571, 1068548.632) ( 83d 7'56.55"W, 32d41'16.79"N) Band 1 Block=20x4 Type=Int16, ColorInterp=Undefined Min=0.000 Max=1371.000 NoData Value=-32768 _______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user