技术领域本发明特别涉及地图设计领域，并且特别涉及用于设计等值域图的计算机实现的方法、计算机程序、以及系统。
背景技术：
在地理信息系统(GIS)中经常使用等值域图。其涉及在地图上表示统计变量的问题，地图的每个区域上的统计变量的值由颜色强度表示。正在被显示的统计变量可以具有不同的特性(人口密度、居民汽车数量、降雨量、疾病比例)。在此上下文中，在这种表示中，在被选择的强度集合中进行选择是最初工作。用于这种选择的背景技术的两种主要算法是Jenks/Fisher的自然分解优化(如在Jenks的文章“TheDataModelConceptinStatisticalMapping”，1967，以及Fisher的文章“Ongroupingformaximumhomogeneity”，1958中所描述的)和头-尾分解(如在Lin和Yue的文章“Acomparisonstudyonnaturalandhead/tailbreaksinvolvingdigitalelevationmodels”，2013，以及Jiang和Bing的文章“Head/tailbreaksforvisualizingthefractalorscalingstructureofgeographicfeatures”，2014或“Ht-indexforquantifyingthefractalorscalingstructureofgeographicfeatures”中所描述的)。这些解决方案实际上基于群集分析算法。使用群集分析的等值域图系统的其他示例包括文献US20050278182A1、US8412419B1、以及US8510080B2。等值域图的设计以一种方式涉及量化(即，根据预定距离，由预定的值的集合中的值中最接近的一个值来替换输入值)。实际上，地图的区域被分组到一起，并且针对这样的区域的统计变量的值相应地由该分组的代表值替换。这样的整体框架涉及更一般的群集分析领域。群集分析涉及将一组对象划分成分组(称为群集)的任务，以使得在每个分组中数据是相似的(参考Jain等的文章，“DataClustering:AReview”)。其作为数据挖掘(参考文章Chen等的文章，“Datamining:anoverviewfromadatabaseperspective”)、机器学习(参考Murphy的书“MachineLearning,AProbabilisticPerspective”)、以及大规模搜索(参考Goodrum的文章“ImageInformationRetrieval:AnOverviewofCurrentResearch”)中的核心问题而出现。群集分析是用于量化的重要工具：将中心分配各每个群集，就具有了由将每个点量化成其群集的中心构成的简单量化。K均值聚类问题是群集分析的最著名的问题，并且作为用于脉冲编码调制的技术于1957年在贝尔实验室由StuartLloyd引入。Lloyd算法将p维点的集合作为输入，并且输出对这些点的划分，该划分的目标是使“总失真”最小化。该算法只是试探性的(其不提供最优聚类)。但实际上我们不能期望精确的算法，这是因为K均值聚类问题在非一维的情况下是NP困难的。当前，Lloyd算法仍然被广泛使用。已经提出了多种变型(参考J.A.Hartigan(1975),“Clusteringalgorithms”,JohnWiley&Sons,Inc.”)。一维应用特别重要。用于该问题的最著名的算法中的一个实际上就是在1967年中开发的上文引用的Jenks自然分解优化(参考Jenks的书,“TheDataModelConceptinStatisticalMapping”,国际制图年鉴(InternationalYearbookofCartography))，并且出于制图的目的而被引入，如上文所述。与Lloyd算法一样，其只是试探性的。在2011年，Wang和Song开发了被称为CK均值的精确算法(参考Wang和Song的文章“Optimalk-meansClusteringinOneDimensionbyDynamicProgramming”)。该算法是文献US1543036A的基础。该算法以时间O(K*n^2)运行，其中，K是所请求的群集的数量，并且n是实数的数量。甚至最近的(2013年)，MaartenHilferink已经开发出更加有效的算法，并且提供了该算法的实现。该实现实际上专用于制图，更准确地说，专用于等值域图，然而，该算法的唯一文档是维基百科页面(Fisher的自然分解分类，在优先权日期在下述URL可访问：http://wiki.objectvision.nl/index.php/Fisher％27s_Natural_Breaks_Classification)。然而，所有这些现有的方法都受到限制，因为它们不产生最优的K均值群集，或者太慢。在此上下文中，仍然存在对用于设计等值域图的改进的解决方案的需求。
技术实现要素：
因此提供了一种设计等值域图的计算机实现的方法。该方法包括提供地图和表示地图的相应区域上的统计变量的一定数量的数值的步骤。方法还包括针对预定的群集数量来计算数值的最优K均值聚类的步骤。该计算步骤包括将应用到阶等于数值的数量的正方矩阵的线性时间行最小值搜索算法迭代与预定的群集数量相对应的次数。方法还包括将相应的着色分配给计算出的聚类的每个群集的步骤。方法还包括在提供相应的数值的地图的所有区域处，应用被分配给群集的着色的步骤其中，所述相应的数值属于所述群集。所述方法可以包括以下中的一个或多个：-所述方法还包括提供预定的颜色，并且被分配给相应的群集的着色是预定的颜色的强度，该强度取决于相应的群集的数值；-被分配给相应的群集的着色是预定颜色的强度，该强度取决于相应的群集的中心的值而递增；-数值被相应地分类和索引，并且计算步骤内的迭代包括，在每个相应的迭代等级上，并且对于次于数值的数量的每个相应的索引，根据应用到正方矩阵的线性时间行最小值搜索算法，计算针对索引小于相应的索引的数值的子集能够获取的最小总失真，群集的数量与所述相应的迭代等级相对应；-在每个相应的迭代等级上，并且对于次于数值的数量的每个相应索引，对于每个行索引和每个列索引，矩阵项与在针对该行索引之前的索引的先前的迭代中计算的最小总失真，以及行索引和列索引之间的数值的连续的子集的失真之和相对应；-该方法还包括在每个相应的迭代等级上，存储由所述行最小值搜索算法返回的索引；-该方法还包括在计算步骤中，根据存储的索引来确定最优聚类；和/或-根据存储的索引来确定最优聚类包括迭代地对数值进行划分，从存储的索引中的最后被索引的数值开始，其中，在每个相应的迭代等级上，当前形成的群集的起始数值的索引等于在计算步骤内的迭代期间，在等于所述预定的群集数量减去所述相应的迭代等级的等级的迭代处存储的索引，所述相应的迭代等级针对等于所述当前形成的群集的最后被所索引的数值的索引的行索引。还提供了一种计算机程序，包括用于执行该方法的指令。还提供了一种计算机可读存储介质，其具有记录在其上的计算机程序。还提供了一种系统，包括耦合到存储器的处理器，所述存储器具有记录在其上的计算机程序。还提供了一种能够通过该方法获得的等值域图。还提供了一种数据存储介质，具有记录在其上的等值域图。附图说明现在作为非限制性示例并且参考附图来描述本发明的实施例，在附图中：-图1示出了方法的示例的流程图。-图2示出了系统的示例；以及-图3-4示出了方法。具体实施方式参照图1的流程图，提出了设计(例如，绘制/产生/定义)等值域图的计算机实现的方法。该方法包括提供S10(例如，空白的和/或地理的)地图(例如，至少包括表示地理区域的线的2D视图)以及多个(例如，任何数量)数值的步骤，所述多个数值表示(例如，通过以下取得的值)在该地图的相应区域上的统计变量(即，任何人口/地理/地质/地缘政治/医疗/商业/营销变量，例如，人口密度、居民汽车数量、降雨量、疾病比例)，(即，数值的每个与表示地图的局部地区/区域的地图位置相关联，由此地图可以被提供有针对线和这样的位置的参考，例如，地图是地理区域的简单绘图并且数值与像素相关联，或者地图被提供作为更加复杂的数据段，其包括描述表示地球或地球的一部分的球体在方形、矩形、圆形或椭圆形平面上的投影的数据，所述数值随后被提供有与投影上的位置相对应的坐标)。所述方法还包括针对预定的群集数量，计算S20数值的最优K均值聚类的步骤。计算步骤S20包括将应用到阶等于数值的数量的正方矩阵的线性时间行最小搜索算法迭代与预定的群集数量相对应的次数。并且方法包括将相应的着色分配(S30)给计算出的聚类的每个群集的步骤，以及在提供了相应数值的地图的所有区域上应用(S40)被分配给相应的数值所属于的群集的着色，由此将视觉细节(即，着色)添加到地图上。这样的方法改进了等值域图的设计。值得注意的是，方法允许对统计数据的值进行聚类S20，以设计等值域图，这也是由现有技术已知的并且被提供的。由于方法通过计算S20最优K均值聚类来执行这样的聚类S30，因此方法设计了相对高质量的等值域图，正如地根据理信息系统的技术所公知的，所述地理信息系统通常目标在于进行K均值聚类以产生等值域图。但最重要的是，方法通过将应用到阶等于数值的数量的正方矩阵的线性时间行最小值搜索算法迭代与预定的群集的数量相对应的次数来执行这样的计算S20。由于由方法实现的该具体算法框架，最优K均值聚类的计算被快速地执行，如下文所详述的。方法是计算机实现的。这表示方法的步骤(或者基本上所有步骤)都是由至少一个计算机或者任何相似的系统来执行的。由此，方法的步骤是由计算机执行的，可能是全自动地(例如，除了提供S10之外的所有步骤)，或者，半自动地。在示例中，可以通过用户-计算机交互来执行对方法的至少一些步骤的触发(例如，提供S10)。所需的用户-计算机交互的水平可以取决于预见的自动化水平，并且与实现用户的愿望的需求相平衡。在示例中，该水平可以是用户定义的和/或预定义的。方法的计算机实现的典型示例是利用适合于该目的的系统来执行方法。系统可以包括耦合到存储器的处理器，存储器具有记录在其上的包括用于执行方法的指令的计算机程序。存储器还可以存储适于维护由方法处理的数据的数据库。存储器是适于这样的存储的任何硬件，可能包括多个物理不同的部分(例如，一个用于程序，并且可能一个用于数据库)。图2示出了系统的示例，其中，系统是客户端计算机系统，例如，用户的工作站。示例的客户端计算机包括连接到内部通信总线1000的中央处理单元(CPU)1010、同样连接到该总线的随机存取存储器(RAM)1070。客户端计算机进一步被提供有与连接到总线的视频随机存取存储器1100相关联的图形处理单元(GPU)1110。视频RAM1100在本领域中还被称为帧缓冲器。大容量存储设备控制器1020管理对大容量存储器设备(例如，硬盘驱动器1030)的访问。适合于有形地实现计算机程序指令和数据的大容量存储器设备包括所有形式的非易失性存储器，举例而言，包括诸如EPROM、EEPROM、以及闪速存储器等的半导体存储器设备、诸如内部硬盘和可移动盘等的磁盘；磁光盘；以及CD-ROM盘1040。前述中的任何一项可以由专门设计的ASIC(专用集成电路)来补充或者被集成到专门设计的ASIC当中。网络适配器1050管理对网络1060的访问。客户端计算机还可以包括触觉设备1090，例如，光标控制设备、键盘或其类似物。在客户端计算机中，使用光标控制设备来允许用户选择性地将光标定位在显示器1080上任何期望的位置处。此外，光标控制设备允许用户选择各种命令，以及输入控制信号。光标控制设备包括用于向系统输入控制信号的多个信号生成设备。通常，光标控制设备可以是鼠标，鼠标的按键用于生成信号。可替换地或除此之外，客户端计算机系统可以包括触摸板，和/或触敏屏幕。计算机程序可以包括可由计算机执行的指令，所述指令包括用于使上述系统执行方法的模块。程序在包括系统的存储器的任何数据存储介质上是可记录。程序可以例如在数字电子电路中、或者在计算机硬件、固件、软件或其组合中实现。程序可以被实现为装置，例如，有形地实现在用于由可编程处理器执行的机器可读存储设备中的产品。方法步骤可以由可编程处理器执行，所述可编程处理器通过对输入数据进行操作并且生成输出来执行指令程序，以执行方法的功能。处理器可以由此是可编程的并且被耦合以从数据存储系统、至少一个输入设备、以及至少一个输出设备接收数据和指令，以及将数据和指令发送到所述数据存储系统、至少一个输入设备以及至少一个输出设备。应用程序可以以高级过程或面向对象的编程语言来实现，或者在需要时以汇编或机器语言来实现。在任何情况下，语言都可以是编译的或解释的语言。程序可以是完全安装程序或更新程序。在任何情况下，该程序在系统上的应用产生用于执行方法的指令。方法提出了一种改进的聚类算法，其具体应用到等值域图设计。但是在详述方法的算法解决方案之前，现在详述其背景(即应用)。方法用于设计等值域图。如上文所述以及根据现有技术公知的，这种类型的地图是一种视图(通常为2D，但也可能是3D)，其根据地理区域的子区域处的统计变量/测量的值，利用不同的着色(例如，颜色/纹理)来表示所述地理区域的子区域，每个着色与统计的相应的值或连续范围相关联。在此上下文中，术语着色(例如，颜色/纹理)一般指定分布在所述子区域上的任何视觉信息。通常，着色是其强度为统计变量的增函数的颜色(例如，RGB或灰度)。但是，在棱柱地图(等值域图的特定情况)的情况中，着色也可以是与统计变量成比例的高度。图3示出了可以通过方法获得的等值域图的示例，不同的灰度强度被分配给地图的不同区域。一般而言，等值域图的设计可以基于预定的有序的并且一维的着色值的集合与统计变量的值之间预定的递增关系(由此与一维有序域相关联(如果不是被包含在其中))。在该方法的情况中，这样的地图的设计遵从现有技术的传统流水线。首先，方法在S10处开始于简单的等值线图(通常为空白的)以及与地图的相应区域相关联的统计变量的值。目标是以某种方式将统计的这种数值表示进行变换，以获得简单并且有意义的视觉化。这就是聚类所介入的地方。实际上，等值域图可以在理论上将不同的着色分配给统计变量的每个值，但是这对计算终端而言将不是繁重的，并且这对于用户而言也不是很有意义(地图将太过“繁杂”)。如已知的，从语义视点，K均值聚类是特别有利的，因为其以有意义的方式来聚集统计变量的值。由此，方法利用K均值聚类来对统计变量的值进行聚类(群集的数量取决于地图的查看者和/或设计者，并且，在该方法的上下文中，群集的数量是预定的)，并且随后，与统计变量的值相关联的地图的每个区域(可能是所有区域，取决于在S10中提供的数据)现在与相应的群集(针对该区域的统计变量的值的群集)相关联。现在，由于方法将相应的着色分配S30给计算出的聚类的每个群集，因此方法可以通过S40来产生等值域图。这都是非常传统的做法，并且在上文引用的现有技术中被广泛解释。方法可以实现用于执行着色的任何的传统方式。一种传统的方式是提供预定的颜色，并且随后将该预定颜色的强度分配给相应的群集。例如，如果颜色是以RGB来提供的，则每个相应的强度与应用到RGB权重的因子(例如，在0和1之间)相关联。可以对灰度颜色执行相同的操作。预定颜色的强度可以取决于相应的群集的中心的值而递增。换言之，群集的中心(被分组到该群集中的统计变量的平均值)递增地定义所应用的强度(即，对于可视化区域，强度越大意味着统计变量的值越高)。可以应用很多其他方法，例如，在下面的URL中所描述的那些(在优先权日期是可访问的)：http://en.wikipedia.org/wiki/Choropleth_map#Color_progression.在示例中，方法可以简单地应用0(黑色)与255(白色)之间的不同灰度等级(如图3上所示)。例如，如果聚类的结果提供了群集中心c1,...,cK(c1<c2<..<cK)，则方法可以应用将c1映射到0(白色)并且将cK映射到255(黑色)的仿射变换：所有这些都是非常传统的，并且不需要进一步解释。现在，由于方法找到了具有降低的算法复杂度(相对于现有技术的精确解决方案算法，并且甚至相对于合理试探法)的最优K均值聚类，因此其允许使用该最优K均值聚类来设计等值域图。这与现有技术不同，现有技术或者应用不产生最优K均值聚类的试探方法，或者在实践中是不可应用的(因为过高的复杂度)。值得注意的是，方法发现了对n大于1000(例如，10000左右)的值的最优K均值聚类，尤其是对于K高于50，甚至高于200，例如，256(用于灰度级值的8位编码)的值而言。例如，所述方法起作用。方法产生地理统计数据的最优量化，以使得颜色映射层是最优的(相对于作为标准目标函数的“总失真”而言)，比当前的其他最优量化器更快。实际上，该方法的计算时间与最好的试探方法(其产生非最优的量化)相当(并且通常甚至更快)。就此，已经提供了方法的背景以及方法的许多示例，但是尚未提供有关该方法的核心的细节(即，定义后面的设计S30-S40的计算S20步骤)。这在下文中完成，应当注意，下文提供的所有实现示例可以在上文提供的示例应用的任何一个中应用。如上文指示的，方法包括针对群集的预定数量K，计算S20数值的最优K均值聚类的步骤。但这不是暴力地完成的。实际上，计算步骤S20包括将线性时间行最小值搜索算法迭代与K相对应的次数(在下文讨论的示例中为K-1次)。由于其使用的是公知的行最小值搜索算法范畴的任何预定的线性时间算法，因此计算步骤S20具有低的复杂度。由此，方法实现了K均值聚类问题的新的并且具有算法效率的解决方案。在提供有关计算步骤S20的更多细节之前，现在讨论量化。实际上，在示例中，统计数值可以在S30-S40之前以某种方式被“量化”。实际上，每个区域的统计值可以被视为由通过方法与该区域相关联的群集的中心所替换，因为它们的相对应的区域被应用了该群集的着色。如已知的，标量量化是使用有限集来近似真实值的计算工具，其中，V的元素(称为数字阶梯)是用作近似的值。标量量化器在数学上被定义为映射：这使得x和q(x)之间的距离小，该距离是任何预定距离(其中，距离的概念可以取决于上下文)，例如，欧几里得距离。在实现中，总是通过将划分成区间I1＝[-∞，a1]，I2＝[a1，a2]，…，IK＝[aK-1，∞](a1＜…＜aK-1)来定义标量量化器，并且随后，针对每一个将Ii指示为唯一区间，以使得x∈Ii，我们关联了q(x)＝ci。实数a1，…，aK-1称为“决策边界”。Gray和Neuhoff的文章“Quantization”提供了对量化的全面调查。方法集中于K均值设置，这是公知的，并且在例如MacQueen的文章“SomeMethodsforclassificationandAnalysisofMultivariateObservations”中定义，K均值设置是使用最为广泛的。在这样的设置中，方法在S20中，对于按照升序分类的给定元组并且对于给定的整数K(通常，K＝2b，其中b是可用于对每个xi进行编码的位数)找到实现最小总失真的采用K个数字阶梯的量化器q，总失真被定义为：TD(q)=Σ1≤i≤n(xi-q(xi))2]]>显然，要使该量最小化，方法可以仅处理将每个真实值x映射到其最接近的数字阶梯的量化器。因此，问题正好等价于找到使以下最小化的K个中心c1，…，cK：TD({c1,...,cK