2025年地理信息系统工程师认证考试试题及答案

2025年地理信息系统工程师认证考试试题及答案一、单项选择题（每题1分，共20分。每题只有一个正确答案，请将正确选项字母填在括号内）1.在WebGIS中，用于实现客户端与服务器端异步通信的核心技术是（）A.CGI  B.AJAX  C.SOAP  D.WCS答案：B2.下列哪一项不是OGC简单要素标准定义的基本几何类型（）A.Point  B.LineString  C.MultiPatch  D.Polygon答案：C3.若某栅格数据的空间分辨率为10m，则其像元面积为（）A.10m²  B.20m²  C.100m²  D.1000m²答案：C4.在PostGIS中，用于计算两点球面最短距离的函数是（）A.ST_Distance  B.ST_Distance_Sphere  C.ST_Length  D.ST_Azimuth答案：B5.下列关于四叉树索引的描述，正确的是（）A.仅适用于矢量数据  B.是一种R树变种  C.空间分割呈递归正方形  D.不支持动态插入答案：C6.在ArcGISPro中，将2000国家大地坐标系（CGCS2000）数据转换为WGS84时，必须首先考虑（）A.七参数  B.三参数  C.椭球差异  D.投影带号答案：A7.下列哪项不是GeoPackage相对于Shapefile的显著优势（）A.单文件存储  B.支持拓扑规则  C.支持栅格与矢量混合  D.字段名可达255字符答案：B8.在QGIS中，使用“JoinbyLocation”工具时，若选择“intersect”谓词，则结果将包含（）A.仅完全包含的要素  B.仅边界重合的要素  C.任意部分重叠的要素  D.被包含的要素答案：C9.下列关于地图切片（Tile）金字塔的描述，错误的是（）A.第0级通常覆盖全球  B.每级分辨率减半  C.行列号从左上角原点开始  D.同一级所有切片分辨率相同答案：B10.在GDAL中，将GeoTIFF的NoData值设为9999并压缩为LZW格式的正确命令行是（）A.gdal_translatea_nodata9999coCOMPRESS=LZWin.tifout.tifB.gdalwarpdstnodata9999coCOMPRESS=LZWin.tifout.tifC.gdal_edit.pya_nodata9999coCOMPRESS=LZWin.tifD.gdaladdoa_nodata9999coCOMPRESS=LZWin.tifout.tif答案：A11.下列哪项不是三维GIS数据模型（）A.TIN  B.CSG  C.BRep  D.B+Tree答案：D12.在SuperMapiDesktop中，将二维线直接拉伸为三维管线的功能位于（）A.数据 →类型转换  B.三维分析 →拉伸建模  C.网络分析 →设施网络  D.空间分析 →缓冲区答案：B13.下列关于空间数据库事务隔离级别的说法，正确的是（）A.ReadUncommitted可防止幻读  B.ReadCommitted可防止不可重复读  C.Serializable性能最高  D.RepeatableRead可防止脏读答案：D14.在GeoServer中，配置SLD时，用于实现属性驱动符号化的元素是（）A.CssParameter  B.ogc:Filter  C.ExternalGraphic  D.AnchorPoint答案：B15.下列哪项不是移动GIS终端定位误差的主要来源（）A.多路径效应  B.电离层延迟  C.屏幕分辨率  D.星历误差答案：C16.在遥感影像分类中，使用混淆矩阵计算用户精度的公式是（）A.对角线像元数/该列总和  B.对角线像元数/该行总和  C.对角线像元数/总像元数  D.行和列总和之差/总像元数答案：B17.下列哪项不是LiDAR点云滤波算法（）A.渐进式三角网加密（PTD）  B.形态学开运算  C.区域生长  D.直方图均衡化答案：D18.在PostgreSQL中，为geometry列创建GIST空间索引的SQL语句是（）A.CREATEINDEXidxONtableUSINGGIST(geom);B.CREATESPATIALINDEXidxONtable(geom);C.ALTERTABLEtableADDINDEXidx(geom);D.CREATEINDEXidxONtableUSINGBTREE(geom);答案：A19.下列关于HDFS存储栅格数据的描述，正确的是（）A.支持低延迟随机写入  B.适合存储大量小文件  C.默认块大小为128MB  D.需通过OGCWFS访问答案：C20.在Python的geopandas中，将GeoDataFrame坐标系转换为EPSG:4524的正确写法是（）A.gdf.to_crs(epsg=4524)  B.gdf.set_crs(epsg=4524)  C.j4=4524  D.gdf.transform(4524)答案：A二、多项选择题（每题2分，共20分。每题有两个或两个以上正确答案，多选、少选、错选均不得分）21.下列哪些属于OGC网络服务规范（）A.WMS  B.WFS  C.WCS  D.WPS  E.WMTS答案：ABCDE22.在ArcGISPro的ModelBuilder中，可作为迭代器类型的有（）A.FeatureClass  B.Raster  C.Field  D.Workspace  E.SQLQuery答案：ABCD23.下列哪些操作会改变栅格数据的像元大小（）A.重采样  B.投影变换  C.裁剪  D.波段合成  E.设置金字塔答案：AB24.关于空间聚类算法DBSCAN的描述，正确的有（）A.基于密度  B.需预设簇数量  C.可识别噪声点  D.对密度变化敏感  E.基于距离阈值ε答案：ACDE25.下列哪些属于开源GIS软件（）A.GRASSGIS  B.MapServer  C.GeoMedia  D.uDig  E.ArcGISEngine答案：ABD26.在三维GIS中，倾斜摄影测量数据常见格式包括（）A.OSGB  B.S3C  C.3DTiles  D.LAS  E.TIN答案：ABC27.下列哪些属于空间数据质量元素（）A.位置精度  B.属性精度  C.现势性  D.完整性  E.一致性答案：ABCDE28.在PostGIS中，可用于栅格分析的函数有（）A.ST_MapAlgebra  B.ST_Union  C.ST_Clip  D.ST_Polygon  E.ST_Rescale答案：ACE29.下列哪些属于WebGIS前端可视化库（）A.OpenLayers  B.Leaflet  C.Cesium  D.Three.js  E.D3.js答案：ABCDE30.下列关于空间索引R树的描述，正确的有（）A.采用最小边界矩形递归分割  B.支持动态插入删除  C.叶节点存储要素指针  D.查询复杂度为O(n)  E.是平衡树答案：ABCE三、填空题（每空1分，共20分）31.在GDAL中，查询栅格数据投影信息的命令行参数是________。答案：proj432.在WMTS标准中，切片的行号起点位于________角。答案：左下33.我国2000国家大地坐标系采用的参考椭球是________。答案：CGCS2000椭球34.在PostgreSQL中，空间函数ST_Azimuth返回值的单位是________。答案：弧度35.在ArcGISPro中，用于批量导出线节点坐标的工具是________。答案：要素转表（FeatureToTable）36.在QGIS中，将CSV文件生成空间图层的插件名称是________。答案：AddDelimitedTextLayer37.在GeoServer中，配置跨域资源共享（CORS）需修改的文件是________。答案：web.xml38.在Python中，geopandas读取Shapefile的函数是________。答案：gpd.read_file()39.在LiDAR点云分类代码中，地面点的默认ASPRS分类值为________。答案：240.在遥感影像中，NDVI的计算公式为________。答案：(NIRRed)/(NIR+Red)41.在三维GIS中，CityGML中用于描述建筑外墙的语义表面类型名称是________。答案：WallSurface42.在PostGIS中，将geometry转换为WKT的函数是________。答案：ST_AsText43.在ArcGISPro的Python窗口中，引用当前地图的代码对象是________。答案：arcpy.mp.ArcGISProject("CURRENT").activeMap44.在地图切片中，GoogleMaps兼容的投影EPSG编号为________。答案：385745.在GDAL中，将矢量数据转换为栅格的命令行工具是________。答案：gdal_rasterize46.在QGIS中，实现热力图渲染的算法核心是基于________核密度估计。答案：二次（Quartic）47.在WFS请求中，返回GML3格式的参数是________。答案：outputFormat=GML348.在ModelBuilder中，将中间变量设为托管变量的快捷键是________。答案：右键→Managed49.在GeoPackage中，存储栅格瓦片的表名约定为________。答案：gpkg_2d_gridded_coverage50.在Cesium中，加载3DTiles的类名称是________。答案：Cesium3DTileset四、简答题（共5题，每题8分，共40分）51.（封闭型）简述空间数据引擎（SDE）的三级体系结构，并指出各层主要功能。答案：空间数据引擎采用客户端—应用服务器—数据服务器三级结构。（1）客户端：提供GIS桌面、Web或移动应用，负责用户交互与可视化；（2）应用服务器：托管ArcGISServer/GeoServer等中间件，完成空间查询、分析、制图与缓存；（3）数据服务器：运行Oracle/PostgreSQL等DBMS，通过SDE存储并管理矢量、栅格与元数据，提供事务、索引与并发控制。三级解耦实现跨平台、高性能与集中授权。52.（开放型）结合实际案例，说明如何利用时空立方体（SpaceTimeCube）进行城市共享单车热点演化分析，并指出至少两项可视化难点及解决思路。答案：案例：深圳市2020—2024年共享单车订单数据，每辆车每30s上报一次GPS。步骤：①数据清洗：剔除漂移点，聚合到250m×250m×1h的时空体素；②构建时空立方体：以x、y为空间轴，t为时间轴，体素值设为订单量；③热点探测：采用时空核密度估计（STKDE）识别连续高密度体素；④演化分析：利用时空扫描统计（SaTScan）检测异常增长簇，结合LandSat夜间灯光数据验证城市功能变化。可视化难点：（1）视觉遮挡：高密度体素叠加导致内部信息丢失，解决思路是采用体渲染透明度传递函数与交互式切片；（2）时间维度感知弱：用户难以关联颜色与时间，解决思路是引入动画时间轴与3D路径流线，同时提供2D热度图联动刷选。53.（封闭型）写出利用PostGIS完成“查询北京市六环内所有海拔高于100m的坡度大于15°区域”的完整SQL，并说明所需数据层。答案：数据层：dem（栅格，SRID=4547，单位米）、beijing_ring6（矢量，六环面，SRID=4547）。SQL：WITHslopeAS(SELECT(ST_DumpAsPolygons(ST_Slope(dem.rast,1,'32BF',111120,true))).geomASgeom,(ST_DumpAsPolygons(ST_Slope(dem.rast,1,'32BF',111120,true))).valASslope_degFROMdemJOINbeijing_ring6ONST_Intersects(dem.rast,beijing_ring6.geom)),high_areaAS(SELECT(ST_DumpAsPolygons(ST_Clip(dem.rast,beijing_ring6.geom))).geomASgeom,(ST_DumpAsPolygons(ST_Clip(dem.rast,beijing_ring6.geom))).valASelevFROMdemJOINbeijing_ring6ONST_Intersects(dem.rast,beijing_ring6.geom))SELECTST_Union(ST_Intersection(slope.geom,high_area.geom))ASresult_geomFROMslopeJOINhigh_areaONST_Intersects(slope.geom,high_area.geom)WHEREslope_deg>15ANDelev>100;54.（开放型）说明如何利用深度学习（DeepLabv3+）从高分辨率影像提取建筑物轮廓，并评价其与传统面向对象影像分析（OBIA）在精度与效率上的差异。答案：流程：①样本制作：采集0.2m分辨率航空正射影像，人工勾绘建筑物矢量，栅格化后按512×512切片，按8:1:1划分训练、验证、测试集；②模型训练：在PyTorch框架下采用ResNet101骨干，学习率0.001，批量16，迭代80k次，使用DiceLoss解决类别不平衡；③后处理：采用形态学闭运算填充空洞，利用OpenCV轮廓跟踪生成矢量，简化道格拉斯普克算法容差0.5m；④精度评价：IoU=0.89，边界F1=0.85，单张GPU（RTX3080）耗时0.8s/km²。与OBIA对比：精度：OBIA在多尺度分割参数80/100时IoU仅0.76，树冠阴影导致过分割；DeepLabv3+通过空洞卷积捕获多尺度上下文，误检率下降35%。效率：OBIA需人工试错分割与规则库，单人处理10km²需2d；深度学习端到端自动化，同等硬件10km²仅需15min，效率提升约64倍。55.（封闭型）解释克里金（Kriging）插值中的“块金效应（NuggetEffect）”与“基台值（Sill）”的物理意义，并给出半方差函数公式。答案：半方差函数：γ(h)=C₀+C[1−e^(−h/a)]，其中C₀为块金效应，C₀+C为基台值。物理意义：块金效应C₀：当采样间距h→0时半方差不为零的值，反映小于最小取样尺度的空间变异与测量误差；基台值C₀+C：半方差达到平稳时的上限，代表区域化变量总空间变异大小；比值C₀/(C₀+C)称块金系数，>0.75表明空间自相关弱。五、应用题（共3题，共70分）56.（综合类，25分）某市计划新建5座垃圾转运站，要求：①距现有学校≥500m；②距河流≥200m；③位于建设用地内；④服务半径2km覆盖人口最大化；⑤避开地质灾害高易发区。已有数据：school（点，CGCS2000）、river（线，CGCS2000）、landuse（栅格，分辨率为30m，类别：1建设用地/2绿地/3水体）、pop（栅格，人/像元）、hazard（矢量，高易发区）。请给出完整技术路线、关键工具及评价指标，并写出ArcGISPro的ModelBuilder流程图文字描述（含工具名称、参数设置、连接关系）。答案：技术路线：（1）数据预处理：统一投影至CGCS2000_3_Degree_GK_Zone_38，栅格重采样至30m；（2）约束因子生成： a.学校500m缓冲区（Erase）→禁止区A； b.河流200m缓冲区（Erase）→禁止区B； c.hazard面（Erase）→禁止区C； d.landuse提取建设用地（Reclassify=1）→允许区D；（3）候选区：允许区D与禁止区A/B/C取反后求交（Intersect），得到候选栅格掩膜；（4）人口聚合：用FocalStatistics（圆形2km，SUM）计算每个30m像元作为转运站时的服务人口；（5）选址优化：在候选区内按服务人口降序排序，用MaximumCoverage模型选前5个且两两间距≥1.5km（LocationAllocation）；（6）评价指标：覆盖率=服务人口/总人口，平均距离=Σ(人口×最近站距离)/服务人口，公平性=基尼系数。ModelBuilder描述：输入变量：school、river、landuse、pop、hazard。工具链：①Buffer（school，500m）→输出school_buf；②Buffer（river，200m）→输出river_buf；③Union（school_buf,river_buf,hazard）→输出forbid；④Reclassify（landuse，重映射1=1，其余NoData）→输出build；⑤Erase（build，forbid）→输出candidate；⑥FocalStatistics（pop，邻域圆2km，SUM）→输出pop_sum；⑦ExtractbyMask（pop_sum，candidate）→输出pop_candidate；⑧LocateRegions（输入pop_candidate，计数5，最小间距1500m）→输出sites；⑨连接顺序：①②→③，③+④→⑤，⑤+⑥→⑦，⑦→⑧；⑩添加派生输出sites作为参数，设置sites为模型参数P。57.（计算类，20分）已知某区域DEM栅格，分辨率为5m，高程基准为1985国家高程基准，坐标系为CGCS2000_3_Degree_GK_Zone_40。（1）请用GDALPython脚本计算该区域平均坡度（百分比）与坡向标准差；（2）若需将坡度重分类为≤3%、3–8%、8–15%、>15%四档，并统计各档面积（公顷），给出完整代码与结果保留两位小数。答案：```pythonimportgdal,numpyasnpfromosgeoimportgdalconstds=gdal.Open('dem.tif')gt=ds.GetGeoTransform()cols=ds.RasterXSize;rows=ds.RasterYSizeband=ds.GetRasterBand(1)nd=band.GetNoDataValue()dem=band.ReadAsArray().astype(np.float32)dem[dem==nd]=np.nan计算坡度百分比fromgdalimportDEMProcessinggdal.DEMProcessing('slope.tif',ds,'slope',scale=111120,slopeFormat='percent')slope_ds=gdal.Open('slope.tif')slope=slope_ds.ReadAsArray().astype(np.float32)slope[slope==slope_ds.GetRasterBand(1).GetNoDataValue()]=np.nanmean_slope=np.nanmean(slope)计算坡向标准差gdal.DEMProcessing('aspect.tif',ds,'aspect')aspect_ds=gdal.Open('aspect.tif')aspect=aspect_ds.ReadAsArray().astype(np.float32)aspect[aspect==aspect_ds.GetRasterBand(1).GetNoDataValue()]=np.nanaspect_std=np.nanstd(aspect)print(f"平均坡度={mean_slope:.2f}%")print(f"坡向标准差={aspect_std:.2f}°")重分类与面积统计bins=[0,3,8,15,100]classified=np.digitize(slope,bins)1pixel_area=55/10000公顷areas=[]foriinrange(4):count=np.sum(classified==i)areas.append(countpixel_area)print("各档面积(公顷):",[round(a,2)forainareas])```运行结果示例：平均坡度=12.45%坡向标准差=68.32°各档面积(公顷):[45.60,123.40,278.90,152.10]58.（分析类，25分）某省拟评估高铁线路（300km/h）对周边土地价格的时空影响。已知：（1）2015–2023年逐年住宅土地交易点（含成交单价、时间）；（2）高铁线路矢量，2020–12–30通车；（3）控制变量：到市中心距离、到最近商圈距离、人口密度、GDP栅格。请构建双重差分（DID）空间计量模型，说明实验组与对照组划分、空间权重矩阵构建、变量构造、估计方法与稳健性检验，并给出Python伪代码框架。答案：（1）实验组：通车后距高铁站点2km内交易点；对照组：距高铁>10km且位于同一条高铁走廊县域、预测趋势相似（采用PSM1:1匹配）。（2）空间权重矩阵：采用k=10最近邻，距离衰减函数wij=exp(−dij/5)，标准化行和为1。（3）变量构造： 被解释变量：log(成交单价)； 核心解释变量：Treat×Post，Treat=1为实验组，Post=1为2021–2023； 控制变量：X=[log(到市中心),log(到商圈),pop_density,log(GDP)]。（4）模型： Yit=α+β1Treati×Postt+γXit+ρWYit+μi+λt+εit 其中ρ为空间自回归系数，μi为个体固定效应，λt为年份固定效应。（5）估计：采用MaximumLikelihood估计，使用PySAL的spreg.Panel_FE_Lag。（6）稳健性： a.置换检验：随机虚构通车时间1000次； b.改变带宽1–3km； c.加入区县×年份