[GRASS-user] Re: GRASS-user] Help with reprojection

Hamish hamish_b at yahoo.com
Sat Jan 3 16:25:31 EST 2009


Hamish wrote:
> - create a simple xy location
>     (lat/lon location will not allow north > 90, and your image while
>      still not geo-referenced will go to 8100)
> - run r.region to set n,s,e,w bounds to 90,-90,180,-180
> - run g.setproj to rejig the location into a lat/lon one.

umm, that might not work -- r.info will still know the map is XY even
if the location is changed to lat/lon.

> - check resolution is correct (nicely "0:01:20") with r.info.
> - zoom to area of interest. you probably do not want to reproject
>   entire planet. after zooming check resolution is preserved,
>   (g.region -p, maybe with "g.region res=0:01:20 -a" after zoom to fix)
>   and run "r.mapcalc cropmap=fullmap" to perform the crop.
> 
> then from the lambert location run r.proj to pull the
> cropped image across.
> 
> there are some examples of this process in the GRASS wiki,
> look at the "Global datasets" page.

the above issue should be covered there; also check the mailing list
archives.

so your "easiest" solution is to create a "world file". see the GDAL
JPG or GeoTiff format import page, or do a web search for instructions.


Hamish



      



More information about the grass-user mailing list