PIE-engine 教程 ——坡度、坡向、山阴、粗糙度计算和显示案例分析

//示例区域影像加载
//加载全国省级行政区矢量数据
var china = pie.FeatureCollection("NGCC/CHINA_PROVINCE_BOUNDARY");
//通过数据属性中的code字段过滤出福建省数据    
var fj = china.filter(pie.Filter.eq("code", "350000"));      
//加载显示福建省图层         
Map.addLayer(fj, { color: "ffff00ff", fillColor: "00000000" }, "fj");  
//加载计算用的栅格数据，选择B1波段
var image = pie.Image('user/15010217868/edu/FJDEM0519').select("B1"); 
Map.addLayer(image)
print(image)


//影像断层信息增强处理
//山体阴影提取
//调用hillShade方法计算山体阴影，传入影像、太阳高度角、太阳方位角、垂直方向因子
var imageShade = pie.Terrain.hillShade(image, 45, 315, 8);  
//计算输出山体阴影，计算结果取绝对值。 
var image_out = pie.Image(255).subtract(imageShade).abs();   
//显示加载山体阴影图层
Map.addLayer(image_out, { color: "ffff00ff", fillColor: "00000000" }, "山体阴影");  


//坡度提取
var imageSlope = pie.Terrain.slope(image).multiply(180 / 3.1415926);  
//调用slope方法计算坡度，计算结果需要弧度角度转换
Map.addLayer(imageSlope, { color: "ffff00ff", fillColor: "00000000" }, "坡度");  
//显示加载坡度图层
//粗糙度计算提取 M=(AC*AB)/(AC*AC)=1/cos a
var num1 = pie.Image(1);    //构建一个像元值为1的影像
var imagecucaodu = num1.divide((imageSlope).cos());     

//计算粗糙度
var visParam = {  //构建显示样式style的json对象
    min: -150,
    max: 150,
    opacity: 3,
    palette: 'FF0000,00FF00,FFFF00'
};
Map.addLayer(imagecucaodu, visParam, "粗糙度");   
//显示加载粗糙度图层


//坡向提取
Map.setCenter(116, 24, 4);    //设置地图显示中心点以及显示级别
var visParam = {   //构建显示样式style的json对象
    min: 0,
    max: 360,
    opacity: 60,
    palette: 'FF0000,00FFFF,00FF00,FF00FF,0000FF,FFFF00,000FFF'
};
var imageAspect = pie.Terrain.aspect(image);    

//调用aspect方法计算坡向
var image_out2 = pie.Image(255).subtract(imageAspect).abs();    

//计算输出坡向，计算结果取绝对值
Map.addLayer(image_out2, visParam, "坡向");     

//显示加载坡向图层
//断层提取实例展示
var featureCollection0 = pie.FeatureCollection('user/15010217868/edu/fujian');
Map.addLayer(featureCollection0, { color: "ff000000", fillColor: "00000000" }, "断层");

// 相关图例设置 
var data = {
    title: "坡向",
    colors: ['#FF0000', '#00FFFF', '#00FF00', '#FF00FF', '#0000FF', '#FFFF00'],
    labels: ["60°", "300°"],
    step: 1
};

//图例样式位置的设定
var style = {
    top: "55%",
    left: "30%",
    width: "300px",
    height: "70px"
};

var legend = ui.Legend(data, style);    

//生成图例对象，传入图例的组成样式数据、在地图上的位置数据，数据为对象。

//加载显示图例，传入图例对象
Map.addUI(legend);                      

var data1 = {
    title: "坡度",
    colors: ['#000000', '#FFFFFF'],
    labels: ["0°", "90°"],
    step: 2
};
//加载第二个图例的位置和大小
var style1 = {
    top: "25%",
    left: "30%",
    width: "300px",
    height: "70px"
};
//生成图例对象，传入图例的组成样式数据、在地图上的位置数据，数据为对象。
var legend1 = ui.Legend(data1, style1);   
//加载显示图例，传入图例对象
Map.addUI(legend1);