CN107424210A - The processing method and processing device that mountain area mima type microrelief boundary profile smoothly transits - Google Patents

The processing method and processing device that mountain area mima type microrelief boundary profile smoothly transits Download PDF

Info

Publication number
CN107424210A
CN107424210A CN201710459339.1A CN201710459339A CN107424210A CN 107424210 A CN107424210 A CN 107424210A CN 201710459339 A CN201710459339 A CN 201710459339A CN 107424210 A CN107424210 A CN 107424210A
Authority
CN
China
Prior art keywords
point
mountain area
mima type
type microrelief
boundary profile
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201710459339.1A
Other languages
Chinese (zh)
Inventor
张宏杰
李峰
杨风利
韩军科
党会学
党会博
廖志军
郭勇
潘峰
沈建国
张飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Energy Source In China Construction Group Zhejiang Province Power Design Institute Co Ltd
STATE GRID JIANGXI ELECTRIC POWER Co
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
Energy Source In China Construction Group Zhejiang Province Power Design Institute Co Ltd
STATE GRID JIANGXI ELECTRIC POWER Co
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Energy Source In China Construction Group Zhejiang Province Power Design Institute Co Ltd, STATE GRID JIANGXI ELECTRIC POWER Co, State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI filed Critical Energy Source In China Construction Group Zhejiang Province Power Design Institute Co Ltd
Priority to CN201710459339.1A priority Critical patent/CN107424210A/en
Publication of CN107424210A publication Critical patent/CN107424210A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05Geographic models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Remote Sensing (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides the processing method and processing device that a kind of mountain area mima type microrelief boundary profile smoothly transits.Wherein, this method comprises the following steps:Selecting step, multiple points are chosen on the boundary profile of the pending mountain area mima type microrelief model built in advance;Difference step, single hop cubic spline interpolation is carried out respectively using the functional expression pre-established and to each point;Step is determined, determines the cubic spline difference result of each point along the boundary condition of bar difference three times according to the single hop of default each point.Method provided by the present invention is realized using cubic spline and pending mountain area mima type microrelief boundary profile is smoothly transitted, and then it ensure that in transient process as far as possible so that air-flow is smooth to flow through changeover portion, it largely reducing influence of the mutation to incoming, so as to more closing to reality situation, the precision of Simulation of Wind Field analysis result is improved.

Description

The processing method and processing device that mountain area mima type microrelief boundary profile smoothly transits
Technical field
The present invention relates to mountain area mima type microrelief wind field characteristic technical field, in particular to a kind of mountain area mima type microrelief border The processing method and processing device of contour smoothing transition.
Background technology
When distinguished and admirable area with a varied topography through mountain region, hills etc., wind speed and wind profile in subaerial atmosphere Significant changes can occur.According to the building and structures of flat ground Specification Design, such as transmission line of electricity, the near-earth in mountain region It is easier to wreck in wind field.It is one of feasible way of clearly this potential risk by mima type microrelief flow field simulation, micro-ly The modeling of shape earth's surface is to influence the key of Simulation of Wind Field analysis result.
Traditional earth's surface modeling uses arc section transition, and this processing method can guarantee that the profile with mima type microrelief on circular arc top Line is concordant, but then without the characteristics of smoothly transitting, is grown suddenly equivalent on ground with ground intersection, arc section The protuberance distortion face that one height increases rapidly, has obvious perturbation action to incoming, the model established in this case is simultaneously Actual conditions can not be reflected, and then reduce the precision of Simulation of Wind Field analysis result.
The content of the invention
In consideration of it, the present invention proposes the processing method and processing device that a kind of mountain area mima type microrelief boundary profile smoothly transits, purport Solving the problems, such as that current earth's surface modeling is low using Simulation of Wind Field analysis result precision caused by arc transition.
On one side, the present invention proposes the processing method that a kind of mountain area mima type microrelief boundary profile smoothly transits.This method Comprise the following steps:Selecting step, multiple points are chosen on the boundary profile of the pending mountain area mima type microrelief model built in advance; Difference step, single hop cubic spline interpolation is carried out respectively to each point using the functional expression pre-established;Step is determined, according to default The single hop of each point the cubic spline difference result of each point is determined along the boundary condition of bar difference three times.
Further, in the processing method that above-mentioned mountain area mima type microrelief boundary profile smoothly transits, pending mountain area mima type microrelief The structure of model comprises the following steps:Obtaining step, obtain the altitude data of pending mountain area mima type microrelief;Construction step, according to Altitude data builds pending mountain area mima type microrelief model.
Further, in the processing method that above-mentioned mountain area mima type microrelief boundary profile smoothly transits, functional expression is f (xi)= axi 3+bxi 2+cxi+d;In above formula, xiFor i-th point of abscissa variable on boundary profile, a, b, c and d are default system Number.
Further, in the processing method that above-mentioned mountain area mima type microrelief boundary profile smoothly transits, boundary condition is f (xi) |Xi=0=Zi, f (xi)|Xi=X=0, f ' (xi)|Xi=X=0;In above formula, ZiFor i-th point of height variable on boundary profile.
Further, in the processing method that above-mentioned mountain area mima type microrelief boundary profile smoothly transits, according to default each point Single hop the cubic spline difference result of each point is determined along the boundary condition of bar difference three times, including:Establishment step, Establish the equation group of boundary profile along the boundary condition of bar difference three times according to the single hop of functional expression and each point, equation group isIn above formula, xiFor i-th point of abscissa variable on boundary profile, a, B, c and d is predetermined coefficient;ZiFor i-th point of height variable on boundary profile;Solution procedure, to solving equations, with Determine the coefficient value of each point.
The present invention does not carry out list by the functional expression pre-established to each point on pending mountain area mima type microrelief boundary profile Section cubic spline interpolation, and the boundary condition of the single hop cubic spline difference according to default each point, may finally be determined The cubic spline difference result of each point, realizes and pending mountain area mima type microrelief boundary profile is smoothly transitted, and then ensure that Cause air-flow is smooth to flow through changeover portion as far as possible in transient process, influence of the mutation to incoming is largely reducing, so as to more Closing to reality situation, improve the precision of Simulation of Wind Field analysis result.
On the other hand, the invention also provides the processing unit that a kind of mountain area mima type microrelief boundary profile smoothly transits.The dress Put including:Choose module, difference block and setting module.Wherein, module is chosen for micro- in the pending mountain area built in advance Multiple points are chosen on the boundary profile of relief model;Difference block is used to carry out list respectively to each point using pre-establishing functional expression Section cubic spline interpolation;Determining module is used to determine each point along the boundary condition of bar difference three times according to the single hop of default each point Cubic spline difference result.
Further, in the processing unit that above-mentioned mountain area mima type microrelief boundary profile smoothly transits, choose in module, it is pending The structure of mountain area mima type microrelief model includes:Acquisition module, for obtaining the altitude data of pending mountain area mima type microrelief;Build mould Block, for building pending mountain area mima type microrelief model according to altitude data.
Further, in the processing unit that above-mentioned mountain area mima type microrelief boundary profile smoothly transits, in difference block, functional expression For f (xi)=axi 3+bxi 2+cxi+d;In above formula, xiFor i-th point of abscissa variable on boundary profile, a, b, c and d are equal For predetermined coefficient.
Further, in the processing unit that above-mentioned mountain area mima type microrelief boundary profile smoothly transits, boundary condition is f (xi) |Xi=0=Zi, f (xi)|Xi=X=0, f ' (xi)|Xi=X=0;In above formula, ZiFor i-th point of height variable on boundary profile.
Further, in the processing unit that above-mentioned mountain area mima type microrelief boundary profile smoothly transits, determining module includes:Establish Submodule, for establishing the equation group of boundary profile along the boundary condition of bar difference three times according to the single hop of functional expression and each point, Equation group isIn above formula, xiFor i-th point of abscissa on boundary profile Variable, a, b, c and d are predetermined coefficient;ZiFor i-th point of height variable on boundary profile;Solve submodule, for pair Solving equations, to determine the coefficient value of each point.
The present invention does not carry out list by the functional expression pre-established to each point on pending mountain area mima type microrelief boundary profile Section cubic spline interpolation, and the boundary condition of the single hop cubic spline difference according to default each point, may finally be determined The cubic spline difference result of each point, realizes and pending mountain area mima type microrelief boundary profile is smoothly transitted, and then ensure that Cause air-flow is smooth to flow through changeover portion as far as possible in transient process, influence of the mutation to incoming is largely reducing, so as to more Closing to reality situation, improve the precision of Simulation of Wind Field analysis result.
Brief description of the drawings
By reading the detailed description of hereafter preferred embodiment, it is various other the advantages of and benefit it is common for this area Technical staff will be clear understanding.Accompanying drawing is only used for showing the purpose of preferred embodiment, and is not considered as to the present invention Limitation.And in whole accompanying drawing, identical part is denoted by the same reference numerals.In the accompanying drawings:
Fig. 1 is the flow chart for the processing method that mima type microrelief boundary profile in mountain area provided in an embodiment of the present invention smoothly transits;
It is pending in the processing method that Fig. 2 smoothly transits for mima type microrelief boundary profile in mountain area provided in an embodiment of the present invention The model schematic in some longitude and latitude section of mountain area mima type microrelief;
It is pending in the processing method that Fig. 3 smoothly transits for mima type microrelief boundary profile in mountain area provided in an embodiment of the present invention The result schematic diagram that the cubic spline in some longitude and latitude section on mountain area mima type microrelief border smoothly transits;
It is pending in the processing method that Fig. 4 smoothly transits for mima type microrelief boundary profile in mountain area provided in an embodiment of the present invention The flow chart of the structure of mountain area mima type microrelief model;
Fig. 5 is in the processing method that mima type microrelief boundary profile in mountain area provided in an embodiment of the present invention smoothly transits, it is determined that step Rapid flow chart;
Fig. 6 is the structural frames for the processing unit that mima type microrelief boundary profile in mountain area provided in an embodiment of the present invention smoothly transits Figure;
It is pending in the processing unit that Fig. 7 smoothly transits for mima type microrelief boundary profile in mountain area provided in an embodiment of the present invention The structured flowchart of the structure of mountain area mima type microrelief model;
Fig. 8 is in the processing unit that mima type microrelief boundary profile in mountain area provided in an embodiment of the present invention smoothly transits, and determines mould The structured flowchart of block.
Embodiment
The exemplary embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although the disclosure is shown in accompanying drawing Exemplary embodiment, it being understood, however, that may be realized in various forms the disclosure without should be by embodiments set forth here Limited.On the contrary, these embodiments are provided to facilitate a more thoroughly understanding of the present invention, and can be by the scope of the present disclosure Completely it is communicated to those skilled in the art.It should be noted that in the case where not conflicting, embodiment in the present invention and Feature in embodiment can be mutually combined.Describe the present invention in detail below with reference to the accompanying drawings and in conjunction with the embodiments.
Embodiment of the method:
Referring to Fig. 1, Fig. 1 is the flow for the processing method that the mountain area mima type microrelief boundary profile that the present embodiment provides smoothly transits Figure.As illustrated, this method comprises the following steps:
Selecting step S110, multiple points are chosen on the boundary profile of the pending mountain area mima type microrelief model built in advance.
Specifically, multiple points, of point are chosen on the boundary profile of the pending mountain area mima type microrelief model built in advance Number is relevant with precision, and precision if desired is higher, then needs the number of point chosen more.The pending mountain area built in advance Mima type microrelief model can be the model of some longitude and latitude of pending mountain area mima type microrelief, specifically may refer to Fig. 2.
Difference step S120, single hop cubic spline interpolation is carried out respectively to each point using the functional expression pre-established.
Specifically, for the rugged elevation problem of pending mountain area mima type microrelief boundary contour, using pre-establishing Cubic spline difference functions formula smoothly transitting for independent single hop cubic spline difference is carried out respectively to each point of selection.
Step S130 is determined, determines each point three times along the boundary condition of bar difference three times according to the single hop of default each point Batten difference result.
Specifically, the knot of the cubic spline difference of each point is determined according to the functional expression and default boundary condition that pre-establish Fruit, to carry out smooth transition to pending mountain area mima type microrelief boundary profile, for certain pending mountain area in longitude and latitude section micro-ly The cubic spline on the shape border result that smoothly transits may refer to Fig. 3, the pending mountain area mima type microrelief border in other longitude and latitude sections It can determine that step S140 method is smoothly transitted using selecting step S110-.
As can be seen that the present embodiment by the functional expression that pre-establishes to each on pending mountain area mima type microrelief boundary profile Point not carry out single hop cubic spline interpolation, and according to default boundary condition, the cubic spline of each point may finally be determined Difference result, realize and pending mountain area mima type microrelief boundary profile is smoothly transitted, and then ensure that in transient process to the greatest extent It may cause air-flow is smooth to flow through changeover portion, largely reducing influence of the mutation to incoming, so that more closing to reality situation, is carried The high precision of Simulation of Wind Field analysis result.
Referring to Fig. 4, in one embodiment of the invention, the structure of pending mountain area mima type microrelief model includes following step Suddenly:
Obtaining step S410, obtain the altitude data of pending mountain area mima type microrelief.
Specifically, the altitude data of pending mountain area mima type microrelief can be come from satellite surveying and mapping data, height value of map etc. Obtain, but in view of satellite surveying and mapping data the real space size of space it is smaller, therefore can be based on ASTER GDEM V2 come To the altitude data of pending mountain area mima type microrelief.If the pending mountain area mima type microrelief model built in advance is some longitude and latitude Model, then only obtain the altitude data of the longitude and latitude.
Construction step S420, pending mountain area mima type microrelief model is built according to altitude data.
Specifically, pending mountain area mima type microrelief model is built according to the altitude data of the pending mountain area mima type microrelief of acquisition. It should be noted that the method for pending mountain area mima type microrelief model is built as originally according to the altitude data of pending mountain area mima type microrelief Well known to the personnel of field, here is omitted.
As can be seen that according to the altitude data that obtains pending mountain area mima type microrelief to build pending mountain area micro- for the present embodiment Shape model, the accuracy of structure model is ensure that, is based particularly on ASTER GDEM V2 to obtain pending mountain area mima type microrelief Altitude data, it further ensure that the accuracy of structure model.
In above-described embodiment, in difference step S120, functional expression is f (xi)=axi 3+bxi 2+cxi+ d, wherein, xiFor side I-th point of abscissa variable on boundary's profile, a, b, c and d are predetermined coefficient.xiInterval can be [0, X], xi =0 represents the point on boundary contour, xi=X represents this on the ground border of wind-tunnel wall or numerical model.
In above-described embodiment, boundary condition is f (xi)|Xi=0=Zi, f (xi)|Xi=X=0, f ' (xi)|Xi=X=0, wherein, Zi For i-th point of height variable on boundary profile.It is this because the derivative of mima type microrelief profile boundary position has uncertainty Only consider the condition of continuity, i.e. f (xi)|Xi=0=Zi, and xI=0The first derivative values of point are taken as 0.
Referring to Fig. 5, in one embodiment of the invention, determine that step S140 may further include:
Establishment step S510, boundary profile is established along the boundary condition of bar difference according to the single hop of functional expression and each point three times Equation group, equation group isIn above formula, xiFor on the boundary profile The abscissa variable of i point, a, b, c and d are predetermined coefficient;ZiFor i-th point of height variable on the boundary profile.
Solution procedure S520, to solving equations, to determine the coefficient value of each point.
Specifically, to solving equations, then the coefficient value for obtaining each point is:It is micro- to obtain pending mountain area The cubic spline difference result of each point on landform boundary profile.
As can be seen that the present embodiment by the functional expression and the single hop of each point that pre-establish three times along the perimeter strip of bar difference Part establishes the equation group of boundary profile, and equation group is passed in and out and solved, and can obtain on pending mountain area mima type microrelief boundary profile The cubic spline difference result of each point, and then realize and pending mountain area mima type microrelief boundary profile is smoothly transitted, algorithm is simple Easy, computational efficiency is high.
To sum up, the present embodiment is other to each point on pending mountain area mima type microrelief boundary profile by the functional expression pre-established Single hop cubic spline interpolation, and the boundary condition of the single hop cubic spline difference according to default each point are carried out, may finally The cubic spline difference result of each point is determined, realizes and pending mountain area mima type microrelief boundary profile is smoothly transitted, and then It ensure that in transient process as far as possible so that air-flow is smooth to flow through changeover portion, largely reducing influence of the mutation to incoming, from And more closing to reality situation, improve the precision of Simulation of Wind Field analysis result.
Device embodiment:
Referring to Fig. 6, Fig. 6 is the structure for the processing unit that the mountain area mima type microrelief boundary profile that the present embodiment provides smoothly transits Block diagram.As illustrated, the device includes:Choose module 100, difference block 200 and determining module 300.
Wherein, module 100 is chosen to be used to choose on the boundary profile of the pending mountain area mima type microrelief model built in advance Multiple points.The functional expression that difference block 200 is used to pre-establish carries out single hop cubic spline interpolation respectively to each point.Determining module The 300 cubic spline difference result for determining each point along the boundary condition of bar difference three times according to the single hop of default each point. Wherein, the specific implementation process of the device will not be repeated here referring to the explanation in above method embodiment, the present embodiment.
As can be seen that the present embodiment by the functional expression that pre-establishes to each on pending mountain area mima type microrelief boundary profile Point not carry out single hop cubic spline interpolation, and the boundary condition of the single hop cubic spline difference according to default each point, finally The cubic spline difference result of each point can be determined, realizes and pending mountain area mima type microrelief boundary profile is smoothly transitted, And then ensure that in transient process as far as possible so that air-flow is smooth to flow through changeover portion, it largely reducing shadow of the mutation to incoming Ring, so as to more closing to reality situation, improve the precision of Simulation of Wind Field analysis result.
Referring to Fig. 7, in above-described embodiment, in shown selection module 100, the structure of pending mountain area mima type microrelief model can be with Including:Acquisition module 110 and structure module 120.Wherein, acquisition module 110 is used for the elevation for obtaining pending mountain area mima type microrelief Data.Module 120 is built to be used to build pending mountain area mima type microrelief model according to altitude data.Wherein, mould is obtained in the device The specific implementation process of block 110 and structure module 120 is referring to the explanation in above method embodiment, and the present embodiment is herein not Repeat again.
As can be seen that according to the altitude data that obtains pending mountain area mima type microrelief to build pending mountain area micro- for the present embodiment Shape model, the accuracy of structure model is ensure that, is based particularly on ASTER GDEM V2 to obtain pending mountain area mima type microrelief Altitude data, it further ensure that the accuracy of structure model.
In above-described embodiment, in difference block 200, functional expression is f (xi)=axi 3+bxi 2+cxi+ d, wherein, xiFor border I-th point of abscissa variable on profile, a, b, c and d are predetermined coefficient.xiInterval can be taken as [0, X], xi =0 means the point on boundary contour, xi=X means this on the ground border of wind-tunnel wall or numerical model On.
In above-described embodiment, boundary condition is f (xi)|Xi=0=Zi, f (xi)|Xi=X=0, f ' (xi)|Xi=X=0, wherein, Zi For i-th point of height variable on boundary profile.It is this because the derivative of mima type microrelief profile boundary position has uncertainty Only consider the condition of continuity, i.e. f (xi)|Xi=0=Zi, and xI=0The first derivative values of point are taken as 0.
Referring to Fig. 8, Fig. 8 is the structured flowchart for the determining module 300 that the present embodiment provides.As illustrated, determining module 300 can With including:Setting up submodule 310 and solution submodule 320.Wherein, setting up submodule 310 is used for the list according to functional expression and each point Boundary conditions of the Duan Sanci along bar difference establishes the equation group of boundary profile, and equation group is In above formula, xiFor i-th point of abscissa variable on the boundary profile, a, b, c and d are predetermined coefficient;ZiFor the side I-th point of height variable on boundary's profile.Submodule 320 is solved to be used for solving equations.To determine the coefficient value of each point. Wherein, in the device specific implementation process of determining module 300 referring to the explanation in above method embodiment, the present embodiment It will not be repeated here.
As can be seen that the present embodiment by the functional expression and the single hop of each point that pre-establish three times along the perimeter strip of bar difference Part establishes the equation group of boundary profile, and equation group is passed in and out and solved, and can obtain on pending mountain area mima type microrelief boundary profile The cubic spline difference result of each point, and then realize and pending mountain area mima type microrelief boundary profile is smoothly transitted, algorithm is simple Easy, computational efficiency is high.
To sum up, the present embodiment the present embodiment by the functional expression that pre-establishes on pending mountain area mima type microrelief boundary profile Each point do not carry out single hop cubic spline interpolation, and the boundary condition of the single hop cubic spline difference according to default each point, The cubic spline difference result of each point may finally be determined, realizes the smooth mistake to pending mountain area mima type microrelief boundary profile Cross, and then ensure that in transient process as far as possible so that air-flow is smooth to flow through changeover portion, largely reducing mutation to incoming Influence, so as to more closing to reality situation, improve the precision of Simulation of Wind Field analysis result.
It should be noted that processing method that the mountain area mima type microrelief boundary profile in the present invention smoothly transits and mountain area are micro-ly The processing unit principle that shape boundary profile smoothly transits is similar, and related part can be with cross-referenced.
Obviously, those skilled in the art can carry out the essence of various changes and modification without departing from the present invention to the present invention God and scope.So, if these modifications and variations of the present invention belong to the scope of the claims in the present invention and its equivalent technologies Within, then the present invention is also intended to comprising including these changes and modification.

Claims (10)

1. the processing method that a kind of mountain area mima type microrelief boundary profile smoothly transits, it is characterised in that comprise the following steps:
Multiple points are chosen on the boundary profile of the pending mountain area mima type microrelief model built in advance;
Single hop cubic spline interpolation is carried out respectively to each point using the functional expression pre-established;
Determine the cubic spline difference of each point along the boundary condition of bar difference three times according to the single hop of default each point As a result.
2. the processing method that mima type microrelief boundary profile in mountain area according to claim 1 smoothly transits, it is characterised in that described The structure of pending mountain area mima type microrelief model comprises the following steps:
Obtain the altitude data of the pending mountain area mima type microrelief;
The pending mountain area mima type microrelief model is built according to the altitude data.
3. the processing method that mima type microrelief boundary profile in mountain area according to claim 1 smoothly transits, it is characterised in that
The functional expression is f (xi)=axi 3+bxi 2+cxi+d;
In above formula, xiFor i-th point of abscissa variable on the boundary profile, a, b, c and d are coefficient to be determined.
4. the processing method that mima type microrelief boundary profile in mountain area according to claim 3 smoothly transits, it is characterised in that
The boundary condition is f (xi)|Xi=0=Zi, f (xi)|Xi=X=0, f ' (xi)|Xi=X=0;
In above formula, ZiFor i-th point of height variable on the boundary profile.
5. the processing method that mima type microrelief boundary profile in mountain area according to claim 4 smoothly transits, it is characterised in that described Determine the cubic spline difference result of each point along the boundary condition of bar difference three times according to the single hop of default each point, Including:
Establish the side of the boundary profile along the boundary condition of bar difference three times according to the single hop of the functional expression and each point Journey group, the equation group are
In above formula, xiFor i-th point of abscissa variable on the boundary profile, a, b, c and d are coefficient to be determined;Zi For i-th point of height variable on the boundary profile;
To the solving equations, to determine the coefficient value of each point.
A kind of 6. processing unit that mountain area mima type microrelief boundary profile smoothly transits, it is characterised in that including:
Module is chosen, for choosing multiple points on the boundary profile of the pending mountain area mima type microrelief model built in advance;
Difference block, for carrying out single hop cubic spline interpolation respectively to each point using the functional expression pre-established;
Determining module, for determining each point along the boundary condition of bar difference three times according to the single hop of default each point Cubic spline difference result.
7. the processing unit that mima type microrelief boundary profile in mountain area according to claim 6 smoothly transits, it is characterised in that enter one Step includes:
Acquisition module, for before choosing multiple points on the boundary profile of the pending mountain area mima type microrelief model built in advance, Obtain the altitude data of the pending mountain area mima type microrelief;
Module is built, for building the pending mountain area mima type microrelief model according to the altitude data.
8. the processing unit that mima type microrelief boundary profile in mountain area according to claim 6 smoothly transits, it is characterised in that described Difference block is used to carry out single hop cubic spline interpolation respectively to each point using following functional expressions:
The functional expression is f (xi)=axi 3+bxi 2+cxi+d;
In above formula, xiFor i-th point of abscissa variable on the boundary profile, a, b, c and d are coefficient to be determined.
9. the processing unit that mima type microrelief boundary profile in mountain area according to claim 8 smoothly transits, it is characterised in that
The boundary condition is f (xi)|Xi=0=Zi, f (xi)|Xi=X=0, f ' (xi)|Xi=X=0;
In above formula, ZiFor i-th point of height variable on the boundary profile.
10. the processing unit that mima type microrelief boundary profile in mountain area according to claim 9 smoothly transits, it is characterised in that institute Stating determining module includes:
Setting up submodule, for establishing institute along the boundary condition of bar difference three times according to the single hop of the functional expression and each point The equation group of boundary profile is stated, the equation group is
In above formula, xiFor i-th point of abscissa variable on the boundary profile, a, b, c and d are coefficient to be determined;Zi For i-th point of height variable on the boundary profile;
Submodule is solved, for the solving equations, to determine the coefficient value of each point.
CN201710459339.1A 2017-06-16 2017-06-16 The processing method and processing device that mountain area mima type microrelief boundary profile smoothly transits Pending CN107424210A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710459339.1A CN107424210A (en) 2017-06-16 2017-06-16 The processing method and processing device that mountain area mima type microrelief boundary profile smoothly transits

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710459339.1A CN107424210A (en) 2017-06-16 2017-06-16 The processing method and processing device that mountain area mima type microrelief boundary profile smoothly transits

Publications (1)

Publication Number Publication Date
CN107424210A true CN107424210A (en) 2017-12-01

Family

ID=60428237

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710459339.1A Pending CN107424210A (en) 2017-06-16 2017-06-16 The processing method and processing device that mountain area mima type microrelief boundary profile smoothly transits

Country Status (1)

Country Link
CN (1) CN107424210A (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008025138A (en) * 2006-07-19 2008-02-07 Public Works Research Institute Safety degree evaluating method for rock-bed slope
CN104008572A (en) * 2014-06-17 2014-08-27 武汉大学 Visual measurement and three-dimensional reconstruction method of model terrain in muddy water
CN104700151A (en) * 2014-05-26 2015-06-10 国网辽宁省电力有限公司 Wind power assessment method based on cubic spline interpolation curve-fitting
CN105321204A (en) * 2014-12-18 2016-02-10 马鞍山城智信息技术有限公司 Three-dimensional geological raster body model construction method
CN105868443A (en) * 2016-03-22 2016-08-17 国网安徽省电力公司电力科学研究院 Construction method for microtopographical micrometeorological element field

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008025138A (en) * 2006-07-19 2008-02-07 Public Works Research Institute Safety degree evaluating method for rock-bed slope
CN104700151A (en) * 2014-05-26 2015-06-10 国网辽宁省电力有限公司 Wind power assessment method based on cubic spline interpolation curve-fitting
CN104008572A (en) * 2014-06-17 2014-08-27 武汉大学 Visual measurement and three-dimensional reconstruction method of model terrain in muddy water
CN105321204A (en) * 2014-12-18 2016-02-10 马鞍山城智信息技术有限公司 Three-dimensional geological raster body model construction method
CN105868443A (en) * 2016-03-22 2016-08-17 国网安徽省电力公司电力科学研究院 Construction method for microtopographical micrometeorological element field

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
宋志明: "实时三维地形建模与仿真研究", 《第十一届中国体视学与图像分析学术会议论文集》 *
张雪洪 等: "《生命科学数据处理与MATLAB应用》", 30 April 2013, 上海交通大学出版社 *

Similar Documents

Publication Publication Date Title
CN108648271B (en) Interpolation method for generating complex terrain grid model based on GIS data
CN108921943B (en) Road three-dimensional model modeling method based on lane-level high-precision map
Sester Generalization based on least squares adjustment
CN104252489B (en) A method of position character description information is quickly obtained according to longitude and latitude data
CN107063197B (en) Reservoir characteristic curve extraction method based on spatial information technology
CN105136054B (en) The fine deformation monitoring method of structures and system based on Three Dimensional Ground laser scanning
CN105513133B (en) A kind of city wind environment numerical map makes and display methods
CN105677899B (en) A kind of enhanced vector numerical map production method comprising road direction of travel
CN106021342B (en) City space annual ring cartography and analysis method
CN105787226A (en) Reconstruction of parameterization model of quadrilateral finite element grid model
CN105258704A (en) Multi-scale space-time hot point path detection method based on rapid road network modeling
CN109584364A (en) A kind of generation method and device of the 2.5 dimension electronic maps based on inclination threedimensional model
CN105279793A (en) Modeling method and system based on DEM real three-dimensional map and greenway
CN104548597B (en) The automatic generation method and device of navigation grid
CN107833282A (en) A kind of terrain modeling and mess generation method and device
CN111508015B (en) Building height extraction method and device based on three-dimensional live-action data
CN107291223A (en) A kind of super large data volume virtual reality space Information Visualization System and method
CN106934093A (en) The Trple grid multi-level finite element modeling method of simulation THREE DIMENSIONAL GROUNDWATER FLOW motion
CN105354882A (en) Method for constructing big data architecture based three-dimensional panoramic display platform for large-spatial-range electricity transmission
CN115471634B (en) Modeling method and device for urban green plant twins
CN106875330B (en) Method for rotating plane model into spherical model
CN106569260B (en) For high-dip structure area two-dimension earthquake speed change sky school drawing methods
CN107169080A (en) A kind of geospatial analysis system being combined based on GIS and spatial database
CN107391794A (en) A kind of typhoon continuous stereo Wind-field Retrieval method
CN109710994A (en) Airport obstacle limitation surface based on digital earth transfinites analysis method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20171201

RJ01 Rejection of invention patent application after publication