一种基于随机切片的社区尺度城市形态数据批量生成方法技术

本发明专利技术公开了一种基于随机切片的社区尺度城市形态数据批量生成方法，包括以下步骤：S1：采用遥感数据、人工选定或数据校核的方式获取城市下垫面的原始数据；S2：将原始数据裁剪成规则的形状，并转换成矩阵形式文件，根据随机数生成切片子矩阵；S3：获取子矩阵中的坐标值，将城市下垫面的数据对应插值在坐标网格中，输出切片文件；S4：根据切片文件中城市下垫面的形态类型，筛选符合要求的切片文件。本发明专利技术能够在质量控制的基础上，实现任意指定尺度的城市形态随机切片和生成，不仅实现了下垫面要素“随机”组合的要求，而且通过程序控制，实现了数据批量生成和处理，能够满足进一步开展城市形态定量分析的需要。城市形态定量分析的需要。城市形态定量分析的需要。

【技术实现步骤摘要】
一种基于随机切片的社区尺度城市形态数据批量生成方法


[0001]本专利技术涉及城市下垫面数据分析领域，具体涉及一种基于随机切片的社区尺度城市形态数据批量生成方法。

技术介绍

[0002]城市下垫面主要由房屋建筑、道路、绿地、广场、水域等组成，由人工构筑物与自然地势复合构建而形成。由于城市地表构筑物类型复杂，由复杂城市下垫面构成的城市形态复杂程度远高于流域尺度。在进行城市形态分析、城市内涝数值模拟等研究时，常需要对不同尺度下的城市形态进行随机提取和定量分析，以获取更多城市形态定量分析信息。采用原始城市下垫面遥感数据下载或人工获取等方式，通常难以达到“随机”和“批量”的要求，为此，提出一种城市形态数据随机生成方法，对于科学统计和分析城市形态特征，具有十分重要的意义。

技术实现思路

[0003]针对现有技术的上述不足，本专利技术提供了一种基于随机切片的社区尺度城市形态数据批量生成方法，能够在质量控制的基础上，实现任意指定尺度(例如社区尺度)的城市形态随机切片和生成。
[0004]为达到上述专利技术目的，本专利技术所采用的技术方案为：
[0005]提出一种基于随机切片的社区尺度城市形态数据批量生成方法，其包括以下步骤：
[0006]S1：采用遥感数据、人工选定或数据校核的方式获取城市下垫面的原始数据；
[0007]S2：将原始数据裁剪成规则的形状，并转换成矩阵形式的文件，生成原始数据矩阵，获取原始数据矩阵中行数和列数的范围，并在行数和列数的范围内生成随机数，根据随机数生成切片子矩阵；
[0008]S3：获取切片子矩阵中的坐标值，设置关于坐标值的变量X和Y，生成基于X和Y的坐标网格，坐标网格按矩阵的行堆叠成一个二维矩阵，将城市下垫面的数据对应插值在坐标网格中，输出切片文件；
[0009]S4：根据切片文件中城市下垫面的形态类型，筛选符合要求的切片文件。
[0010]进一步地，步骤S2包括：
[0011]S21：将原始数据裁剪成规则的形状，并转换成矩阵形式的文件，生成原始数据矩阵；
[0012]S22：设定子矩阵的行数X
b
和列数Y
b
，根据原始数据矩阵的列数Y
a
和行数X
a
，分别计算能够随机生成的子矩阵的最大行数X
max
和最大列数Y
max
：
[0013]X
max
＝X
a
‑
X
b
；
[0014]Y
max
＝Y
a
‑
Y
b
；
[0015]S23：确定关于行坐标的随机数的取值范围(0，X
max
)、列坐标的随机数的取值范围
(0，Y
max
)；
[0016]S24：在范围(0，X
max
)和范围(0，Y
max
)内分别随机生成一个行坐标值X1和列坐标值Y1；
[0017]S25：将随机坐标(X1，Y1)作为随机生成的子矩阵的左上角顶点，再以左上角顶点为起点向右侧取列数Y
b
、向下侧取行数X
b
，得到切片子矩阵。
[0018]S26：重复步骤S24
‑
S25，在范围(0，X
max
)和范围(0，Y
max
)内分别随机生成不同的行坐标值和列坐标值，得到若干个切片子矩阵。
[0019]进一步地，步骤S3包括：
[0020]S31：根据切片子矩阵中的每个坐标参数，根据城市下垫面的原始数据加载每个坐标参数对应的原始数据，原始数据包括高程数据和用地类型数据，并将原始数据赋值、保存在对应的坐标参数上，原始数据作为坐标参数对应的变量；
[0021]S32：设定单元网格的大小，将切片子矩阵划分成x
×
y个网格，研究区切片的分辨率取决于每个网格的大小，将切片子矩阵按照网格生成一维数组，保存一维数组中对应的变量；
[0022]S33：基于X、Y生成的空间坐标，将一维数组按行列堆叠成一个二维数组，得到新的坐标数组G
xy
：
[0023][0024]其中，G为坐标数组G
xy
的坐标；
[0025]S34：提取切片子矩阵中任一坐标参数作为样本坐标(c,d)，遍历坐标数组G
xy
中的每个坐标作为插值坐标(a,b)，计算样本坐标与每个插值坐标之间的距离S：
[0026][0027]S35：将样本坐标与每个插值坐标之间的距离值进行比较，筛选出距离最小值S
min
对应的样本坐标，并将距离最小值S
min
与设定的距离阈值S
阈值
进行比较：
[0028]若S
min
≤S
阈值
，则将该样本坐标对应的原始数据插值到最近的插值坐标上，执行步骤S36；
[0029]否则，执行步骤S37；
[0030]S36：判断插值的原始数据中的高程数据，在城市形态解析中考虑道路、城市或其他三种用地类型，并进行插值判断分析是否存在异常值的情况：
[0031]若存在空值，即插值为
‑
9999，表示缺乏数据或者无效数据，通常对该插值进行特殊处理，判断为异常值，并进行数据清洗或排除；
[0032]删除该插值坐标上的原始数据，提取步骤S35中比较出的第二小的距离值对应的样本坐标，返回执行步骤S35；
[0033]否则，插值坐标插值成功；
[0034]S37：返回步骤S34，在切片子矩阵中重新提取一个样本坐标，执行步骤S34
‑
S36；
[0035]S38：直到坐标数组G
xy
中的每个坐标均插值成功后，输出城市形态数据的切片文件。
[0036]进一步地，步骤S4包括：
[0037]利用GIS软件对城市形态数据切片进行栅格统计，分别统计出城市形态数据切片的数据中要素的占比，根据各项要素的占比大小进行筛选，筛选的规则包括：对道路面积占比进行上下限控制、对房屋面积占比进行上下限控制、对道路和房屋占比之和进行上下限控制；在质量控制基础上，可以进一步筛选出出符合要求的城市形态数据切片，完成基本数据的获取。
[0038]进一步地，步骤S1包括：
[0039]采用遥感数据、人工选定或数据校核的方式获取城市下垫面的原始数据：
[0040]遥感数据：通过下载公开的遥感数据进行下垫面的解析或者直接下载已经解析过的土地覆盖数据进行实验数据准备；通过遥感影像识别城市道路等要素的轮廓，根据实验不同需求选定不同城市下垫面条件；
[0041]人工选定：通过制定的规则选择特定的城市下垫面研究区域；
[0042]数据校核：对已有数据进行校核，核实数据可用性，确保本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于随机切片的社区尺度城市形态数据批量生成方法，其特征在于，包括以下步骤：S1：采用遥感数据、人工选定或数据校核的方式获取城市下垫面的原始数据；S2：将原始数据裁剪成规则的形状，并转换成矩阵形式的文件，生成原始数据矩阵，获取原始数据矩阵中行数和列数的范围，并在行数和列数的范围内生成随机数，根据随机数生成切片子矩阵；S3：获取切片子矩阵中的坐标值，设置关于坐标值的变量X和Y，生成基于X和Y的坐标网格，坐标网格按矩阵的行堆叠成一个二维矩阵，将城市下垫面的数据对应插值在坐标网格中，输出切片文件；S4：根据切片文件中城市下垫面的形态类型，筛选符合要求的切片文件。2.根据权利要求1所述的基于随机切片的社区尺度城市形态数据批量生成方法，其特征在于，步骤S2包括：S21：将原始数据裁剪成规则的形状，并转换成矩阵形式的文件，生成原始数据矩阵；S22：设定子矩阵的行数X
b
和列数Y
b
，根据原始数据矩阵的列数Y
a
和行数X
a
，分别计算能够随机生成的子矩阵的最大行数X
max
和最大列数Y
max
：X
max
＝X
a
‑
X
b
；Y
max
＝Y
a
‑
Y
b
；S23：确定关于行坐标的随机数的取值范围(0，X
max
)、列坐标的随机数的取值范围(0，Y
max
)；S24：在范围(0，X
max
)和范围(0，Y
max
)内分别随机生成一个行坐标值X1和列坐标值Y1；S25：将随机坐标(X1，Y1)作为随机生成的子矩阵的左上角顶点，再以左上角顶点为起点向右侧取列数Y
b
、向下侧取行数X
b
，得到切片子矩阵。S26：重复步骤S24
‑
S25，在范围(0，X
max
)和范围(0，Y
max
)内分别随机生成不同的行坐标值和列坐标值，得到若干个切片子矩阵。3.根据权利要求1所述的基于随机切片的社区尺度城市形态数据批量生成方法，其特征在于，步骤S3包括：S31：根据切片子矩阵中的每个坐标参数，根据城市下垫面的原始数据加载每个坐标参数对应的原始数据，原始数据包括高程数据和用地类型数据，并将原始数据赋值、保存在对应的坐标参数上，原始数据作为坐标参数对应的变量；S32：设定单元网格的大小，将切片子矩阵划分成x
×
...

【专利技术属性】
技术研发人员：梅超，石虹远，刘家宏，高希超，王佳，宋天旭，董莉榕，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：