专利名称：房屋变动判断方法及房屋变动判断程序的制作方法
技术领域：
本发明涉及房屋变动判断方法及房屋变动判断程序。
背景技术：
在下述专利文献中提出了利用空中测量技术进行固定资产的变动判断等的方法。以往，提出了从飞机等对不同的两个时期取得高分辨率的图像及DSM(Digital Surface Model 数字表面模型)数据、将这两个时期的数据按照像素比较、来判断该两个时期间的 房屋的变动的方法。具体而言，在使用从飞机摄影的图像的房屋变动判断中，进行过去得到的图像与 新得到的图像的对位(registration，配准)，对于对应的位置的像素求出新旧的亮度值等 的像素数据的差，将该差与规定的阈值比较来判断变化的有无。并且，基于变化的像素群提 取图像的变化区域，判断房屋的消失及新建的房屋的变动。此外，在使用DSM的房屋变动判断中，例如将飞行位置不同的多个航空图像进行 对照，通过基于同一地物的该多个航空图像中的观察方式的差异求出该地物的高度而生成 DSM。并且，与上述新旧图像的比较同样，按照像素比较新旧的DSM的数据，提取DSM的变化 区域。专利文献1 日本特开2004-117245号公报专利文献2 日本特开2007-3244号公报像素单位的解析处理由于数据的噪声或配准的偏差等的影响，有难以适当地设定 阈值、容易发生对变化的漏检测及误检测的问题。特别是，在高分辨率图像的对比中，有难 以以高精度进行配准、难以进行高精度的变化检测的问题。此外，不容易将像素单位的变化有无的判断结果以房屋单位汇总。例如，即使利用 网孔单位的解析、或房屋多边形等的辅助数据，但由于有增改建，所以房屋单位下的变动判 断并不容易。此外，在如改建的情况下图像数据等的变化较小，难以提取变化区域，所以有 变动判断较困难的问题。
发明内容
有关本发明的房屋变动判断方法，具有凸部区域提取步骤，根据对判断对象区域 内的各点从上空取得的、包含地物的地面表层的标高数据，对与作为判断对象的房屋相应 的大小的上述地面表层的每一个凸部，决定包含该凸部的凸部区域；轮廓提取步骤，根据上 述标高数据，从上述各凸部区域中提取上述房屋的轮廓；目标定义步骤，根据所提取的上述 轮廓及该轮廓内的上述标高数据，定义具有上述房屋的三维信息的房屋目标；以及变动判 断步骤，对关于不同的两个时期分别得到的上述房屋目标进行对比，检测上述两个时期间 的变化，并根据该变化判断房屋变动。可以构成为，上述凸部区域提取步骤具有归一化处理步骤，该归一化处理步骤从 上述地面表层的标高减去不包含上述地物的地表的标高，并生成上述地物的高度数据，上述凸部区域提取步骤根据上述地物的高度数据，决定上述凸部区域。在本发明的优选方式中可以构成为，上述凸部区域提取步骤通过基于分水岭法的 区域分割，求出上述凸部区域。此外，在本发明的优选方式中可以构成为，面积比较步骤，根据上述轮廓计算上述 房屋目标的面积，判断对比对象的上述房屋目标相互的上述面积的异同；轮廓比较步骤，判 断对比对象的上述房屋目标相互的上述轮廓的异同；平均高度比较步骤，根据上述三维信 息计算上述房屋目标的平均高度，并判断对比对象的上述房屋目标相互的上述平均高度的 异同；高度直方图比较步骤，根据上述三维信息求出上述房屋目标的高度直方图，并判断对 比对象的上述房屋目标相互的上述高度直方图的异同；以及判断步骤，在上述各比较步骤 中都判断为没有上述两个时期间的差异的情况下判断为无房屋变动，另一方面，在某一个 比较步骤中判断为有差异的情况下判断为有房屋变动。有关本发明的房屋变动判断程序，使计算机执行凸部区域提取步骤，根据对判断 对象区域内的各点从上空取得的、包含地物的地面表层的标高数据，对与作为判断对象的 房屋相应的大小的上述地面表层的每一个凸部，决定包含该凸部的凸部区域；轮廓提取步 骤，根据上述标高数据，从上述各凸部区域中提取上述房屋的轮廓；目标定义步骤，根据所 提取的上述轮廓及该轮廓内的上述标高数据，定义具有上述房屋的三维信息的房屋目标； 以及变动判断步骤，对关于不同的两个时期分别得到的上述房屋目标进行对比，检测上述 两个时期间的变化，并根据该变化判断房屋变动。可以构成为，上述凸部区域提取步骤具有归一化处理步骤，该归一化处理步骤从 上述地面表层的标高减去不包含上述地物的地表的标高，并生成上述地物的高度数据，上 述凸部区域提取步骤根据上述地物的高度数据，决定上述凸部区域。在本发明的优选方式中可以构成为，上述凸部区域提取步骤基于通过分水岭法的 区域分割，求出上述凸部区域。此外，在本发明的优选方式中可以构成为，上述变动判断步骤具有面积比较步 骤，根据上述轮廓计算上述房屋目标的面积，并判断对比对象的上述房屋目标相互的上述 面积的异同；轮廓比较步骤，判断对比对象的上述房屋目标相互的上述轮廓的异同；平均 高度比较步骤，根据上述三维信息计算上述房屋目标的平均高度，并判断对比对象的上述 房屋目标相互的上述平均高度的异同；高度直方图比较步骤，根据上述三维信息求出上述 房屋目标的高度直方图，并判断对比对象的上述房屋目标相互的上述高度直方图的异同； 以及判断步骤，在上述各比较步骤中都判断为没有上述两个时期间的差异的情况下判断为 无房屋变动，另一方面，在某一个比较步骤中判断为有差异的情况下判断为有房屋变动。根据本发明，由于基于在DSM中出现的起伏提取对应于各个房屋的房屋目标，以 房屋目标单位进行两个时期相互的对比，所以房屋移动的判断变得容易，并且精度提高。进 而，房屋目标除了具有房屋所位置的区域的信息以外，还具有基于DSM得到的高度的信息， 由于三维地掌握房屋，所以如改建那样的以往难以判别的变化也能够高精度地检测。
图1是表示作为本发明的实施方式的房屋变动判断方法的概况的处理流程的流 程图。
图2是房屋目标提取处理的概况的流程图。图3是说明对于DSM的标准化处理的示意图。图4是说明区域分割处理的示意性的俯视图。图5是表示变化检测处理的概况的流程图。图6是表示比较判断处理的一例的流程图。
具体实施例方式以下，基于附图对本发明的实施方式(以下称作实施方式)进行说明。图1是表示作为实施方式的房屋变动判断方法的概况的处理的流程的流程图。对 关于判断对象区域在不同的时期Tl、T2得到的DSM(DigitalSurface Model 数字表面模 型)数据文件2、4分别进行房屋目标提取处理S6、S8。房屋目标提取处理S6、S8基于DSM 所表示的地面表层面的起伏，决定包含相应于房屋的凸部的凸部区域，从各凸部区域提取 房屋的轮廓，定义具有该轮廓与该轮廓内的各点的高度的信息的房屋目标。关于各时期T1、 T2得到的房屋目标的提取结果存储在房屋目标文件10、12中。对于时期Tl、T2各自的房 屋目标的信息，进行变化检测处理S14。变化检测处理S14将时期T1的房屋目标与时期T2 的房屋目标进行对照、对比，检测两时刻之间的房屋目标的变化，并将其结果作为变动判断 结果文件16输出。这里，在已经进行了房屋目标提取处理、且房屋目标文件已经存在的情况下，能够 将其用在变化检测处理S14中而省略房屋目标提取处理。图2是房屋目标提取处理S6、S8的概况的流程图。DSM被输入到归一化处理S20 中，在归一化处理S20中，进行将包含在DSM中的DTM的影响除去的归一化处理，图3是说 明归一化处理的示意图。图3(a)是地面的示意性的垂直剖视图。DSM30表示房屋32、树 木34等的地物表面的水平坐标-标高、和在地物间露出的地表面36的水平坐标-标高。 这里，DSM30中的房屋的标高为其下方的地表面的标高与房屋的高度的合计。归一化处理 S20是从该DSM30的值减去DTM (Digital Terrain Model 数字地形模型)40的值、并生成 NDSM(Normalized Digital Surface Model 归一化数字表面模型)的处理。图3 (b)是示 意地表示对应于图3(a)的地面而得到的NDSM42的垂直剖视图。在该NDSM42中，地物(房 屋32a、树木34a)具有距地表的高度信息，另一方面，地表36a基本上为对应于高度0的高 度。例如，作为DTM40可以使用根据DSM通过形态滤波等各种滤波处理推测的DTM。或者也 可以使用有关判断对象区域的已有的DTM。在NDSM中，房屋以外的地物也具有有意义的高度，特别是，树木34a等的植被能 够成为相当于房屋的高度，所以在房屋判断时成为噪声的主要原因。植被除去处理S44是 从NDSM42除去有关植被的部分的处理，判断植被的区域，并通过将该区域的NDSM数据改写 为0而进行遮蔽。房屋的屋檐及屋顶遍及较大的水平距离展开着平滑的面，相对于此，在植 被区域中，以比较小的空间周期发生高度变动。例如，可以利用该差异判断植被区域，并从 NDSM42中除去。具体而言，可以将NDSM数据的相邻点间的梯度的分散超过规定的阈值的区 域判断为植被区域。此外，也可以基于因高度的变动而在从飞机等摄影的图像中产生的纹 理(texture)进行判断。进而，也可以使用多光谱图像，利用树木的叶子强烈地反射近红外 光的性质来检测植被区域。例如，可以将使用红(R)成分的像素值％的近红外(IR)成分的像素值 Dik 由下式计算出的 NDVI (Normalized DifferenceVegetation Index 归一化植 被指数)为规定的阈值以上的区域判断为植被区域。NDVI = (Die-De)/(Die+De)接着，将用除去了植被的NDSM表示的判断对象区域分割为分别各包含1个在该 NDSM中出现的凸部的凸部区域(区域分割处理S46)。在从上述DSM向NDSM的变换处理 中，以比房屋的展宽大的范围表现的地表的凹凸基本上被除去，留下相应于房屋的大小的 凸部。由此，能够对每一个房屋32a设定凸部区域。图4是说明区域分割处理的示意性的 俯视图。图4(a)是将植被除去处理S44后的判断对象区域中的NDSM的凹凸简单地进行二 值化表示的图，涂黑的区域50对应于凸部的位置。区域分割处理S46将判断对象区域分割 为多个凸部区域52，以使各凸部相互分离。作为该区域分割的方法，例如可以使用分水岭 (Watershed)法。以上，对通过归一化处理S20、植被除去处理S44、区域分割处理S46提取的凸部区 域52，进行提取房屋的轮廓的轮廓提取处理S60。在轮廓提取处理中，例如基于用亮度或颜 色显示从飞机等取得的地面表层的标高的图像(标高图像)62的变化，提取能够构成边缘 的点群。另外，也可以代替标高图像而使用从飞机等摄影的地面表层的图像来进行轮廓提 取。边缘提取可以使用各种边缘滤波器进行。对于构成边缘的点群进行提取线段而生成线 段图像的处理。例如，执行进行跟踪及分割处理的滤波，将该点群分割为能够分别作为1个 汇集而捕捉的单位的点的集合，并且去除没有构成足够大的汇集而能够掌握为噪声的点。 然后，例如通过Hough (哈夫)变换来提取线圈。接着，进行根据线段图像提取房屋的轮廓 的处理。例如，可以对线段图像使用Perceptual Grouping(知觉编组)技术提取房屋的轮 廓。对提取的每一个轮廓定义房屋目标(目标定义处理S64)。房屋目标将轮廓及该轮 廓内的NDSM数据的信息建立对应。以上使用图2等说明的处理作为房屋目标提取处理S6、S8对时期Tl、T2各自的 DSM进行。并且，对时期T1、T2各自的房屋目标的信息进行变化检测处理S14。图5是表示变化检测处理S14的概况的流程图。设从时期T1的判断对象区域提 取的房屋目标的个数是m个，将该房屋目标用01 (i) (1 ^ i ^ m)表示，设从时期T2的判断 对象区域提取的房屋目标的个数是n个，将该房屋目标用02 (j) (1 ^ j ^ n)表示。例如，首先对先于时期T2的时期T1的任意的房屋目标01 (i)，搜索作为比较对象 的02 (j)(步骤S72 S80)。这里，将处于与该01 (i)基本上相同的位置的02 (j)作为比较 对象。如果有作为比较对象的0(j) (S76中“Y”的情况)，则将表示比较对象的检测的有无 的标志F设定为1(S82)，进行与该02(j)的比较判断(S84)，将其结果向判断结果文件输出 并记录(S86)。在通过搜索在所有的02 (j)之中没有找到作为比较对象的02 (j)的情况下 (S88中标志F = 0的情况下)，向判断结果文件记录对应于01⑴的房屋在时期T2中不存 在、即消失了内容的变动判断结果(S90)。对所有的01 (i)尝试该比较判断(S92，S94)，然后，这次与对上述01(i)进行的 比较判断同样地进行对02 (j)的比较判断。具体而言，首先1个个地选择房屋目标02(j) (S100, S124)，搜索作为该02(j)的比较对象的01(i)(步骤S102 S110)。如果有与该 02 (j)相同位置的01(i)(在S106中是“Y”的情况)，则进行比较判断(S114，S116)。另一方面，在没有作为比较对象的01 (i)的情况下(在S118中标志F = 0的情况下)，向判断 结果文件记录对应于02 (j)的房屋在时期T1中不存在、S卩，将02 (j)是新建筑等内容的变 动判断结果(S120)。另外，这里比较判断处理S114和比较判断处理S84可以为基本上相 同的处理。在此情况下，对于作为比较对象的01 (i)及02 (j)的组合，在两个比较判断处理 S84、S114中处理重复，所以例如可以构成为，在比较判断处理S114中省略具体的处理、在 向判断结果文件的输出处理S116中、根据需要而复制比较判断处理S84中的对应的判断结 果，作为对于该02 (j)的结果。图6是表示比较判断处理S84、S114的一例的流程图，对时期T1、T2各自的房屋目 标中作为比较对象的一对房屋目标01、02判断两个时期间的变动的有无。首先，基于房屋目标具有的轮廓信息，计算01、02各自的面积，比较它们的大小 (S130)。如果在两者的面积上存在有意义的差，则判断在两目标表示的房屋之间有变动 (S140)。在面积上没有差的情况下，计算01、02相互的轮廓的一致度，如果该一致度小于 规定的基准值(S132)，则判断为有变动(S140)。在轮廓的一致度是基准值以上的情况下，基于房屋目标具有的、轮廓内的NDSM的 值进行关于房屋的高度的比较。作为该高度信息的比较，例如首先计算01、02各自的平均 高度，将它们进行比较(S134)。如果在两者的平均高度上存在有意义的差，则判断为有变动 (S140)。在平均高度上没有差的情况下，进一步对01、02各自在轮廓内的多个点具有的 NDSM数据进行详细的对比。例如，求出01、02各自的高度(NDSM数据)的直方图，如果它 们的相似度小于规定的基准值(S136)，则判断为有变动(S140)。此外，也可以将轮廓内的 各点的01、02的高度的差(A)的平方和、或A的绝对值超过规定的阈值的像素数等作为 指标来判断变动的有无。在如上所述的各种判断基准中都没有在01、02中发现差异的情况下，判断为两者 没有变动(S150)。另外，在图6所示的比较判断处理中，主要以变动的有无的检测为目的，如果在中 途的判断基准中检测到差异，则对于其后的判断基准省略判断。相对于此，也可以构成为对 所有的判断基准进行判断、基于它们的判断结果综合性地判断变动的有无、并且将它们的 判断结果也记录到判断结果文件中。此外，也可以根据多个判断基准的判断结果判别变动 的种类，将其记录到判断结果文件中。例如，作为变动的种类，除了房屋的单纯的消失、新建 以外，可以基于房屋目标具有的高度信息判断相同地方的改建、或将一层建筑的上方整体 或部分地做成二层建筑的增建、或相反的减建。上述房屋变动判断方法可以作为在计算机 上执行的程序实现。计算机通过执行该程序，基于由飞机等取得的DSM数据文件2、4，执行 房屋目标提取处理S6、S8，生成房屋目标文件10、12，此外，对于保存在房屋目标文件10、12 中的两个时期的房屋目标执行变化检测处理S14，生成变动判断结果文件16。此外，在上述实施方式中，说明了使用从飞机摄影的图像及由该图像生成的DSM 的方式，但也可以使用从高分辨率卫星取得的图像及DSM。此外，DSM也可以是从飞机等向 地面照射激光而取得的。工业实用性
使用有关本发明的房屋变动判断方法及房屋变动判断程序得到的房屋变动的信 息在固定资产税的评价、及灾害场所的确定等方面成为有用的信息。
权利要求
一种房屋变动判断方法，其特征在于，具有凸部区域提取步骤，根据对于判断对象区域内的各点从上空取得的、包含地物的地面表层的标高数据，对与作为判断对象的房屋相应的大小的上述地面表层的每一个凸部，决定包含该凸部的凸部区域；轮廓提取步骤，根据上述标高数据，从上述各凸部区域中提取上述房屋的轮廓；目标定义步骤，根据所提取的上述轮廓及该轮廓内的上述标高数据，定义具有上述房屋的三维信息的房屋目标；以及变动判断步骤，对关于不同的两个时期分别得到的上述房屋目标进行对比，检测上述两个时期间的变化，并根据该变化判断房屋变动。
2.如权利要求1所述的房屋变动判断方法，其特征在于，上述凸部区域提取步骤具有归一化处理步骤，该归一化处理步骤从上述地面表层的标 高减去不包含上述地物的地表的标高，并生成上述地物的高度数据，上述凸部区域提取步骤根据上述地物的高度数据，决定上述凸部区域。
3.如权利要求1所述的房屋变动判断方法，其特征在于，上述凸部区域提取步骤通过基于分水岭法的区域分割，求出上述凸部区域。
4.如权利要求1所述的房屋变动判断方法，其特征在于， 上述变动判断步骤具有面积比较步骤，根据上述轮廓计算上述房屋目标的面积，并判断对比对象的上述房屋 目标相互的上述面积的异同；轮廓比较步骤，判断对比对象的上述房屋目标相互的上述轮廓的异同； 平均高度比较步骤，根据上述三维信息计算上述房屋目标的平均高度，并判断对比对 象的上述房屋目标相互的上述平均高度的异同；高度直方图比较步骤，根据上述三维信息求出上述房屋目标的高度直方图，并判断对 比对象的上述房屋目标相互的上述高度直方图的异同；以及判断步骤，在上述各比较步骤中都判断为没有上述两个时期间的差异的情况下判断为 无房屋变动，另一方面，在某一个比较步骤中判断为有差异的情况下判断为有房屋变动。
5.—种房屋变动判断程序，其特征在于，使计算机执行凸部区域提取步骤，根据对于判断对象区域内的各点从上空取得的、包含地物的地面 表层的标高数据，对与作为判断对象的房屋相应的大小的上述地面表层的每一个凸部，决 定包含该凸部的凸部区域；轮廓提取步骤，根据上述标高数据，从上述各凸部区域中提取上述房屋的轮廓； 目标定义步骤，根据所提取的上述轮廓及该轮廓内的上述标高数据，定义具有上述房 屋的三维信息的房屋目标；以及变动判断步骤，对关于不同的两个时期分别得到的上述房屋目标进行对比，检测上述 两个时期间的变化，并根据该变化判断房屋变动。
6.如权利要求5所述的房屋变动判断程序，其特征在于，上述凸部区域提取步骤具有归一化处理步骤，该归一化处理步骤从上述地面表层的标 高减去不包含上述地物的地表的标高，并生成上述地物的高度数据，上述凸部区域提取步骤根据上述地物的高度数据，决定上述凸部区域。
7.如权利要求5所述的房屋变动判断程序，其特征在于，上述凸部区域提取步骤通过基于分水岭法的区域分割，求出上述凸部区域。
8.如权利要求5所述的房屋变动判断程序，其特征在于， 上述变动判断步骤具有面积比较步骤，根据上述轮廓计算上述房屋目标的面积，并判断对比对象的上述房屋 目标相互的上述面积的异同；轮廓比较步骤，判断对比对象的上述房屋目标相互的上述轮廓的异同； 平均高度比较步骤，根据上述三维信息计算上述房屋目标的平均高度，并判断对比对 象的上述房屋目标相互的上述平均高度的异同；高度直方图比较步骤，根据上述三维信息求出上述房屋目标的高度直方图，并判断对 比对象的上述房屋目标相互的上述高度直方图的异同；以及判断步骤，在上述各比较步骤中都判断为没有上述两个时期间的差异的情况下判断为 无房屋变动，另一方面，在某一个比较步骤中判断为有差异的情况下判断为有房屋变动。
全文摘要
提高基于由飞机等取得的图像等进行的房屋的变动判断的精度。从由飞机等取得的数字表面模型(DSM)的标高值减去地形的标高，生成归一化DSM(NDSM)。对在NDSM中出现的、相应于房屋的大小的每一个凸部，将判断对象区域分割为多个凸部区域。从各凸部区域中提取房屋的轮廓，根据该轮廓和轮廓内的NDSM数据，定义具有上述房屋的三维信息的房屋目标。对关于不同的两个时期分别得到的房屋目标进行对比，检测该两个时期间的变化，根据该变化判断房屋变动。
文档编号G01C11/00GK101828093SQ20088011191
公开日2010年9月8日 申请日期2008年10月20日 优先权日2007年10月19日
发明者岛村秀树, 朱林 申请人:株式会社博思科