|
@@ -1,397 +0,0 @@
|
|
|
-import { v4 as uuidv4 } from "uuid";
|
|
|
-import moment from "moment";
|
|
|
-import {
|
|
|
- cartesian3ToWGS84,
|
|
|
- mapQuery,
|
|
|
- flatten,
|
|
|
- mercator2lonLat,
|
|
|
- undergroundMode,
|
|
|
-} from "@/utils/MapHelper/MapHelper.js";
|
|
|
-import {
|
|
|
- area,
|
|
|
- intersect,
|
|
|
- polygon,
|
|
|
- point,
|
|
|
- midpoint,
|
|
|
- difference,
|
|
|
-} from "@turf/turf";
|
|
|
-const BenchmarkLandPrice = {
|
|
|
- /**
|
|
|
- * 基准地价分析
|
|
|
- * @param {*} BenchmarkLandPrice {
|
|
|
- ProjectScope: [],分析范围
|
|
|
- ProjectName: "",项目名称
|
|
|
- ConstructionUnit: "",建设单位
|
|
|
- imageBase64: [],输出图片
|
|
|
- }
|
|
|
- * @param {*} isOutputImg 是否输出范围图片 true false 如需输出图片,请确保要输出的范围已在屏幕正确认为
|
|
|
- * @param {*} Callback 回调
|
|
|
- */
|
|
|
- async calculateLandPrice(BenchmarkLandPrice, isOutputImg = false, Callback) {
|
|
|
- let GHDKList = await this.getGHDKs(BenchmarkLandPrice);
|
|
|
- await this.GetJZDJ(BenchmarkLandPrice, GHDKList);
|
|
|
- let data = await this.CalculateJZDJ(BenchmarkLandPrice, GHDKList);
|
|
|
- if (isOutputImg) {
|
|
|
- if (data && data.picBase64List) {
|
|
|
- var promise = scene.outputSceneToFile();
|
|
|
- Cesium.when(promise, function (base64data) {
|
|
|
- data.picBase64List.push(base64data);
|
|
|
- if (Callback) Callback(data);
|
|
|
- });
|
|
|
- } else {
|
|
|
- if (Callback) Callback(false);
|
|
|
- }
|
|
|
- } else {
|
|
|
- if (data) {
|
|
|
- if (Callback) Callback(data);
|
|
|
- } else {
|
|
|
- if (Callback) Callback(false);
|
|
|
- }
|
|
|
- }
|
|
|
- },
|
|
|
-
|
|
|
- /**
|
|
|
- * 查询区域内规划地块
|
|
|
- * @param {*} BenchmarkLandPrice 基准地价信息
|
|
|
- * @returns {*} 图层返回信息
|
|
|
- * {
|
|
|
- id: uuidv4(),
|
|
|
- layerInfo: layer,
|
|
|
- data: [],
|
|
|
- geometry: feature.geometry,
|
|
|
- }
|
|
|
- */
|
|
|
- async getGHDKs(BenchmarkLandPrice) {
|
|
|
- let featureDataList = [];
|
|
|
- //查询规划地块
|
|
|
- let layer = flatten(window.layerTree).find((item) => item.core == "003006");
|
|
|
- let points = BenchmarkLandPrice.ProjectScope;
|
|
|
- let geo = {
|
|
|
- id: 0,
|
|
|
- style: null,
|
|
|
- parts: [points.length],
|
|
|
- points: points,
|
|
|
- type: "REGION",
|
|
|
- prjCoordSys: {
|
|
|
- epsgCode: null,
|
|
|
- },
|
|
|
- };
|
|
|
- let queryByGeometryParameters = {
|
|
|
- getFeatureMode: "SPATIAL",
|
|
|
- datasetNames: [
|
|
|
- layer.date_server.dataSourceName + ":" + layer.date_server.datasetName,
|
|
|
- ],
|
|
|
- geometry: geo,
|
|
|
- spatialQueryMode: "INTERSECT",
|
|
|
- hasGeometry: "true",
|
|
|
- };
|
|
|
- let e = await mapQuery(
|
|
|
- layer.date_server.url + "/featureResults.json?returnContent=true",
|
|
|
- queryByGeometryParameters
|
|
|
- );
|
|
|
- if (e && e.totalCount > 0) {
|
|
|
- e.features.forEach((feature) => {
|
|
|
- let featureData = {
|
|
|
- id: uuidv4(),
|
|
|
- layerInfo: layer,
|
|
|
- data: [],
|
|
|
- geometry: feature.geometry,
|
|
|
- };
|
|
|
-
|
|
|
- feature.fieldNames.forEach((fieldName, i) => {
|
|
|
- let Field = e.datasetInfos
|
|
|
- ? e.datasetInfos[0].fieldInfos.find(
|
|
|
- (c) =>
|
|
|
- c.name &&
|
|
|
- c.name.toUpperCase() == fieldName.toUpperCase() &&
|
|
|
- c.name.toUpperCase().indexOf("SM") == -1
|
|
|
- )
|
|
|
- : null;
|
|
|
- if (Field) {
|
|
|
- featureData.data.push({
|
|
|
- label: fieldName,
|
|
|
- labelCN: Field ? Field.caption : fieldName,
|
|
|
- value: feature.fieldValues[i],
|
|
|
- });
|
|
|
- }
|
|
|
- });
|
|
|
- featureDataList.push(featureData);
|
|
|
- });
|
|
|
- }
|
|
|
- return featureDataList;
|
|
|
- },
|
|
|
-
|
|
|
- /**
|
|
|
- * 获取基准地价
|
|
|
- * @param {*} BenchmarkLandPrice
|
|
|
- * @param {*} GHDKList
|
|
|
- * @returns
|
|
|
- */
|
|
|
- async GetJZDJ(BenchmarkLandPrice, GHDKList) {
|
|
|
- for (let GHDKindex = 0; GHDKindex < GHDKList.length; GHDKindex++) {
|
|
|
- let GHDK = GHDKList[GHDKindex];
|
|
|
- GHDK.JZDJList = [];
|
|
|
- let KZXXGYDDM = GHDK.data.find((c) => c.label == "KZXXGYDDM");
|
|
|
- let landUseNatureRelation = flatten(window.landUseNatureRelation).find(
|
|
|
- (c) => {
|
|
|
- let GHDKType = c.GHDKType.find((a) =>
|
|
|
- KZXXGYDDM.value.startsWith(a.ydxzdl)
|
|
|
- );
|
|
|
-
|
|
|
- return GHDKType != null;
|
|
|
- }
|
|
|
- );
|
|
|
- if (landUseNatureRelation) {
|
|
|
- let layerinfos = landUseNatureRelation.layerinfo;
|
|
|
- for (let index = 0; index < layerinfos.length; index++) {
|
|
|
- let layerinfo = layerinfos[index];
|
|
|
- let layer = flatten(window.layerTree).find(
|
|
|
- (item) => item.core == layerinfo.layerCore
|
|
|
- );
|
|
|
- if (layer && layer.date_server && layer.date_server.url) {
|
|
|
- let points = BenchmarkLandPrice.ProjectScope;
|
|
|
- let geo = {
|
|
|
- id: 0,
|
|
|
- style: null,
|
|
|
- parts: [points.length],
|
|
|
- points: points,
|
|
|
- type: "REGION",
|
|
|
- prjCoordSys: {
|
|
|
- epsgCode: null,
|
|
|
- },
|
|
|
- };
|
|
|
- let queryByGeometryParameters = {
|
|
|
- getFeatureMode: "SPATIAL",
|
|
|
- datasetNames: [
|
|
|
- layer.date_server.dataSourceName +
|
|
|
- ":" +
|
|
|
- layer.date_server.datasetName,
|
|
|
- ],
|
|
|
- geometry: geo,
|
|
|
- spatialQueryMode: "INTERSECT",
|
|
|
- hasGeometry: "true",
|
|
|
- };
|
|
|
- let e = await mapQuery(
|
|
|
- layer.date_server.url + "/featureResults.json?returnContent=true",
|
|
|
- queryByGeometryParameters
|
|
|
- );
|
|
|
- if (e && e.totalCount > 0) {
|
|
|
- e.features.forEach((feature) => {
|
|
|
- let layerData = {
|
|
|
- id: uuidv4(),
|
|
|
- layerInfo: layerinfo,
|
|
|
- data: [],
|
|
|
- geometry: feature.geometry,
|
|
|
- };
|
|
|
- feature.fieldNames.forEach((fieldName, i) => {
|
|
|
- let Field = e.datasetInfos
|
|
|
- ? e.datasetInfos[0].fieldInfos.find(
|
|
|
- (c) =>
|
|
|
- c.name &&
|
|
|
- c.name.toUpperCase() == fieldName.toUpperCase() &&
|
|
|
- c.name.toUpperCase().indexOf("SM") == -1
|
|
|
- )
|
|
|
- : null;
|
|
|
- if (Field) {
|
|
|
- layerData.data.push({
|
|
|
- label: fieldName,
|
|
|
- labelCN: Field ? Field.caption : fieldName,
|
|
|
- value: feature.fieldValues[i],
|
|
|
- });
|
|
|
- }
|
|
|
- });
|
|
|
- GHDK.JZDJList.push(layerData);
|
|
|
- });
|
|
|
- }
|
|
|
- } else {
|
|
|
- console.log(
|
|
|
- "服务编码【" + layerinfo.layerCore + "】配置不正确,请检查后重试"
|
|
|
- );
|
|
|
- }
|
|
|
- }
|
|
|
- } else {
|
|
|
- console.log(KZXXGYDDM + "未配置当前土地用途,请联系管理员添加");
|
|
|
- }
|
|
|
- }
|
|
|
- },
|
|
|
-
|
|
|
- /**
|
|
|
- * 计算基准地价
|
|
|
- */
|
|
|
- async CalculateJZDJ(BenchmarkLandPrice, GHDKList) {
|
|
|
- //转换成超图格式面
|
|
|
- let HZgeometry = {
|
|
|
- partTopo: [1],
|
|
|
- parts: [BenchmarkLandPrice.ProjectScope.length],
|
|
|
- points: BenchmarkLandPrice.ProjectScope,
|
|
|
- };
|
|
|
- // 循环规划地块区域信息
|
|
|
- GHDKList.forEach((GHDK) => {
|
|
|
- GHDK.fromData = [];
|
|
|
- GHDK.JZDJList.forEach((JZDJ) => {
|
|
|
- let area = this.calculateIntersectArea(
|
|
|
- GHDK.geometry,
|
|
|
- JZDJ.geometry,
|
|
|
- HZgeometry
|
|
|
- );
|
|
|
- if (Number(area) > 0) {
|
|
|
- let fromData = { data: JZDJ.data, area: Number(area) / 666.66 };
|
|
|
- //土地级别
|
|
|
- let tdjb = JZDJ.data.find(
|
|
|
- (c) => c.label == JZDJ.layerInfo.tdjbField
|
|
|
- ).value;
|
|
|
- fromData.tdjb = tdjb;
|
|
|
- //单价
|
|
|
- let dj = JZDJ.data.find(
|
|
|
- (c) => c.label == JZDJ.layerInfo.ydlxPriceField
|
|
|
- ).value;
|
|
|
- fromData.dj = (Number(dj) * 666.66) / 10000;
|
|
|
- //占用面积地价
|
|
|
- fromData.zdj = Number(area * dj) / 10000;
|
|
|
- fromData.geometry = JZDJ.geometry;
|
|
|
- GHDK.fromData.push(fromData);
|
|
|
- }
|
|
|
- });
|
|
|
- });
|
|
|
-
|
|
|
- //计算分析面积
|
|
|
- let fxzmj = this.CalculateAnalysisArea(BenchmarkLandPrice);
|
|
|
- //总地价
|
|
|
- let zdj = 0;
|
|
|
- //参与地价计算的面积(规划地块+绘制范围+地价地块范围三者交集)
|
|
|
- let jszmj = 0;
|
|
|
- //规划地块详细地价
|
|
|
- let ghdkDetailedDjList = [];
|
|
|
- GHDKList.forEach((GHDK) => {
|
|
|
- let ghdkDetailed = {};
|
|
|
- ghdkDetailed.id = GHDK.id;
|
|
|
- ghdkDetailed.geometry = GHDK.geometry;
|
|
|
- ghdkDetailed.data = GHDK.data;
|
|
|
- ghdkDetailed.ghyt = GHDK.data.find((c) => c.label == "KZXXGYDMC").value;
|
|
|
- ghdkDetailed.dkbm = GHDK.data.find((c) => c.label == "DKBM").value;
|
|
|
- //占用土地
|
|
|
- ghdkDetailed.zytdList = GHDK.fromData;
|
|
|
- ghdkDetailedDjList.push(ghdkDetailed);
|
|
|
- GHDK.fromData.forEach((fromData) => {
|
|
|
- zdj += fromData.zdj;
|
|
|
- jszmj += fromData.area;
|
|
|
- });
|
|
|
- });
|
|
|
- ghdkDetailedDjList = ghdkDetailedDjList.filter(
|
|
|
- (c) => c.zytdList && c.zytdList.length > 0
|
|
|
- );
|
|
|
-
|
|
|
- return {
|
|
|
- projectName: BenchmarkLandPrice.ProjectName, //项目名称
|
|
|
- analysisDate: moment(new Date()).format("YYYY-MM-DD HH:mm:ss"), //分析时间
|
|
|
- // projectType: BenchmarkLandPrice.ProjectType, //项目类型
|
|
|
- buildUnit: BenchmarkLandPrice.ConstructionUnit, //建设单位
|
|
|
- analysisArea: fxzmj, //分析面积(亩)
|
|
|
- ParticipationArea: jszmj, //参与地价计算的面积(亩)
|
|
|
- totalLandPrice: zdj, //总地价(万元) 总地价预估合计
|
|
|
- ghdkDetailedDjList: ghdkDetailedDjList, //规划地块详细地价
|
|
|
- picBase64List: [], //图片集 分析范围
|
|
|
- };
|
|
|
- },
|
|
|
-
|
|
|
- /**
|
|
|
- * 计算分析面积
|
|
|
- * @param {*} BenchmarkLandPrice 基准地价信息
|
|
|
- * @returns {number} 分析总面积
|
|
|
- */
|
|
|
- CalculateAnalysisArea(BenchmarkLandPrice) {
|
|
|
- let fxzmj = 0;
|
|
|
- let points = BenchmarkLandPrice.ProjectScope;
|
|
|
- let DataPs = points.map((item) => [item.x, item.y]);
|
|
|
- if (DataPs.length > 0) {
|
|
|
- if (
|
|
|
- JSON.stringify(DataPs[0]) != JSON.stringify(DataPs[DataPs.length - 1])
|
|
|
- ) {
|
|
|
- DataPs.push(DataPs[0]);
|
|
|
- }
|
|
|
- debugger;
|
|
|
- let polygonPs = polygon([DataPs]);
|
|
|
- let sdsd = area(polygonPs);
|
|
|
- let areaPs = area(polygonPs) * window.earthRadius;
|
|
|
- fxzmj = Number(areaPs / 666.66).toFixed(2);
|
|
|
- }
|
|
|
- return fxzmj;
|
|
|
- },
|
|
|
-
|
|
|
- /**
|
|
|
- * 计算两个面的交集面积
|
|
|
- * @param Points1
|
|
|
- * @param Points2
|
|
|
- * @param Points3
|
|
|
- */
|
|
|
- calculateIntersectArea(Points1, Points2, Points3) {
|
|
|
- var geometry1;
|
|
|
- var geometry2;
|
|
|
- var geometry3;
|
|
|
- if (Points1) {
|
|
|
- geometry1 = this.AssemblySurface(Points1);
|
|
|
- }
|
|
|
- if (Points2) {
|
|
|
- geometry2 = this.AssemblySurface(Points2);
|
|
|
- }
|
|
|
- if (Points3) {
|
|
|
- geometry3 = this.AssemblySurface(Points3);
|
|
|
- }
|
|
|
- let areaPs;
|
|
|
- if (geometry1 && geometry2) {
|
|
|
- areaPs = intersect(geometry1, geometry2);
|
|
|
- }
|
|
|
- if (areaPs && geometry3) {
|
|
|
- areaPs = intersect(areaPs, geometry3);
|
|
|
- }
|
|
|
- if (areaPs) {
|
|
|
- let areadifference = area(areaPs) * window.earthRadius;
|
|
|
- return Number(areadifference.toFixed(2));
|
|
|
- } else {
|
|
|
- return 0;
|
|
|
- }
|
|
|
- },
|
|
|
- /**
|
|
|
- * 组装带洞和不带洞的面
|
|
|
- * @param geometry 超图返回的图形数组
|
|
|
- */
|
|
|
- AssemblySurface(geometry) {
|
|
|
- let point3ds = [];
|
|
|
- let pointholes = [];
|
|
|
- let startindex = 0;
|
|
|
- for (let index = 0; index < geometry.parts.length; index++) {
|
|
|
- let thisps = [];
|
|
|
- let endindex = startindex + geometry.parts[index];
|
|
|
- let geometryPoints = geometry.points.slice(startindex, endindex);
|
|
|
- for (
|
|
|
- let pointindex = 0;
|
|
|
- pointindex < geometryPoints.length;
|
|
|
- pointindex++
|
|
|
- ) {
|
|
|
- let point = geometryPoints[pointindex];
|
|
|
- thisps.push([point.x, point.y]);
|
|
|
- }
|
|
|
- if (geometry.partTopo[index] === 1) {
|
|
|
- point3ds.push(thisps);
|
|
|
- } else {
|
|
|
- pointholes.push(thisps);
|
|
|
- }
|
|
|
-
|
|
|
- startindex = endindex;
|
|
|
- }
|
|
|
-
|
|
|
- // 主多边形
|
|
|
- var mainPolygon = polygon(point3ds);
|
|
|
-
|
|
|
- // 岛洞多边形
|
|
|
- var holePolygon = polygon(pointholes);
|
|
|
-
|
|
|
- // 计算差集 得到有导洞的图形
|
|
|
- var getdifference = difference(mainPolygon, holePolygon);
|
|
|
- return getdifference;
|
|
|
- },
|
|
|
-};
|
|
|
-/**
|
|
|
- * 基准地价分析
|
|
|
- */
|
|
|
-export default BenchmarkLandPrice;
|