CN106323670B - Soil sampling point determines method and device - Google Patents
Soil sampling point determines method and device Download PDFInfo
- Publication number
- CN106323670B CN106323670B CN201610829349.5A CN201610829349A CN106323670B CN 106323670 B CN106323670 B CN 106323670B CN 201610829349 A CN201610829349 A CN 201610829349A CN 106323670 B CN106323670 B CN 106323670B
- Authority
- CN
- China
- Prior art keywords
- grid
- soil sampling
- sampling area
- count
- sampled point
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/02—Devices for withdrawing samples
- G01N1/04—Devices for withdrawing samples in the solid state, e.g. by cutting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Sampling And Sample Adjustment (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The application discloses a kind of soil sampling point and determines method and device, this method comprises: determining soil sampling area；According to the area in soil sampling area and specified sampled point quantity, the initial soil sampled point spacing in soil sampling area is determined；According to initial soil sampled point spacing, soil sampling area is converted into raster data, and determine effective grid quantity of current soil sample region grid；According to effective grid quantity of current soil sample region grid, the size of current soil sample region grid is adjusted, so that the difference of effective grid quantity of soil sampling area adjusted grid grid quantity equal or effective with sampled point quantity and sampled point quantity is less than a preset difference threshold；And according to soil sampling area adjusted grid, determine each soil sampling point in soil sampling area.This method can be realized the soil sampling point that specified quantity is automatically determined in specified range, and guarantee that the spacing between each soil sampling dot grid is strictly equal, improve the efficiency that soil sampling is layouted.
Description
Technical field
The present invention relates to soil science technologies, determine method and device in particular to a kind of soil sampling point.
Background technique
Determine that soil sampling point is the important technical for obtaining soil physics and chemical characteristic.Currently, generalling use net
The method of formatting determines soil sampling point, more common in space interpolation analysis and other statistical methods.
In soil sampling practice, it is often necessary to be adopted according to the soil that sampling funds design specified quantity within the specified range
Sampling point, but due to the irregular characteristic of sampling area shape, lead to be difficult the spacing for accurately determining sampling network lattice point.Currently, logical
Sampling network lattice point is tentatively often generated by rough estimate mesh point spacing mode, then according to sample region shape and is falling into sampling
The quantity that dot grid is sampled in area carries out manual setting fractional-sample dot grid position, supplement or deletes fractional-sample dot grid,
To generate the sampling dot grid of specified quantity or approximate specified quantity.But this method low efficiency needs a large amount of manual operationss, raw
At sampling network lattice point spacing can not be essentially equal, cannot achieve sampling dot grid automation generate.
Above- mentioned information are only used for reinforcing the understanding to background of the invention, therefore it disclosed in the background technology part
It may include the information not constituted to the prior art known to persons of ordinary skill in the art.
Summary of the invention
In view of this, the present invention, which provides a kind of soil sampling point, determines method and device, can be realized in specified range from
The dynamic soil sampling point for determining specified quantity, and guarantee that the spacing between each soil sampling dot grid is strictly equal, improve soil
Sample the efficiency layouted.
Other characteristics and advantages of the invention will be apparent from by the following detailed description, or partially by the present invention
Practice and acquistion.
According to an aspect of the present invention, it provides a kind of soil sampling point and determines method, comprising: determine soil sampling area；
According to the area in the soil sampling area and specified sampled point quantity, the initial soil sampled point in the soil sampling area is determined
Spacing；According to the initial soil sampled point spacing, the soil sampling area is converted into raster data, and determination is presently described
Effective grid quantity of soil sampling area grid；According to effective grid quantity of presently described soil sampling area grid, adjustment is worked as
The size of preceding soil sampling area grid so that effective grid quantity of soil sampling area grid adjusted with it is described
Sampled point quantity is equal or the difference of effective grid quantity and the sampled point quantity is less than a preset difference threshold；
And according to the soil sampling area grid adjusted, determine each soil sampling point in the soil sampling area.
According to an embodiment of the present invention, the soil sampling area includes at least one piece of irregular sampling area.
According to an embodiment of the present invention, determine that soil sampling area includes: according to GIS-Geographic Information System acquisition
The vectorial geographical information data in soil sampling area.
According to an embodiment of the present invention, the initial soil sampled point spacing is calculated according to following formula:
Wherein, S_Area is the calculated soil sampling of vectorial geographical information data according to the soil sampling area
The area in area, P_Count are the sampled point quantity, and D_CellSize is the initial soil sampled point spacing.
According to an embodiment of the present invention, according to effective grid quantity of presently described soil sampling area grid, adjustment
The size of presently described soil sampling area grid, so that effective grid quantity of soil sampling area grid adjusted and institute
The difference of equal sampled point quantity or effective grid quantity and the sampled point quantity is stated less than a preset difference threshold
Value includes: step a, the effective grid quantity and the sampled point quantity of more presently described soil sampling area grid, if worked as
Preceding effective grid quantity is equal with the sampled point quantity or presently described effective grid quantity and the sampled point quantity
Difference be less than the difference threshold, then follow the steps b；Otherwise, step c is executed；Step b determines presently described soil sampling
Area's grid is soil sampling area grid adjusted；Step c, if presently described effective grid quantity is greater than the sampling
Point quantity or presently described effective grid quantity be greater than the sampled point quantity and presently described effective grid quantity with it is described
The difference of sampled point quantity is greater than the difference threshold, then increases the size of presently described soil sampling area grid；If current
Effective grid quantity is less than the sampled point quantity or presently described effective grid quantity is less than the sampled point quantity
And the difference of presently described effective grid quantity and the sampled point quantity is greater than the difference threshold, then reduces presently described soil
The size of earth sample region grid；The adjustment size of the size for increasing or reducing presently described soil sampling area grid is under
State formula calculating:
CellSize_change=sign (T_Count-P_Count) * Step_Change；
Wherein, CellSize_change is the adjustment size, and P_Count is the sampled point quantity, and CellSize is
The size of presently described soil sampling area grid, T_Count are presently described effective grid quantity, and Step_Change is one pre-
If adjusting step；And
Step d determines new soil sampling area grid according to the size of soil sampling area grid adjusted
Lattice, and it is presently described soil sampling area grid that the new soil sampling area grid, which is arranged, executes step b again；
According to an embodiment of the present invention, according to the soil sampling area grid adjusted, the soil is determined
Each soil sampling point of earth sample region includes: that the soil sampling area grid adjusted is converted to vectorial geographical information
Data, wherein each soil sampling point that the center position coordinates of each grid are the soil sampling area is believed in the vectorial geographical
Cease the coordinate in data.
According to another aspect of the present invention, a kind of soil sampling point determining device is provided, comprising: sample region determining module,
For determining soil sampling area；Spacing determining module, for the area and specified sampling number according to the soil sampling area
Amount, determines the initial soil sampled point spacing in the soil sampling area；Initial determining module, for according to the initial soil
The soil sampling area is converted to raster data, and determines having for presently described soil sampling area grid by earth sampled point spacing
Imitate grid quantity；Grid adjusts module, and for effective grid quantity according to presently described soil sampling area grid, adjustment is current
The size of soil sampling area grid, so that effective grid quantity of soil sampling area grid adjusted is adopted with described
Sampling point quantity is equal or the difference of effective grid quantity and the sampled point quantity is less than a preset difference threshold；With
And sampled point determining module, for determining the soil sampling area according to the soil sampling area grid adjusted
Each soil sampling point.
According to an embodiment of the present invention, the soil sampling area includes at least one piece of irregular sampling area.
According to an embodiment of the present invention, the grid adjustment module includes: quantity Comparative sub-module, is worked as comparing
The effective grid quantity and the sampled point quantity of preceding soil sampling area grid；Grid determines submodule, for when described
The more the currently active grid quantity of quantity Comparative sub-module is equal with the sampled point quantity or presently described effective grid quantity
When being less than the difference threshold with the difference of the sampled point quantity, determine that presently described soil sampling area grid is adjusted
Soil sampling area grid；Spacing adjusting submodule, for working as the more presently described effective grid of the quantity Comparative sub-module
Lattice quantity is greater than the sampled point quantity or presently described effective grid quantity and is greater than the sampled point quantity and presently described
When the difference of effective grid quantity and the sampled point quantity is greater than the difference threshold, increase presently described soil sampling area grid
The size of lattice；Or work as the more presently described effective grid quantity of the quantity Comparative sub-module and be less than the sampled point quantity,
Or presently described effective grid quantity is less than the sampled point quantity and presently described effective grid quantity and the sampled point
When the difference of quantity is greater than the difference threshold, reduce the size of presently described soil sampling area grid；It is described to increase or reduce
The adjustment size of the size of presently described soil sampling area grid is calculated according to following formula:
CellSize_change=sign (T_Count-P_Count) * Step_Change；
Wherein, CellSize_change is the adjustment size, and P_Count is the sampled point quantity, and CellSize is
The size of presently described soil sampling area grid, T_Count are presently described effective grid quantity, and Step_Change is one pre-
If adjusting step；And
New grid determines submodule, for the size according to soil sampling area grid adjusted, determines new institute
Soil sampling area grid is stated, and it is presently described soil sampling area grid that the new soil sampling area grid, which is arranged,.
According to an embodiment of the present invention, the sampled point determining module includes: raster data transform subblock, is used for
The soil sampling area grid adjusted is converted into vectorial geographical information data, wherein the center of each grid
Coordinate is coordinate of each soil sampling point in the soil sampling area in the vectorial geographical information data
Soil sampling point according to the present invention determines method, may be implemented to automatically generate specified quantity within the specified range
Soil sampling dot grid, and the spacing between each grid is strictly equal, substantially increases the efficiency that soil sampling is layouted.In addition,
Soil sampling point provided by the invention determines that method applies also for irregular sampling area, solves in the prior art for irregular
Specified sampled point is arranged and the unequal problem of spacing between a large amount of manual operationss of needs generated and each sampled points in sample region.
It should be understood that the above general description and the following detailed description are merely exemplary, this can not be limited
Invention.
Detailed description of the invention
Its example embodiment is described in detail by referring to accompanying drawing, above and other target of the invention, feature and advantage will
It becomes more fully apparent.
Fig. 1 is the flow chart that a kind of soil sampling point shown according to an illustrative embodiments determines method.
Fig. 2 is the flow chart that another soil sampling point shown according to an illustrative embodiments determines method.
Fig. 3 A is the schematic diagram according to the soil sampling area shown in an example.
Fig. 3 B is the distribution schematic diagram according to the soil sampling area initial samples point grid shown in an example.
Fig. 4 is a kind of block diagram of soil sampling point determining device shown according to an illustrative embodiments.
Fig. 5 is the block diagram of another soil sampling point determining device shown according to an illustrative embodiments.
Specific embodiment
Example embodiment is described more fully with reference to the drawings.However, example embodiment can be with a variety of shapes
Formula is implemented, and is not understood as limited to example set forth herein；On the contrary, thesing embodiments are provided so that the present invention will more
Fully and completely, and by the design of example embodiment comprehensively it is communicated to those skilled in the art.Attached drawing is only the present invention
Schematic illustrations, be not necessarily drawn to scale.Identical appended drawing reference indicates same or similar part in figure, thus
Repetition thereof will be omitted.
In addition, described feature, structure or characteristic can be incorporated in one or more implementations in any suitable manner
In mode.In the following description, many details are provided to provide and fully understand to embodiments of the present invention.So
And it will be appreciated by persons skilled in the art that technical solution of the present invention can be practiced and omit one in the specific detail
Or more, or can be using other methods, constituent element, device, step etc..In other cases, it is not shown in detail or describes
Known features, method, apparatus, realization or operation are to avoid a presumptuous guest usurps the role of the host and each aspect of the present invention is made to thicken.
Fig. 1 is the flow chart that a kind of soil sampling point shown according to an illustrative embodiments determines method.Such as Fig. 1 institute
Show, soil sampling point determines that method 10 includes:
In step s 102, soil sampling area is determined.
The soil sampling area may include at least one piece of irregular sampling area.
In step S104, according to the area in soil sampling area and specified sampled point quantity, determine that initial soil samples
Point spacing.
In step s 106, to specify lattice dimensions size, soil sampling area is converted into raster data.
When first operating procedure S106, using the initial soil sampled point spacing that step S104 is determined as meaning in this step
" specified lattice dimensions size ".
In step S108, effective grid quantity of current soil sample region grid is determined.
In step s 110, according to effective grid quantity of current soil sample region grid, it is big to adjust current grid size
It is small, so that the effective grid quantity for the sample region raster data that lattice dimensions generate after adjustment and specified sampled point quantity are more
It is close or equal.
Wherein close to the difference threshold for for example referring to that the difference of effective grid quantity and sampled point quantity is specified less than one.
In step S112, according to adjustment post-sampling point grid, the soil sampling point in the soil sampling area is determined.
Step S106, step S108 and step S110 may need to be repeated as many times, until actual arrangement sampling number with
Specified sampled point quantity is consistent or the difference between the two is less than above-mentioned specified difference threshold.
Soil sampling point provided by the invention determines method, may be implemented to automatically generate specified quantity within the specified range
Soil sampling dot grid, and the spacing between each grid is strictly equal, substantially increases the efficiency that soil sampling is layouted.In addition,
Soil sampling point provided by the invention determines that method applies also for irregular sampling area, solves in the prior art for irregular
Specified sampled point is arranged and the unequal problem of spacing between a large amount of manual operationss of needs generated and each sampled points in sample region.
It will be clearly understood that the present disclosure describe how being formed and using particular example, but the principle of the present invention is not limited to
These exemplary any details.On the contrary, the introduction based on present disclosure, these principles can be applied to many other
Embodiment.
Fig. 2 is the flow chart that another soil sampling point shown according to an illustrative embodiments determines method.Such as Fig. 2
Shown, soil sampling point determines that method 20 includes:
In step S202, according to the geography information that GIS-Geographic Information System obtains, soil sampling area is determined.
In the present embodiment, can according to the geography information in the data acquisition soil sampling area in GIS-Geographic Information System,
In the geography information in soil sampling area that gets be vector data.
Fig. 3 A is the schematic diagram according to the soil sampling area shown in an example.As shown in Figure 3A, soil sampling area A is (in figure
Dash area) it include at least one piece of irregular sampling area A1-A3.
In step S204, the area in soil sampling area is determined.
According to the vectorial geographical information in the soil sampling area obtained from GIS-Geographic Information System, the soil sampling area is calculated
Area.
In step S206, according to the area in determining soil sampling area and specified sampled point quantity, initial soil is determined
The spacing of earth sampled point.
For example, it is assumed that the area in the soil sampling area determined is S_Area, specified sampled point quantity is P_Count, then
It can determine the spacing of initial soil sampled point:
In step S208, according to specified lattice dimensions value (being denoted as CellSize), by the vector number in soil sampling area
The raster data that grid side length is CellSize is converted to according to geography information, so that it is determined that initial samples within the scope of soil sampling area
The distribution of point grid.
It should be noted that when vector data is converted to raster data, it should be ensured that grid central point is fallen in vector data
Portion.When first time operating procedure S208 CellSize value be step S206 in determine initial soil sampled point space D _
CellSize。
Fig. 3 B is the distribution schematic diagram according to the soil sampling area initial samples point grid shown in an example.Such as Fig. 3 B institute
Show, wherein the grid comprising point " " is the raster data that vector data conversion in soil sampling area generates, wherein " " is sampling
The central point of area's grid, other are then non-sampled area's grid.Non-sampled area's grid is typically expressed as " no data grid " or invalid
Grid.Therefore the effective grid quantity generated by sampled point spacing CellSize by the conversion of sample region vector data, as spacing
The quantity layouted is sampled for the grid of CellSize, the position of " " is the position of grid sampled point.
In step S210, the quantity of the effective grid generated in statistic procedure S208 is denoted as T_Count, which is
For initial samples point grid point quantity.
In step S212, according to the T_Count being calculated in step S210, current soil sample region grid is adjusted
Size, so that be more nearly before the quantity of effective grid is compared to adjustment after adjustment with specified sampled point quantity, or
Person makes the quantity of effectively grid after adjusting be equal to specified sampled point quantity.
For example, the quantity of more the currently active grid and specified sampled point quantity, if the quantity of the currently active grid
Greater than specified sampled point quantity, then increase the size of current grid, if the quantity of the currently active grid is less than specified adopt
Sampling point quantity then reduces the size of current grid.The adjustment size CellSize_change of lattice dimensions for example can be with are as follows:
CellSize_change=sign (T_Count-P_Count) * Step_Change
The grid sampling number that wherein T_Count, which is the lattice dimensions obtained in step S210, generates when being CellSize
Amount；Sign (T_Count-P_Count) is a discriminant function, and when T_Count-P_Count is greater than zero, which is 1,
When T_Count-P_Count is less than zero, which is -1, when T_Count-P_Count is equal to zero, the function value
It is 0；Step_Change is a preset adjusting step, and length can be configured according to actual needs, and the present invention is not with this
It is limited.
Lattice dimensions CellSize adjusted is equal to current grid size and adds CellSize_change.
S208 to step S212 is repeated, until the quantity of the currently active grid is equal to specified sampled point quantity, or
The difference of the quantity of the currently active grid of person and specified sampled point quantity is less than preset difference threshold, stops above-mentioned adjustment operation,
And generate the sample region raster data for meeting quantitative requirement.
In step S214, according to the sample region raster data for meeting quantitative requirement, the soil in the soil sampling area is determined
Sampled point.
Stop repeat S208 to step S212 after, the soil sampling area raster data of generation, the raster data it is effective
Grid quantity is equal with specified sampled point quantity or the two difference specified difference threshold above-mentioned less than one.
Further according to the raster data, each soil sampling point in final soil sampling area is determined.
Specifically, above-mentioned raster data is converted into vector point data again, wherein the center position coordinates of each grid
Corresponding point is each soil sampling point in soil sampling area.
It will be appreciated by those skilled in the art that realizing that all or part of the steps of above embodiment is implemented as being held by CPU
Capable computer program.When the computer program is executed by CPU, execute above-mentioned defined by the above method provided by the invention
Function.The program can store in a kind of computer readable storage medium, which can be read-only memory,
Disk or CD etc..
Further, it should be noted that above-mentioned attached drawing is only according to the present invention included by the method for illustrative embodiments
Processing schematically illustrates, rather than limits purpose.It can be readily appreciated that above-mentioned processing shown in the drawings does not indicate or limits these
The time sequencing of processing.In addition, being also easy to understand, these processing, which can be, for example either synchronously or asynchronously to be executed in multiple modules
's.
Following is apparatus of the present invention embodiment, can be used for executing embodiment of the present invention method.For apparatus of the present invention reality
Undisclosed details in example is applied, embodiment of the present invention method is please referred to.
Fig. 4 is a kind of block diagram of soil sampling point determining device shown according to an illustrative embodiments.Such as Fig. 4 institute
Show, soil sampling point determining device 30 includes: that sample region determining module 302, spacing determining module 304, initial determine mould
Block 306, grid adjustment module 308 and sampled point determining module 310.
Wherein, sample region determining module 302 is for determining soil sampling area.
The soil sampling area for example may include at least one piece of irregular sampling area as shown in Figure 3A.
Spacing determining module 304 is used for area and specified sampled point quantity according to the soil sampling area, determines institute
State the initial soil sampled point spacing in soil sampling area.
The initial soil sampled point spacing can be calculated according to following formula:
Wherein, S_Area is the calculated soil sampling of vectorial geographical information data according to the soil sampling area
The area in area, P_Count are the sampled point quantity, and D_CellSize is the initial soil sampled point spacing.
Initial determining module 306 is used to be turned the soil sampling area according to the initial soil sampled point spacing
It is changed to raster data, and determines effective grid quantity of presently described soil sampling area grid.
Grid adjusts module 308 for effective grid quantity according to presently described soil sampling area grid, and adjustment is current
The size of soil sampling area grid, so that effective grid quantity of soil sampling area grid adjusted is adopted with described
Sampling point quantity is equal or the difference of effective grid quantity and the sampled point quantity is less than a preset difference threshold.
Sampled point determining module 310 is used to determine the soil according to the soil sampling area grid adjusted
Each soil sampling point of sample region.
Soil sampling point determining device provided by the invention, may be implemented to automatically generate specified quantity within the specified range
Soil sampling point grid, and the spacing between each grid is strictly equal, substantially increases the efficiency that soil sampling is layouted.In addition,
Soil sampling point provided by the invention determines that method applies also for irregular sampling area, solves in the prior art for irregular
Specified sampled point is arranged and the unequal problem of spacing between a large amount of manual operationss of needs generated and each sampled points in sample region.
Fig. 5 is the block diagram of another soil sampling point determining device shown according to an illustrative embodiments.Such as Fig. 5 institute
Show, soil sampling point determining device 40 and soil sampling point determining device 30 shown in Fig. 4 the difference is that:
Sample region determining module 302 may include: that vector data acquisition submodule 3022 is used for according to a geography information system
System obtains the vectorial geographical information data in the soil sampling area.
Grid adjustment module 308 may include: that quantity Comparative sub-module 3084, grid determine submodule 3086, spacing tune
Whole submodule 3088 and new grid determine submodule 3090.
Wherein, quantity Comparative sub-module 3084 for more presently described soil sampling area grid effective grid quantity and
The sampled point quantity；
Grid determine submodule 3086 for when the more the currently active grid quantity of the quantity Comparative sub-module with it is described
When sampled point quantity is equal or the difference of presently described effective grid quantity and the sampled point quantity is less than the difference threshold,
Determine that presently described soil sampling area grid is soil sampling area grid adjusted.
Spacing adjusting submodule 3088 is used for when the more presently described effective grid quantity of quantity Comparative sub-module 3084 is big
It is greater than the sampled point quantity and presently described effective grid in the sampled point quantity or presently described effective grid quantity
When the difference of quantity and the sampled point quantity is greater than the difference threshold, increase the ruler of presently described soil sampling area grid
It is very little；Perhaps when the more presently described effective grid quantity of the quantity Comparative sub-module is less than the sampled point quantity or works as
Preceding effective grid quantity is less than the sampled point quantity and presently described effective grid quantity and the sampled point quantity
When difference is greater than the difference threshold, reduce the size of presently described soil sampling area grid；It is described to increase or reduce current institute
The adjustment size for stating the size of soil sampling area grid is calculated according to following formula:
CellSize_change=sign (T_Count-P_Count) * Step_Change
Wherein CellSize_change is the adjustment size, and P_Count is the sampled point quantity, and CellSize is
The size of presently described soil sampling area grid, T_Count are presently described effective grid quantity, and Step_Change is one pre-
If adjusting step.New grid determines that submodule 3090 is used for the size according to soil sampling area grid adjusted, really
Fixed new soil sampling area grid, and it is presently described soil sampling area that the new soil sampling area grid, which is arranged,
Grid.
Sampled point determining module 310 may include: raster data transform subblock 3102, be used for the institute adjusted
It states soil sampling area grid and is converted to vectorial geographical information data, wherein the center position coordinates of each grid are that the soil is adopted
Coordinate of each soil sampling point in sample area in the vectorial geographical information data.
It should be noted that above-mentioned block diagram shown in the drawings is functional entity, not necessarily must with physically or logically
Independent entity is corresponding.Can realize these functional entitys using software form, or in one or more hardware modules or
These functional entitys are realized in integrated circuit, or are realized in heterogeneous networks and/or processor device and/or microcontroller device
These functional entitys.
Through the above description of the embodiments, those skilled in the art is it can be readily appreciated that example described herein is implemented
Mode can also be realized by software realization in such a way that software is in conjunction with necessary hardware.Therefore, according to the present invention
The technical solution of embodiment can be embodied in the form of software products, which can store non-volatile at one
Property storage medium (can be CD-ROM, USB flash disk, mobile hard disk etc.) in or network on, including some instructions are so that a calculating
Equipment (can be personal computer, server, mobile terminal or network equipment etc.) executes embodiment according to the present invention
Method.
It is particularly shown and described exemplary embodiments of the present invention above.It should be appreciated that the present invention is unlimited
In detailed construction described herein, set-up mode or implementation method；On the contrary, it is intended to cover included in appended claims
Spirit and scope in various modifications and equivalence setting.
Claims (10)
1. a kind of soil sampling point determines method characterized by comprising
Determine soil sampling area；
According to the area in the soil sampling area and specified sampled point quantity, determine that the initial soil in the soil sampling area is adopted
Sampling point spacing；
According to the initial soil sampled point spacing, the soil sampling area is converted into raster data, and determination is presently described
Effective grid quantity of soil sampling area grid；
According to effective grid quantity of presently described soil sampling area grid, the size of presently described grid is adjusted, so that adjustment
Effective grid quantity of soil sampling area grid afterwards is equal with the sampled point quantity or effective grid quantity
Difference with the sampled point quantity is less than a preset difference threshold；And
According to the soil sampling area grid adjusted, each soil sampling point in the soil sampling area is determined.
2. the method according to claim 1, wherein the soil sampling area includes at least one piece of irregular sampling
Area.
3. method according to claim 1 or 2, which is characterized in that determine that soil sampling area includes: according to a geography information
System obtains the vectorial geographical information data in the soil sampling area.
4. according to the method described in claim 3, it is characterized in that, the initial soil sampled point spacing is according to following formula meters
It calculates:
Wherein, S_Area is the calculated soil sampling area of vectorial geographical information data according to the soil sampling area
Area, P_Count are the sampled point quantity, and D_CellSize is the initial soil sampled point spacing.
5. according to the method described in claim 3, it is characterized in that, according to effective grid of presently described soil sampling area grid
Quantity adjusts the size of presently described soil sampling area grid, so that effective grid of soil sampling area grid adjusted
Lattice quantity is equal with the sampled point quantity or the difference of effective grid quantity and the sampled point quantity is pre- less than one
If difference threshold include:
Step a, the effective grid quantity and the sampled point quantity of more presently described soil sampling area grid, if current institute
State that effective grid quantity is equal with the sampled point quantity or the difference of presently described effective grid quantity and the sampled point quantity
Value is less than the difference threshold, thens follow the steps b；Otherwise, step c is executed；
Step b determines that presently described soil sampling area grid is soil sampling area grid adjusted；
Step c, if presently described effective grid quantity is greater than the sampled point quantity or presently described effective grid quantity
It is greater than the difference threshold greater than the difference of the sampled point quantity and presently described effective grid quantity and the sampled point quantity
Value, then increase the size of presently described soil sampling area grid；If presently described effective grid quantity is less than the sampled point
Quantity or presently described effective grid quantity are less than the sampled point quantity and presently described effective grid quantity is adopted with described
The difference of sampling point quantity is greater than the difference threshold, then reduces the size of presently described soil sampling area grid；It is described increase or
The adjustment size for reducing the size of presently described soil sampling area grid is calculated according to following formula:
CellSize_change=sign (T_Count-P_Count) * Step_Change；
Wherein, CellSize_change is the adjustment size, and P_Count is the sampled point quantity, and CellSize is current
The size of soil sampling area grid, T_Count are presently described effective grid quantity, and Step_Change is one preset
Adjusting step, sign (T_Count-P_Count) they are a discriminant function, when T_Count-P_Count is greater than zero, sign
(T_Count-P_Count) value is 1, when T_Count-P_Count is less than zero, sign (T_Count-P_Count) value
It is -1, when T_Count-P_Count is equal to zero, sign (T_Count-P_Count) value is 0；And
Step d determines new soil sampling area grid according to the size of soil sampling area grid adjusted, and
It is presently described soil sampling area grid that the new soil sampling area grid, which is arranged, executes step b again.
6. according to the method described in claim 3, it is characterized in that, according to the soil sampling area grid adjusted,
Each soil sampling point for determining the soil sampling area includes: that the soil sampling area grid adjusted is converted to arrow
Geographic information data is measured, wherein the center position coordinates of each grid are each soil sampling point in the soil sampling area described
Coordinate in vectorial geographical information data.
7. a kind of soil sampling point determining device characterized by comprising
Sample region determining module, for determining soil sampling area；
Spacing determining module, for according to the soil sampling area area and specified sampled point quantity, determine the soil
The initial soil sampled point spacing of sample region；
Initial determining module, for according to the initial soil sampled point spacing, the soil sampling area to be converted to grid
Lattice data, and determine effective grid quantity of presently described soil sampling area grid；
Grid adjusts module and adjusts presently described soil for effective grid quantity according to presently described soil sampling area grid
The size of earth sample region grid, so that effective grid quantity of soil sampling area grid adjusted and the sampling number
The difference of equal or described effective grid quantity and the sampled point quantity is measured less than a preset difference threshold；And
Sampled point determining module, for determining the soil sampling area according to the soil sampling area grid adjusted
Each soil sampling point.
8. device according to claim 7, which is characterized in that the soil sampling area includes at least one piece of irregular sampling
Area.
9. device according to claim 7 or 8, which is characterized in that the grid adjusts module and includes:
Quantity Comparative sub-module, for effective grid quantity of more presently described soil sampling area grid and the sampling number
Amount；
Grid determines submodule, for when the more the currently active grid quantity of the quantity Comparative sub-module and the sampling number
When measuring the difference of equal or presently described effective grid quantity and the sampled point quantity less than the difference threshold, determine current
Soil sampling area grid is soil sampling area grid adjusted；
Spacing adjusting submodule, for being greater than described adopt when the more presently described effective grid quantity of the quantity Comparative sub-module
Sampling point quantity or presently described effective grid quantity are greater than the sampled point quantity and presently described effective grid quantity and institute
When stating the difference of sampled point quantity greater than the difference threshold, increase the size of presently described soil sampling area grid；Or work as
The more presently described effective grid quantity of quantity Comparative sub-module is less than the sampled point quantity or presently described effective
Grid quantity is less than the sampled point quantity and the difference of presently described effective grid quantity and the sampled point quantity is greater than institute
When stating difference threshold, reduce the size of presently described soil sampling area grid；It is described to increase or reduce presently described soil sampling
The adjustment size of the size of area's grid is calculated according to following formula:
CellSize_change=sign (T_Count-P_Count) * Step_Change；
Wherein, CellSize_change is the adjustment size, and P_Count is the sampled point quantity, and CellSize is current
The size of soil sampling area grid, T_Count are presently described effective grid quantity, and Step_Change is one preset
Adjusting step, sign (T_Count-P_Count) they are a discriminant function, when T_Count-P_Count is greater than zero, sign
(T_Count-P_Count) value is 1, when T_Count-P_Count is less than zero, sign (T_Count-P_Count) value
It is -1, when T_Count-P_Count is equal to zero, sign (T_Count-P_Count) value is 0；And
New grid determines submodule, for the size according to soil sampling area grid adjusted, determines the new soil
Earth sample region grid, and it is presently described soil sampling area grid that the new soil sampling area grid, which is arranged,.
10. device according to claim 7 or 8, which is characterized in that the sampled point determining module includes: raster data
Transform subblock, for the soil sampling area grid adjusted to be converted to vectorial geographical information data, wherein often
The center position coordinates of a grid are each soil sampling point in the soil sampling area in the vectorial geographical information data
Coordinate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610829349.5A CN106323670B (en) | 2016-09-18 | 2016-09-18 | Soil sampling point determines method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610829349.5A CN106323670B (en) | 2016-09-18 | 2016-09-18 | Soil sampling point determines method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106323670A CN106323670A (en) | 2017-01-11 |
| CN106323670B true CN106323670B (en) | 2019-08-09 |
Family
ID=57787282
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610829349.5A Expired - Fee Related CN106323670B (en) | 2016-09-18 | 2016-09-18 | Soil sampling point determines method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106323670B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110411778A (en) * | 2019-06-17 | 2019-11-05 | 淮北矿业股份有限公司 | A kind of coal slime method of sampling |
| CN110457422B (en) * | 2019-08-20 | 2023-06-06 | 重庆壤科农业数据服务有限公司 | Automatic soil sample collection and distribution system and method |
| CN110623555B (en) * | 2019-09-19 | 2022-02-08 | 杭州九阳小家电有限公司 | Food material size calculation method in oven, size identification device and oven |
| CN112200356A (en) * | 2020-09-30 | 2021-01-08 | 杭州鲁尔物联科技有限公司 | Landslide prediction method, device, equipment and storage medium |
| CN112200361A (en) * | 2020-09-30 | 2021-01-08 | 杭州鲁尔物联科技有限公司 | Landslide prediction method, device, equipment and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102645350A (en) * | 2012-03-08 | 2012-08-22 | 浙江大学 | Soil-sampling method based on high-resolution satellite remote sensing data |
| CN103196698A (en) * | 2013-03-20 | 2013-07-10 | 浙江大学 | Soil sampling method based on near-earth sensor technology |
-
2016
- 2016-09-18 CN CN201610829349.5A patent/CN106323670B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102645350A (en) * | 2012-03-08 | 2012-08-22 | 浙江大学 | Soil-sampling method based on high-resolution satellite remote sensing data |
| CN103196698A (en) * | 2013-03-20 | 2013-07-10 | 浙江大学 | Soil sampling method based on near-earth sensor technology |
Non-Patent Citations (2)
| Title |
|---|
| Determining optimal density of grid soil-sampl ing;Wang Hongbin等;《农业工程学报》;20060831;第22卷(第8期);第145-148页 * |
| 土壤采样密度及栅格尺度对处方图结果的影响;王秀等;《中墨农业工程学会2005年学术年会论文集》;20051231;第296-298页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106323670A (en) | 2017-01-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106323670B (en) | Soil sampling point determines method and device | |
| CN106407682B (en) | Soil sampling point determines method and device | |
| Christen et al. | A general purpose sampling algorithm for continuous distributions (the t-walk) | |
| Kalus et al. | Constraints on anisotropic cosmic expansion from supernovae | |
| CN106408604A (en) | Filtering method and device for point cloud data | |
| CN105469446A (en) | Point cloud mesh simplification system and method | |
| US10019649B2 (en) | Point cloud simplification | |
| US7747401B2 (en) | Fast intrinsic mode decomposition of time series data with sawtooth transform | |
| US20110213587A1 (en) | Method and computer program product for finding statistical bounds, corresponding parameter corners, and a probability density function of a performance target for a circuit | |
| CN111047627A (en) | Smooth constraint unscented Kalman filtering method and target tracking method | |
| Cea et al. | Event based sampling in non-linear filtering | |
| CN109142896A (en) | Lightning Warning method based on three-dimensional atmospheric electric field and MEMD | |
| Dang et al. | An EM-based algorithm for clustering data streams in sliding windows | |
| Jin et al. | Stochastic mirror descent method for linear ill-posed problems in Banach spaces | |
| EP2863319A1 (en) | Data sampling method and data sampling device | |
| Bai et al. | A new approach to maximum likelihood estimation of the three‐parameter gamma and Weibull distributions | |
| CN111192302A (en) | Feature matching method based on motion smoothness and RANSAC algorithm | |
| Grall-Maes et al. | Assessing the number of clusters in a mixture model with side-information | |
| JP5391635B2 (en) | Analysis device, data storage method, and data storage program | |
| CN110120058A (en) | A kind of method that elevation scatterplot generates compact outer boundary | |
| JP5801242B2 (en) | Estimated interest score database generation apparatus, method, and program | |
| Foralewski et al. | On the moduli and characteristic of monotonicity in Orlicz-Lorentz function spaces | |
| Aydiner et al. | A sparse data fast Fourier transform (SDFFT) | |
| Xu et al. | ExtremeCast: Boosting Extreme Value Prediction for Global Weather Forecast | |
| CN116522781B (en) | Sample data generation method, model training method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190809 Termination date: 20200918 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |