Landscape planning-oriented image acquisition method for spatial accurate positioning

Abstract

The invention relates to the technical field of urban planning, in particular to a landscape planning-oriented image acquisition method for spatial accurate positioning, which comprises the following steps of: 1, determining a working area and a measuring point; step 2, determining the position of a reference station and parameter conversion; step 3, determining an operation path; step 4, image acquisition; step 5, link processing; and step 6, data storage and management. According to the method, the position of each point can be accurately defined by defining a coordinate point calculation method, and the uniqueness of each point in a track is ensured through an ID (Identity) system, so that a traceable space structure is provided, a parameterization mode is introduced to describe the track, the coordinates are allowed to be dynamically adjusted on the basis of time or progress, and the tracking accuracy is improved. The method adapts to the change in the actual moving process, and achieves the beneficial effects of more comprehensive consideration and no deviation during the calculation of the activity characteristics.

Description

Technical Field

[0001] The invention relates to the technical field of urban planning, and in particular to an image acquisition method for spatial precise positioning of landscape planning. Background Art

[0002] The route navigation instruments that tourists currently carry in scenic spots are based on LBS base station positioning on mobile phones and GPS positioning based on instruments. They are mainly used for location query and route navigation, as well as marking locations when acquiring images. Some scenic spots are displayed through some photos of some scenic spots, but they cannot fully reflect the overall situation. The problem of how tourists can choose routes based on real-life photos has the disadvantages of strong subjectivity and difficulty in operation. It is necessary to provide an image acquisition technology for precise spatial positioning for landscape planning.

[0003] In order to solve the defects of the prior art, a spatial precise positioning image acquisition method for landscape planning is disclosed with publication number CN103335645A. The method is based on the existing open scenic spots that are difficult to accurately locate, and proposes an image acquisition technology based on GPRS fixed points and GPS spatial precise positioning as described in the present invention. It promotes the comprehensive display of high-quality tourist landscapes, provides tourists with the convenience of completely independently choosing routes, has good social, economic and ecological benefits, provides a rational and scientific basis for urban cruise line planning, urban tourism planning and related construction, and makes more efficient and reasonable use of precious urban tourism resources.

[0004] However, the aforementioned technical solution sets the straight line group interval m based on specific tour activity characteristics, but does not take into account the impact of factors such as different terrains and weather conditions on tour activities, resulting in deviations in the final calculation results. At the same time, the aforementioned technical solution mentions the settings of the base station and the mobile station, but does not clearly explain how to ensure the measurement accuracy of the GPS receiver, which leads to errors in the final results. Summary of the invention

[0005] 1. Technical issues to be solved

[0006] In view of the deficiencies of the prior art, the present invention provides an image acquisition method for spatial precise positioning of landscape planning, which has the advantages of taking into account more comprehensiveness when calculating activity characteristics without deviation, and introducing an error correction mechanism to improve the reliability of measurement, thereby solving the problems of the above-mentioned technologies.

[0007] (II) Technical solution

[0008] To achieve the above object, the present invention provides the following technical solution: an image acquisition method for spatial precise positioning of landscape planning, comprising the following steps:

[0009] Step 1. Determination of the working area and measuring points: Obtain complete and accurate land vector graphics and coordinate data of the working area and the surrounding area. On the land vector graphics within the working area, draw straight line groups along the due north-south and due east-west directions, with an interval of m, covering the entire working area. According to the real-time monitored environmental data: terrain, vegetation density, weather conditions, dynamically adjust the value of the interval m: when walking, the interval range is 80 meters to 120 meters; when boating, the interval range is 200 meters to 300 meters; draw points at the intersection of the straight line groups, determine the coordinate data of each point, and number them to form a coordinate grid within the working area;

[0010] Step 2: Determine the position of the base station and parameter conversion: Place a base station GPS receiver at known coordinate point A, and inject the coordinates and elevation of the base station into the receiver through a data recording notebook. Set up another GPS receiver at known coordinate point B to receive satellite signals and data transmitted by the base station, and inject the coordinates and elevation of the point. Calculate the conversion parameters relative to the WGS-84 ellipsoid coordinate system, including coordinate difference, elevation difference, scale factor, rotation angle, introduce differential GPS technology, correct the measurement data, improve measurement accuracy, and place the GPS receiver with the obtained parameters on the survey pole, connected to the notebook, and set as a mobile station.

[0011] Step 3: Determine the operation path: Set the operation trajectory and initial measurement point that can pass through each numbered coordinate point in the coordinate grid, and use the handheld mobile station to accurately reach each numbered coordinate point along the operation trajectory;

[0012] Step 4: Image acquisition: Arrive at each numbered coordinate point, use the coordinate point as the observation point, take a 360° panoramic photo, and assign a corresponding panoramic photo number to each coordinate point;

[0013] Step 5: Linking: Linking the panoramic photo with the corresponding number to form a touch screen and display screen display system;

[0014] Step 6. Data storage and management: Use cloud storage or local database to store photos and coordinate data to ensure data security and convenience;

[0015] In the step 1, the value of m is determined according to the real-time monitored environmental data:

[0016] When walking: 80,m<m<120,m

[0017] When on a boat tour: 200m<m<300m.

[0018] Preferably, the number of coordinate points in step 1 is calculated as:

[0019]

[0020]

[0021] Where: N x is the number of points along the X-axis; N y is the number of points along the Y axis.

[0022] Preferably, the position of each coordinate point in step 1 is calculated as:

[0023] P(i, j)=(X1+j·m, Y1+i·m)

[0024] Where: (i = 0, 1, 2, ..., N y -1), (j=0,1,2,…,N x -1);

[0025] The final coordinate point number is:

[0026] Each coordinate point is given an ID, which is calculated as follows:

[0027] ID(i, j) = i·N x +j

[0028] Where: ID(i, j) is the coordinate point number of the i-th row and j-th column.

[0029] Preferably, the coordinate difference calculated in step 2 is:

[0030] ΔX=X B -X A

[0031] ΔY=Y B -Y A

[0032] ΔZ=Z B -Z A

[0033] Where: A(X A , Y A , Z A ) are the three-dimensional coordinates X, Y, Z of the known coordinate point A; B(X B , Y B , Z B ) are the three-dimensional coordinates X, Y, and Z of the known coordinate point B; ΔX is the coordinate difference in the X direction; ΔY is the coordinate difference in the Y direction; and ΔZ is the elevation difference in the Z direction.

[0034] Preferably, the proportionality coefficient in step 2 is calculated as:

[0035]

[0036] Among them: the proportional coefficient S is obtained by calculating the distance between two coordinate points, D is an arbitrary standard distance, which is taken as the distance between known points. If the actual distance between the points is known, it is used directly.

[0037] Preferably, the step 2 further comprises calculating the rotation angle:

[0038] The rotation angle θ is obtained by calculating the direction vectors of the two coordinate points to get the angle between them. There are two three-dimensional vectors:

[0039] V A =(X A , Y A , Z A )

[0040] V B =(X B , Y B , Z B )

[0041] Calculate the direction vector:

[0042] D=V B -V A =(X B -X A , Y B -Y A , Z B -Z A )

[0043] Where: D is composed of V B and V A The difference between point A and point B is calculated, indicating the change or difference from point A to point B; V A Indicates point A, V B Indicates point B;

[0044] Compute the magnitude of a vector:

[0045]

[0046] Use the inverse tangent function to calculate the rotation angle:

[0047] Preferably, the operation path of step 3 is determined and calculated as:

[0048] Coordinate system setting: Two-dimensional coordinate system, where each numbered coordinate point is represented by (x i ,y i ), where i is a serial number;

[0049] Trajectory description, describes the trajectory in a parameterized form, using parameter t to represent time or progress:

[0050] P(t)=x(t),y(t)

[0051] Where: P(t) is a vector containing two components: x(t) and y(t); it means that the coordinates of the point are x(t), y(t) as a function of time or parameter t;

[0052] The path is calculated as:

[0053]

[0054] Where: Δt is the movement from point P i To P i+1 the time required;

[0055] In order to ensure that the mobile station reaches each point accurately along the trajectory, an objective function is set to minimize the path deviation:

[0056]

[0057] Among them: |·| represents the norm of the vector, t i is the time corresponding to each coordinate point.

[0058] Preferably, the image acquisition calculation in step 4 is:

[0059] I i =Capture(x i ,y i )

[0060] Among them, Capture(x i ,y i ) means that at the coordinate (x i ,y i ) Panoramic photos taken at the location i ;

[0061] The photo number is calculated as:

[0062] P i =AssignPhotoNumber(I i )

[0063] Where: AssignPhotoNumber(I i ) is a function that converts the panoramic photo I i The corresponding number P i association.

[0064] Preferably, the data structure in step 5 is defined as:

[0065] D=(P1,I1),(P2,I2),…,(P N , IN )

[0066] Where: D is a set, each element of which is a two-tuple containing the panoramic photo number P i and the corresponding panoramic photo I i , the value range of i is from 1 to N.

[0067] Preferably, the data security and convenience evaluation calculation in step 6 is:

[0068] Evaluate Security:S f (D, S)

[0069] Evaluate Convenience:C f (D, S)

[0070] Where: D is a collection of panoramic photos and coordinate data, S is a suitable storage method selected according to the needs; S f (D, S) is a security assessment function that returns the security level of the storage method.

[0071] Compared with the prior art, the present invention provides an image acquisition method for spatial precise positioning of landscape planning, which has the following beneficial effects:

[0072] 1. The present invention clarifies the method of using coordinate point calculation, can accurately define the position of each point, and ensure the uniqueness of each point in the trajectory through a numbering (ID) system, thereby providing a traceable spatial structure, calculation of differences and proportional coefficients: by calculating the coordinate differences ΔX, ΔY and Δz, and introducing the proportional coefficient S, the spatial relationship between two known points is understood in detail, which helps to accurately reflect the change of position during the movement process and reduce the deviation caused by inaccurate position. By calculating the angle between the two coordinate points, the change of the direction vector can be better understood, thereby ensuring the accuracy of the desired direction in the definition of the motion trajectory, further improving the accuracy of data calculation, introducing a parametric method to describe the trajectory, allowing dynamic adjustment of coordinates based on time or progress, adapting to changes in the actual movement process, flexibly responding to uncertainties in the scene, ensuring real-time optimization of the path, setting an objective function that minimizes the path deviation, ensuring that during the movement process, the system can automatically adjust the path to minimize the deviation from the ideal trajectory, and ensuring the fit and accuracy of the activity in a dynamic environment. The beneficial effect of considering more comprehensively when calculating activity characteristics without deviation is achieved.

[0073] 2. The present invention calculates the vector D=V B - VA, the material can directly reflect the position change from point A to point B, providing a basis for subsequent error analysis, monitoring the deviation between the actual movement between the two points and the theoretical expected value, through Such calculations can quantify the true distance between two points and provide a specific measure for correcting errors. The analysis can identify the directional errors that occur during the movement, judge and correct the polarity and tilt of the mobile station, ensure the accuracy of the activity direction, describe the movement trajectory in a parameterized way, and apply the formula The moving trajectory can be dynamically updated and adjusted according to the actual moving progress, thereby reducing the error caused by irregular movement, by minimizing the objective function of path deviation. The system is able to correct and optimize the path at each step, and promptly correct deviations caused by environmental changes or perception errors. During the image acquisition process, the panoramic photos are associated with their corresponding coordinate numbers, which helps to ensure that the data is consistent with the required position during the acquisition process, thereby reducing errors caused by shooting angles or time differences, achieving the beneficial effect of introducing an error correction mechanism to improve the reliability of measurements. BRIEF DESCRIPTION OF THE DRAWINGS

[0074] Figure 1 It is a schematic diagram of the steps of the present invention. DETAILED DESCRIPTION

[0075] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0076] See also Figure 1 , an image acquisition method for spatial precise positioning of landscape planning, comprising the following steps:

[0077] Step 1. Determination of the working area and measuring points: Obtain complete and accurate land vector graphics and coordinate data of the working area and the surrounding area. On the land vector graphics within the working area, draw straight line groups along the due north-south and due east-west directions, with an interval of m, covering the entire working area. According to the real-time monitored environmental data: terrain, vegetation density, weather conditions, dynamically adjust the value of the interval m: when walking, the interval range is 80 meters to 120 meters; when boating, the interval range is 200 meters to 300 meters; draw points at the intersection of the straight line groups, determine the coordinate data of each point, and number them to form a coordinate grid within the working area;

[0078] Step 2: Determine the position of the base station and parameter conversion: Place a base station GPS receiver at known coordinate point A, and inject the coordinates and elevation of the base station into the receiver through a data recording notebook. Set up another GPS receiver at known coordinate point B to receive satellite signals and data transmitted by the base station, and inject the coordinates and elevation of the point. Calculate the conversion parameters relative to the WGS-84 ellipsoid coordinate system, including coordinate difference, elevation difference, scale factor, rotation angle, introduce differential GPS technology, correct the measurement data, improve measurement accuracy, and place the GPS receiver with the obtained parameters on the survey pole, connected to the notebook, and set as a mobile station.

[0079] Step 3: Determine the operation path: Set the operation trajectory and initial measurement point that can pass through each numbered coordinate point in the coordinate grid, and use the handheld mobile station to accurately reach each numbered coordinate point along the operation trajectory;

[0080] Step 4: Image acquisition: Arrive at each numbered coordinate point, use the coordinate point as the observation point, take a 360° panoramic photo, and assign a corresponding panoramic photo number to each coordinate point;

[0081] Step 5: Linking: Linking the panoramic photo with the corresponding number to form a touch screen and display screen display system;

[0082] Step 6. Data storage and management: Use cloud storage or local database to store photos and coordinate data to ensure data security and convenience;

[0083] In step 1, the value of m is determined based on the real-time monitored environmental data:

[0084] When walking: 80,m<m<120,m

[0085] When on a boat tour: 200m<m<300m.

[0086] Specifically, the number of coordinate points in step 1 is calculated as:

[0087]

[0088]

[0089] Where: N x is the number of points along the X-axis; N y is the number of points along the Y axis.

[0090] Specifically, the position of each coordinate point in step 1 is calculated as:

[0091] P(i, j)=(X1+j·m, Y1+i·m)

[0092] Where: (i = 0, 1, 2, ..., N y -1), (j=0,1,2,…,N x -1);

[0093] The final coordinate point number is:

[0094] Each coordinate point is given an ID, which is calculated as follows:

[0095] ID(i, j) = i·N x +j

[0096] Where: ID(i, j) is the coordinate point number of the i-th row and j-th column.

[0097] Specifically, the coordinate difference calculated in step 2 is:

[0098] ΔX=X B -X A

[0099] ΔY=Y B -Y A

[0100] ΔZ=Z B -Z A

[0101] Where: A(X A , Y A , Z A ) are the three-dimensional coordinates X, Y, Z of the known coordinate point A; B(X B , Y B , Z B ) are the three-dimensional coordinates X, Y, and Z of the known coordinate point B; ΔX is the coordinate difference in the X direction; ΔY is the coordinate difference in the Y direction; and ΔZ is the elevation difference in the Z direction.

[0102] Specifically, the proportional coefficient in step 2 is calculated as:

[0103]

[0104] Among them: the proportional coefficient S is obtained by calculating the distance between two coordinate points, D is an arbitrary standard distance, which is taken as the distance between known points. If the actual distance between the points is known, it is used directly.

[0105] Specifically, step 2 also includes calculating the rotation angle:

[0106] The rotation angle θ is obtained by calculating the direction vectors of the two coordinate points to get the angle between them. There are two three-dimensional vectors:

[0107] V A =(X A , Y A, Z A )

[0108] V B =(X B , Y B , Z B )

[0109] Calculate the direction vector:

[0110] D=V B -V A =(X B -X A , Y B -Y A , Z B -Z A )

[0111] Where: D is composed of V B and V A The difference between point A and point B is calculated, indicating the change or difference from point A to point B; V A Indicates point A, V B Indicates point B;

[0112] Compute the magnitude of a vector:

[0113]

[0114] Use the inverse tangent function to calculate the rotation angle:

[0115] Specifically, the calculation of the operation path in step 3 is determined as follows:

[0116] Coordinate system setting: Two-dimensional coordinate system, where each numbered coordinate point is represented by (x i ,y i ), where i is a serial number;

[0117] Trajectory description, describes the trajectory in a parameterized form, using parameter t to represent time or progress:

[0118] P(t)=x(t),y(t)

[0119] Where: P(t) is a vector containing two components: x(t) and y(t); it means that the coordinates of the point are x(t), y(t) as a function of time or parameter t;

[0120] The path is calculated as:

[0121]

[0122] Where: Δt is the movement from point P i To P i+1 the time required;

[0123] In order to ensure that the mobile station reaches each point accurately along the trajectory, an objective function is set to minimize the path deviation:

[0124]

[0125] Among them: |·| represents the norm of the vector, t i is the time corresponding to each coordinate point.

[0126] Specifically, the image acquisition calculation in step 4 is:

[0127] I i =Capture(x i ,y i )

[0128] Among them, Capture(x i ,y i ) means that at the coordinate (x i ,y i ) Panoramic photos taken at the location i ;

[0129] The photo number is calculated as:

[0130] P i =AssignPhotoNumber(I i )

[0131] Where: AssignPhotoNumber(I i ) is a function that converts the panoramic photo I i The corresponding number P i association.

[0132] Specifically, the data structure in step 5 is defined as:

[0133] D=(P1,I1),(P2,I2),…,(P N , I N )

[0134] Where: D is a set, each element of which is a two-tuple containing the panoramic photo number P i and the corresponding panoramic photo I i , the value range of i is from 1 to N.

[0135] Specifically, the data security and convenience evaluation calculation in step 6 is:

[0136] Evaluate Security:S f (D, S)

[0137] Evaluate Convenience:C f (D, S)

[0138] Where: D is a collection of panoramic photos and coordinate data, S is a suitable storage method selected according to the needs; S f (D, S) is a security assessment function that returns the security level of the storage method.

[0139] Embodiment 1

[0140] Calculate the working area and measuring points:

[0141] The coordinates of the lower left corner of the working area (A(0,0))

[0142] Coordinates of the upper right corner of the work area (B(500, 300))

[0143] Select walking tour, set (m=100)

[0144] The calculation steps are as follows:

[0145]

[0146]

[0147] Coordinate point calculation:

[0148] P(0,0)=(0+0·100,0+0·100)=(0,0);

[0149] P(1,0)=(0+0·100,0+1·100)=(0,100);

[0150] And so on, until P(3, 5).

[0151] Through the above steps, a complete coordinate grid is formed in the working area.

[0152] Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the appended claims and their equivalents.

Claims

1. An image acquisition method for spatial precise positioning of landscape planning, characterized in that: The following steps are involved: Step 1. Determination of the working area and measuring points: Obtain complete and accurate land vector graphics and coordinate data of the working area and the surrounding area. On the land vector graphics within the working area, draw straight line groups along the due north-south and due east-west directions, with an interval of m, covering the entire working area. According to the real-time monitored environmental data: terrain, vegetation density, weather conditions, dynamically adjust the value of the interval m: when walking, the interval range is 80 meters to 120 meters; when boating, the interval range is 200 meters to 300 meters; draw points at the intersection of the straight line groups, determine the coordinate data of each point, and number them to form a coordinate grid within the working area; Step 2: Determine the position of the base station and parameter conversion: Place a base station GPS receiver at known coordinate point A, and inject the coordinates and elevation of the base station into the receiver through a data recording notebook. Set up another GPS receiver at known coordinate point B to receive satellite signals and data transmitted by the base station, and inject the coordinates and elevation of the point. Calculate the conversion parameters relative to the WGS-84 ellipsoid coordinate system, including coordinate difference, elevation difference, scale factor, rotation angle, introduce differential GPS technology, correct the measurement data, improve measurement accuracy, and place the GPS receiver with the obtained parameters on the survey pole, connected to the notebook, and set as a mobile station. Step 3: Determine the operation path: Set the operation trajectory and initial measurement point that can pass through each numbered coordinate point in the coordinate grid, and use the handheld mobile station to accurately reach each numbered coordinate point along the operation trajectory; Step 4: Image acquisition: Arrive at each numbered coordinate point, use the coordinate point as the observation point, take a 360° panoramic photo, and assign a corresponding panoramic photo number to each coordinate point; Step 5: Linking: Linking the panoramic photo with the corresponding number to form a touch screen and display screen display system; Step 6. Data storage and management: Use cloud storage or local database to store photos and coordinate data to ensure data security and convenience; In the step 1, the value of m is determined according to the real-time monitored environmental data: When walking: 80 m <m<120,m When cruising: 200m <m<300,m。 2. The image acquisition method for spatial precise positioning of landscape planning according to claim 1, characterized in that: The number of coordinate points in step 1 is calculated as: Where: N x is the number of points along the X-axis; N y is the number of points along the Y axis.

3. The image acquisition method for spatial precise positioning of landscape planning according to claim 1, characterized in that: The position of each coordinate point in step 1 is calculated as: P(i,j)=(X1+j·m,Y1+i·m) Where: (i=0,1,2,…,N y -1), (j=0,1,2,…,N x -1); The final coordinate point number is: Each coordinate point is given an ID, which is calculated as follows: ID(i,j)=i·N x +j Where: ID(i,j) is the coordinate point number of the i-th row and j-th column.

4. The image acquisition method for spatial precise positioning of landscape planning according to claim 1, characterized in that: The coordinate difference calculated in step 2 is: ΔX=X B -X A ΔY=Y B -Y A ΔZ=Z B -Z A Where: A(X A ,Y A ,Z A ) are the three-dimensional coordinates X, Y, Z of the known coordinate point A; B(X B ,Y B ,Z B ) are the three-dimensional coordinates X, Y, and Z of the known coordinate point B; ΔX is the coordinate difference in the X direction; ΔY is the coordinate difference in the Y direction; and ΔZ is the elevation difference in the Z direction.

5. The image acquisition method for spatial precise positioning of landscape planning according to claim 4, characterized in that: The proportionality coefficient in step 2 is calculated as: Among them: the proportional coefficient S is obtained by calculating the distance between two coordinate points, D is an arbitrary standard distance, which is taken as the distance between known points. If the actual distance between the points is known, it is used directly.

6. The image acquisition method for spatial precise positioning of landscape planning according to claim 5, characterized in that: The second step also includes calculating the rotation angle: The rotation angle θ is obtained by calculating the direction vectors of the two coordinate points to get the angle between them. There are two three-dimensional vectors: V A =(X A ,Y A ,Z A ) V B =(X B ,Y B ,Z B ) Calculate the direction vector: D=V B -V A =(X B -X A ,Y B -Y A ,Z B -Z A ) Where: D is composed of V B and V A The difference between point A and point B is calculated, indicating the change or difference from point A to point B; V A Indicates point A, V B Indicates point B; Compute the magnitude of a vector: Use the inverse tangent function to calculate the rotation angle:

7. The image acquisition method for spatial precise positioning of landscape planning according to claim 6, characterized in that: The operation path of step 3 is determined and calculated as follows: Coordinate system setting: Two-dimensional coordinate system, where each numbered coordinate point is represented by (x i ,y i ), where i is a serial number; Trajectory description, describes the trajectory in a parameterized form, using parameter t to represent time or progress: P(t)=x(t),y(t) Where: P(t) is a vector containing two components: x(t) and y(t); it means that the coordinates of the point are x(t), y(t) as a function of time or parameter t; The path is calculated as: Where: Δt is the movement from point P i To P i+1 the time required; In order to ensure that the mobile station reaches each point accurately along the trajectory, an objective function is set to minimize the path deviation: Among them: |·| represents the norm of the vector, t i is the time corresponding to each coordinate point; MinimizeJ is the objective function for minimizing the path deviation.

8. The image acquisition method for spatial precise positioning of landscape planning according to claim 1, characterized in that: The image acquisition calculation in step 4 is: I i =Capture(x i ,y i ) Among them, Capture(x i ,y i ) indicates that at coordinate (x i ,y i ) Panoramic photos taken at the location i ; The photo number is calculated as: P i =AssignPhotoNumber(I i ) Where: AssignPhotoNumber(I i ) is a function that converts the panoramic photo I i The corresponding number P i association.

9. The image acquisition method for spatial precise positioning of landscape planning according to claim 1, characterized in that: The data structure in step 5 is defined as: D=(P1,I1),(P2,I2),…,(P N ,I N ) Where: D is a set, each element of which is a two-tuple containing the panoramic photo number P i and the corresponding panoramic photo I i , the value range of i is from 1 to N.

10. The image acquisition method for spatial precise positioning of landscape planning according to claim 1, characterized in that: The data security and convenience evaluation calculation in step 6 is: Evaluate Security:S f (D,S) Evaluate Convenience:C f (D,S) Where: D is a collection of panoramic photos and coordinate data, S is a suitable storage method selected according to the needs; S f (D,S) is a security assessment function that returns the security level of the storage method.

Images