I have some 10 m digital elevation model rasters in UTM, and I need to generate 30 m resolution slope rasters in Albers.

Which should I do first: reproject the rasters, or calculate slope?

In case it matters, I am planning to use gdalwarp and gdaldem to do this.

  • 1
    i would reproject first – ziggy Dec 22 '16 at 18:52
  • What is the rationale? – user1521655 Dec 22 '16 at 20:08
  • 2
    as a general rule you should always have the correct datums and coordinates systems for your layers before you proceed with any type of geoprocessing – ziggy Dec 23 '16 at 4:22
  • 3
    I posted a discussion of this issue at gis.stackexchange.com/a/40464/664. Although it focuses on global projections, the points it makes apply to computing slopes for any raster data. If you are working in a small region and use a conformal projection designed for that region (which includes the UTM coordinate system when used appropriately), then you will be perfectly fine computing the slopes without further ado. Otherwise you need to be more careful and follow the procedures I outlined. – whuber Dec 23 '16 at 23:03

Reproject should be the first operation.

This will convert the data, and probably resize the pixels, so that any geometric operations that you would then do, would have values in reference to that Coordinate System.

If you ran Slope First, those value would be in the context of the source Coordinate System, and would not match with the values in the target Coordinate System.

  • 1
    As stated in such an unqualified manner, this answer looks incorrect. It is likely that UTM already is one of the best possible coordinate systems for computing slopes, assuming this raster does not lie outside its UTM zone. If the reprojection is into some other projection--especially a nonconformal one--then anybody following your advice will be getting worse results than otherwise. – whuber Dec 23 '16 at 23:07
  • @whuber, I'm interested in hearing more of your thoughts on this. One thing I was wondering about in posing my question is what relevant data might be lost through the resampling or interpolation necessary for reprojection. Is this considered to be an issue at all? I'd be interested in an answer that addresses both this issue and the conformal/nonconformal issue you mention here. You're more than welcome to provide an answer. – user1521655 Dec 28 '16 at 20:04
  • 2
    Resampling can be an issue. A thorough answer would take a long time to formulate. Maybe it would suffice here to remark that one should aim to minimize reprojections, but that usually a reprojection into a coordinate system that is suitable for one's analysis should be preferred. For instance, if your data were originally in a noticeably non-conformal projection (which you nevertheless needed to employ for other analyses), it would be wise to reproject to a conformal system, compute slopes, and project the results back. In your present case, it looks like no reprojection is needed. – whuber Dec 28 '16 at 20:11

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

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