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

[Michael Barton]

On May 25, 2010, at 1:29 AM, grass-user-request at lists.osgeo.org wrote:

> Date: Tue, 25 May 2010 08:50:10 +0100
> From: Kim Besson <kimbesson1981 at gmail.com>
> Subject: [GRASS-user] Get current location projection in a Python
> 	Script
> To: grass-user at lists.osgeo.org
> Message-ID:
> 	<AANLkTinVGIHZ3X8oStE4MWkdynHjPZyAo_Ve4KcvZs2U at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
> 
> Greetings
> 
> 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?
> 
> Thanks
> 
> Best regards,
> Kim

Use the region() command from the GRASS Python library. <http://download.osgeo.org/grass/grass6_progman/pythonlib.html>

Michael