GDAL山体阴影伪影

我正在使用gdal创建不同类型的图层，如彩色浮雕和山体阴影，并使用Mapnik将它们组合成一幅图像，用作从单个.hgt文件获得的3D dem模型的纹理。
假设我是gdal新手，我在hillshade层上遇到了一些问题。
使用gdal命令：

gdal\u translate N44E007.hgt N44E007.tif

我得到了N44E007.tif文件，它在IrfanView中看起来像这样

使用以下gdal命令：

gdaldem hillshade-of PNG。\tif\N44E007.tif。\hillshade\u PNG\N44E007\u hillshade.PNG

我得到的N44E007_hillshade.png文件如下所示

如何防止gdal在hillshade.png中创建这些工件？
我正在使用Windows 7和cmd。
更新1

这是我用

gdalwarp-t\srs EPSG:32632-r双线性N44E006.hgt N44E006.tif替换
gdal\u translate N44E007.hgt N44E007.tif
得到的图像



问题是重新投影的图像稍微旋转和拉伸。如何获得一个正方形的直线图像作为3d平面的纹理？
第一次尝试中使用人工制品的原因是光栅水平距离单位以度为单位，垂直距离单位以米为单位。您可以使用标准化水平和垂直距离单位，例如，尝试：

gdaldem hillshade -s 111120 -compute_edges -of PNG N44E007.hgt N44E007_hs.png


第二次尝试（更新1）重新投影到，这是一个以9°E子午线为中心的横向墨卡托投影，接近以7.5°E为中心的SRTM光栅。由于两条子午线不同，预计光栅将旋转。由于南北方向和东西方向的真实度数距离不相等（赤道除外），因此它会被拉伸。
第一次尝试使用人工制品的原因是光栅水平距离单位以度数为单位，垂直距离单位以米为单位。您可以使用标准化水平和垂直距离单位，例如，尝试：

gdaldem hillshade -s 111120 -compute_edges -of PNG N44E007.hgt N44E007_hs.png


第二次尝试（更新1）重新投影到，这是一个以9°E子午线为中心的横向墨卡托投影，接近以7.5°E为中心的SRTM光栅。由于两条子午线不同，预计光栅将旋转。而且，它会被拉伸，因为除赤道外，N-S和E-W方向的真实度数距离并不相等。
这些瑕疵可能是您输入文件的结果，因为它们显然只出现在山谷/平坦区域。也许你可以和我们共享.hgt文件。此外，这个问题可能更适合在。这是文件：。我试着把它发到GIS StackExchange上，但我还没有5点声誉。更糟糕的是，我用
gdalwarp-t\u srs EPSG:32632-r双线性N44E007.hgt N44E007.tif替换了
gdal\u translate N44E007.hgt N44E007.tif
，现在它看起来还不错，没有瑕疵。问题是，我需要一个“直”版本来用作平面纹理，EPSG重投影不会给我一个正方形，而是更像一个梯形。这些伪影可能是您输入文件的结果，因为它们显然只出现在山谷/平坦区域。也许你可以和我们共享.hgt文件。此外，这个问题可能更适合在。这是文件：。我试着把它发到GIS StackExchange上，但我还没有5点声誉。更糟糕的是，我用
gdalwarp-t\u srs EPSG:32632-r双线性N44E007.hgt N44E007.tif替换了
gdal\u translate N44E007.hgt N44E007.tif
，现在它看起来还不错，没有瑕疵。问题是，我需要一个“直”版本来用作平面纹理，EPSG重投影不会给我一个正方形，而是更像一个梯形。我知道了，非常感谢！是否有可能知道您是如何计算gdaldem命令的比例因子的？它实际上在我链接的帮助页面上，但它是从赤道周长40075016.68 m/360°=111319.49 m/°推导出来的。有趣的是，这和推荐的很接近，但不同，我不知道为什么。@MikeT，因为地球是椭球体，你的公式是一个完美的球体。111319.49可能会工作得很好，因为NA/Europe大约在曲面的一半，处于椭圆的两个极端之间。我知道了，非常感谢！是否有可能知道您是如何计算gdaldem命令的比例因子的？它实际上在我链接的帮助页面上，但它是从赤道周长40075016.68 m/360°=111319.49 m/°推导出来的。有趣的是，这和推荐的很接近，但不同，我不知道为什么。@MikeT，因为地球是椭球体，你的公式是一个完美的球体。111319.49可能会运行良好，因为NA/欧洲大约位于表面的一半，处于椭圆极限之间。