<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7654.12">
<TITLE>get projection zone and datum using gdal functions</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->
<P><FONT SIZE=2>Hi list,<BR>
<BR>
I use gdal2tiles.py to generate KML for Google Earth.<BR>
<BR>
gdal2tiles.py requires gdal_translate and gdalwarp to convert images to EPSG:4326<BR>
If an optical image is in UTM, then "gdalwarp -s_srs '+proj=utm +zone=## ...'<BR>
<BR>
Has anybody extracted projection info using Gdal Python function? I can "img_ds.GetProjection()" and then parse the projection string, but this is tedious.<BR>
Are there Python functions I can call to get the UTM zone and datum?<BR>
<BR>
Thanks,<BR>
Shawn<BR>
<BR>
<BR>
</FONT>
</P>
</BODY>
</HTML>