时空演变格局(土地利用遥感龙口市)
20世纪以来，世界格局发生了哪3次变化
1。20世纪初期，国际中心位于欧洲，第1次世界大战之后，欧洲的地位有所下降，美国一步步上升。 2。两级格局：雅尔塔会议分治世界。二战之后，美国大发战争财，急剧膨胀，欧洲国家基本成为战争废墟，急需恢复，苏联经过几个五年计划，实力大大增强。 3。一超多强，多极化趋势明显。东欧剧变，苏联解体，两极中的一极倒塌，两极格局不复存在。美国为超级大国，俄罗斯，欧洲诸国，中国，日本等国为多个强国。
1991年到到现在世界格局的转换时间
四个阶段。1945~1946年雅尔塔体系形成。1946~1955年以美苏对抗为核心，社会主义阵营与资本主义阵营的两极对峙格局形成。20世纪50年代中期90年代初美国、苏联争夺世界霸权。1991年到现在“一超多强”向多极化格局演变。中国地位，第1阶段，中国是世界反法西斯同盟的重要成员，为世界反法西斯战争胜利做出重要贡献，国际地位大为提高，成为联合国安理会的常任理事国。第2阶段，中国实行“一边倒”条文，是社会主义阵营的重要成员。为反对美帝国主义侵略，维护世界和平做出重要贡献。第3阶段，20世纪50年代末，中苏联系破裂，社会主义阵营不复存在。中国团结亚非拉广大新兴国家，推动世界反帝反殖反霸斗争。70年代初中美两国开始关系正常化，以共同对付苏联威胁。第4阶段，中国是世界向多极化格局演变的重要力量，在反对霸权主义，维护和平、促进发展方面起到重要作用。咱们国家的发展强大将推动世界向多极化发展。
时空分异和时空演变一样吗
不一样。时空分异就是时间和空间上的区别规律。时空演变从古到今， 关于时间与空间的探究从未间歇过。哲学和物理学都比较注重研究时间和空间，但两者在这样一个话题上显现了分歧。
未来世界格局的演变?
中国会逐渐崛起，欧洲会更加统一，战争还会继续，美俄中欧四国鼎立。 到20二十四中国会有本人的一支航母编队（可作战的，不是花架子了），甚至两支。紧接着中美间的利益冲突会日益增大，仅需要一个导火索，战争不可避开，也许钓鱼岛，也许老调重提的，甚至或许是南海的岛礁归属，中美间较量一旦从经济升格为战争，就会牵动整个世界，俄，欧盟，日均会参加其中，当然免不了朝韩，东南亚等国这几个跳梁小丑的表演。
湖泊水体时空演化规律和水质有关系吗
湖泊水体时空演化规律和水质有关系。例如利用ERDAS IMAGINE 9、2处理1957。1974。1989。1990、2001。2006。2011和2012年8个年代的航空和卫星影像，ArcMap 10处理剑湖演变矢量数据和栅格数据，Fragstats 3、3计算景观指数值，Excel统计景观指数值，采取图、文、表相结合的方式方法，剖析了1957~2012年剑湖水域面积和形状的时空演变规律，同时从驱 动因素、驱动方向、驱动速度和驱动过程4个方面剖析了剑湖水域演变的驱动机制。结果表明:(1)面积和形状时空演变规律。在1957~2012年的55年 间，剑湖水域面积呈随着时间变化而逐年减少的趋势，其形状呈随着时间变化而日渐复杂的趋势。(2)驱动机制。驱动因素对剑湖萎缩的效果按其作用大小从大到 小的排序，在1983年以前是疏浚河道→泥沙淤积→湖滩开发→湖水污染，在1983年之后是泥沙淤积→湖滩开发→湖水污染;剑湖入湖河流入湖口和出湖河流 出湖口周围是剑湖萎缩的主要驱动方向，而带着大量泥沙入湖的金龙河入湖口周围是剑湖萎缩的最主要驱动方向。
基于RS和GIS技术的龙口市土地利用时空变化监测与剖析
徐秋晓1 于明洋2
（1、山东省地矿工程勘察院，济南250014；2、山东建筑大学土木工程学院，济南250101）
作者简单介绍：徐秋晓（1979—），女，助理设计师，主要从事遥感、水文地质、环境地质勘查工作。
摘要：遥感技术与地理信息系统技术的发展，为研究全球变化和可持续发展提供了信息源和技术手段。而土地利用作为地球表层系统最突出的景观标志，其变化是最近几年以来全球变化钻石的重要范畴。本文以山东省龙口市作为研究区，应用基于遥感影像综合理解模型的龙口市土地利用/土地覆盖分类方法，建立地学规律知识库，提取不同时期的土地利用类型。最后利用地学信息图谱监测与剖析土地利用的时空变化。
关键词：土地利用；遥感影像；地学辅助信息；信息图谱
遥感技术与地理信息系统技术的发展，为研究全球变化和可持续发展提供了信息源和技术手段。而土地利用作为地球表层系统最突出的景观标志，其变化是最近几年以来全球变化钻石的重要范畴。只有对土地利用时空变化进行监测与剖析，更佳地了解土地利用变化的过程和机制，并且通过调整人类社会经济活动，促使土地利用更趋合理，保证国家宏观战略决策的针对性、有效性，才能达到土地资源可持续利用的意图。龙口市作为我国沿海对外开放较早的城市，其土地利用变化具有表现性。
1 研究区概况
龙口市位于胶东半岛西北部，东与蓬莱市接壤，南与栖霞、招远市毗连，西、北濒临渤海；全境东西最大横距46、08km，南北最大纵距37、43km；土地总面积（含桑岛、依岛）893、32km2。
本次研究区范围以1∶10000 地形图矢量化生成的研究区边界对遥感影像进行裁剪所得，作为遥感数据信息，面积共89217、45hm2（不含桑岛、依岛）。
2 基于遥感影像综合理解模型的龙口市土地利用/土地覆盖分类
影像理解是研究通过计算机系统来解释图像，从而实现类似人类视觉系统理解外部世界的一门学科。在影像理解系统中，存在两项根本的任务：从输入图像中提取出与模型相适应的图像结构或线索，而后完成输入图像中图像结构与模型中目标的正确影射（周成虎，1999）。影像理解不同于模式识别，模式识别通常来讲按预先规定的测量集对对象作简单容易的分类，而影像理解则要对影像作出描述和解释，需要涉及不同处理层次实体间的互相作用（王润生，1994）。
针对土地利用/土地覆盖的分类特点，本文构筑了基于GIS信息的遥感影像综合理解模型。模型分为两个过程，即遥感影像理解过程和地理信息系统处理过程，两个过程的具体内容为：
（一）遥感影像理解过程，主要完成遥感影像的前期理解过程，包括以下几个过程：
1）遥感影像预处理：包括图像格式转换，图像纠正和图像增强变换等。图像格式转换完成遥感影像到REDAS软件系统处理格式的转换（﹡。img）；图像纠正完成遥感影像的大气校正、几何纠正、辐射增强以及遥感影像的匹配、镶嵌等；遥感影像信息增强和变换处理可以突出有关的专题信息，本次使用方法主要有线性拉伸、K-T空间变换、边界增强等。
2）与地学辅助信息的配准：主要完成遥感影像与地学辅助信息（各类专题GIS数据着重是坡度和高度）的坐标、投影系统转换，使得遥感影像与所采用的辅助地学信息纳入到统一坐标与投影系统下。
3）“训练区”选择与计算：经过对遥感影像信息特征的初步理解，同时结合地学辅助信息以及实地考察，确定样本“训练区”。“训练区”要有典型性与可分性，确定“训练区”后，对“训练区”数据进行计算，确定样本的统计信息（均值、最大最小值、方差矩阵、协方差矩阵等）。
本次提取的土地利用类型为六大种类型，即建筑用地、耕地、水域、园地、林地和未利用地，采用AOI扩展方式进行训练区选择，运用此种方法进行训练区选择时，初始种子（seed pixel）和光谱距离的阈值十分重要，任一训练区初始种子和光谱距离的阈值经数次实验后才能确定。不同地类这两个参数是不同的，一般在选取了2～3个训练区后，观看报警掩膜的情形，即时对这两个参数进行修改和调整。增添训练区时，及时将报警掩膜和以前的报警掩膜叠加，判断样本数据的质量变化情况，并做出调整，直至报警掩膜和实际地类比较符合时，即认为这一类的训练区选择完毕。一种地类训练完毕，需要对最后的训练区作一个整体的剖析，保证光谱的纯度。
（二）地理信息系统处理过程，在GIS系统支持下，完成地学辅助信息的处理，主要包括以下过程：
1）专题信息导入与预处理：通过地面调查或专业人士知识经验，收集地学辅助信息（包括各类专题信息、统计资料等），导入到GIS系统中，完成辅助信息的前期预处理，包括各类矢量数据的数字化、编辑、拓扑关系的建立等，统计数据的整理与地学编码、统计数据的空间化等工作。
2）辅助数据的生成：利用前期处理好的辅助数据，进行各类数据的格式转换，如矢量数据的栅格化，点状统计数据插值（IDW/Kriging等方法）生成面状数据；最后统一坐标与投影系统，达到与遥感数据的配准。
3）建立各类地学辅助因子数据库：GIS软件支持下，基于前面生成的各类辅助数据，建立专题数据库，形成地学辅助因子数据库。
（三）知识库的生成过程，建立专业人士知识库，着重是地学规律知识库，包括以下过程：
1）知识获取：经过对照遥感影像，进行野外考察，针对各类有表现性的影像特征对该地区所有地貌条件下的土地利用/土地覆盖状况、植被分布、生态环境条件进行实地考察，获取各类实践知识。
2）知识库生成：对获取的知识进行整理，从总体上归纳出各类规律，形成知识规则，最后经过对“训练区”的数据不断训练，修改和调试知识规则，形成地学规律知识库。
土地覆盖/土地利用类型对高程有明显的依靠关系。对研究区已有的土地利用现状图剖析，建筑用地、园地主要分布于海拔300m以下；耕地绝多数分布于海拔150m以下；林地的分布范围较广，不同海拔皆有分布，分布于10m以下的林地着重是沿海防护林和公路绿化林地。350m以上的海拔，只有林地分布，没有其他地类。
坡度数据可以 使用于区别某些土地覆盖/土地利用类型。依据实地考察结果和对地形图、已有土地利用现状图的剖析，建筑用地、水浇地、园地主要分布于坡度小于20 °区域，坡度大于10 °时很少有水浇地。因 此，遥感数据的光谱特征同样表现为绿色植被，难以判断是水浇地还是林地时，坡度数据是一个有价值的参数。下面是几种地类与高程和坡度的详细关系。
在本次的分类知识库，专业人士规则采用下面的基本形式来表达：
IF（条件）THEN（结论），Confidence（结论信任度）
其中信任度的值域为［0，1］，值为0 时，完全排除当前像元为所给类别的可能性；当取值为1时，维持像元原始信任度，且预示当前结论永久成立。信任度可以依据地学经验或专业人士打分等方法来核实确定。
知识的预示与知识库的构造要结合地学问题的研究特点。通过不断修改和调试知识库，使影像解译结果基本达到人工目视解译的作用。当信任度的值为0时，则排除了当前像元为（结论）所给出的类别的可能性；而当取值为1 时，已有信任度值不该变。这种预示方法，不仅考虑到了遥感影像解译的特征，而且明显地减少了知识库中规则的数量。这对于大数据量的遥感数据处理是极为重要的。下面给出规则库中的规则：
IF VALUE ＝1 DEM ＜300 AND SLOPE ＜20 THEN 建筑用地 CF ＝1
IF VALUE ＝2 DEM ＜50 AND SLOPE ＜10 THEN 水浇地 CF ＝1
I VALUE ＝3 THEN 水域 CF ＝1
IF VALUE ＝4 DEM ＜300 AND SLOPE ＜20 THEN 旱地 CF ＝1
IF VALUE ＝5 DEM ＜300 AND SLOPE ＜20 THEN 园地 CF ＝1
IF VALUE ＝6 THEN 林地 CF ＝1
ELSE IF DEM＞＝300 OR SLOPE＞＝20 THEN 林地 CF ＝1
ELSE IF 50＝＜DEM＜＝250 AND 10＝＜SLOPE＜20 THEN 园地 CF＝0。5
ELSE 林地 CF ＝1
本次分类结果精度评价采用分层随机采样法，主要参考龙口市土地利用现状图，同时结合目视判读的结果以及现场验证，对两个时期的分类结果进行精度评价。各时期的遥感分类结果直接参照同年的土地利用现状图进行精度验证。因为遥感分类体系与土地利用现状的分类体系有一定差别，于是，在进行随机采样之前，将遥感分类图像的土地利用类型和土地利用现状的分类进行适当统一。紧接着分别对每期遥感分类结果选取300个样本点，并保证每类有10个以上的样点，用基于误差矩阵的精度评价方法，对龙口市1989年和2003年的分类结果进行评价。实用Kappa系数计算表明1989年和2003年龙口市土地利用TM遥感分类结果的总体精度和使用者精度皆在75%以上，Kappa系数也皆在0。8以上，达到最低允许判别精度0。7的要求。这几个表明了龙口市两期图像的土地利用遥感分类结果均比较理想，各地类的分类精度也较高。
3 龙口市土地利用时空变化剖析
随着地球信息科学的兴起与发展，人们可获取的资源的极大丰富以及信息处理技术的极大提高，特别是动态可视化技术获得新的冲破。在此需求与技术背景下，陈述彭先生倡导在传统地学图谱的基础上开展地学信息图谱的探讨与研究。地学信息图谱是地学图谱在地球信息科学基础上的自然延伸，是依照一定指标递变规律或分类体系排列的一组能够反映地学空间信息规律的数字地图、图表、曲线或图像。地学信息图谱是“图”与“谱”的结合，兼有图形与谱系的双重特征。
本文所采用土地利用图谱剖析模型包括两部分内容：①转移矩阵，从中可以看出各个时序单元土地利用变化的主要类型以及各地类的补给来源。②不同时序单元内土地利用图谱剖析，考察图谱单元的空间组合与时空位移。
3、1 转移矩阵
转移矩阵对于剖析土地利用类型之间的流向具有重要作用，它不但可以定量说明土地利用类型之间的互相转化状况，而且可以揭示不同景观类型间的转移概率，从而可以更佳地了解土地利用的时空演化过程。转移矩阵包括转移面积矩阵、概率矩阵。
表1 1989～2003年土地利用转移矩阵（单位：hm2）
注：R为各土地利用类型的转移比例（%）。
从表1可以看出，在1989年至2003年，建筑用地发生用途流转面积为487、61hm2，占初始建筑用地面积的比率为3、62%，不存在明显流向。耕地发生用途流转面积为17632、40hm2，占初始耕地面积的比率为 49、85%，主要流向是园地，该流向13835、82hm2，占初始耕地面积的比率为39、11%，其次是建筑用地，该流向3508、77hm2，占初始耕地面积的比率为9、92%。水域发生用途流转面积为492、92hm2，占初始水域面积的比率为8、72%，主要流向是园地和林地，共占初始水域面积的比率为6、4%。园地发生用途流转面积为1316、62hm2，占初始园地面积的比率为7、89%，主要流向是建筑用地，该流向 759、25hm2，占初始园地面积的比率为 4、55%。林地发生用途流转面积为1405、45hm2，占初始林地面积的比率为 10、45%，主要流向是建筑用地，该流向6二十四。64hm2，占初始林地面积的比率为 4、66%。未利用地发生用途流转面积为2158、74hm2，占初始未利用地面积的比率为 47、二十四%，主要流向是园地，该流向1306、85hm2，占初始未利用地面积的比率为28、60%，其次是流向林地和建筑用地，共占初始未利用地面积的比率为12、89%。
3、2 土地利用信息图谱
在图谱中，共有36类图谱单元，即土地利用变化类型，其中有30类显示为土地利用类型发生了变化，占研究区总面积的26、34%。为了更简单并且明了地读取土地利用类型流转的主要方向，认识土地利用变化的主要特点，将该30类变化的图谱单元依照面积大小进行排序，计算各类图谱单元的转换面积百分率和累计转换百分率，统计其中涵盖变化总面积的92、04%的10类图谱单元，得到1989～2003年土地利用主要图谱单元类型的面积排序表（表2）。
表2 1989～2003年土地利用主要图谱单元类型的面积排序表
从表2可以看出，龙口市在1989～2003年间土地利用变化最显著的结构特征是耕地向园地的转化，该流向的耕地面积共13835、82hm2，占到整个变化面积的58、89%。其次是耕地向建筑用地的转化，面积为3508、77hm2，占到整个变化面积的14、93%。再次是未利用地向园地的转化，面积为1306、85 hm2，占到整个变化面积的5、56%。可以看出，研究时段内龙口市的土地利用类型着重是流向建筑用地和园地。
4 结论
（一）本研究利用陆地卫星资料ETM＋对龙口市土地利用进行了时空监测与剖析，取得了较好的作用；不过ETM＋的影像分辨率较低，它着重是反映一些综合的地类信息，对于一些土地利用图斑较为破碎，用地类型交错复杂的地区，其地物提取就有一定的难度，所以在土地利用变化的详细监测上就有所欠缺。在今后的研究中，应结合高分辨率卫星影像进行研究，积极探求新的动态监测方法以充分利用这几个高分辨率的遥感数据来获得更可靠准确的区域土地利用变化信息。
（二）怎样充分地利用地理信息系统空间数据库提供的富饶的地理辅助数据，进而自动发现知识，并融合多尺度、多时相的高分辨率的遥感数据，建立灵活高效的推理机制，从而完成遥感影像的专题信息自动提取，是需要进一步的研究方向。
（三）“地学信息图谱”是一新兴的学术思想，目前还处于认识阶段，对它的理解尚不很成熟，需要更加的多的学者和更加的多的研究工作来完善对其的认识。本文利用遥感技术开展土地利用演化与发展的信息图谱研究，来反演时空变化，进而认识客观世界，揭示和再现过去，是一先进可行的技术途径。
（四）研究结果表明龙口市在1998～2003年间，土地利用方式互相转化较为频繁，变化强度大，矛盾非常突出。主要表现为耕地的不断减少，园地和建筑用地持续增长；其主要流向是从耕地流向园地和建筑用地。这几个变化造成了龙口市耕地质量下降和利用程度加强，给耕地保护带过来了巨大压力。于是，土地管理部门应该从宏观决策上给予注重和重视，处理好社会经济发展与土地后备资源储备的关系实施可持续发展。
参考书籍
陈述彭。2001、地学信息图谱探索研究。商务印书馆
李小英，彭望碌，曹彤。2002、土地利用演化信息图谱的研究——的北边京市顺义县为例。地球信息科学，2：（58～63）
齐清文，池天河。2001、地学信息图谱的论理和方法。地理学报，56（增刊）：8～18
术洪磊，毛赞猷。1997、GIS辅助下的基于知识的遥感影像分类方法研究——以土地覆盖土地利用类型为例。测绘学报，11（四）
叶庆华。2001、黄河三角洲土地利用/土地覆盖变化的时空复合模式研究。中国科学院地理科学与资源研究所博士学位论文
周成虎，骆剑承，杨存建，刘庆生。2001、遥感影像地学理解与剖析。科学出版社，180～二十四0
Lucas I F J，Frans J M，Wel V D。1994、Accuracy assessment of satellite derived land-cover data：a review。Photogrammetric Engineering and Remote Sensing，60（四）：410～432