[GRASS-user] Get current location projection in a Python Script

[Michael Barton]

On May 25, 2010, at 8:25 AM, Nikos Alexandris wrote:

> Kim Besson:
>>> I'm building a GRASS Python Script to automatically import a set of data
>>> to my location. At GRASS book (chapter 3.3.3) it uses g.proj -wf to to
>>> check current location projection. Can I use the same command in Python
>>> to retrieve the same?
>>> 
>>> 
>>> And how about calling gdalwarp in a Python-GRASS Script? can it be called
>>> and used by a Python_GRASS Script? if so, is there any Python example
>>> with that?
>>> 
>>> Regarding "g.proj -wf" command. I get a lot of information. I suppose
>>> GDALWARP can read this right ? If not, How can I obtain current Location
>>> EPSG code?
> 
> Michael Barton: 
>> Use the region() command from the GRASS Python library.
>> <http://download.osgeo.org/grass/grass6_progman/pythonlib.html>
> 
> Just for the records, a thread not _directly_ about it but close (enough?):
> <http://osgeo-org.1803224.n2.nabble.com/Getting-rows-cols-of-a-region-in-a-
> script-tc2787474.html#a2788062>
> 
> Nikos

The thread is for BASH and does not work for Python.

Michael