|
Hi Ken On 07/17/2016 02:06 PM, Ken Mankoff
wrote:
Well, shame on me for missing that.Hi Micha, Sorry if my terminology was incorrect and confused the issue.On 2016-07-17 at 05:45, Micha Silver <[email protected]> wrote:First, I have never heard of UTM 33X. Normally UTM zones are either north or south. So you might have UTM 33N, or UTM 33S. All UTM north zones that are based on the WGS84 datum are numbered: 326zz, where zz is the zone. So if your data is in UTM 33N, the EPSG code will be 32633. All UTM south projections (again based on WGS84) are numbered 327zz, with the zz=zone again.33X exists. It covers Svalbard. See https://en.wikipedia.org/wiki/Universal_Transverse_Mercator_coordinate_system#Exceptions I like to take part in these mail lists because I learn something new all the time! Helmut already pointed to trying gdalinfo on the file.grid. If that does not return a full projection description then either get the full projection info from the data provider, or try to manually setup a GRASS location with the parameters for UTM 33X, then import the original into that correctly defined location with the -o option to r.in.gdal . At that point you will be able to re-project to other projections. And if I'm in another projection, and I try to import with =r.proj=, it does not work:The r.proj module does not import. Rather it does reprojection between one CRS and another.Yes, I mean reproject, although I think of it as importing from an existing location to a new location.So: I think you should setup a GRASS LOCATION and MAPSET (using the wizard) based on the EPSG code 32633 (or 32733 if your data is UTM 33S), then, within a GRASS session in that LOCATION just do: r.in.gdal -c input=file.grd output=fileI've let grass define the location using the "-c" flag to the grass command, then imported it with r.in.gdal. But because it is 33X, the location isn't fully defined, and I cannot reproject. If I follow your advice, it does not work: $ grass70 -c epsg:32633 tmp GRASS 7.0.3 (tmp): > r.in.gdal input=file.grid output=file ERROR: Projection of dataset does not appear to match current location. -k. |
_______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
