mirror of
https://github.com/jiawanlong/Cesium-Examples.git
synced 2025-07-04 15:17:36 +00:00
坡度坡向分析；新增捐赠
This commit is contained in:
jiawanlong
parent
73a886eefd
commit
11fd999a4a
23
README.md
23
README.md
@ -56,17 +56,18 @@
|
|||||||
|
|
||||||
## 打赏列表:
|
## 打赏列表:
|
||||||
|
|
||||||
| 捐赠人 | 额度 | 时间 |
|
| 捐赠人 | 额度 | 时间 |
|
||||||
| ---------- | ------- | -----------|
|
| ---------- | ------- | -----------|
|
||||||
| 不浪 | 88 | 2025.02.17 |
|
| 不浪 | 88 | 2025.02.17 |
|
||||||
| 一心一意 | 5 | 2025.03.07 |
|
| 一心一意 | 5 | 2025.03.07 |
|
||||||
| HAN_韩 | 100 | 2025.03.13 |
|
| HAN_韩 | 100 | 2025.03.13 |
|
||||||
| 阿强同学 | 5 | 2025.03.22 |
|
| 阿强同学 | 5 | 2025.03.22 |
|
||||||
| 🎾 🎾 🎾 | 20 | 2025.03.24 |
|
| 🎾 🎾 🎾 | 20 | 2025.03.24 |
|
||||||
| 雫咲poi | 20 | 2025.03.24 |
|
| 雫咲poi | 20 | 2025.03.24 |
|
||||||
| 着急 | 20 | 2025.03.25 |
|
| 着急 | 20 | 2025.03.25 |
|
||||||
| 孟庆祥 | 100 | 2025.03.25 |
|
| 孟庆祥 | 100 | 2025.03.25 |
|
||||||
| 熊重钢 | 20 | 2025.04.09 |
|
| 熊重钢 | 20 | 2025.04.09 |
|
||||||
|
| 一直游到海水变蓝 | 66 | 2025.04.09 |
|
||||||
|
|
||||||
### 关于打赏
|
### 关于打赏
|
||||||
- 您的每一笔打赏，我都会原封不动以**您的名义**捐赠到腾讯公益。
|
- 您的每一笔打赏，我都会原封不动以**您的名义**捐赠到腾讯公益。
|
||||||
|
|||||||
BIN
donate/159861999120250409EWDDWQPYJ8.jpg
Normal file
BIN
donate/159861999120250409EWDDWQPYJ8.jpg
Normal file
Binary file not shown.
|
After Width: | Height: | Size: 399 KiB |
@ -24,6 +24,7 @@
|
|||||||
<span>【着急】同学选择不公示图片</span>
|
<span>【着急】同学选择不公示图片</span>
|
||||||
<img src="./159297686120250326APBPJIJFRW.jpg">
|
<img src="./159297686120250326APBPJIJFRW.jpg">
|
||||||
<img src="./159297686120250409EWI1O7V3LF.jpg">
|
<img src="./159297686120250409EWI1O7V3LF.jpg">
|
||||||
|
<img src="./159861999120250409EWDDWQPYJ8.jpg">
|
||||||
</body>
|
</body>
|
||||||
|
|
||||||
</html>
|
</html>
|
||||||
BIN
donate/list1.jpg
BIN
donate/list1.jpg
Binary file not shown.
|
Before Width: | Height: | Size: 287 KiB After Width: | Height: | Size: 289 KiB |
217
examples/cesiumEx/4.1.12、坡度分析.html
Normal file
217
examples/cesiumEx/4.1.12、坡度分析.html
Normal file
@ -0,0 +1,217 @@
|
|||||||
|
<!--********************************************************************
|
||||||
|
* by jiawanlong
|
||||||
|
*********************************************************************-->
|
||||||
|
<!DOCTYPE html>
|
||||||
|
<html>
|
||||||
|
|
||||||
|
<head>
|
||||||
|
<meta charset="UTF-8" />
|
||||||
|
<link rel="stylesheet" href="./../../libs/cesium/Cesium1.98/Widgets/widgets.css" />
|
||||||
|
<script type="text/javascript" src="./../../libs/cesium/Cesium1.98/Cesium.js"></script>
|
||||||
|
<script src="https://unpkg.com/axios/dist/axios.min.js"></script>
|
||||||
|
<script src="./turf.min.js"></script>
|
||||||
|
</head>
|
||||||
|
|
||||||
|
<body style="
|
||||||
|
margin: 0;
|
||||||
|
overflow: hidden;
|
||||||
|
">
|
||||||
|
<div id="map" style="margin: 0 auto; width: 100%; height: 100vh"></div>
|
||||||
|
|
||||||
|
<script type="text/javascript">
|
||||||
|
Cesium.Ion.defaultAccessToken =
|
||||||
|
"eyJhbGciOiJIUzI1NiIsInR5cCI6IkpXVCJ9.eyJqdGkiOiI3ZjQ5ZGUzNC1jNWYwLTQ1ZTMtYmNjYS05YTY4ZTVmN2I2MDkiLCJpZCI6MTE3MTM4LCJpYXQiOjE2NzY0NDUyODB9.ZaNSBIfc1sGLhQd_xqhiSsc0yr8oS0wt1hAo9gbke6M";
|
||||||
|
const viewer = new Cesium.Viewer("map", {});
|
||||||
|
viewer.scene.debugShowFramesPerSecond = true;
|
||||||
|
viewer.scene.globe.depthTestAgainstTerrain = true;
|
||||||
|
|
||||||
|
viewer.terrainProvider = new Cesium.CesiumTerrainProvider({
|
||||||
|
url: 'http://data.mars3d.cn/terrain'
|
||||||
|
});
|
||||||
|
|
||||||
|
|
||||||
|
var polygonPos;
|
||||||
|
var slopes;
|
||||||
|
const delta = 0.002;
|
||||||
|
var width;
|
||||||
|
var height;
|
||||||
|
function addSlopeCanvas(extent) {
|
||||||
|
polygonPos = turf.getCoord(extent);
|
||||||
|
const rectangle = Cesium.Rectangle.fromDegrees(...polygonPos);
|
||||||
|
const positions = [];
|
||||||
|
width = Math.ceil((polygonPos[2] - polygonPos[0]) / delta);
|
||||||
|
height = Math.ceil((polygonPos[3] - polygonPos[1]) / delta);
|
||||||
|
radio = width / height;
|
||||||
|
const ddx =
|
||||||
|
Cesium.Cartesian3.distance(
|
||||||
|
Cesium.Cartesian3.fromDegrees(polygonPos[0], polygonPos[1]),
|
||||||
|
Cesium.Cartesian3.fromDegrees(polygonPos[2], polygonPos[1])
|
||||||
|
) / width;
|
||||||
|
const ddy =
|
||||||
|
Cesium.Cartesian3.distance(
|
||||||
|
Cesium.Cartesian3.fromDegrees(polygonPos[0], polygonPos[1]),
|
||||||
|
Cesium.Cartesian3.fromDegrees(polygonPos[0], polygonPos[3])
|
||||||
|
) / height;
|
||||||
|
for (let y = 0; y < height; y++) {
|
||||||
|
for (let x = 0; x < width; x++) {
|
||||||
|
const longitude = Cesium.Math.toDegrees(
|
||||||
|
Cesium.Math.lerp(rectangle.west, rectangle.east, x / (width - 1))
|
||||||
|
);
|
||||||
|
const latitude = Cesium.Math.toDegrees(
|
||||||
|
Cesium.Math.lerp(rectangle.north, rectangle.south, y / (height - 1))
|
||||||
|
);
|
||||||
|
positions.push(
|
||||||
|
Cesium.Cartographic.fromDegrees(longitude - delta / 2, latitude)
|
||||||
|
);
|
||||||
|
positions.push(
|
||||||
|
Cesium.Cartographic.fromDegrees(
|
||||||
|
longitude - delta / 2,
|
||||||
|
latitude - delta / 2
|
||||||
|
)
|
||||||
|
);
|
||||||
|
positions.push(
|
||||||
|
Cesium.Cartographic.fromDegrees(longitude, latitude + delta / 2)
|
||||||
|
);
|
||||||
|
positions.push(
|
||||||
|
Cesium.Cartographic.fromDegrees(
|
||||||
|
longitude - delta / 2,
|
||||||
|
latitude + delta / 2
|
||||||
|
)
|
||||||
|
);
|
||||||
|
positions.push(
|
||||||
|
Cesium.Cartographic.fromDegrees(longitude + delta / 2, latitude)
|
||||||
|
);
|
||||||
|
positions.push(
|
||||||
|
Cesium.Cartographic.fromDegrees(
|
||||||
|
longitude + delta / 2,
|
||||||
|
latitude + delta / 2
|
||||||
|
)
|
||||||
|
);
|
||||||
|
positions.push(
|
||||||
|
Cesium.Cartographic.fromDegrees(longitude, latitude - delta / 2)
|
||||||
|
);
|
||||||
|
positions.push(
|
||||||
|
Cesium.Cartographic.fromDegrees(
|
||||||
|
longitude + delta / 2,
|
||||||
|
latitude - delta / 2
|
||||||
|
)
|
||||||
|
);
|
||||||
|
positions.push(Cesium.Cartographic.fromDegrees(longitude, latitude)); // mid
|
||||||
|
}
|
||||||
|
}
|
||||||
|
slopes = [];
|
||||||
|
Cesium.sampleTerrainMostDetailed(viewer.terrainProvider, positions).then(
|
||||||
|
(updatedPositions) => {
|
||||||
|
for (let i = 0; i < updatedPositions.length; i += 9) {
|
||||||
|
let center = updatedPositions[i + 8];
|
||||||
|
const westHeight = updatedPositions[i + 0].height;
|
||||||
|
let maxHeight = westHeight;
|
||||||
|
let maxIndex = i + 0;
|
||||||
|
const westSouthHeight = updatedPositions[i + 1].height;
|
||||||
|
if (westSouthHeight > maxHeight) {
|
||||||
|
maxHeight = westSouthHeight;
|
||||||
|
maxIndex = i + 1;
|
||||||
|
}
|
||||||
|
const northHeight = updatedPositions[i + 2].height;
|
||||||
|
if (northHeight > maxHeight) {
|
||||||
|
maxHeight = northHeight;
|
||||||
|
maxIndex = i + 2;
|
||||||
|
}
|
||||||
|
const westNorthHeight = updatedPositions[i + 3].height;
|
||||||
|
if (westNorthHeight > maxHeight) {
|
||||||
|
maxHeight = northHeight;
|
||||||
|
maxIndex = i + 3;
|
||||||
|
}
|
||||||
|
const eastHeight = updatedPositions[i + 4].height;
|
||||||
|
if (eastHeight > maxHeight) {
|
||||||
|
maxHeight = eastHeight;
|
||||||
|
maxIndex = i + 4;
|
||||||
|
}
|
||||||
|
const eastNorthHeight = updatedPositions[i + 5].height;
|
||||||
|
if (eastNorthHeight > maxHeight) {
|
||||||
|
maxHeight = eastNorthHeight;
|
||||||
|
maxIndex = i + 5;
|
||||||
|
}
|
||||||
|
const southHeight = updatedPositions[i + 6].height;
|
||||||
|
if (southHeight > maxHeight) {
|
||||||
|
maxHeight = southHeight;
|
||||||
|
maxIndex = i + 6;
|
||||||
|
}
|
||||||
|
const eastSouthHeight = updatedPositions[i + 7].height;
|
||||||
|
if (eastSouthHeight > maxHeight) {
|
||||||
|
maxHeight = eastSouthHeight;
|
||||||
|
maxIndex = i + 7;
|
||||||
|
}
|
||||||
|
const fx =
|
||||||
|
(westSouthHeight +
|
||||||
|
2 * southHeight +
|
||||||
|
eastSouthHeight -
|
||||||
|
(westNorthHeight + 2 * northHeight + eastNorthHeight)) /
|
||||||
|
(8 * ddx);
|
||||||
|
const fy =
|
||||||
|
(eastNorthHeight +
|
||||||
|
2 * eastHeight +
|
||||||
|
eastSouthHeight -
|
||||||
|
(westNorthHeight + 2 * westHeight + westSouthHeight)) /
|
||||||
|
(8 * ddy);
|
||||||
|
const slope = Math.atan(Math.sqrt(fx ** 2 + fy ** 2));
|
||||||
|
slopes.push(slope);
|
||||||
|
}
|
||||||
|
createMaterial();
|
||||||
|
}
|
||||||
|
);
|
||||||
|
}
|
||||||
|
function createMaterial() {
|
||||||
|
const canvas = toCanvas();
|
||||||
|
const positions = [
|
||||||
|
Cesium.Cartesian3.fromDegrees(polygonPos[0], polygonPos[1]),
|
||||||
|
Cesium.Cartesian3.fromDegrees(polygonPos[0], polygonPos[3]),
|
||||||
|
Cesium.Cartesian3.fromDegrees(polygonPos[2], polygonPos[3]),
|
||||||
|
Cesium.Cartesian3.fromDegrees(polygonPos[2], polygonPos[1]),
|
||||||
|
];
|
||||||
|
const polygon = viewer.entities.add({
|
||||||
|
polygon: {
|
||||||
|
hierarchy: new Cesium.PolygonHierarchy(positions),
|
||||||
|
material: new Cesium.ImageMaterialProperty({
|
||||||
|
image: canvas,
|
||||||
|
}),
|
||||||
|
classificationType: Cesium.ClassificationType.BOTH,
|
||||||
|
},
|
||||||
|
});
|
||||||
|
}
|
||||||
|
function toCanvas() {
|
||||||
|
const w = width;
|
||||||
|
const h = height;
|
||||||
|
const canvas = document.createElement("canvas");
|
||||||
|
canvas.height = h;
|
||||||
|
canvas.width = w;
|
||||||
|
const ctx = canvas.getContext("2d");
|
||||||
|
const bitmap = new Uint8ClampedArray(w * h * 4); // 各像素点rgba
|
||||||
|
console.log(slopes);
|
||||||
|
for (let y = 0; y < h; y++) {
|
||||||
|
for (let x = 0; x < w; x++) {
|
||||||
|
const slope = slopes[y * w + x];
|
||||||
|
const hue = (1 - slope / (Math.PI / 2)) / 2;
|
||||||
|
const color = Cesium.Color.fromHsl(hue, 1, 0.6).withAlpha(0.7);
|
||||||
|
const bitmapIndex = (y * w + x) * 4;
|
||||||
|
bitmap[bitmapIndex + 0] = color.red * 255;
|
||||||
|
bitmap[bitmapIndex + 1] = color.green * 255;
|
||||||
|
bitmap[bitmapIndex + 2] = color.blue * 255;
|
||||||
|
bitmap[bitmapIndex + 3] = color.alpha * 255;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
const imageData = new ImageData(bitmap, w, h);
|
||||||
|
ctx.putImageData(imageData, 0, 0);
|
||||||
|
return canvas;
|
||||||
|
}
|
||||||
|
|
||||||
|
addSlopeCanvas([78.64, 33.72, 78.79, 34.02]);
|
||||||
|
|
||||||
|
viewer.camera.setView({
|
||||||
|
destination: Cesium.Cartesian3.fromDegrees(78.7, 33.85, 30000),
|
||||||
|
});
|
||||||
|
|
||||||
|
</script>
|
||||||
|
</body>
|
||||||
|
|
||||||
|
</html>
|
||||||
292
examples/cesiumEx/4.1.13、坡向分析.html
Normal file
292
examples/cesiumEx/4.1.13、坡向分析.html
Normal file
@ -0,0 +1,292 @@
|
|||||||
|
<!--********************************************************************
|
||||||
|
* by jiawanlong
|
||||||
|
*********************************************************************-->
|
||||||
|
<!DOCTYPE html>
|
||||||
|
<html>
|
||||||
|
|
||||||
|
<head>
|
||||||
|
<meta charset="UTF-8" />
|
||||||
|
<link rel="stylesheet" href="./../../libs/cesium/Cesium1.98/Widgets/widgets.css" />
|
||||||
|
<script type="text/javascript" src="./../../libs/cesium/Cesium1.98/Cesium.js"></script>
|
||||||
|
<script src="https://unpkg.com/axios/dist/axios.min.js"></script>
|
||||||
|
<script src="./turf.min.js"></script>
|
||||||
|
</head>
|
||||||
|
|
||||||
|
<body style="
|
||||||
|
margin: 0;
|
||||||
|
overflow: hidden;
|
||||||
|
">
|
||||||
|
<div id="map" style="margin: 0 auto; width: 100%; height: 100vh"></div>
|
||||||
|
|
||||||
|
<script type="text/javascript">
|
||||||
|
Cesium.Ion.defaultAccessToken =
|
||||||
|
"eyJhbGciOiJIUzI1NiIsInR5cCI6IkpXVCJ9.eyJqdGkiOiI3ZjQ5ZGUzNC1jNWYwLTQ1ZTMtYmNjYS05YTY4ZTVmN2I2MDkiLCJpZCI6MTE3MTM4LCJpYXQiOjE2NzY0NDUyODB9.ZaNSBIfc1sGLhQd_xqhiSsc0yr8oS0wt1hAo9gbke6M";
|
||||||
|
const viewer = new Cesium.Viewer("map", {});
|
||||||
|
viewer.scene.debugShowFramesPerSecond = true;
|
||||||
|
viewer.scene.globe.depthTestAgainstTerrain = true;
|
||||||
|
|
||||||
|
viewer.terrainProvider = new Cesium.CesiumTerrainProvider({
|
||||||
|
url: 'http://data.mars3d.cn/terrain'
|
||||||
|
});
|
||||||
|
|
||||||
|
|
||||||
|
//坡向分析
|
||||||
|
class falling {
|
||||||
|
constructor(viewer, config) {
|
||||||
|
this.viewer = viewer;
|
||||||
|
this.config = config;
|
||||||
|
this.loadResult();
|
||||||
|
}
|
||||||
|
|
||||||
|
loadResult() {
|
||||||
|
var _this = this;
|
||||||
|
|
||||||
|
var startx1 = _this.config.positons[0];
|
||||||
|
var starty1 = _this.config.positons[1];
|
||||||
|
var startx2 = _this.config.positons[2];
|
||||||
|
var starty2 = _this.config.positons[3];
|
||||||
|
var startx3 = _this.config.positons[4];
|
||||||
|
var starty3 = _this.config.positons[5];
|
||||||
|
var startx4 = _this.config.positons[6];
|
||||||
|
var starty4 = _this.config.positons[7];
|
||||||
|
|
||||||
|
//坡向显示
|
||||||
|
function analysis(startx1, starty1, startz1, startx2, starty2, startz2, startx3, starty3, startz3, startx4, starty4, startz4) {
|
||||||
|
//高度z全为0
|
||||||
|
var count = 100;
|
||||||
|
var slopelineposition = [];
|
||||||
|
var hireacys = [];
|
||||||
|
var hireacysdistance = []
|
||||||
|
for (let j = 0; j < 100; j++) {
|
||||||
|
for (var i = 0; i < 100; i++) {
|
||||||
|
var hireacy = [];
|
||||||
|
//分割成小面，切分经纬度
|
||||||
|
hireacy.push(new Cesium.Cartesian3(startx1 + (startx2 - startx1) / count * i + (startx4 + (startx3 - startx4) / count * (i) - startx1 - (startx2 - startx1) / count * i) / count * j,
|
||||||
|
starty1 + (starty2 - starty1) / count * i + (starty4 + (starty3 - starty4) / count * (i) - starty1 - (starty2 - starty1) / count * i) / count * j,
|
||||||
|
startz1 + (startz2 - startz1) / count * i + (startz4 + (startz3 - startz4) / count * (i) - startz1 - (startz2 - startz1) / count * i) / count * j))
|
||||||
|
hireacy.push(new Cesium.Cartesian3(startx1 + (startx2 - startx1) / count * (i + 1) + (startx4 + (startx3 - startx4) / count * (i + 1) - startx1 - (startx2 - startx1) / count * (i + 1)) / count * j,
|
||||||
|
starty1 + (starty2 - starty1) / count * (i + 1) + (starty4 + (starty3 - starty4) / count * (i + 1) - starty1 - (starty2 - starty1) / count * (i + 1)) / count * j,
|
||||||
|
startz1 + (startz2 - startz1) / count * (i + 1) + (startz4 + (startz3 - startz4) / count * (i + 1) - startz1 - (startz2 - startz1) / count * (i + 1)) / count * j))
|
||||||
|
hireacy.push(new Cesium.Cartesian3(startx4 + (startx3 - startx4) / count * (i + 1) - (startx4 + (startx3 - startx4) / count * (i + 1) - startx1 - (startx2 - startx1) / count * (i + 1)) / count * (count - 1 - j),
|
||||||
|
starty4 + (starty3 - starty4) / count * (i + 1) - (starty4 + (starty3 - starty4) / count * (i + 1) - starty1 - (starty2 - starty1) / count * (i + 1)) / count * (count - 1 - j),
|
||||||
|
startz4 + (startz3 - startz4) / count * (i + 1) - (startz4 + (startz3 - startz4) / count * (i + 1) - startz1 - (startz2 - startz1) / count * (i + 1)) / count * (count - 1 - j)))
|
||||||
|
hireacy.push(new Cesium.Cartesian3(startx4 + (startx3 - startx4) / count * i - (startx4 + (startx3 - startx4) / count * (i) - startx1 - (startx2 - startx1) / count * i) / count * (count - 1 - j),
|
||||||
|
starty4 + (starty3 - starty4) / count * i - (starty4 + (starty3 - starty4) / count * (i) - starty1 - (starty2 - starty1) / count * i) / count * (count - 1 - j),
|
||||||
|
startz4 + (startz3 - startz4) / count * i - (startz4 + (startz3 - startz4) / count * (i) - startz1 - (startz2 - startz1) / count * i) / count * (count - 1 - j)))
|
||||||
|
hireacys.push(hireacy);
|
||||||
|
//取出面的8个点坐标，拿点坐标去求高度值
|
||||||
|
slopelineposition.push(Cesium.Cartographic.fromDegrees(hireacy[0].x, hireacy[0].y));
|
||||||
|
slopelineposition.push(Cesium.Cartographic.fromDegrees((hireacy[0].x + hireacy[1].x) / 2, (hireacy[0].y + hireacy[1].y) / 2));
|
||||||
|
slopelineposition.push(Cesium.Cartographic.fromDegrees(hireacy[1].x, hireacy[1].y));
|
||||||
|
slopelineposition.push(Cesium.Cartographic.fromDegrees((hireacy[1].x + hireacy[2].x) / 2, (hireacy[1].y + hireacy[2].y) / 2));
|
||||||
|
slopelineposition.push(Cesium.Cartographic.fromDegrees(hireacy[2].x, hireacy[2].y));
|
||||||
|
slopelineposition.push(Cesium.Cartographic.fromDegrees((hireacy[2].x + hireacy[3].x) / 2, (hireacy[2].y + hireacy[3].y) / 2));
|
||||||
|
slopelineposition.push(Cesium.Cartographic.fromDegrees(hireacy[3].x, hireacy[3].y));
|
||||||
|
slopelineposition.push(Cesium.Cartographic.fromDegrees((hireacy[3].x + hireacy[0].x) / 2, (hireacy[3].y + hireacy[0].y) / 2));
|
||||||
|
|
||||||
|
}
|
||||||
|
}
|
||||||
|
var promise = Cesium.sampleTerrainMostDetailed(_this.config.terrainProvider, slopelineposition);
|
||||||
|
promise.then(function (updatedPositions) {
|
||||||
|
|
||||||
|
//拿到所有的高度数据
|
||||||
|
var heighmm = [];
|
||||||
|
var m = 0
|
||||||
|
//计算坡度比的次数
|
||||||
|
var countcolor1 = 0;
|
||||||
|
var countcolor2 = 0;
|
||||||
|
var countcolor3 = 0;
|
||||||
|
var countcolor4 = 0;
|
||||||
|
var countcolor5 = 0;
|
||||||
|
var countcolor6 = 0;
|
||||||
|
var countcolor7 = 0;
|
||||||
|
for (var k = 0; k < updatedPositions.length / 8; k++) {
|
||||||
|
//第一个点与第五个点的坡度
|
||||||
|
var slope1 = (updatedPositions[m].height - updatedPositions[m + 4].height) / (Cesium.Cartesian3.distance(
|
||||||
|
Cesium.Cartesian3.fromDegrees(updatedPositions[m].latitude, updatedPositions[m].longitude, updatedPositions[m].height),
|
||||||
|
Cesium.Cartesian3.fromDegrees(updatedPositions[m + 4].latitude, updatedPositions[m + 4].longitude, updatedPositions[m + 4].height)))
|
||||||
|
//第二个点与第六个点的坡度
|
||||||
|
var slope2 = (updatedPositions[m + 1].height - updatedPositions[m + 5].height) / (Cesium.Cartesian3.distance(
|
||||||
|
Cesium.Cartesian3.fromDegrees(updatedPositions[m + 1].latitude, updatedPositions[m + 1].longitude, updatedPositions[m + 1].height),
|
||||||
|
Cesium.Cartesian3.fromDegrees(updatedPositions[m + 5].latitude, updatedPositions[m + 5].longitude, updatedPositions[m + 5].height)))
|
||||||
|
//第三个点与第七个点的坡度
|
||||||
|
var slope3 = (updatedPositions[m + 2].height - updatedPositions[m + 6].height) / (Cesium.Cartesian3.distance(
|
||||||
|
Cesium.Cartesian3.fromDegrees(updatedPositions[m + 2].latitude, updatedPositions[m + 2].longitude, updatedPositions[m + 2].height),
|
||||||
|
Cesium.Cartesian3.fromDegrees(updatedPositions[m + 6].latitude, updatedPositions[m + 6].longitude, updatedPositions[m + 6].height)))
|
||||||
|
//第四个点与第八个点的坡度
|
||||||
|
var slope4 = (updatedPositions[m + 3].height - updatedPositions[m + 7].height) / (Cesium.Cartesian3.distance(
|
||||||
|
Cesium.Cartesian3.fromDegrees(updatedPositions[m + 3].latitude, updatedPositions[m + 3].longitude, updatedPositions[m + 3].height),
|
||||||
|
Cesium.Cartesian3.fromDegrees(updatedPositions[m + 7].latitude, updatedPositions[m + 7].longitude, updatedPositions[m + 7].height)))
|
||||||
|
// console.log("slope1:"+slope1+";slope2:"+slope2+";slope3:"+slope3+";slope4:"+slope4);
|
||||||
|
var arrposition = [Math.abs(slope1), Math.abs(slope2), Math.abs(slope3), Math.abs(slope4)];//取绝对值
|
||||||
|
arrposition.sort();
|
||||||
|
var slope = arrposition[3]; // 拿到最大的坡度值
|
||||||
|
var lineposition = [];//画方向线的坐标
|
||||||
|
if (slope == Math.abs(slope1)) {
|
||||||
|
if (slope1 > 0) {
|
||||||
|
lineposition.push(Cesium.Math.toDegrees(updatedPositions[m].longitude), Cesium.Math.toDegrees(updatedPositions[m].latitude),
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 4].longitude), Cesium.Math.toDegrees(updatedPositions[m + 4].latitude));
|
||||||
|
} else {
|
||||||
|
lineposition.push(
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 4].longitude), Cesium.Math.toDegrees(updatedPositions[m + 4].latitude),
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m].longitude), Cesium.Math.toDegrees(updatedPositions[m].latitude));
|
||||||
|
|
||||||
|
}
|
||||||
|
} else if (slope == Math.abs(slope2)) {
|
||||||
|
if (slope2 > 0) {
|
||||||
|
lineposition.push(Cesium.Math.toDegrees(updatedPositions[m + 1].longitude), Cesium.Math.toDegrees(updatedPositions[m + 1].latitude),
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 5].longitude), Cesium.Math.toDegrees(updatedPositions[m + 5].latitude));
|
||||||
|
} else {
|
||||||
|
lineposition.push(
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 5].longitude), Cesium.Math.toDegrees(updatedPositions[m + 5].latitude),
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 1].longitude), Cesium.Math.toDegrees(updatedPositions[m + 1].latitude));
|
||||||
|
}
|
||||||
|
} else if (slope == Math.abs(slope3)) {
|
||||||
|
if (slope3 > 0) {
|
||||||
|
lineposition.push(Cesium.Math.toDegrees(updatedPositions[m + 2].longitude), Cesium.Math.toDegrees(updatedPositions[m + 2].latitude),
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 6].longitude), Cesium.Math.toDegrees(updatedPositions[m + 6].latitude));
|
||||||
|
} else {
|
||||||
|
lineposition.push(
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 6].longitude), Cesium.Math.toDegrees(updatedPositions[m + 6].latitude),
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 2].longitude), Cesium.Math.toDegrees(updatedPositions[m + 2].latitude));
|
||||||
|
}
|
||||||
|
} else if (slope == Math.abs(slope4)) {
|
||||||
|
if (slope4 > 0) {
|
||||||
|
lineposition.push(Cesium.Math.toDegrees(updatedPositions[m + 3].longitude), Cesium.Math.toDegrees(updatedPositions[m + 3].latitude),
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 7].longitude), Cesium.Math.toDegrees(updatedPositions[m + 7].latitude));
|
||||||
|
} else {
|
||||||
|
lineposition.push(
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 7].longitude), Cesium.Math.toDegrees(updatedPositions[m + 7].latitude),
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 3].longitude), Cesium.Math.toDegrees(updatedPositions[m + 3].latitude));
|
||||||
|
}
|
||||||
|
}
|
||||||
|
slope = (Math.abs(slope1) + Math.abs(slope2) + Math.abs(slope3) + Math.abs(slope4)) / 4; //四个坡度值大小有的差值特别大，这里取的平均值用来配置颜色
|
||||||
|
//console.log(slope);
|
||||||
|
let slopecolor;
|
||||||
|
if (0 <= slope && slope < 0.29) {
|
||||||
|
slopecolor = '#ff9090'
|
||||||
|
countcolor1++;
|
||||||
|
} else if (0.29 <= slope && slope < 0.5) {
|
||||||
|
slopecolor = '#ff8080'
|
||||||
|
countcolor2++;
|
||||||
|
} else if (0.5 <= slope && slope < Math.sqrt(2) / 2) {
|
||||||
|
slopecolor = '#ff7070'
|
||||||
|
countcolor3++;
|
||||||
|
} else if (Math.sqrt(2) / 2 <= slope && slope < 0.87) {
|
||||||
|
slopecolor = '#ff6060'
|
||||||
|
countcolor4++;
|
||||||
|
} else if (0.87 <= slope && slope < 0.91) {
|
||||||
|
slopecolor = '#ff5050'
|
||||||
|
countcolor5++;
|
||||||
|
} else if (0.91 <= slope && slope < 0.95) {
|
||||||
|
slopecolor = '#ff4040'
|
||||||
|
countcolor6++;
|
||||||
|
} else {
|
||||||
|
slopecolor = '#ff3030'
|
||||||
|
countcolor7++;
|
||||||
|
}
|
||||||
|
_this.viewer.entities.add({
|
||||||
|
//type: 'drawSloperectange',
|
||||||
|
rectangle: {
|
||||||
|
coordinates: Cesium.Rectangle.fromDegrees(
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m].longitude), Cesium.Math.toDegrees(updatedPositions[m].latitude),
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 4].longitude), Cesium.Math.toDegrees(updatedPositions[m + 4].latitude))
|
||||||
|
,
|
||||||
|
material: Cesium.Color.fromCssColorString(slopecolor)
|
||||||
|
|
||||||
|
},
|
||||||
|
// polyline: {
|
||||||
|
// clampToGround: true,
|
||||||
|
// positions: Cesium.Cartesian3.fromDegreesArray(lineposition),
|
||||||
|
// material: new Cesium.PolylineArrowMaterialProperty(Cesium.Color.fromCssColorString("#ff9090")),// Cesium.Color.BLUE.withAlpha(0.5) ,
|
||||||
|
// width: 8,
|
||||||
|
// distanceDisplayCondition: new Cesium.DistanceDisplayCondition(0, updatedPositions[m + 4].height + 4500)
|
||||||
|
// },
|
||||||
|
}
|
||||||
|
);
|
||||||
|
var tt = _this.viewer.entities.add({
|
||||||
|
name: "Red polygon on surface",
|
||||||
|
rectangle: {
|
||||||
|
coordinates: Cesium.Rectangle.fromDegrees(
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m].longitude), Cesium.Math.toDegrees(updatedPositions[m].latitude),
|
||||||
|
Cesium.Math.toDegrees(updatedPositions[m + 4].longitude), Cesium.Math.toDegrees(updatedPositions[m + 4].latitude))
|
||||||
|
,
|
||||||
|
material: Cesium.Color.fromCssColorString(slopecolor)
|
||||||
|
},
|
||||||
|
});
|
||||||
|
|
||||||
|
_this.viewer.zoomTo(tt);
|
||||||
|
|
||||||
|
_this.viewer.entities.add({
|
||||||
|
name: "Red polygon on surface",
|
||||||
|
polyline: {
|
||||||
|
clampToGround: true,
|
||||||
|
positions: Cesium.Cartesian3.fromDegreesArray(lineposition),
|
||||||
|
material: new Cesium.PolylineArrowMaterialProperty(Cesium.Color.fromCssColorString("#ff9090")),// Cesium.Color.BLUE.withAlpha(0.5) ,
|
||||||
|
width: 8,
|
||||||
|
distanceDisplayCondition: new Cesium.DistanceDisplayCondition(0, updatedPositions[m + 4].height + 4500)
|
||||||
|
}
|
||||||
|
});
|
||||||
|
|
||||||
|
m += 8;
|
||||||
|
}
|
||||||
|
// var contents = "<div class='slope_layer'><div> <div style='height:19px;background-color:#ff9090; width: 64px;float: left;'> </div><span >" + (countcolor1 / (updatedPositions.length / 8) * 100).toFixed(2) + "% (<0.29)</span></div>";
|
||||||
|
// contents += " <div style='clear: both;'><div style='height:19px;background-color:#ff8080 ; width: 64px;float: left;'></div><span > " + (countcolor2 / (updatedPositions.length / 8) * 100).toFixed(2) + "%(<0.5)</span></div>";
|
||||||
|
// contents += "<div style='clear: both;'><div style='height:19px;background-color:#ff7070;width: 64px;float: left;'> </div><span >" + (countcolor3 / (updatedPositions.length / 8) * 100).toFixed(2) + "%(<0.7)</span></div>";
|
||||||
|
// contents += " <div style='clear: both;'><div style='height:19px;background-color:#ff6060;width: 64px;float: left;'> </div><span >" + (countcolor4 / (updatedPositions.length / 8) * 100).toFixed(2) + "%(<0.87)</span></div>";
|
||||||
|
// contents += " <div style='clear: both;'><div style='height:19px;background-color:#ff5050;width: 64px;float: left;'> </div><span > " + (countcolor5 / (updatedPositions.length / 8) * 100).toFixed(2) + "%(<0.91)</span></div>";
|
||||||
|
// contents += " <div style='clear: both;'><div style='height:19px;background-color:#ff4040;width: 64px;float: left;'> </div><span >" + (countcolor6 / (updatedPositions.length / 8) * 100).toFixed(2) + "%(<0.95)</span></div>";
|
||||||
|
// contents += " <div style='clear: both;'><div style='height:19px;background-color:#ff3030;width: 64px;float: left;'> </div><span > " + (countcolor7 / (updatedPositions.length / 8) * 100).toFixed(2) + "%(<1)</span></div></div>";
|
||||||
|
// layer.closeAll();
|
||||||
|
|
||||||
|
|
||||||
|
// layer.open({
|
||||||
|
// type: 1,
|
||||||
|
// offset: ["500px", "1200px"],
|
||||||
|
// area: ['260px', '179px'],
|
||||||
|
// title: "坡度分析信息"
|
||||||
|
// ,
|
||||||
|
// content: contents
|
||||||
|
// ,
|
||||||
|
// btnAlign: 'c' //按钮居中
|
||||||
|
// ,
|
||||||
|
// shade: 0 //不显示遮罩
|
||||||
|
// ,
|
||||||
|
// cancel: function (index) {
|
||||||
|
// layer.close(index);
|
||||||
|
// }
|
||||||
|
// });
|
||||||
|
|
||||||
|
|
||||||
|
});
|
||||||
|
// Cesium.when(promise,
|
||||||
|
// function (updatedPositions) {})
|
||||||
|
}
|
||||||
|
|
||||||
|
analysis(startx1, starty1, 0, startx2, starty2, 0, startx3, starty3, 0, startx4, starty4, 0)
|
||||||
|
}
|
||||||
|
}
|
||||||
|
let terrain = viewer.terrainProvider;
|
||||||
|
viewer.camera.flyTo({
|
||||||
|
destination: Cesium.Cartesian3.fromDegrees(
|
||||||
|
114.20597076416016,
|
||||||
|
31.15141366186558
|
||||||
|
, 5000.0)
|
||||||
|
});
|
||||||
|
var contour = new falling(viewer, {
|
||||||
|
positons: [
|
||||||
|
114.20331001281738,
|
||||||
|
31.148034709323344,
|
||||||
|
114.21009063720703,
|
||||||
|
31.148034709323344,
|
||||||
|
114.21009063720703,
|
||||||
|
31.154204879640236,
|
||||||
|
114.20331001281738,
|
||||||
|
31.154204879640236,
|
||||||
|
],
|
||||||
|
terrainProvider: terrain
|
||||||
|
});
|
||||||
|
</script>
|
||||||
|
</body>
|
||||||
|
|
||||||
|
</html>
|
||||||
Loading…
Reference in New Issue
Block a user