[GRASS-user] UTM 33x

Helmut Kudrnovsky hellik at web.de
Sun Jul 17 05:46:33 PDT 2016


>> please post the output of g.proj -p
>> (https://grass.osgeo.org/grass73/manuals/g.proj.html)
>
>XY location (unprojected)
>
>(that is from =g.proj -p= run in the GRASS location where the GRASS CLI
creates the location via the "-c" flag and the file I then import.
>
>> and also post the output of: gdalinfo file.grid
>
>Driver: GS7BG/Golden Software 7 Binary Grid (.grd)
>Files: file.grd
>Size is 101, 166
>Coordinate System is `'
>Origin = (509450.000000000000000,8562350.000000000000000)
>Pixel Size = (100.000000000000000,-100.000000000000000)
>Corner Coordinates:
>Upper Left  (  509450.000, 8562350.000)
>Lower Left  (  509450.000, 8545750.000)
>Upper Right (  519550.000, 8562350.000)
>Lower Right (  519550.000, 8545750.000)
>Center      (  514500.000, 8554050.000)
>Band 1 Block=101x1 Type=Float64, ColorInterp=Undefined
>  Min=-2.316 Max=385.699
>  NoData Value=1.70141000000000007e+38
> 
>> or try: testepsg file.grid (see attached the file
>> gdal_testepsg_espsg_32633.txt
>
>ERROR 1: Error occured translating file.grd. 

and from your original report in this thread:

>WARNING: <PROJ_INFO> file not found for location <foo>
>ERROR: Unable to get projection info of input map

as you can see your data has no projection information, thus imported to an
xy location (that means a location without any projection info).

GRASS GIS (or any other GIS) can't re-project such data to another
projection.

some work arounds:

- ask the data provider for data _with_ projection informations
- have a look at e.g. [1], assign projection information to your raster by
e.g. gdal_translate 

HTH

[1] https://lists.osgeo.org/pipermail/gdal-dev/2009-November/022607.html



-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/UTM-33x-tp5276714p5276752.html
Sent from the Grass - Users mailing list archive at Nabble.com.


More information about the grass-user mailing list