I have this problem accessing to the coordinate reference system(crs) of a raster layer map.
I wonder if such thing is possible?
I have used this xml from GDAL open map service to load in a raster layer(crs id 900913) to my stand-alone pyqgis application (http://www.gdal.org/frmt_wms_googlemaps_tms.xml)
It displays the map just fine, but when I tried to draw an icon at a specific lat long position (crs id 4326), the icon appeared weirdly.
I used the QgsCoordinateTransform class to convert the coordinates.
Please advise me what I should do to get it right? I am very new to pyqgis.
Thanks a ton!
[GIS] How to get CRS of a raster layer in pyqgis
coordinate systempyqgis
Best Answer
I am not entirely sure if i understood your question right, but you flagged it with pyqgis and it is quite simple to get the CRS here.
Using gdal it would look similar to this