1

SRTM data is given in Lat/Lon coordinates, while Landsat is given in UTM coordinates. Is there a way to reproject the SRTM data into the Landsat CRS in GDAL (Python)?

Improve this question
3

1 Answer 1

Reset to default
2

The tool in gdal that can reproject your data is gdalwarp.

With gdalwarp, you must provide the source and target coordinate systems with -s_srs and -t_srs, respectively. In your case, it would be

gdalwarp -s_srs EPSG:4326 -t_srs EPSG 326xx your_source _file_name your_dest_file_name, where xx are the numbers of the UTM zone.

Note that you can get the coordinate system automatically if your have images everywhere in the world (see, e.g., GDAL/Python: How do I get coordinate system name from SpatialReference?)

Improve this answer
1
  • The code reported in the link you provided does not work in my case (it seems that srs.IsProjected() is returning False). Is there any way to automate this within a Python script?
    – fmonegaglia
    Commented May 3, 2017 at 7:38

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.