Patent Document ID: 8065124
Application ID: 12176387
Patent Flag: 1

Claim One:
1. A method for determining relative elevation values of an allotment by using a computer, the allotment comprising a surface having a length L, a width W, and a plurality of nodes i, wherein the total number of said nodes is n; i is a first integral variable and is in the range of [1,n]; said nodes have a row spacing dy and a column spacing dx; and each of said nodes has a relative elevation value relative to said surface; said relative elevation values have a mean z and a standard deviation S d ; the method comprising: a) generating an initial data set of said relative elevation values comprising z 1 0 , z 2 0 ,. .. , and z n 0 by using Monte-Carlo simulation by: a1) entering into the computer values for allotment parameters, the allotment parameters being the length L, the width W, the row spacing dy, the column spacing dx, the mean z , and the standard deviation S d , and setting i equal to 1; a2) by means of the computer, calculating n, wherein n=(L/dy+1)(W/dx+1); a3) by means of the computer, randomly generating an evenly distributed number r, wherein r is between [0,1]; a4) by means of the computer, determining z from the following equation: ∫ - ∞ z ⁢ 1 S d ⁢ 2 ⁢ ⁢ π ⁢ exp ⁡ [ - 1 2 ⁢ ( z - z _ S d ) 2 ] ⁢ ⅆ z = r ; a5) by means of the computer, determining whether z falls into the numerical range of [ z −3S d , z +3S d ]; when z falls into the numerical range of [ z −3S d , z +3S d ], obtaining an initial relative elevation value z i 0 for the i th node through z i 0 =z and continuing to a6); and when z does not fall into the numerical range of [ z −3S d , z +3S d ], returning to a3); and a6) by means of the computer, determining whether i=n; when i=n, continuing to b); when i is not equal to n, increasing i by 1 and returning to a3); b) by means of the computer, adjusting the initial data set to establish a first adjusted data set of said relative elevation values having spatial dependence and comprising z 1 1 , z 2 1 ,. .. , and z n 1 by: b1) determining an allotment type of the allotment from the length L and the width W, wherein the allotment type is a strip type, a narrow type, or a wide type; the strip type has a ratio of the length L to the width W larger than 3 and an absolute value of the width W smaller than 10 m; the narrow type has a ratio of the length L to the width W larger than 3 and an absolute value of the width W larger than or equal to 10 m; and the wide type has a ratio of the length L to the width W smaller than 3 and an absolute value of the width W larger than or equal to 10 m; b2) calculating spatial variation parameters C 0 , C 0 +C, and R from the following empirical formulas 
 for the strip type, C 0 =0.21 S d 2 , ( C 0 +C )= S d 2 , and R= 0.18 L+ 1.53; 
 for the narrow type, C 0 =0.34 S d 2 , ( C 0 +C )= S d 2 and R= 0.21 L− 4.11; and 
 for the wide type, C 0 =0.32 S d 2 , ( C 0 +C )= S d 2 , and R= 16.69( L*W )+5.26; wherein S d 2 is the square of S d ; C 0 is a nugget, C 0 +C is a sill, and R is a range; and calculating a spatial structure function γ(h) for said relative elevation values from the following formula: Y ⁡ ( h ) = { 0 h = 0 C 0 + C · ( 3 2 · h R - 1 2 · h 3 R 3 ) 0 < h ≤ R C 0 + C h > R. wherein h is a spacing distance between two of said nodes; b3) calculating spatial-position weight coefficients of said nodes using Kriging interpolation, wherein the weight coefficient of the j th node is λ j , j is a second integral variable, and j is in the range of [1,n]; and b4) calculating the first adjusted data set from the following equation: 
 z j 1 =Σ k=1 M λ k z k 0 , wherein k is a third integral variable, k is in a range of [1,n], M is an amount of surrounding nodes of said nodes and each of said surrounding nodes has a distance from the j th node smaller than or equal to R; z j 1 is a relative elevation value of said first adjusted data set and z k 0 is a relative elevation value of said initial data set; c) by means of the computer, making statistical parameter corrections to the first adjusted data set to obtain a final data set of said relative elevation values comprising z 1 3 , z 2 3 ,. .. , and z n 3 by: c1) making a mean correction of the first adjusted data set to form a second adjusted data set comprising z 1 2 , z 2 2 ,. .. , and z n 2 by using the following equation: z j 2 = z _ z 2 _ · z j 1 , wherein z 1 is a mean of the first adjusted data set; and c2) making a standard deviation correction to the second adjusted data set to form the final data set by using the following equation: z j 3 = ( z j 2 - z _ ) · s d s d 2 + z _ , wherein S d 2 is a standard deviation of the second adjusted data set; d) saving the final data set of said relative elevation values; and e) determining elevation distribution of the allotment based on said final data set of said relative elevation values saved in step d) and fixing the elevation distribution of the allotment on a physical medium.