vector hillshade 制作

http://igis.me/vector-hillshade-zhi-zuo/
1 首先下载Dem数据，在NASA SRTM DEM下载 http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp
2 使用GDAL工具将数据投影到3857

gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 -r bilinear srtm.tif srtm-3857.tif

3 生成hillshade

gdaldem hillshade -z 5 srtm-3857.tif srtm-3857-hs.tif

4 在ArcGIS中打开srtm-3857-hs.tif，并使用remap/reclass工具进行重分类，一般分两类：shadow和hightlight
5 在ArcGIS中使用RasterToPolygon工具将重分类结果转成多边形srtm-3857-hs-polygon
6 使用行政区划边界裁剪结果srtm-3857-hs-polygon_clip
7 使用属性工具选择面积小于一定数量的多边形并删除，最后得到结果
另外，还可以使用下面命令生成等高线：

gdal_contour -a elev -i 20 srtm-3857.tif $srtm-3857-contour.shp

放一张效果图


查看原文：http://www.giser.net/?p=1281