<div dir="ltr">Thanks for the suggestions. I opened the image up in python using gdal, looked at the geotransform and found that the image is ever-so-slightly shifted west (-180.00000000012497). Hence the negative W-E center. Maybe this is why GRASS thinks it has a 180E origin also, though I don't know. <br>
<br>I checked the region - same problem - then set it manually and re-imported. Same results. In the end I was able to assign correct extents using r.region. Though the boundary isn't exactly at 180W, I think nine 0s after the decimal point is a bit much. :-)<br>
<br>Thanks again.<br><br>-Jamie<br><br><div class="gmail_quote">On Mon, Oct 13, 2008 at 1:35 PM, Nikos Alexandris <span dir="ltr"><<a href="mailto:nikos.alexandris@felis.uni-freiburg.de">nikos.alexandris@felis.uni-freiburg.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div><div></div><div class="Wj3C7c">On Mon, 2008-10-13 at 12:41 -0700, Jamie Adams wrote:<br>
> Hello all,<br>
><br>
> When I import a global elevation dataset with extents 180w 90n 180e<br>
> 90s using r.in.gdal, the GRASS raster shows up with w-e extents both<br>
> showing 180E. While this doesn't seem to be affecting raster<br>
> calculations, afaik, it's causing issues with QGIS. It displays it<br>
> with it's origin at 180E and extends to 540E. Has anyone run into<br>
> this problem?<br>
><br>
> gdalinfo:<br>
><br>
> Driver: GTiff/GeoTIFF<br>
> Files: 2_00-res/out-2_00-nozero-res.tif<br>
> Size is 43200, 21600<br>
> Coordinate System is:<br>
> GEOGCS["WGS 84",<br>
> DATUM["WGS_1984",<br>
> SPHEROID["WGS 84",6378137,298.2572235630016,<br>
> AUTHORITY["EPSG","7030"]],<br>
> AUTHORITY["EPSG","6326"]],<br>
> PRIMEM["Greenwich",0],<br>
> UNIT["degree",0.0174532925199433],<br>
> AUTHORITY["EPSG","4326"]]<br>
> Origin = (-180.000000000124970,90.000000000000000)<br>
> Pixel Size = (0.008333333333333,-0.008333333333333)<br>
> Metadata:<br>
> AREA_OR_POINT=Area<br>
> Image Structure Metadata:<br>
> INTERLEAVE=BAND<br>
> Corner Coordinates:<br>
> Upper Left (-180.0000000, 90.0000000) (180d 0'0.00"W, 90d 0'0.00"N)<br>
> Lower Left (-180.0000000, -90.0000000) (180d 0'0.00"W, 90d 0'0.00"S)<br>
> Upper Right ( 180.0000000, 90.0000000) (180d 0'0.00"E, 90d 0'0.00"N)<br>
> Lower Right ( 180.0000000, -90.0000000) (180d 0'0.00"E, 90d 0'0.00"S)<br>
> Center ( -0.0000000, 0.0000000) ( 0d 0'0.00"W, 0d 0'0.00"N)<br>
> Band 1 Block=43200x1 Type=Float32, ColorInterp=Gray<br>
> NoData Value=10<br>
><br>
> <a href="http://r.info" target="_blank">r.info</a>:<br>
><br>
> +----------------------------------------------------------------------------+<br>
> | Layer: 2_00_res Date: Fri Oct 10 17:45:37<br>
> 2008<br>
> | Mapset: PERMANENT Login of Creator: jamie<br>
><br>
> | Location: Bathy<br>
><br>
> | DataBase: /work/grass<br>
><br>
> | Title: ( 2_00_res )<br>
><br>
> | Timestamp: none<br>
><br>
> |----------------------------------------------------------------------------<br>
> |<br>
><br>
> | Type of Map: raster Number of Categories: 255<br>
><br>
> | Data Type: FCELL<br>
><br>
> | Rows: 21600<br>
><br>
> | Columns: 43200<br>
><br>
> | Total Cells: 933120000<br>
><br>
> | Projection: Latitude-Longitude<br>
><br>
> | N: 90N S: 90S Res: 0:00:30<br>
><br>
> | E: 180E W: 180E Res: 0:00:30<br>
><br>
> | Range of data: min = -10814.000000 max = -0.010000<br>
><br>
> |<br>
><br>
> | Data Description:<br>
><br>
> | generated by r.in.gdal<br>
><br>
> |<br>
><br>
> | Comments:<br>
><br>
> | r.in.gdal input="out-2_00-nozero-res.tif" output="2_00_res"<br>
><br>
> |<br>
><br>
> +----------------------------------------------------------------------------+<br>
><br>
<br>
</div></div>1. I am not 100% sure but I think the line "Center ( -0.0000000,<br>
0.0000000) ( 0d 0'0.00"W, 0d 0'0.00"N)" in the geo-metadata should<br>
read enter ( 0.0000000, -0.0000000) ( 0d 0'0.00"W, 0d<br>
0'0.00"N)<br>
<br>
(Could it be that the file is not ok?)<br>
<br>
2. What doeas g.region -p say? If you get the same value for east and<br>
west it's wrong I guess. west should be, well... west: 180W and east:<br>
180E.<br>
<br>
3. Try to set "g.region e=180 s=-90 w=-180 n=90 -p". Re-import again and<br>
check for errors.<br>
<br>
Regards, Nikos<br>
<br>
<br>
</blockquote></div><br></div>