云南省测绘学会学术年会论文集

李喆1 马润赓2

（1.北京大学遥感与地理信息系统研究所，北京，100871；2.中国地质大学地球科学与资源学院，北京，100083）

摘要：利用1999年2月和2001年4月2幅TM影像资料和其他辅助资料为资料源，对主成分分析处理后的遥感影像进行计算机自动识别，采用人机交互式解译方法获得各种土地利用类型分布信息，并使用地理信息系统的空间分析和数理统计功能分析丽江坝子地区土地利用变化的数量变化，土地利用动态度等特征，以及各类型之间的转化情况。结果表明：本区域耕地面积增加了，城乡用地面积增加了，水体面积减少了，林地和草地基本保持不变；各类型的转化突出表现为水域向城乡用地，未利用土地向城乡用地和耕地转换。

关键词：土地利用变化；马尔柯夫矩阵；云南丽江地区；遥感和 GIS

区域土地利用/覆被变化（LUCC）是全球环境变化研究中的重要领域之一。土地利用变化数据可以提供研究区域土地利用状况的现实状况，保证国家及时、准确地掌握土地利用变化情况，为制定国民经济发展规划、计划和宏观决策提供科学依据［1］。土地利用变化包括时间变化、空间变化和质量变化等三个方面，其中空间变化反映土地利用变化的空间类型、变化类型的空间分布以及区域差异，是土地管理和规划关注的焦点问题［2］。这一问题的解决关键在于如何提取土地利用变化数据以及如何对获取的信息进行科学的分析处理。遥感技术以其快速、准确、周期短等优点在大中尺度的土地利用/覆盖变化的监测中具有明显的优势。本文利用遥感和 GIS 相结合的方法，通过数理统计分析，描述了丽江坝子地区1999～2001年间土地利用的数量变化和空间变化特征，为该区域土地管理决策、生态环境保护、资源合理开发等奠定了基础，对于区域土地可持续发展和利用具有重要意义。

1 研究区概况

丽江位于云南省西北部［3］，金沙江中游，地理坐标为：99°23′E～101°31′E，25°59′N～27°56′N。四周分别与迪庆州、怒江州、大理州、楚雄州和四川省攀枝花市、凉山州毗邻，辖丽江纳西族自治县、宁蒗彝族自治县、永胜县和华坪县，国土总面积有。由于地处青藏高原和云贵高原的结合部，属横断山区，地形复杂，山高谷深，长期自然环境相对封闭，形成具有典型特点的垂直自然带和立体气候，动植物区系极为复杂，珍稀野生动植物丰富。全区总人口109 万，包括汉、纳西、彝、傈僳、藏、白、普米等23个民族，其中少数民族人口占57%［3］。

2 研究方法

研究资料

分别拍摄于1999年2月和2001年4月2幅 LANDSAT－TM 图像，地面分辨率为30m，在研究区域内没有云层覆盖，质量较好；另有1∶1 万地形图和1∶5 万土地利用图，用于辅助训练样区的选择和目视解译；还有其他统计资料（气象、水文、人口、土壤、社会经济）等。

主要工作流程

使用遥感图像信息进行变化检测的方法主要有三种：分类后比较法，多时相图像直接求交法和多时相图像分类法［4］。多时相图像直接求交法要求两个时相比较接近，检测变化信息一般使用图像差值、比值或主成分分析，操作简单，但是得不到变化的具体地物类型。本文所使用遥感图像时间差异接近3个月，不宜使用直接求交法。多时相图像分类法在使用时要使用静态类型和动态类型，动态类型的训练样本一般不易确定。本文主要使用分类后比较法。

本研究的主要工作流程如图1 所示：首先使用 PCI 软件分别对2 幅 TM 图像进行图像预处理（如几何纠正等），并分别对以上两个时相的影像做主成分分析，然后充分结合各种知识（如地形图，植被图等）进行监督分类［3］。监督分类时采用《中国土地利用现状调查技术规程》规定，将土地利用类型划分为：耕地，林地，草地，水域，城乡工矿居民地和未利用土地共6 类。采用人机交互式方法判读分类结果，并将解译结果以栅格形式转入地理信息系统软件 ARC GIS 中进行数据处理，对二期数据进行空间叠置分析，获得土地利用变化数据，以此为基础进行土地利用动态变化分析。

图1 项目主要工作流程图

几个关键技术的处理

遥感图像的几何配准

两个时相遥感图像的几何配准是动态变化研究的基础。为了便于利用地学信息进行辅助分析，可以将遥感图像配准到大地坐标系中。本文选择1∶1 万地形图为基准，在遥感图像上均匀选择12个地面控制点，使用二次多项式拟合，灰度采样方式为双三次卷积，进行几何配准和灰度重采样处理。结果总体均方误差为，配准误差小于1个像元，满足土地利用动态监测的要求。

遥感图像的主成分分析

主成分分析（又称主分量变换）是一种尽可能不丢失信息用几个综合性指标汇集多个变量的测量值而进行描述的方法［5］。在多光谱图像中，由于各波段的数据间存在相关的情况很多，通过主成分分析就可以把图像的大部分信息用少数波段表示出来，使得信息几乎不丢失但数据量可以减少。本文分别对于二期 LANDSAT－TM 图像的除热红外波段的6个波段采用主成分分析，把获得的第一到第三主成分进行彩色合成，得到处理后的结果图像。

结果数据的生成

根据各种辅助资料，在经过主成分分析处理后的二期遥感图像上选择6个土地利用类型的训练样本，采用最大释然法进行计算机自动识别。由于遥感影像上存在同物异谱、异物同谱以及混合像元等情况，计算机自动识别在某些类别的区分上不太理想，例如未利用土地和城镇用地极易发生混淆。实际工作还需要辅助各种实测资料进行目视修正，以得到两个时相的最终分类结果，然后将分类结果以栅格形式转入地理信息系统软件ARC GIS中，对二期数据进行空间叠置分析，获得土地利用变化数据，并进行必要的数理统计，以便结果分析使用。

3 结果分析

土地利用变化幅度分析

区域土地利用变化幅度主要体现在不同土地利用类型的面积总量变化上，可以提供区域土地利用变化总态势和土地利用结构变化信息。

利用二期遥感资料，对丽江坝子地区两期土地利用数据进行统计分析，得到结果见表1。

表1 1999～2001年土地利用变化表

表1说明，3年来本区域耕地面积增加了，城乡用地面积增加了，水体面积减少了，林地和草地则基本保持不变。

土地利用动态度分析

土地利用动态度可以简单采用单一土地利用动态度和综合土地利用动态度来描述［6］。其中单一土地利用动态度可定量描述区域一定时间范围内某种土地利用类型变化的速度，可提供土地利用变化的区域差异和预测未来土地利用变化趋势；综合土地利用动态度用于表征区域土地利用变化的速度。

单一土地利用动态度公式表达式为：

土地信息技术的创新与土地科学技术发展：2006年中国土地学会学术年会论文集

式中，Ua，Ub分别表示研究初期及研究末期某一种土地利用类型的数量；T表示研究期时段长。

综合土地利用动态度表达式为：

土地信息技术的创新与土地科学技术发展：2006年中国土地学会学术年会论文集

式中，LUi 表示研究初期第 i 地物类研究末期转为其他土地利用的面积；ΔLUi 表示第i地物类研究末期转为其他土地利用类型的面积；T表示研究时段。

根据公式（1）和（2）计算丽江坝子地区土地利用六种类型的年变化率。结果说明，1999～2001年丽江坝子地区土地利用的年变化率为，其中城乡用地和未利用地变化速度最大，年变化率分别达到和，林地和草地基本保持不变，耕地和水域的年变化率分别为和。

土地利用区域差异分析

由于地形地貌和气候等自然条件差异，经济发展和人口增长速度不同，使得土地利用区域差异显著。对于土地利用变化的区域差异，可用某一特定土地利用类型相对变化率来表示。相对变化率是一种很好地反映土地利用变化区域差异的方法，其表达式为：

土地信息技术的创新与土地科学技术发展：2006年中国土地学会学术年会论文集

式中，Kb、Ka分别代表某区域某特定土地利用类型研究初期和研究末期的面积；Cb、Ca分别代表全区域某特定土地利用类型研究初期和研究末期的面积。

根据行政区划，将丽江划分成丽江县、宁蒗县、永胜县和华坪县4个区域，根据遥感数据分别计算6 种土地利用类型的相对变化率，其结果见表2。

表2 1999～2001年土地利用区域差异表

从表2中可以看出，土地利用变化存在明显差异：①就耕地而言，丽江县最大，达到，明显大于永胜县，宁蒗县和华坪县也大于永胜县；②林地，草地，水域变化各区域差异不大；③城乡用地的区域差异最为明显，丽江县大于其他3县很多，达到；④未利用地的变化在华坪县较突出，约为，几乎是永胜县的3倍。

土地利用类型转换矩阵分析

土地利用类型之间的相互转化情况，可采用马尔柯夫转移矩阵模型来进一步描述［7］。马尔柯夫链是一种具有“无后效性”的特殊随机过程，它反映的是一系列特定时间间隔下，一个亚稳定系统由n时刻向n＋1 时刻状态转换的一系列过程中，n＋1 时刻的状态只与n时刻的状态有关。由于土地利用类型演变具有马尔柯夫随机过程的性质：①一定区域内，不同土地利用类型之间具有相互可转化性；②土地利用类型之间的相互转化过程包含着较多尚难用函数关系准确描述的事件，故可采用马尔柯夫转移矩阵模型来描述土地利用类型动态转换。

马尔柯夫模型在土地利用类型转化上应用的关键是确定土地利用类型之间相互转化的转移概率矩阵P。若采用斑块相互之间面积的转移概率为矩阵元素，则转移矩阵模型为：

土地信息技术的创新与土地科学技术发展：2006年中国土地学会学术年会论文集

式中，Pij为土地利用类型i转化为土地利用类型j的转移概率。

遥感图像计算机自动识别中得到的土地利用变化转移概率矩阵如表3所示。

表3 1999～2001年土地利用变化转移概率矩阵单位：%

从表3可以看出：

（1）耕地的增加量主要是来自草地和未利用地，分别占和，减少量大部分转化为未利用地和草地，各占和。

（2）林地的增加量主要来自于草地和耕地，各有和，有1%的林地转换为草地。

（3）草地的增加量中有和 来自于耕地和林地，减少量中有 和分别转化为林地和耕地。

（4）水域有转化为城乡用地，有转化为未利用地，有转化为林地，表明水体的减少主要受人为因素影响。

（5）城乡用地的增加主要来自于耕地、未利用地和水域，分别占 、和。

（6）城乡用地、耕地、草地和林地各有 、、和 转化为未利用地，表明本区域土地利用状况不太好，大量的农田遭到占用却未合理开发。

4 结论

土地利用/土地覆被变化（LUCC）研究旨在深入了解土地利用状况及其变化的动态过程、主导原因及演化机理，以提高人们对土地利用变化的预测、管理、决策和调控能力，对于其数据的获取和分析处理十分重要。通过对本次研究获得数据的分析，可以得到如下结论：

（1）1999～2001年间云南丽江坝子地区的土地利用状况发生一定的变化。3年来耕地面积增加了，城乡用地面积增加了，水体面积减少了，林地和草地基本保持不变；各类型的转化突出表现为水域向城乡用地、未利用土地向城乡用地和耕地转换。

（2）土地利用的年变化率为 。其中城乡用地增加速度最快，以平均每年的速度递增，未利用地则以的速度递减，林地和草地基本保持不变。

（3）土地利用类型相互间的转化突出表现为水域向城乡用地，未利用土地向城乡用地和耕地转换。3年来本区域共有 即的水域面积转换为城乡用地，未利用土地中有 和 转化为城乡用地和耕地，分别占未利用土地面积的和。

参考文献

［1］朱会议，李秀彬，何书金等.环渤海地区土地利用的时空变化分析［J］.地理学报，2001，56 （3）：253～260

［2］王思远，刘纪元，张增祥等.近10年中国土地利用格局及其演变［J］.地理学报，2002，57 （5）：523～530

［3］丽江地区地方志办公室.丽江年鉴（1997） ［M］.昆明：云南民族出版社，1997，293～294

［4］李天宏，韩鹏.厦门市土地利用/覆盖动态变化的遥感检测与分析［J］.地理科学，2001，21 （6）：537～544

［5］村井俊治等.遥感精解［M］.北京.测绘出版社，1993，12：194～196

［6］王秀兰，包玉海.土地利用动态变化研究方法探讨［J］.地理科学进展，1999，18 （1）：81～87

［7］李忠峰，王一谋，王建华等.基于RS与GIS的榆林地区土地利用变化分析.水土保持学报，2003，6：97～99

[1]赵军,李东平,朱瑜馨. 基于Map Info的森林资源调查管理信息系统的设计与实现[A]. 中国地理学会地图学与GIS专业委员会.地图学与GIS学术讨论会论文集[C].中国地理学会地图学与GIS专业委员会:,2002:6.[2]赵国虎,高锦章,赵军,索南,孙看军,杨武. 吗啡对B-Z振荡反应的影响[A]. 甘肃省化学会.甘肃省化学会成立六十周年学术报告会暨二十三届年会论文集[C].甘肃省化学会:,2003:3.[3]赵军,武江民,赵生龙. 中国地理信息技术普及教育现状与对策刍议[A]. 中国地理学会地图学与地理信息系统专业委员会、福建师范大学地球信息科学研究中心、福建省空间信息工程研究中心.全国地图学与GIS学术会议论文集[C].中国地理学会地图学与地理信息系统专业委员会、福建师范大学地球信息科学研究中心、福建省空间信息工程研究中心:,2004:4.[4]赵军,陶明娟. 国外在野生动物伯劳鸟保护中对地理信息技术的利用[A]. 中国地理学会地图学与地理信息系统专业委员会、福建师范大学地球信息科学研究中心、福建省空间信息工程研究中心.全国地图学与GIS学术会议论文集[C].中国地理学会地图学与地理信息系统专业委员会、福建师范大学地球信息科学研究中心、福建省空间信息工程研究中心:,2004:5.[5]汪明冲,潘竟虎,赵军,孙晓茹. 基于GIS与RS的土壤侵蚀变化定量监测——以黄土高原水保二期世行贷款庆城项目区为例[A]. .《自然地理学与生态建设》论文集[C].:,2006:7.[6]贾艳红,赵军,南忠仁,赵传燕. 基于格网GIS的草原生态安全评价研究初探——以甘肃牧区为例[A]. .《自然地理学与生态建设》论文集[C].:,2006:7.[7]甄计国,霍静,赵军,潘竟虎,李伟元,何瑞东. 我国西北工业与旅游城市可持续发展与节约集约用地研究——以甘肃省嘉峪关市为例[A]. 中国土地学会.2007年中国土地学会年会论文集[C].中国土地学会:,2007:12.[8]符海月,赵军,李满春. 从Google Maps看我国全球化地理信息服务面临的挑战和对策[J]. 地理与地理信息科学,2006,02:113-115.[9]符海月,李满春,赵军,刘永学. 人口数据格网化模型研究进展综述[J]. 人文地理,2006,03:115-119+114.[10]孙晓茹,赵军. Google Maps API在WEBGIS中的应用[J]. 微计算机信息,2006,19:224-226.[11]贾艳红,赵军,南忠仁,赵传燕,王胜利. 基于熵权法的草原生态安全评价——以甘肃牧区为例[J]. 生态学杂志,2006,08:1003-1008.[12]汪明冲,赵军,李玉琳. 空间数据库引擎及其解决方案分析[J]. 地理信息世界,2006,04:63-66.[13]李红星,赵军. 基于GIS的区域环境噪声评价方法研究[J]. 遥感技术与应用,2006,04:363-367.[14]贾艳红,赵军,南忠仁,赵传燕. 熵权法在草原生态安全评价研究中的应用——以甘肃牧区为例[J]. 干旱区资源与环境,2007,01:17-21.[15]李传华,潘竟虎,赵军. 基于GIS的兰州市普通住宅价格空间格局分析[J]. 云南地理环境研究,2007,02:88-91.[16]赵军,刘琳,李霞. 基于RS与GIS的半干旱区城市绿化三维量测算研究——以兰州市安宁区为例[J]. 地理与地理信息科学,2007,04:37-40.[17]杜怀玉,赵军,冯翠琴. 西北干旱区经济社会发展与土地利用相关分析—以嘉峪关市为例[J]. 干旱区资源与环境,2007,11:90-94.[18]赵军,李旺平,李飞. 黄土高原太阳总辐射气候学计算及特征分析[J]. 干旱区研究,2008,01:53-58.[19]潘竟虎,赵军,董晓峰. Geo-CA与GIS在区域土壤侵蚀动态监测与模拟中的应用[J]. 干旱区研究,2008,01:131-136.[20]李飞,赵军,赵传燕,郝君明,郑佳佳. 中国潜在植被空间分布格局[J]. 生态学报,2008,11:5347-5355.[21]赵军,魏伟,冯翠芹. 天祝草原景观格局分析及景观利用格局优化[J]. 资源科学,2008,02:281-287.[22]魏伟,赵军,王旭峰. 天祝高寒草原区NDVI,DEM与地表覆盖的空间关系[J]. 干旱区研究,2008,03:394-401.[23]李娟,赵军. 基于DEM的西北干旱区典型地貌类型坡度提取分析[J]. 遥感技术与应用,2008,02:214-218.[24]张斌才,赵军. 大气污染扩散的高斯烟羽模型及其GIS集成研究[J]. 环境监测管理与技术,2008,05:17-19+55.[25]甄计国,赵军. 区域积温插值的GIS方法[J]. 冰川冻土,2005,04:591-597.[26]赵军,陶明娟. 兰州市2002年生态足迹计算与可持续发展状况分析[J]. 地域研究与开发,2005,06:113-116+128.[27]刘振波,倪绍祥,赵军. 绿洲生态预警信息系统初步设计[J]. 干旱区地理,2004,01:19-23.[28]赵军,胡秀芳. 区域生态安全与构筑我国21世纪国家安全体系的策略[J]. 干旱区资源与环境,2004,02:1-4.[29]贾艳红,赵军. 白银市区域生态环境质量评价研究[J]. 西北师范大学学报(自然科学版),2004,04:91-95.[30]刘振波,赵军,倪绍祥. 绿洲生态环境质量评价指标体系研究——以张掖市绿洲为例[J]. 干旱区地理,2004,04:580-585.[31]赵军,武江民,赵生龙. 我国地理信息技术普及教育现状与对策刍议[J]. 地理信息世界,2005,01:39-41+50.[32]赵军,刘春雨,潘竟虎,刘英英,杨东辉. 基于MODIS数据的甘南草原区域蒸散发量时空格局分析[J]. 资源科学,2011,02:341-346.[33]李飞,赵军,赵传燕,张小强. 中国干旱半干旱区潜在植被演替[J]. 生态学报,2011,03:689-697.[34]胡秀芳,赵军,查书平,郑雪萍. 草原生态安全模糊评价方法研究——以甘肃天祝高寒草原为例[J]. 干旱区资源与环境,2011,03:71-77.[35]李飞,赵军,赵传燕,王旭峰. 中国西北干旱区潜在植被模拟与动态变化分析[J]. 草业学报,2011,04:42-50.[36]赵军,师银芳,王大为,付鹏. 1961～2008年中国大陆极端气温时空变化分析[J]. 干旱区资源与环境,2012,03:52-56.[37]魏伟,石培基,赵军,王旭峰. 石羊河流域海拔、植被覆盖与景观类型空间关系研究[J]. 干旱区地理,2012,01:91-98.[38]赵军,李霞. 中国农牧交错带研究进展[J]. 草业科学,2009,01:94-99.[39]李海亮,赵军. 草地遥感估产的原理与方法[J]. 草业科学,2009,03:34-38.[40]王世杰,赵军. 甘肃省区域经济时空差异GIS-ESDA分析[J]. 干旱区资源与环境,2009,08:5-8.[41]韦莉,赵军,潘竟虎,李霞. 基于MODIS数据的黄土高原草地净初级生产力的估算研究[J]. 遥感技术与应用,2009,05:660-664+553.[42]魏伟,赵军,王旭峰. GIS、RS支持下的石羊河流域景观利用优化研究[J]. 地理科学,2009,05:750-754.[43]赵军,韦莉,陈姗. 石羊河流域上游生态系统服务价值的变化研究[J]. 干旱区资源与环境,2010,01:36-40.[44]魏伟,赵军,王旭峰. 石羊河流域土地利用类型景观异质性[J]. 生态学杂志,2010,04:760-765.[45]索俊锋,赵军. 基于CLUE-S模型的千岛湖镇土地利用动态模拟分析[J]. 长江流域资源与环境,2010,05:522-528.[46]赵军. 兰州城市边缘区土地利用及其变化研究[J]. 西北师范大学学报(自然科学版),2000,03:67-71.[47]梁会民,赵军. 基于GIS的黄土塬区居民点空间分布研究[J]. 人文地理,2001,06:81-83.[48]梁会民,赵军. 小流域综合治理的生态经济效益评估研究[J]. 生态经济,2001,08:12-14.[49]赵军,符海月. GIS在人口重心迁移研究中的应用[J]. 测绘工程,2001,03:41-43.[50]梁会民,赵军. 地理信息系统在居民点空间分布研究中的应用[J]. 西北师范大学学报(自然科学版),2001,02:76-80.[1]赵军,王雪平,金蓉. ArcGIS在草地资源规划中的应用[J]. 草业科学,2006,03:1-4.[2]赵军,付金霞. 虚拟水理论在河西走廊的应用研究[J]. 人民黄河,2006,02:38-40.[3]甄计国,王义德,赵军. 兴隆山国家级自然保护区的植被指数及其变化特征[J]. 遥感技术与应用,2006,04:294-301.[4]付金霞,赵军. 地理空间数据的空间性分析应用[J]. 地球信息科学,2006,04:65-69.[5]付彩菊,潘竟虎,赵军. 基于RS和GIS的额济纳旗土地利用变化研究[J]. 国土资源科技管理,2006,06:71-75.[6]胡秀芳,赵军,钱鹏,李东平. 草原生态安全理论与评价研究[J]. 干旱区资源与环境,2007,04:93-97.[7]冯翠芹,赵军. GIS在现代景观生态研究中的应用[J]. 安徽农业科学,2007,19:5951-5952.[8]潘竟虎,任皓晨,赵军. 甘肃省土地利用状况的计量地理分析[J]. 农业系统科学与综合研究,2007,03:377-381+384.[9]汪明冲,潘竟虎,赵军. 基于GIS与RS的土壤侵蚀变化定量监测——以黄土高原水保二期世行贷款庆城项目区为例[J]. 干旱地区农业研究,2007,06:116-121.[10]王世杰,赵军. 全球化背景下的地理信息服务发展现状及思考[J]. 矿山测量,2007,04:48-50+4.[11]李娟,赵军. 基于WebGIS的遥感影像数据发布系统的设计与实现[J]. 北京测绘,2008,04:60-63.[12]潘竟虎,赵军. 高师遥感课程实践教学的改革[J]. 理工高教研究,2008,01:118-120.[13]李娟,赵军. 基于DEM的西北干旱区坡度提取分析[J]. 北京测绘,2008,01:42-45.[14]党国锋,赵军. 兰州市城市化水平与耕地面积变化的定量研究[J]. 干旱区资源与环境,2008,04:17-21.[15]潘竟虎,刘普幸,赵军. 黑河下游土地利用与景观格局时空特征分析[J]. 土壤,2008,02:306-311.[16]田凤宾,赵军. 基于Google Maps的同城生活信息平台设计策略研究[J]. 电脑知识与技术,2008,16:1340-1341.[17]张俊,赵军. 基于Web Service的WebGIS应用系统模式及其与其他模式的比较[J]. 科学技术与工程,2008,17:4873-4876.[18]赵军,贾新会,甄计国. 石羊河流域人口经济压力定量分析[J]. 西北人口,2008,06:100-103.[19]孙建国,赵军,甄计国,李东平. 基于DTM的黄土丘陵沟壑区太阳辐射值计算模型及应用研究[J]. 测绘工程,2003,01:28-30.[20]赵军,党国锋. 祁连山草地资源利用面临的问题及治理对策[J]. 草业科学,2003,07:1-4.[21]陶明娟,赵军. 兰州市2002年可持续发展状况的生态足迹分析[J]. 云南地理环境研究,2005,02:51-54.[22]赵军,贾艳红. 国外GIS在地下水管理与研究中的应用和启示[J]. 地下水,2005,03:166-168.[23]刘琳,赵军,付金霞,王雪平. 我国航空航天遥感技术进展初探[J]. 华北航天工业学院学报,2005,03:10-13.[24]赵军,刘琳,闫浩文. 基于GIS的危险废物管理信息系统设计研究——以兰州市危险废物地理信息系统为例[J]. 西北师范大学学报(自然科学版),2005,05:75-78.[25]刘振波,赵军,倪绍祥. 基于GIS的区域发展均衡性测度研究——以甘肃省为例[J]. 地域研究与开发,2003,02:35-38.[26]孙建国,赵军,闫浩文. 一个空间数据挖掘原型系统的设计与实现[J]. 兰州交通大学学报,2004,04:94-97.[27]朱瑜馨,张锦宗,赵军. 基于人工神经网络的森林资源预测模型研究[J]. 干旱区资源与环境,2005,01:101-104.[28]赵军,许兆霞. 中国人口增长率与人口文化素质相关分析[J]. 干旱区资源与环境,2005,02:1-4.[29]赵军,胡自治. 从生态信息图谱的角度看草原综合顺序分类法检索图[J]. 草原与草坪,2005,02:12-14.[30]魏伟,赵军,王旭峰. RS、GIS支持下高寒草原退化过程中土壤流失及治理模式仿真研究[J]. 干旱区地理,2009,01:47-56.[31]赵军,田英. 区域人口经济压力的空间格局演变分析——以甘肃省为例[J]. 干旱区资源与环境,2009,03:15-19.[32]赵军,任皓晨,赵传燕,张勃,孟宝. 黑河流域土壤含水量遥感反演及不同地类土壤水分效应分析[J]. 干旱区资源与环境,2009,08:139-144.[33]赵军,郑珊,胡秀芳. 基于GIS的天祝高寒草原生态安全模糊评价研究[J]. 干旱区资源与环境,2010,04:66-71.[34]李飞,王春,赵军,郑佳佳. 中国陆地多年平均积温空间化研究[J]. 自然资源学报,2010,05:778-784.[35]王世杰,赵军. 基于GIS的甘肃省区域经济时空差异研究[J]. 测绘科学,2010,02:170-172.[36]赵军,杨东辉,潘竟虎. 基于空间化技术和土地利用的兰州市GDP空间格局研究[J]. 西北师范大学学报(自然科学版),2010,05:92-96+102.[37]赵军,杨梅,杨东辉. 基于GIS的甘南高原积温资源空间插值分析[J]. 甘肃科技,2010,13:40-42.[38]杨东辉,赵军,张智慧,付鹏,刘春雨. 近10年甘南牧区草地净初级生产力变化研究[J]. 干旱地区农业研究,2011,01:257-263.[39]赵军,刘喆,赵媛. 生态脆弱区生态旅游开发中的生态保护[J]. 甘肃科技,2011,04:8-10.[40]赵军,陈姗. 基于GIS的石羊河流域森林生态效益估算与空间分布特征研究[J]. 干旱区地理,2011,01:12-19.[41]魏伟,石培基,赵军,王旭峰. 基于GIS和RS的石羊河流域景观功能分区[J]. 干旱区资源与环境,2011,08:47-52.[42]刘春雨,赵军,刘英英,党国锋,师银芳. 25年来额济纳天然绿洲LUCC及景观格局时空变化[J]. 干旱区资源与环境,2011,09:32-38.[43]赵媛,赵军. 基于移动GIS和智能手机的电子导游系统发展及其关键技术[J]. 数字技术与应用,2011,08:51-52.[44]赵军,王小敏,李东成. 基于MODIS的民勤绿洲植被覆盖变化定量分析[J]. 干旱区资源与环境,2012,10:91-96.[45]朱瑜馨,赵军,曹静. 祁连山山地生态系统稳定性评估模型[J]. 干旱区研究,2002,04:33-37.[46]李东平,赵军. WebGIS主要技术及其发展趋势[J]. 现代计算机(专业版),2001,12:23-25.[47]赵军,李东平,朱瑜馨. 可视化编程语言实现GIS软件的二次开发[J]. 测绘通报,2002,S1:68-70.[48]田英,赵军,管信龙. 1990年～2005年甘肃省人口经济压力空间格局及变化分析[J]. 资源环境与发展,2008,01:19-23+27.[49]张勃,石培基,赵军. 甘肃石羊河流域武威绿洲水资源系统分析[J]. 西北师范大学学报(自然科学版),1994,03:73-79.[50]李传华,赵军. 基于GIS的民勤县生态环境脆弱性演化研究[J]. 中国沙漠,2013,01:302-307.