[GRASS-user] Trouble importing raster dataset

[Thomas Adams]

Markus,

I'm not sure what I had done wrong previously when I sent my last email,
but the steps:

http://www.cgiar-csi.org/data/global-aridity-and-pet-database

gdal_translate -a_ullr -180 90 180 -60 w001001.adf w001001.adf.fixed
r.in.gdal input=/home/thomas.adams/data/PET_he_monthly/pet_he_12/w001001.adf.fixed
output=pet_he_12

worked fine. Thank you!

Best,
Tom

On Wed, Jul 6, 2016 at 3:33 PM, Markus Neteler <neteler at osgeo.org> wrote:

> Hi,
>
> On Tue, Jul 5, 2016 at 11:36 PM, Thomas Adams <tea3rd at gmail.com> wrote:
> > Here is more info from error from gdalinfo...
> >
> > gdalinfo w001001.adf
> ...
> > Coordinate System is:
> > GEOGCS["WGS 84",
> >     DATUM["WGS_1984",
> >         SPHEROID["WGS 84",6378137,298.257223563,
> ...
> >         AUTHORITY["EPSG","9108"]],
> >     AUTHORITY["EPSG","4326"]]
> > Origin = (-180.000000000000000,90.000007823109627)
> > Pixel Size = (0.008333333767951,-0.008333333767951)
> > Corner Coordinates:
> > Upper Left  (-180.0000000,  90.0000078) (180d 0' 0.00"W, 90d 0' 0.03"N)
>
> ^-- as you can see it exceeds the north pole (> 90.0 degN)
>
> > Lower Left  (-180.0000000, -60.0000000) (180d 0' 0.00"W, 60d 0' 0.00"S)
> > Upper Right ( 180.0000188,  90.0000078) (180d 0' 0.07"E, 90d 0' 0.03"N)
>
> ... in addition, it exceeds 180W as well. Essentially a broken dataset.
>
> Solutions:
> https://grass.osgeo.org/grass70/manuals/r.in.gdal.html
> * use flag
>   -l Force Lat/Lon maps to fit into geographic coordinates (90N,S; 180E,W)
>
> or see the trick of applying
> gdal_translate -a_ullr ...
>
> which I added here (the WorldClim data suffer from the same issue):
> https://grass.osgeo.org/grass70/manuals/r.in.gdal.html#worldclim.org-data
>
> HTH
> Markus
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20160726/c3b59156/attachment.html>