Map Reprojection: Difference between revisions

From GRASS-Wiki
Jump to navigation Jump to search
No edit summary
m (Retouched, re-wording)
Line 1: Line 1:
'''Q:''' How to change map/image projections, datums, etc in GRASS GIS?
'''Question'''


'''A:'''
How to change map/image projections, datums, etc in GRASS GIS?
 
'''Answer'''


# create an import location in the projection and datum of the map/image you want to import
# create an import location in the projection and datum of the map/image you want to import
Line 8: Line 10:
# working in the destination location, use r.proj or v.proj (depending on whether the map/image is raster or vector) to reproject the map/image from the import location to the destination location.
# working in the destination location, use r.proj or v.proj (depending on whether the map/image is raster or vector) to reproject the map/image from the import location to the destination location.


If your raster data set is too big (i.e, 2 gb or so), r.proj will crash. Instead, export as a geotif, use gdalwarp to change the project and then import the new geotif file. If you have installed gdal, you almost certainly have gdalwarp as well. Directions for gdalwarp are at http://www.remotesensing.org/gdal/gdalwarp.html. The options for gdalwarp are a bit confusing. Here is an example of projecting a geotif based on srtm elevation data to utm 37N. -t_srs is the output file projection. AfricaHornElev.tif is the input file, AfricaHornElev37n.tif is the output file.
If your raster data set is too big (i.e., 2GB or so), {{cmd|r.proj}} will crash. Instead, export as a GeoTIFF, use <code>gdalwarp</code> to change the project and then import the new GeoTIFF file. If you have installed <code>gdal</code>, you almost certainly have <code>gdalwarp</code> as well. Directions for <code>gdalwarp</code> are at http://www.remotesensing.org/gdal/gdalwarp.html. The options for <code>gdalwarp</code> are a bit confusing. In the following example of projecting a GeoTIFF based on SRTM elevation data to UTM 37N, <code>-t_srs</code> is the output file projection, <code>AfricaHornElev.tif</code> is the input file and <code>AfricaHornElev37n.tif</code> is the output file.
 
<source lang="bash">gdalwarp -t_srs '+proj=utm +zone=37 +ellps=WGS84 +datum=WGS84 +units=m +no_defs' -rcs AfricaHornElev.tif AfricaHornElev37n.tif</source>


gdalwarp -t_srs '+proj=utm +zone=37 +ellps=WGS84 +datum=WGS84 +units=m +no_defs' -rcs AfricaHornElev.tif AfricaHornElev37n.tif
Another (similar) strategy frequently needed when maps need to be imported using, e.g. {{cmd|v.in.ogr}} or {{cmd|r.in.ogr}}, is to tell these programs to create a new location (argument "location"). If the imported data contains proper <code>PROJ_INFO</code> files, they can be reprojected into the current mapset, using the the projection in the current location/mapset.


Another (similar) strategy frequently needed when maps need to be imported using, e.g. v.in.ogr or r.in.ogr, is to tell these programs to create a new location (argument "location"). If the imported data contains proper PROJ_INFO files, they can be reprojected into the current mapset, using the the projection in the current location/mapset.
[[Category:FAQ]]
[[Category:FAQ]]

Revision as of 15:13, 5 February 2013

Question

How to change map/image projections, datums, etc in GRASS GIS?

Answer

  1. create an import location in the projection and datum of the map/image you want to import
  2. import the map/image into the location
  3. create a destination location in the projection and datum you want to reproject the map/image into
  4. working in the destination location, use r.proj or v.proj (depending on whether the map/image is raster or vector) to reproject the map/image from the import location to the destination location.

If your raster data set is too big (i.e., 2GB or so), r.proj will crash. Instead, export as a GeoTIFF, use gdalwarp to change the project and then import the new GeoTIFF file. If you have installed gdal, you almost certainly have gdalwarp as well. Directions for gdalwarp are at http://www.remotesensing.org/gdal/gdalwarp.html. The options for gdalwarp are a bit confusing. In the following example of projecting a GeoTIFF based on SRTM elevation data to UTM 37N, -t_srs is the output file projection, AfricaHornElev.tif is the input file and AfricaHornElev37n.tif is the output file.

gdalwarp -t_srs '+proj=utm +zone=37 +ellps=WGS84 +datum=WGS84 +units=m +no_defs' -rcs AfricaHornElev.tif AfricaHornElev37n.tif

Another (similar) strategy frequently needed when maps need to be imported using, e.g. v.in.ogr or r.in.ogr, is to tell these programs to create a new location (argument "location"). If the imported data contains proper PROJ_INFO files, they can be reprojected into the current mapset, using the the projection in the current location/mapset.