李满春的南京大学科技处处长

基于边缘的多光谱遥感图像分割方法[J].遥感学报.2006,10(3):350-356.

基于特征的土地利用空间数据库模型[J].地理与地理信息科学,2006,22(2):55-58.

基于Agent技术的分布式虚拟地理信息系统研究[J].测绘科学,2006,31 (6):121-122,79.

网络地理信息服务构建初步研究[J].遥感信息,2006 .1 :46-49,63.

基于GIS的土地利用规划公众参与研究[J].现代测绘,2006 ,29(3):7-10.

基于Agent技术的分布式虚拟地理信息系统研究[J].测绘科学.2006,31(6):121-122.

区域土地资源集约利用及其评价指标体系构建研究[J].地域研究与开发,2006,25(4): 117-119,124.

基于矢栅混合数据模型的土地适宜性评价研究.长江流域资源与环境,2006,15(3): 320-324.

标准农田规划管理信息系统的设计与实现.遥感信息,2006,2: 63-66.

人口数据格网化模型研究进展综述.人文地理,2006 ,21(3): 115-119,114.

土地利用规划在城市郊区历史遗迹保护中的应用——以无锡市新区为例.北京第二外国语学院学报,2006, 5: 43-46,22.

GIS支持的县域人口迁移空间模型研究--以浙江省临安市为例.长江流域资源与环境,2006 ,15 (3): 281-286.

中国海岸带专用地图投影设计[J].测绘学报.2006,34(4):395-399.

基于矢量数据模型的模糊地理对象建模[J].地理与地理信息科学.2006,22(6):66-70.

海南岛西部土地沙漠化时空模拟及虚拟政策实验[J].地理与地理信息科学.2006,22(6):66-70.

从Google Maps看我国全球化地理信息服务面临的挑战和对策[J].地理与地理信息科学.2006,22(2):1-3.

乡镇土地利用总体规划对农村居民出行可达性的影响研究[J]。地理与地理信息科学.2006,22(2):62-66.

基于SOA的WebGIS框架探索[J].计算机应用.2006,(9):2225-2228.

基于最短路径距离的宗地地价点状因素自动化评价[J].计算机应用.2006,(9):143-145,178.

基于名义尺度栅格数据空间自相关测度初探[J].地理与地理信息科学.2006,(3):53-56.

面向离散点的空间权重矩阵生成算法与实证研究[J] .地理与地理信息科学.2006,(3):53-56.

GIS支持的县域人口迁移空间模型研究——以浙江省临安市为例[J].长江流域资源与环境.2006,15(3):281-286.

GIS支持的常州市区域经济空间集聚特征研究[J].2006,29(2):3-5.

基于Web的虚拟地景信息系统研究[J].科技通报,2005（6）.

一种基于面积紧凑度的二维空间形状指数及其应用[J].地理与地理信息科学,2005（5）.

基于一维Kohonen网络的遥感影像分类方法研究[J].遥感信息,2005（5）.

土地利用规划建库与制图一体化方案研究[J].江西师范大学学报(自然科学版),2005（4）.

基于特征的空间数据库模型设计与实现[M]。《国土整治与土地资源可持续利用》，中国大地出版社，2005年9月。

区域土地利用战略研究[M]。《国土整治与土地资源可持续利用》，2005年9月中国大地出版社，2005年9月。

计算机辅助土地估价报告编制系统设计与实现 [J].微型电脑应用,2005（1）.

设计模式在GeoMedia WebMap地图数据发布中的应用[J].计算机应用研究,2005（4）.

江苏省耕地后备资源的空间分异及开发时序模型研究 [J].地理与地理信息科学,2005（2）.

基于生命机制概念的城市土地利用演化模型研究[J].经济地理,2005 Vol.25 No.4.

AMSR被动微波数据介绍及主要应用研究领域分析[J].遥感信息,2005（3）.

B/S结构信息系统中可重用界面的实现[J].武汉理工大学学报,2005（4）.

基于SuperMap的专题地图制图技术与方法实践[J].现代测绘,2005（3）.

基于地统计学的宁波市区地价空间分布特征研究

先后主持国家自然科学基金项目3项、国家863重点项目课题1项、863项目子课题1项；主要负责国土资源大调查项目、国家行业公益项目等项目。共发表学术论文100余篇，其中SCI论文20余篇、计算机与技术学科类EI论文30余篇、一级学报论文10余篇；参与制定国家行业标准1部；申请国家发明专利17项，其中已授权9项（其中，第一发明人3项）；登记软件著作权7项。先后获国家科学技术进步奖二等奖（2014，4/15），地理信息科技进步奖一等奖（国家测绘地理信息局，2013，4/10）、测绘科技进步奖一等奖（国家测绘局，2009，3/9）、“十一五”测绘地理信息优秀青年贡献奖（国家测绘地理信息局，2011）、测绘地理信息教学成果奖二等奖（国家测绘地理信息局，3/7）等奖项。2012年入选教育部“新世纪优秀人才支持计划”。

以第一作者、通讯作者发表SCI论文18篇

1. Y.X. Liu*, C. Sun, Y.H. Yang, M.X. Zhou, W.F. Zhan, W.Y. Cheng. 2015. Automatic extraction of offshore platforms using time-series Landsat-8 Operational Land Imager data. Remote sensing of environment. 10.1016/j.rse.2015.12.047.

2. Y.X.Liu*, M.X. Zhou, S.S. Zhao, W.F. Zhan, K. Yang, M.C. Li. 2015. Automated extraction of tidal creeks from airborne laser altimetry data. Journal of Hydrology 527, 1006-1020

3. Y.X.Liu*, J.L. Chen, W.Y. Cheng, C. Sun, S.S. Zhao, Y Pu. 2014. Spatiotemporal dynamics of the urban sprawl in a typical urban agglomeration: a case study on Southern Jiangsu, China (1983–2007). Frontiers of Earth Science 8 (4), 490-504

4. Y.X.Liu*, M.C. Li, L. Mao, L. Cheng, K.F. Chen. 2013. Seasonal pattern of tidal-flat topography along the Jiangsu middle coast, China, using HJ-1 optical images. Wetlands 33 (5), 871-886

5. Y.X.Liu*, M.C. Li, M.X. Zhou, K. Yang, L. Mao. 2013. Quantitative Analysis of the Waterline Method for Topographical Mapping of Tidal Flats: A Case Study in the Dongsha Sandbank, China. Remote Sensing 5 (11), 6138-6158

6. Y.X.Liu*, M.C. Li, L. Cheng, F.X. Li, K.F. Chen. 2012. Topographic mapping of offshore sandbank tidal flats using the waterline detection method: A case study on the Dongsha Sandbank of Jiangsu Radial Tidal Sand Ridges, China. Marine Geodesy 35 (4), 362-378

7. Y.X. Liu*, M.C. Li, L. Mao, L. Cheng, F.X. Li. 2012. Toward a method of constructing tidal flat digital elevation models with MODIS and medium-resolution satellite images. Journal of Coastal Research 29 (2), 438-448

8. Y.X. Liu*, M.C. Li, L. Chen, F.X. Li, Y.M. Su. 2010. A DEM Inversion Method for Inter-tidal Zone based on MODIS Dataset: A Case Study in the Dongsha Sandbank of Jiangsu Radial Tidal Sand-Ridges, China. China Ocean Engineering 24 (4), 735-748

9. Y.X. Liu*, M.C. Li, L. Mao, F.F. Xu, S. Huang. 2006. Review of remotely sensed imagery classification patterns based on object-oriented image analysis. Chinese Geographical Science 16 (3), 282-288

10. Y.W. Duan, Y.X. Liu*, M.C. Li, et al. Reefs survey of the Nansha Islands based on Landsat 8 OLI images. Acta Oceanologica Sinica. Accept.

11. S.Y. Zhang, Y.X. Liu*, Y.H. Yang, C. Sun, F.X. Li. Erosion and deposition within Poyang Lake: Evidence from a decade of satellite data. Journal of great lakes research. Accept. 10.1016/j.jglr.2015.12.012

12. C. Sun, Y.X. Liu*, S.S. Zhao, M.X. Zhou, Y.H. Yang, F.X. Li. 2016. Classification mapping and species identification of salt marshes based on a short-time interval NDVI time-series from HJ-1 optical imagery. International Journal of Applied Earth Observation and Geoinformation 45, 27-41

13. Y.H. Yang, Y.X. Liu*, M.X. Zhou, S.Y. Zhang, W.F. Zhan, C. Sun, Y.W. Duan. 2015. Landsat 8 OLI image based terrestrial water extraction from heterogeneous backgrounds using a reflectance homogenization approach. Remote Sensing of Environment 171, 14-32

14. S.S. Zhao, Y.X. Liu*, M.C. Li. 2015. Analysis of Jiangsu Mudflats Reclamation from 1974 to 2012 using Remote Sensing. China Ocean Engineering 29 (1), 143-154

15. Y.H. Yang, Y.X. Liu*, M.X. Zhou, F.X. Li, C. Sun. 2015. Robustness assessment of urban rail transit based on complex network theory: a case study of the Beijing Subway. Safety science 79, 149-162

16. K. Yang, M.C. Li, Y.X. Liu*, L. Cheng, Y.W. Duan, M.X. Zhou. 2014. River delineation from remotely sensed imagery using a multi-scale classification approach. IEEE Selected Topics in Applied Earth Observations and Remote Sensing.

17. J.S. Wang, M.C. Li, Y.X. Liu*, H.X. Zhang, W. Zou, L. Cheng. 2014. Safety assessment of shipping routes in the South China Sea based on the fuzzy analytic hierarchy process. Safety science 62, 46-57

18. Y.C. Wang, Y.X. Liu*, M.C. Li, L. Tan. 2014. The reconstruction of abnormal segments in HJ-1A/B NDVI time series using MODIS: a statistical method. International Journal of Remote Sensing 35 (23), 7991-8007

以第一作者、通讯作者发表中文核心期刊论文34篇

1. 孙超, 刘永学*, 李满春, 赵赛帅, 周旻曦. 近25a来江苏中部沿海盐沼分布时空演变及围垦影响分析. 自然资源学报. 2015, 30(9):1486-1498.

2. 周旻曦, 刘永学*, 李满春, 等. 多目标珊瑚岛礁地貌遥感信息提取方法——以西沙永乐环礁为例. 地理研究, 2015, 34(4): 677-690.

3. 王加胜, 李满春, 刘永学*, 等. 南沙群岛历史事件时空演化的可视化表达. 地理科学, 2015, 35(5): 608-614.

4. 蒋洁, 刘永学*, 钟礼山, 李满春, 孙超. SAR影像时间序列支持下的互花米草盐沼提取——以江苏盐城东部沿海为例. 测绘科学技术学报. 2015, 32(3):271-276.

5. 邹伟, 刘永学*, 李满春, 等. 基于 WebGIS 的南海历史地理情势分析系统设计与实现. 热带地理, 2015, 35(2): 202-207.

6. 赵赛帅, 刘永学*, 李满春, 张荷霞, 孙超. 基于AHP-变异系数法的越占南沙岛礁战略价值评价. 海洋科学. 2015, 39(6):114-121.

7. 蒋洁, 刘永学*, 李满春, 等. 南海岛礁风能资源及风力发电评价——基于 QuikSCAT 风场数据与 Landsat ETM+ 影像. 资源科学, 2014, 36(1): 139-147.

8. 邹伟, 刘永学*, 李满春, 等. 网络新闻中黄岩岛争端事件舆情研究——以新浪网“中菲黄岩岛争端”专题为例. 现代图书情报技术, 2014, 30(2): 72-78.

9. 邹伟, 刘永学*, 李满春, 等. 南海历史地理争端空间分布与关联性研究. 地球信息科学学报, 2014, 16(2): 249-256.

10. 孙超, 刘永学*, 李满春, 等. 近 35 年来热带风暴对我国南海岛礁的影响分析. 国土资源遥感, 2014, 26(3): 135-140.

11. 蒋洁, 刘永学*, 李满春, 等. 基于 ENVISAT ASAR 影像与 QuikSCAT 风场数据的近海风能资源评价——以香港东南海域为例. 地理与地理信息科学, 2014, 30(2).

12. 张荷霞, 刘永学*, 李满春, 等. 基于 AHP 和 EWCM 的部分南沙岛礁战略价值模糊综合评价. 海洋通报, 2014, 4: 003.

13. 汪业成, 刘永学*, 李满春, 等. 基于场强模型的南沙岛礁战略地位评价. 地理研究, 2013, 32(12).

14. 张荷霞, 刘永学*, 李满春, 等. 南海中南部海域油气资源开发战略价值评价. 资源科学, 2013, 35(11).

15. 赵赛帅, 刘永学*, 李满春, 等. 基于 AHP-变异系数法的越占南沙岛礁战略价值评价. 海洋科学, 2015, 39(6): 114-121.

16. 王加胜, 刘永学*, 李满春, 等. 基于 ENVISAT ASAR 的海洋钻井平台遥感检测方法. 地理研究, 2013, 32(11).

17. 成王玉, 刘永学*,李满春, 等. 基于 AHP 与模糊综合评价方法的南沙东部岛礁战略价值评价. 热带地理, 2013, 33(4): 381-386.

18. 张荷霞, 刘永学*,李满春, 等. 基于 JASON-1 资料的南海海域海面风浪场特征分析. 地理与地理信息科学, 2013, 29(5).

19. 洪武扬, 刘永学*, 李满春, 刘敏, 童礼华. 土地利用总体规划实施评估方法研究 . 国土资源科技管理, 2013, 30(5):68-73.

20. 杨康, 李满春, 刘永学, 程亮, 陈焱明.基于累积相似度表面的空间权重矩阵构建方法. 测绘学报, 2012, 41(2):259-265,272

21. 江冲亚, 李满春, 刘永学. 海岸带水体遥感信息全自动提取方法. 测绘学报. 2011, 40(3): 332-337.

22. 杨康, 李满春, 刘永学, 程亮, 江冲亚. 遥感影像道路的多点同时快速行进提取方法.遥感技术与应用. 2011, 26(3):295-302.

23. 胡炜,刘永学*,李满春,陈洁丽,毛鹍. 基于局部相关分析法的ETM+影像修复方法研究. 地理与地理信息科学.2011(5):29-32.

24. 江冲亚, 李满春, 李飞雪, 李雪, 刘永学. 地理信息系统(GIS)科学体验区建设. 地理信息世界. 2010, 4(2):84-89.

25. 陈洁丽,刘永学*,李满春,申成磊,胡炜,蔡文婷.一种基于遥感数据快速提取居民地信息的新方法.地理与地理信息科学. 2010(5):72-75.

26. 陈洁丽,刘永学*,李满春,申成磊,张栋,蔡文婷.基于归一化指数分析的居民地遥感信息提取.测绘科学.2010(2): 204-206.

27. 刘永学,李满春,刘国洪. 地理信息系统安全初探.遥感信息.2007(2):71-76.

28. 刘永学,李满春,毛亮.基于边缘的多光谱遥感图像分割.遥感学报.2006,10(3):350-356.

29. 刘永学,李满春,高月明.基于角色的土地利用规划管理信息系统访问控制模型设计与实现.地理与地理信息科学,2006.(2):4-6.

30. 刘永学,张忍顺,李满春.应用卫星影像系列海图叠合法分析沙洲动态变化——以江苏东沙为例.地理科学.2004, 24(2):199-204.（学报）

31. 刘永学,张忍顺,李满春.质心分析法在小沙洲动态演化分析中的应用——以江苏辐射沙洲亮月沙为例.海洋通报.2004, 23(1):69-75.

32. 刘永学,张忍顺,李满春.江苏淤泥质潮滩遥感影像特征与信息提取方法研究.遥感信息.2004,(1):23-26.

33. 刘永学,李满春,张忍顺.江苏辐射沙洲水边线自动提取方法研究.海洋科学.2004, (6):42-45.

34. 刘永学,张忍顺,李满春.江苏淤泥质潮滩地物信息遥感提取方法研究.海洋科学进展.2004,22(2).210-214.

35. 刘永学,张忍顺,李满春.江苏沿海互花米草盐沼动态变化及影响因素研究.湿地科学. 2004, 2(2).116-121.

36. 刘永学,张忍顺等.粒径趋势分析法在细颗粒潮坪上的应用.南京师大学报:自然科学版.2001, 24(1):105-109,113.

37. 刘永学,陈君等.江苏海岸盐沼植被演替的遥感图像分析.农村生态环境.2001, 17(3):39-41.

楼立明1 冯秀丽2

（1.宁波市国土资源局镇海分局，宁波，3151202.宁波大学建工学院，宁波，315120）

摘要：城市地价是一个具有时空性质的多维概念，在空间分布上具有较强的关联性和特殊性。本文以宁波市中心城区为研究区域，以地价信息为研究对象，探讨了如何基于地统计学和 GIS 对地价的空间分布特征和规律进行研究的原理和方法。

关键词：城市地价；空间分析；地统计学；宁波市区

1 地统计学的基本概念

地统计学（Geostatistics）是由法国著名数学家G·Matheron教授在研究了南非地质工程师 D.G.Krike 等人工作的基础上，于1963年提出并创立的。地统计学是在地质分析和统计分析的基础上形成的一套分析空间相关变量的理论和方法，它以区域化变量理论为基础，以变差函数为主要工具，研究那些在空间分布上既有随机性又有结构性的自然现象的科学。地统计学能最大限度的利用野外调查所能提供的各种信息，例如样本位置、样本值和样本承载大小等；能利用稀疏的或无规律的空间数据。由于地统计学能够较准确地描述区域化变量的随机性和结构性变化，因而越来越受到重视，除成功应用于自然资源方面外，还广泛应用于环境科学、农林科学、水利科学和土地科学中。

在地价研究领域，地统计学的应用主要体现在三个方面：

一是定量区域化变量的空间相关性，二是对调查数据进行空间插值，三是分析空间数据的时空规律性。相对来说，对空间插值应用较多，由于在通常的地价调查中，野外调查所得的数据不能完全覆盖所要求的区域范围，需要应用地统计学的方法进行插值，将离散的采样点数据内插为连续的数据表面。

应用地统计学的最大好处是它能够在空间相关分析的基础上，利用稀疏的、无规律的调查数据，最大限度的揭示这些数据所能提供的空间信息。但是，地统计学应用在地价研究中还刚刚开始，还存在着诸如空间与时间的协同分析、样点数目、取样位置、方向、大小的设计等等，这些都值得进一步完善和改进。

2 地统计学分析的基本函数

在地统计学上，用于空间相关分析的函数主要是半方差函数（Variograms Function）、协方差函数（Covariance Function）、相关函数（Correlation），其中半方差函数是地统计学最常用的工具。除此之外，还有一般相对方差函数、交叉变差函数、成对相对方差函数、对数方差函数、广义方差函数、特征化方差函数（或指示方差函数）及散点图等，但应用较少，一般不适合对地价进行空间分析。

半方差函数定义为：区域化变量z （xi）和z （xi ＋h）的增量平方的数学期望，即区域化变量增量的方差。半方差函数既是距离h的函数，又是方向α的函数。其计算公式如下：

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

式中，γ（h）为半方差函数值，半方差函数曲线图（Semivariogram）是半方差函数γ（h）对距离 h 的坐标图形。N （h）是被及分隔的数据对的数量，z （xi）和z （xi ＋h）分别是在点xi 和 xi ＋ h处样本的测量值，h是两分隔样点的距离。

对于一个典型的空间聚集分布，半方差函数一般随着距离的增大而增大，亦即区域化变量的空间变异愈来愈大，空间相关性逐渐减小，但增加至某一值时，半方差函数不再增加而是保持稳定，这表示样点间已不存在空间相关关系。将半方差函数值不再增加时的距离称为空间依赖范围（Range of Spatial Dependence），简称变程或相关程（Range），用 a表示。此时的半方差函数值称为基台值（Sill），用C0 ＋C表示。半方差函数曲线在 Y 轴上的截距称为区域不连续值，亦称块金（Nugget）系数或者核方差，用 C0表示。C0的大小可以反映区域化变量的局部随机性大小。（基台一块金）/基台（即 C/（C0 ＋C））的大小可以反映空间变异在总变异中所占的比例，或用随机程度（块金/基台，即C0/（C0 ＋C））的大小反映研究范围内不是由地价的空间自相关引起的那部分变异在总变异中所占的比率，也就是地价随机性和结构性所占成分。

3 宁波市区地价的地统计学分析

3.1 宁波市区地价的地统计学分析范围及样点分布

本次宁波市区城镇土地定级范围包含了宁波的六个区，从土地利用方式和土地市场发育水平来看，都存在着较大的差异，特别是山区，土地交易现象极少，土地价格样点稀缺。从地统计学对样点的要求而言，虽然不要求地价样点规则取样，但是大面积内样点稀缺会对分析结果的可靠性造成很大影响，同时，考虑到宁波市区土地交易主要集中在以市三江片为核心的一个辐射圈内，因此，本次宁波市区地价的地统计学分析的范围定为宁波市三江片向外扩展的一个区域。分析范围及范围内地价样点情况如图1～图3 所示。

图1 分析范围内商业地价样点分布图

图2 分析范围内住宅地价样点分布图

图3 分析范围内工业地价样点分布图

3.2 地价在各向异性条件下的变异分析

地价作为一种区域化变量，在各个方向上都有变化。一个区域化变量如果在不同方向上都有变化，那么当变异函数r （h）在各个方向上的变化都相同时称为各向同性，反之称为各向异性。图4～图6为不同地价类型在0°，45°，90°和135°四个方向上的变异曲线图。

（1）工业地价没有表现出各向异性结构特征，在不同的方向上，不同距离的半方差函数值均不能用合适的模型来拟合出半方差曲线，这说明宁波市区工业用地的发展轴向不明确，同时工业用地价格的政策性因素也较大，造成地价规律性较差。

（2）住宅、商业两种地价均表现出一定的各向异性结构特征，在不同的方向上，块金值、基台值和变程均不相同，具有带状各向异性特征。在135°方向（西北—东南）上地价不同距离上的半方差函数值的拟合效果较好，说明近年来鄞州区中心区的建设对宁波城市商业和居住功能的分布格局产生了显著的影响。宁波市三江片西部和南部大量的居住区建设，形成了好又多、麦德龙等新的商业中心。

3.3 地价在各向同性条件下的变异分析

为了不同地价间的对比和地价扩散情况的分析，往往需要将各向异性结构通过线性变换和矩阵变换转化为各向同性结构。其原理是通过改变不同方向上的距离h，使 γ（h）在各个方向上具有相同的变化情况。地统计学软件 GS＋提供了这种工具，使得具有各向异性的区域化变量可以转化为各向同性结构下进行研究。图7～图9 为不同地价类型在各向同性下的变异曲线图，图10～图12 为不同地价类经过kriging空间插值后的平面和相应的三维曲面图。表1 为不同地价类型各向同性下的变异曲线的模拟公式参数。

表1 不同地价类型各向同性下的变异曲线的模拟公式参数

图4 工业地价在0°，45°，90°，135°四个方向上的变异曲线图（指数模型）

图5 住宅地价在0°，45°，90°，135°四个方向上的变异曲线图（指数模型）

图6 工业地价在0°，45°，90°，135°四个方向上的变异曲线图（球状模型）

工业地价变异曲线模拟方程为：

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

图7 工业地价变异曲线图（球状模型）

住宅地价变异曲线模拟方程为：

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

商业地价变异曲线模拟方程为：

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

图8 住宅地价变异曲线图（球状模型）

图9 商业地价变异曲线图（球状模型）

图10 工业地价经过kriging 空间插值后的平面和相应的三维曲面图

图11 住宅地价经过kriging 空间插值后的平面和相应的三维曲面图

图12 商业地价经过kriging 空间插值后的平面和相应的三维曲面图

4 结论

（1）三种地价在空间一定范围内均存在着空间相关性，它们的空间相关距离为1810～3925 m。工业地价的空间相关距离最大，为3925 m；住宅次之，为2914 m；商业最小，为1810 m。说明土地价格空间变化的梯度是商业用地大于住宅用地，住宅用地大于工业用地。

（2）在三种地价空间变异的总方差中，均是结构方差（C）所占的比例要大于块金效应（C0）所占比例。这说明确定性因素（交通状况、基础设施、环境状况等）对地价的影响要大于随机因素引起的地价差异，地价的构成还是比较合理的。

（3）块金效应（C0），是住宅地价＞ 商业地价＞ 工业地价，这说明三种地价中住宅地价最容易受不确定性因素影响，价格变动最大，工业地价相对最稳定。这与宁波市房地产市场中住宅价格明显提高，近年来政府不时出台宏观调控政策的情况比较符合。

图13 宁波市区商业用地价格分布图

图14 宁波市区住宅用地价格分布图

（4）空间变异系数 C/C0 ＋C，商业地价达到 0.659，住宅为 0.807，工业为 0.874，这说明工业地价的空间变异性最强，受周围地价的影响最大。而商业地价和住宅地价受到宁波市近年来城市规划调整的影响，随着新规划的城市中心（如东部新城）、副中心（鄞州中心区）的建设，在空间上出现了不连续的、突变性特点。

（5）从通过空间插值得到的宁波市中心城区（三江片）的地价分布图，加入道路和河流等控制性基础因素（图13～图15），可以看出，宁波市商业地价的地域分异规律明显，不仅原来的中心城区老的城市建成区地价较高，而且整个城市往东发展和江北区往北发展，海曙区往西发展、鄞州区往西南发展带来的商业用地地价变化也很明显。宁波市区住宅用地价格和工业用地价格变化规律也得到了非常直观的反映。

图15 宁波市区工业用地价格分布图

参考文献

刘卫东等.宁波三江片城市土地价格调查.北京：科学出版社，2002

王政权.地统计学及在生态学中的应用.北京：科学出版社，1999

王瑞萍.GIS 空间分析技术应用研究.油气田地面工程，2003，22 （6）

陈浮，李满春等.城市地价空间分布因式的地统计学分析.南京大学学报（自然科学），1999，35 （6）

单卫东，包浩生.城市中空间扩散的各向异性研究.南京大学学报（自然科学），1996，32 （3）

许晓晖.上海市商品住宅价格空间分布特征分析.经济地理，1997，17 （3）