A gis-based township land space planning base processing method and system

Abstract

The application provides a GIS-based township land space planning base processing method and system, which comprises the following steps: extracting the past land class code from the three-adjustment plot vector to obtain direct conversion vectors and indirect conversion vectors, and updating the direct conversion vectors into first updated three-adjustment plot vectors; updating the indirect conversion vectors into second updated three-adjustment plot vectors and to-be-planned vectors through three types of data; updating the to-be-planned vectors into third updated three-adjustment plot vectors, and then obtaining final three-adjustment data. Compared with the manual processing method, the direct conversion vectors are converted into the first updated three-adjustment plot vectors, which avoids repeated calculation of most stable land blocks, and the coding of the indirect conversion vectors can be automatically assigned, which improves the base processing efficiency. Moreover, the three types of data with dynamic and real-time characteristics are introduced to make up for the tardiness of the three-adjustment data and improve the accuracy.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of data processing, in particular to a GIS-based township land space planning base processing method and system. BACKGROUND

[0002] Land space planning is an arrangement of land space development and protection in space and time in a certain area, which is a guide for spatial development and a blueprint for sustainable development. It is the basic basis for various development, protection and construction activities. Land space planning aims to realize "multi-plan integration", and integrates main function zone planning, land use planning, urban and rural planning and other spatial planning into unified land space planning.

[0003] Township-level land space planning is the most basic "level" of the "five-level three-type" land space planning compilation system, and is a "chessboard" for realizing the accurate landing of various land space elements. Its purpose is to regulate the use of specific plots.

[0004] The current township land space planning base processing relies on manual comparison of three survey data and planning land classification standards, but using only three survey data, the data is out of date, the dynamic is insufficient, the accuracy is low, and the manual comparison method is extremely time-consuming and laborious, and the data processing efficiency is low. SUMMARY

[0005] In view of the deficiencies of the prior art, the purpose of the present application is to provide a GIS-based township land space planning base processing method and system, which aims to solve the technical problems of insufficient dynamic, low accuracy and low data processing efficiency in the prior art by manually comparing three survey data with planning land classification standards.

[0006] In order to achieve the above-mentioned purpose, in a first aspect, the present application provides a GIS-based township land space planning base processing method, comprising the following steps:

[0007] The three survey data is divided into a plurality of three survey polygon vectors, the previous land class code is extracted from the three survey polygon vector, the plurality of three survey polygon vectors are divided into direct conversion vectors and indirect conversion vectors based on the previous land class code, the direct conversion vectors are assigned with a first update code to obtain a first update three survey polygon vector;

[0008] Obtain land supply data, POI data and remote sensing image data corresponding to the three survey data, and align the three survey data, the land supply data, the POI data and the remote sensing image data to obtain three survey conversion data, land supply conversion data, POI conversion data and image conversion data;

[0009] The indirect conversion vector is given a second update code and a blank identification through the land supply conversion data, the POI conversion data and the image conversion data, so as to obtain a second update three-adjustment polygon vector and a vector to be planned respectively;

[0010] The vector to be planned is grid divided to obtain a plurality of sub-grid vectors, the sub-grid vectors are given a development suggestion identification to form an update sub-grid vector, a plurality of the update sub-grid vectors are combined into a third update three-adjustment polygon vector, and the first update three-adjustment polygon vector, the second update three-adjustment polygon vector and the third update three-adjustment polygon vector are combined into final three-adjustment data.

[0011] Further, the step of dividing a plurality of the three-adjustment polygon vectors into direct conversion vectors and indirect conversion vectors based on the past land class code and giving the direct conversion vectors a first update code comprises:

[0012] The past land class code is compared with a preset dictionary, and the preset dictionary comprises a plurality of stable land class codes;

[0013] The three-adjustment polygon vector corresponding to the past land class code matched with the stable land class code is selected as the direct conversion vector, and the stable land class code is selected as the first update code, and the first update code is associated with the direct conversion vector;

[0014] The three-adjustment polygon vector corresponding to the past land class code not matched with any stable land class code is selected as the indirect conversion vector.

[0015] Further, the step of aligning the three-adjustment data, the land supply data, the POI data and the remote sensing image data to obtain three-adjustment conversion data, land supply conversion data, POI conversion data and image conversion data comprises:

[0016] The three-adjustment data is mapped to a geodetic coordinate system to obtain three-adjustment conversion data;

[0017] A plurality of reference points are selected in the three-adjustment conversion data, and reference points corresponding to the reference points are selected in the land supply data, the POI data and the remote sensing image respectively;

[0018] A transformation matrix is constructed through the reference points and the reference points, and the land supply data, the POI data and the remote sensing image data are mapped into land supply conversion data, POI conversion data and image conversion data through the transformation matrix.

[0019] Further, after the step of constructing a transformation matrix by the reference point and the reference point, mapping the land supply data, the POI data and the remote sensing image data into land supply conversion data, POI conversion data and image conversion data by the transformation matrix, further comprising:

[0020] Selecting a calibration point conversion coordinate in the land supply conversion data, the POI conversion data and the image conversion data respectively, and obtaining a calibration point real coordinate corresponding to the calibration point conversion coordinate;

[0021] Judging whether the alignment processing is accurate by the calibration point conversion coordinate and the calibration point real coordinate.

[0022] Further, the land supply conversion data includes a plurality of land property vectors, the POI conversion data includes a plurality of POI vectors, and the image conversion data includes a plurality of image vectors. The step of assigning a second update code and a blank identifier to the indirect conversion vector by the land supply conversion data, the POI conversion data and the image conversion data includes:

[0023] Obtaining a land use property code corresponding to the land property vector, selecting an indirect conversion vector corresponding to the land property vector as a first indirect vector, and selecting the remaining indirect conversion vector as a first to-be-identified vector, selecting the land use property code as a second update code and associating it to the first indirect vector;

[0024] Selecting a first to-be-identified vector with a POI vector as a second indirect vector, and selecting the remaining first to-be-identified vector as a second to-be-identified vector, obtaining a land use type code by the POI vector, selecting the land use type code as a second update code and associating it to the second indirect vector;

[0025] Generating a landform type code and a blank identifier by the landform information of the image vector corresponding to the second to-be-identified vector, selecting the landform type code as a second update code, and associating the second update code and the blank identifier to the second to-be-identified vector respectively.

[0026] Further, the step of obtaining a land use type code by the POI vector includes:

[0027] Obtaining the vector type of different POI vectors in the second indirect vector, classifying a plurality of the vector types to obtain a plurality of category groups;

[0028] Obtaining the category proportion of a plurality of the category groups respectively, selecting the category group with the highest category proportion as a target category, and generating a land use type code by the target category.

[0029] Further, the formula for obtaining the category proportion is:

[0030] ,

[0031] wherein, represents the category proportion of the i-th category group, represents the number of vector types in the i-th category group, represents the total number of vector types.

[0032] Further, the step of grid dividing the to-be-planned vector to obtain a plurality of sub-grid vectors comprises:

[0033] obtaining a minimum circumscribed rectangular frame corresponding to the to-be-planned vector;

[0034] respectively setting a plurality of evenly distributed horizontal division lines and a plurality of evenly distributed vertical division lines based on the width and height of the minimum circumscribed rectangular frame;

[0035] separating the to-be-planned vector into a plurality of sub-grid vectors through the horizontal division lines and the vertical division lines.

[0036] Further, the development suggestion identifier comprises a high suitability identifier, a medium suitability identifier and a restricted development identifier, and the step of assigning a development suggestion identifier to the sub-grid vector comprises:

[0037] obtaining a flat area and a stable foundation area corresponding to the sub-grid vector, and obtaining a suitability index based on the flat area and the stable foundation area;

[0038] if the suitability index is greater than or equal to a first index threshold, a high suitability identifier is generated, if the suitability index is less than or equal to a second index threshold, a restricted development identifier is generated, and if the suitability index is between the first index threshold and the second index threshold, a medium suitability identifier is generated.

[0039] In a second aspect, the embodiments of the present application provide a GIS-based township land space planning base number processing system, which is applied to the GIS-based township land space planning base number processing method as described in the first aspect above, and the system comprises:

[0040] a first processing module, configured to separate three-adjustment data into a plurality of three-adjustment polygon vectors, extract a past land class code from the three-adjustment polygon vectors, separate a plurality of the three-adjustment polygon vectors into direct conversion vectors and indirect conversion vectors based on the past land class code, assign a first update code to the direct conversion vectors to obtain first update three-adjustment polygon vectors;

[0041] The mapping module is configured to acquire land supply data, POI data and remote sensing image data corresponding to the three-inventory data, and perform alignment processing on the three-inventory data, the land supply data, the POI data and the remote sensing image data to obtain three-inventory conversion data, land supply conversion data, POI conversion data and image conversion data.

[0042] The second processing module is configured to assign a second update code and a blank identifier to the indirect conversion vector through the land supply conversion data, the POI conversion data and the image conversion data, to respectively acquire a second update three-inventory polygon vector and a to-be-planned vector.

[0043] The third processing module is configured to perform grid division on the to-be-planned vector to acquire a plurality of sub-grid vectors, assign a development suggestion identifier to the sub-grid vectors to form update sub-grid vectors, combine a plurality of the update sub-grid vectors into a third update three-inventory polygon vector, and combine the first update three-inventory polygon vector, the second update three-inventory polygon vector and the third update three-inventory polygon vector into final three-inventory data.

[0044] In a third aspect, an embodiment of the present application provides a computer, including a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the GIS-based township land space planning base processing method according to the first aspect when executing the computer program.

[0045] In a fourth aspect, an embodiment of the present application provides a storage medium having a computer program stored thereon, and the computer program is executable by a processor to implement the GIS-based township land space planning base processing method according to the first aspect.

[0046] Compared with the prior art, the beneficial effects of the present application are that: by introducing the land supply data, the POI data and the remote sensing image data to the second update coding of the indirect conversion vector, the dynamic and real-time of the three types of data compensate for the lag of the three adjustment data, and the accuracy of the base number processing is improved; by performing the alignment processing, the consistency of the three types of data introduced and the three adjustment data in space is ensured, and the accuracy of the base number processing is further improved; by assigning the development suggestion identifier to the to-be-planned vector, a reasonable development suggestion is provided for the missing area of the three adjustment data, and the rationality of the subsequent development based on the final three adjustment data after the base number processing is ensured; compared with the manual processing method, by assigning the three adjustment data area to the direct conversion vector and the indirect conversion vector, and assigning the first update coding to the direct conversion vector, repeated calculation of most stable plots is avoided, the efficiency of the base number processing is improved, and the coding of the indirect conversion vector can be automatically assigned, and then the final three adjustment data is obtained, the processing efficiency is further improved, and data errors caused by manual processing can be avoided. BRIEF DESCRIPTION OF DRAWINGS

[0047] Figure 1 The flowchart of the GIS-based township land space planning base number processing method in the first embodiment of the present application;

[0048] Figure 2 The structural block diagram of the GIS-based township land space planning base number processing system in the second embodiment of the present application;

[0049] The following specific embodiments will further illustrate the present application in conjunction with the above drawings. DETAILED DESCRIPTION

[0050] In order to facilitate the understanding of the present application, the present application will be described more fully below with reference to the related drawings. The drawings show several embodiments of the present application. However, the present application can be realized in many different forms and is not limited to the embodiments described herein. On the contrary, the purpose of providing these embodiments is to make the disclosure of the present application more thorough and comprehensive.

[0051] It should be noted that when an element is referred to as being "fixedly attached" to another element, it can be directly on the other element or there can be an intervening element. When an element is referred to as being "connected" to another element, it can be directly connected to the other element or intervening elements can be present. The terms "vertical", "horizontal", "left", "right" and similar expressions used herein are for illustrative purposes only.

[0052] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used in the description herein is for describing particular embodiments only and is not intended to be limiting of the application. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.

[0053] Referring to Figure 1 The first embodiment of the application provides a GIS-based township land space planning base processing method, comprising the following steps:

[0054] S10: separating the triennial data into a plurality of triennial polygon vectors, extracting the past land class code from the triennial polygon vectors, separating the plurality of triennial polygon vectors into direct conversion vectors and indirect conversion vectors based on the past land class code, and assigning a first update code to the direct conversion vectors to obtain a first update triennial polygon vector;

[0055] The triennial data is composed of a plurality of triennial polygon vectors, and the triennial polygon vector contains a past polygon code, a past land class code, and a past town and village attribute.

[0056] The step S10 comprises:

[0057] S110: comparing the past land class code with a preset dictionary, wherein the preset dictionary comprises a plurality of stable land class codes;

[0058] The stable land class code corresponds to stable land classes such as farmland and forest land, and its purpose is single and long-term unchanged, while the land class code corresponding to dynamic land classes such as industry and commerce is not in the preset dictionary.

[0059] S120: selecting the triennial polygon vector corresponding to the stable land class code matching the past land class code as a direct conversion vector, and selecting the stable land class code as a first update code, and associating the first update code with the direct conversion vector;

[0060] S130: selecting the triennial polygon vector corresponding to the past land class code which does not match any stable land class code as an indirect conversion vector;

[0061] If the past land class code is 0101 (farmland), the stable land class code 0101 corresponding to the past land class code exists in the preset dictionary, in which case, the three-adjustment patch vector corresponding to the past land class code is selected as the direct conversion vector, and the code 0101 is assigned to the direct conversion vector again. It should be noted that after the assignment of the first updated code is completed, the direct conversion vector is associated with the past land class code and the first updated code at the same time, which is not overwritten, but a new association key is generated. If the past land class code is 0601 (garden), the stable land class code corresponding to the past land class code does not exist in the preset dictionary, and the three-adjustment patch vector is selected as an indirect conversion vector.

[0062] S20: Obtain land supply data, POI data and remote sensing image data corresponding to the three-adjustment data, and perform alignment processing on the three-adjustment data, the land supply data, the POI data and the remote sensing image data to obtain three-adjustment conversion data, land supply conversion data, POI conversion data and image conversion data.

[0063] The land supply conversion data includes a plurality of land property vectors, the POI conversion data includes a plurality of POI vectors, and the image conversion data includes a plurality of image vectors. It can be understood that the land supply data, the POI data and the remote sensing image data also respectively include the foregoing content, only the coordinate system is different. The three-adjustment conversion data still contains the indirect conversion vector, only the spatial coordinates are different.

[0064] Specifically, the step S20 includes:

[0065] S210: Map the three-adjustment data to a geodetic coordinate system to obtain three-adjustment conversion data;

[0066] S220: Select a plurality of reference points in the three-adjustment conversion data, and select reference points corresponding to the reference points in the land supply data, the POI data and the remote sensing image respectively;

[0067] The number of reference points is greater than or equal to three. In this embodiment, the reference points are the vertices of the three-adjustment patch vectors corresponding to the road intersections. The reference points are targeted, and the points at the same positions in the other three types of data are selected as the reference points.

[0068] S230: Construct a transformation matrix through the reference points and the reference points, and map the land supply data, the POI data and the remote sensing image data to land supply conversion data, POI conversion data and image conversion data through the transformation matrix;

[0069] Specifically, the reference points and the reference points correspond to each other to form a matching point pair, and by substituting the horizontal and vertical coordinates of the matching point pair into the initial affine transformation formula, six transformation parameters are obtained, and then by the determined six transformation parameters, a final affine transformation formula is formed to complete data conversion according to the final affine transformation formula. It should be noted that three final affine transformation formulas are generated for three types of data. Affine transformation is now widely used and will not be described here.

[0070] Preferably, the step S20 further comprises:

[0071] S240: selecting a calibration point conversion coordinate in the land supply conversion data, the POI conversion data and the image conversion data respectively, and obtaining a calibration point real coordinate corresponding to the calibration point conversion coordinate;

[0072] In this embodiment, the calibration point conversion coordinate is the position coordinate of the concrete well lid at the road intersection.

[0073] S250: determining whether the alignment processing is accurate by the calibration point conversion coordinate and the calibration point real coordinate;

[0074] Specifically, the position residual between the calibration point conversion coordinate and the calibration point real coordinate is calculated, and the calculation formula of the position residual is:

[0075]

[0076] wherein, represents the position residual, , represents the horizontal coordinate and the vertical coordinate of the calibration point conversion coordinate, respectively, , represents the horizontal coordinate and the vertical coordinate of the calibration point real coordinate, respectively.

[0077] After obtaining the position residual, the position residual is compared with the distance threshold value. If the position residual is greater than the distance threshold value, it is determined that the alignment processing is not accurate, and the alignment processing needs to be performed again. If the position residual is less than the distance threshold value, it is determined that the alignment processing is accurate. In some embodiments, by selecting a plurality of calibration point conversion coordinates, a more accurate verification effect is achieved. At this time, there are a plurality of position residuals, and the position residual mean value is obtained by averaging a plurality of position residuals, and then the alignment processing is verified. In this embodiment, the distance threshold value is 0.3m, which ensures the accuracy requirement (0.5m) of the three adjustment data.

[0078] ​S30: assigning a second update code and a blank mark to the indirect conversion vector by the land supply conversion data, the POI conversion data and the image conversion data, to obtain a second update three-adjustment polygon vector and a to-be-planned vector respectively;

[0079] The step 30 comprises:

[0080] S310: obtaining a land use property code corresponding to the land property vector, selecting an indirect conversion vector corresponding to the land property vector as a first indirect vector and selecting the remaining indirect conversion vectors as first to-be-identified vectors, selecting the land use property code as a second update code and associating the second update code to the first indirect vector;

[0081] The land property vector represents the land use property of a batch, which has the highest priority, and in the subsequent development process, the land use property needs to be consistent with the batch, so when the land property vector corresponding to the indirect conversion vector exists, the land use property code is directly generated according to the land use property.

[0082] S320: selecting a first to-be-identified vector with a POI vector as a second indirect vector and selecting the remaining first to-be-identified vectors as second to-be-identified vectors, obtaining a land use type code by the POI vector, selecting the land use type code as a second update code and associating the second update code to the second indirect vector;

[0083] It should be noted that the POI vector is essentially a POI point, which represents the building type in the region, and in the case where the land property vector is missing, it is identified whether the POI point exists in the first to-be-identified vector, and if the POI point exists, the building types corresponding to all the POI points are counted. Specifically,

[0084] The vector types of different POI vectors in the second indirect vector are obtained, a plurality of the vector types are classified to obtain a plurality of category groups, the category proportions of the plurality of category groups are obtained respectively, the category group with the highest category proportion is selected as a target category group, and a land use type code is generated by the target category group;

[0085] The formula for obtaining the category proportion is:

[0086] ,

[0087] Among them, represents the category proportion of the i-th category group, represents the number of vector types in the i-th category group, represents the total number of vector types. Assuming that the category group with the highest category proportion is a commercial facility category, a land use type code corresponding to the commercial facility category is generated.

[0088] S330: generating a landform type code and a blank mark through the landform information of the image vector corresponding to the second to-be-identified vector, selecting the landform type code as a second update code, and respectively associating the second update code and the blank mark to the second to-be-identified vector;

[0089] It can be understood that according to the image vector, the landform image corresponding to the second to-be-identified vector can be obtained, and then it is determined whether it is a landform type such as a construction site or a road, and then the landform type code is assigned, and if the image vector represents no building, the blank mark is generated. It should be noted that the previous land type code, the first update code, the second update code, the land use property code, the land use type code, and the landform type code all represent the same code, such as 0101 (farmland) as described above, that is, the substantive meaning is the same, and the name is only used to distinguish that it exists in different processing nodes. And the second to-be-identified vector can only be associated with the second update code or the blank mark, and there is no case of simultaneous association of the second update code and the blank mark. By specifying the assignment rules of the three types of data and the indirect conversion vector, a conflict resolution mechanism after the introduction of multi-source data is provided, ensuring the rationality of subsequent development.

[0090] S40: performing grid division on the to-be-planned vector to obtain a plurality of sub-grid vectors, assigning a development suggestion mark to the sub-grid vector to form an updated sub-grid vector, combining a plurality of updated sub-grid vectors into a third update tri-adjustment polygon vector, and combining the first update tri-adjustment polygon vector, the second update tri-adjustment polygon vector and the third update tri-adjustment polygon vector into final tri-adjustment data;

[0091] The development suggestion mark includes a high suitability mark, a medium suitability mark, and a restricted development mark.

[0092] The step S40 includes:

[0093] S410: obtaining a minimum circumscribed rectangle corresponding to the to-be-planned vector;

[0094] S420: setting a plurality of interval-distributed horizontal division lines and a plurality of interval-distributed vertical division lines based on the width and height of the minimum circumscribed rectangle;

[0095] S430: separating the to-be-planned vector into a plurality of sub-grid vectors through the horizontal division lines and the vertical division lines;

[0096] In this embodiment, the size of the sub-grid vector is 100m*100m.

[0097] S440: Obtain a flat area and a stable foundation area corresponding to the sub-grid vector, and obtain a suitability index based on the flat area and the stable foundation area;

[0098] In the sub-grid vector, the area of the region with a slope less than 15% is extracted to obtain the flat area, and the areas of regions with different slopes are obtained. In this embodiment, the sub-grid vector is converted into a slope raster vector corresponding to a pixel point by using the ArcGIS Slope tool. The region with a slope less than 15% is marked as 1, and the remaining regions are marked as 0. Then, the total number of regions marked as 1 is counted to obtain the flat area. The area of the region with a PH value of 6-8 is extracted to obtain the stable foundation area. In this embodiment, the PH values of different regions in the sub-grid vector are defined by introducing soil survey data.

[0099] S450: If the suitability index is greater than or equal to a first index threshold, a high suitability identifier is generated; if the suitability index is less than or equal to a second index threshold, a restricted development identifier is generated; and if the suitability index is between the first index threshold and the second index threshold, a medium suitability identifier is generated.

[0100] The suitability index = (flat area + stable foundation area) / total area of the sub-grid vector. It should be noted that if the sum of the flat area and the stable foundation area is greater than the total area of the sub-grid vector, the suitability index is forcibly defined as 100%. In this embodiment, the first index threshold is 70%, and the second index threshold is 40%.

[0101] By introducing the land supply data, the POI data and the remote sensing image data, the second update coding of the indirect conversion vector is given, the dynamic and real-time of the three types of data is used to compensate for the lag of the three types of data, and the accuracy of the base number processing is improved. By performing the alignment processing, the consistency of the three types of data introduced and the three types of data in space is ensured, and the accuracy of the base number processing is further improved. By giving the development suggestion identifier to the to-be-planned vector, a reasonable development suggestion is provided for the missing area of the three types of data, and the rationality of the subsequent development based on the final three types of data after the base number processing is ensured. Compared with the manual processing method, by converting the three types of data area into the direct conversion vector and the indirect conversion vector, and giving the first update coding to the direct conversion vector, repeated calculation of most stable land blocks is avoided, the efficiency of the base number processing is improved, and the coding of the indirect conversion vector can be completely automated, and then the final three types of data are obtained, the processing efficiency is further improved, and the data error caused by manual processing is avoided.

[0102] Please refer to Figure 2The second embodiment of the present application provides a GIS-based township land space planning base processing system. The system is applied to the GIS-based township land space planning base processing method in the above embodiment, which has been described above and will not be repeated here. As used below, the terms "module", "unit", "sub-unit", and the like can be a combination of software and / or hardware that implements a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, hardware, or a combination of software and hardware can also be implemented and conceived.

[0103] The system comprises:

[0104] The first processing module 10 is configured to separate the three-level data into a plurality of three-level patch vectors, extract a past land class code from the three-level patch vectors, divide the plurality of three-level patch vectors into direct conversion vectors and indirect conversion vectors based on the past land class code, and assign a first update code to the direct conversion vectors to obtain a first update three-level patch vector.

[0105] The first processing module 10 comprises:

[0106] The first unit is configured to compare the past land class code with a preset dictionary, and the preset dictionary comprises a plurality of stable land class codes.

[0107] The second unit is configured to select a three-level patch vector corresponding to a past land class code matching a stable land class code as a direct conversion vector, and select the stable land class code as a first update code, and associate the first update code with the direct conversion vector.

[0108] The third unit is configured to select a three-level patch vector corresponding to a past land class code that does not match any stable land class code as an indirect conversion vector.

[0109] The mapping module 20 is configured to obtain land supply data, POI data, and remote sensing image data corresponding to the three-level data, and perform alignment processing on the three-level data, the land supply data, the POI data, and the remote sensing image data to obtain three-level conversion data, land supply conversion data, POI conversion data, and image conversion data.

[0110] The mapping module 20 comprises:

[0111] The fourth unit is configured to map the three-level data to a geodetic coordinate system to obtain three-level conversion data.

[0112] The fifth unit is configured to select a plurality of reference points in the three-level conversion data, and select reference points corresponding to the reference points in the land supply data, the POI data, and the remote sensing image, respectively.

[0113] a sixth unit configured to construct a transformation matrix by using the reference point and the reference point, and to map the land supply data, the POI data and the remote sensing image data into land supply conversion data, POI conversion data and image conversion data by using the transformation matrix;

[0114] Preferably, the mapping module 20 further comprises:

[0115] a seventh unit configured to select a calibration point conversion coordinate from the land supply conversion data, the POI conversion data and the image conversion data respectively, and to obtain a calibration point real coordinate corresponding to the calibration point conversion coordinate;

[0116] an eighth unit configured to determine whether the alignment processing is accurate by using the calibration point conversion coordinate and the calibration point real coordinate;

[0117] a second processing module 30 configured to assign a second update code and a blank identifier to the indirect conversion vector by using the land supply conversion data, the POI conversion data and the image conversion data, so as to obtain a second update three-regulation polygon vector and a to-be-planned vector respectively;

[0118] The second processing module 30 comprises:

[0119] a ninth unit configured to obtain a land use property code corresponding to the land property vector, select an indirect conversion vector corresponding to the land property vector as a first indirect vector, select the remaining indirect conversion vectors as first to-be-identified vectors, select the land use property code as a second update code, and associate the second update code to the first indirect vector;

[0120] a tenth unit configured to select a first to-be-identified vector with a POI vector as a second indirect vector, select the remaining first to-be-identified vectors as second to-be-identified vectors, obtain a land use type code by using the POI vector, select the land use type code as a second update code, and associate the second update code to the second indirect vector;

[0121] The tenth unit is further configured to obtain a vector type of different POI vectors in the second indirect vector, perform classification processing on a plurality of the vector types, so as to obtain a plurality of category groups, obtain a category proportion of the plurality of category groups respectively, select a category group with a highest category proportion as a target category group, and generate a land use type code by using the target category group;

[0122] an eleventh unit configured to generate a landform type code and a blank identifier by using landform information of an image vector corresponding to the second to-be-identified vector, select the landform type code as a second update code, and associate the second update code and the blank identifier to the second to-be-identified vector respectively.

[0123] The third processing module 40 is configured to perform grid division on the to-be-planned vector to obtain a plurality of sub-grid vectors, assign development suggestion identifiers to the sub-grid vectors to form updated sub-grid vectors, combine a plurality of the updated sub-grid vectors into a third updated tri-adjustment polygon vector, and combine the first updated tri-adjustment polygon vector, the second updated tri-adjustment polygon vector and the third updated tri-adjustment polygon vector into final tri-adjustment data.

[0124] The third processing module 40 comprises:

[0125] The twelfth unit is configured to obtain a minimum circumscribed rectangle corresponding to the to-be-planned vector.

[0126] The thirteenth unit is configured to set a plurality of evenly distributed horizontal division lines and a plurality of evenly distributed vertical division lines based on the width and height of the minimum circumscribed rectangle.

[0127] The fourteenth unit is configured to separate the to-be-planned vector into a plurality of sub-grid vectors by the horizontal division lines and the vertical division lines.

[0128] The fifteenth unit is configured to obtain a flat area and a stable foundation area corresponding to the sub-grid vector, and obtain a suitability index based on the flat area and the stable foundation area.

[0129] The sixteenth unit is configured to generate a high suitability identifier if the suitability index is greater than or equal to a first index threshold, generate a restricted development identifier if the suitability index is less than or equal to a second index threshold, and generate a medium suitability identifier if the suitability index is between the first index threshold and the second index threshold.

[0130] The application further provides a computer comprising a memory, a processor and a computer program stored on the memory and executable on the processor, and the processor implements the GIS-based township land space planning base processing method as described in the above technical solution when executing the computer program.

[0131] The application further provides a storage medium having a computer program stored thereon, and the computer program is executed by a processor to implement the GIS-based township land space planning base processing method as described in the above technical solution.

[0132] In the description of the present specification, the description of the terms "one embodiment", "some embodiments", "an example", "a specific example", or "some examples" and the like means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present application. In the present specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

[0133] The above-described embodiments only express several implementation manners of the present application, which are described in a more specific and detailed manner, but cannot be understood as a limitation on the patent scope of the present application. It should be noted that, for those of ordinary skill in the art, several modifications and improvements can be made without departing from the concept of the present application, which are all within the protection scope of the present application. Therefore, the protection scope of the patent of the present application should be subject to the appended claims.

Claims

1. A GIS-based township land space planning base processing method, characterized in that, The method comprises the following steps: The three-grade data is divided into a plurality of three-grade patch vectors, a previous land class code is extracted from the three-grade patch vectors, a plurality of the three-grade patch vectors are divided into direct conversion vectors and indirect conversion vectors based on the previous land class code, and the direct conversion vectors are assigned a first update code to obtain first update three-grade patch vectors; The step of dividing a plurality of the three-grade patch vectors into direct conversion vectors and indirect conversion vectors based on the previous land class code and assigning the direct conversion vectors a first update code comprises: The previous land class code is compared with a preset dictionary, and the preset dictionary comprises a plurality of stable land class codes; The three-grade patch vector corresponding to the previous land class code matched with the stable land class code is selected as the direct conversion vector, and the stable land class code is selected as the first update code, and the first update code is associated with the direct conversion vector; The three-grade patch vector corresponding to the previous land class code not matched with any stable land class code is selected as the indirect conversion vector; Land supply data, POI data and remote sensing image data corresponding to the three-grade data are obtained, and the three-grade data, the land supply data, the POI data and the remote sensing image data are aligned to obtain three-grade conversion data, land supply conversion data, POI conversion data and image conversion data; The indirect conversion vector is assigned a second update code and a blank identifier through the land supply conversion data, the POI conversion data and the image conversion data to obtain a second update three-grade patch vector and a to-be-planned vector respectively; The to-be-planned vector is grid divided to obtain a plurality of sub-grid vectors, the sub-grid vectors are assigned a development suggestion identifier to form update sub-grid vectors, a plurality of the update sub-grid vectors are combined into a third update three-grade patch vector, and the first update three-grade patch vector, the second update three-grade patch vector and the third update three-grade patch vector are combined into final three-grade data.

2. The GIS-based township land space planning base processing method according to claim 1, characterized in that, The step of aligning the three-grade data, the land supply data, the POI data and the remote sensing image data to obtain three-grade conversion data, land supply conversion data, POI conversion data and image conversion data comprises: The three-grade data is mapped to a geodetic coordinate system to obtain three-grade conversion data; A plurality of reference points are selected in the three-grade conversion data, and reference points corresponding to the reference points are selected in the land supply data, the POI data and the remote sensing image respectively; A transformation matrix is constructed through the reference points and the reference points, and the land supply data, the POI data and the remote sensing image data are mapped into land supply conversion data, POI conversion data and image conversion data through the transformation matrix.

3. The GIS-based township land space planning base processing method according to claim 2, characterized in that, After the step of constructing a transformation matrix through the reference points and the reference points, and mapping the land supply data, the POI data and the remote sensing image data into land supply conversion data, POI conversion data and image conversion data through the transformation matrix, the method further comprises: Respective in the land supply conversion data, the POI conversion data and the image conversion data, select a calibration point conversion coordinate, and obtain a calibration point real coordinate corresponding to the calibration point conversion coordinate; Determine whether the alignment processing is accurate through the calibration point conversion coordinate and the calibration point real coordinate.

4. The GIS-based township land space planning base processing method according to claim 1, characterized in that, The land supply conversion data includes a plurality of land property vectors, the POI conversion data includes a plurality of POI vectors, and the image conversion data includes a plurality of image vectors. The step of assigning a second update code and a blank identifier to the indirect conversion vector through the land supply conversion data, the POI conversion data and the image conversion data includes: Obtain a land use property code corresponding to the land property vector, select an indirect conversion vector corresponding to the land property vector as a first indirect vector, and select the remaining indirect conversion vectors as first to-be-identified vectors. Select the land use property code as a second update code and associate it to the first indirect vector; Select the first to-be-identified vector with a POI vector as a second indirect vector, and select the remaining first to-be-identified vectors as second to-be-identified vectors. Obtain a land use type code through the POI vector, select the land use type code as a second update code, and associate it to the second indirect vector; Generate a landform type code and a blank identifier through the landform information of the image vector corresponding to the second to-be-identified vector. Select the landform type code as a second update code, and associate the second update code and the blank identifier to the second to-be-identified vector respectively.

5. The GIS-based township land space planning base processing method according to claim 4, characterized in that, The step of obtaining a land use type code through the POI vector includes: Obtain the vector types of different POI vectors in the second indirect vector, and perform classification processing on a plurality of the vector types to obtain a plurality of category groups; Obtain the category proportion of a plurality of the category groups respectively, select the category group with the highest category proportion as a target category, and generate a land use type code through the target category. 6.The GIS-based township land space planning base processing method according to claim 5, characterized in that, The formula for obtaining the category proportion is: ， wherein, represents the proportion of the class of the i-th class group, represents the number of vector types in the i-th class group, represents the total number of vector types.

7. The GIS-based township land space planning base processing method according to claim 1, characterized in that, The step of performing grid division on the to-be-planned vector to obtain a plurality of sub-grid vectors includes: Obtain a minimum circumscribed rectangle corresponding to the to-be-planned vector; Set a plurality of horizontally distributed division lines and a plurality of vertically distributed division lines based on the width and height of the minimum circumscribed rectangle; Separate the to-be-planned vector into a plurality of sub-grid vectors through the horizontally distributed division lines and the vertically distributed division lines. 8.The GIS-based township land space planning base processing method according to claim 1, characterized in that, The development suggestion identifier includes a high suitability identifier, a medium suitability identifier and a restricted development identifier. The step of assigning a development suggestion identifier to the sub-grid vector includes: Obtain a flat area and a stable foundation area corresponding to the sub-grid vector, and obtain a suitability index based on the flat area and the stable foundation area; If the suitability index is greater than or equal to a first index threshold, generate a high suitability identifier. If the suitability index is less than or equal to a second index threshold, generate a restricted development identifier. If the suitability index is between the first index threshold and the second index threshold, generate a medium suitability identifier.

9. A GIS-based township land space planning base data processing system applied to the GIS-based township land space planning base data processing method according to any one of claims 1-8, characterized in that, The system comprises: A first processing module is configured to separate the three-level data into a plurality of three-level patch vectors, extract a past land class code from the three-level patch vectors, divide the three-level patch vectors into direct conversion vectors and indirect conversion vectors based on the past land class code, and assign a first updated code to the direct conversion vectors to obtain first updated three-level patch vectors; A mapping module is configured to obtain land supply data, POI data and remote sensing image data corresponding to the three-level data, and perform alignment processing on the three-level data, the land supply data, the POI data and the remote sensing image data to obtain three-level conversion data, land supply conversion data, POI conversion data and image conversion data; A second processing module is configured to assign a second updated code and a blank identification to the indirect conversion vectors through the land supply conversion data, the POI conversion data and the image conversion data to obtain second updated three-level patch vectors and to-be-planned vectors respectively; A third processing module is configured to perform grid division on the to-be-planned vectors to obtain a plurality of sub-grid vectors, assign a development suggestion identification to the sub-grid vectors to form updated sub-grid vectors, combine a plurality of the updated sub-grid vectors into third updated three-level patch vectors, and combine the first updated three-level patch vectors, the second updated three-level patch vectors and the third updated three-level patch vectors into final three-level data.

Images