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

2025年地理信息系统工程师认证考试试题及答案一、单项选择题（每题1分，共30分。每题只有一个最符合题意的选项，错选、多选均不得分）1.在Web墨卡托投影（EPSG:3857）中，赤道处1像素对应地面实际距离约为：A.0.1mB.0.3mC.0.6mD.1.2m答案：B解析：Web墨卡托在赤道处1像素对应0.28m左右，四舍五入后最接近0.3m。2.下列哪项不是GeoPackage标准规定的必选表：A.gpkg_contentsB.gpkg_geometry_columnsC.gpkg_spatial_ref_sysD.gpkg_metadata答案：D解析：gpkg_metadata为可选表，其余三项为强制存在。3.在PostGIS中，执行ST_Union与ST_Collect的主要区别是：A.前者返回多部件几何，后者返回单部件B.前者会溶解相邻边界，后者仅简单收集C.前者仅支持2D，后者支持3DD.两者无区别，可互换答案：B解析：ST_Union做几何合并并溶解公共边，ST_Collect仅返回GeometryCollection。4.当使用ArcGISPro构建空间权重矩阵时，若选择“K最近邻”且K=8，则每个要素的邻居数：A.一定等于8B.可能小于8，若要素位于边界C.可能大于8，若存在距离相等D.取决于投影坐标系答案：A解析：K最近邻强制固定邻居数，边界要素亦会跨区寻找第8邻。5.在GDAL中，下列命令行可将GeoTIFF的NoData值从9999改为0并生成新文件的是：A.gdal_translatea_nodata0in.tifout.tifB.gdalwarpdstnodata0in.tifout.tifC.gdal_edit.pya_nodata0in.tifD.gdal_calc.pyAin.tifoutfile=out.tifcalc="A(A!=9999)"答案：A解析：a_nodata参数直接重写元数据，无需重采样。6.下列关于H3空间索引的说法正确的是：A.分辨率0的六边形面积约为1,100km²B.相邻分辨率面积比精确为1:7C.索引字符串可唯一表示全球任意点D.支持跨分辨率无损合并答案：A解析：H3分辨率0的六边形平均面积约1,100km²，面积比约1:7.07，非精确1:7。7.在QGIS中，使用“几何生成器”符号化时，表达式buffer($geometry,100)的单位：A.总是米B.总是图层单位C.总是地图单位D.取决于项目CRS答案：B解析：几何生成器使用图层坐标系单位，与项目CRS无关。8.下列哪种插值方法对“孤岛”效应最敏感：A.IDWB.普通克里金C.自然邻域D.TIN答案：A解析：IDW无法外推超出采样范围的值，易产生“牛眼”或孤岛。9.在GeoServer中，若需发布PostGIS视图，必须：A.视图含唯一整型列B.视图含几何列C.视图含主键D.视图必须物化答案：C解析：GeoServer需可唯一标识要素，主键或唯一索引即可，不限定物化。10.下列关于SLD的命名空间前缀，错误的是：A.sld:NamedLayerB.ogc:FilterC.se:StrokeD.gml:PointSymbolizer答案：D解析：PointSymbolizer属于sld/se命名空间，非gml。11.在PostGIS中，执行ST_Subdivide(geom,256)的作用是：A.将几何对象切分为顶点数不超过256的片段B.将几何对象切分为面积不超过256m²的片段C.将几何对象切分为256份等面积片段D.将几何对象降维为256个点答案：A解析：ST_Subdivide限制顶点数，提升查询效率。12.下列哪项不是OGCAPIFeatures核心操作：A./collectionsB./collections/{cid}/itemsC./collections/{cid}/coverageD./api答案：C解析：coverage属于CoverageAPI，非Features核心。13.在GDAL中，创建COG（CloudOptimizedGeoTIFF）必须包含的块大小是：A.256×256B.512×512C.任意，但需内嵌概览D.与原始块一致答案：C解析：COG无强制块大小，但需内嵌概览且按瓦块组织。14.下列关于空间数据库RTree索引的描述，正确的是：A.仅支持二维外包矩形B.可加速空间关系运算C.不支持并发写入D.必须定期重建答案：B解析：RTree加速空间查询，支持高并发，无需频繁重建。15.在ArcGISPro中，使用“时空立方体”工具时，时间步长单位不支持：A.秒B.小时C.周D.财年答案：D解析：财年非固定长度，工具未内置。16.下列关于EPSG:4326与EPSG:4490的区别，正确的是：A.两者椭球参数不同B.4490采用CGCS2000，4326采用WGS84C.4490为地心坐标系D.两者单位不同答案：B解析：4490基于CGCS2000椭球，4326基于WGS84，单位均为度。17.在QGIS中，使用“处理模型”时，若需将矢量输出设为临时图层，应：A.输出路径留空B.使用memory:前缀C.勾选“临时输出”D.输出至TEMP环境变量答案：C解析：3.20+版本提供“临时输出”复选框。18.下列关于GeoJSON的坐标精度，正确的是：A.规范要求保留7位小数B.1位小数约对应1kmC.6位小数约对应0.1mD.精度与CRS有关答案：C解析：6位小数≈0.1m，规范未强制位数。19.在PostGIS中，ST_Transform与ST_SetSRID的区别是：A.前者重投影，后者仅改元数据B.前者仅改元数据，后者重投影C.两者均重投影D.两者均不改坐标值答案：A解析：ST_Transform计算新坐标，ST_SetSRID仅改SRID标识。20.下列关于三维Tiles（3DTiles）的说法，错误的是：A.基于glTFB.支持LODC.必须使用四叉树D.支持点云答案：C解析：3DTiles支持四叉树、八叉树及自定义空间索引。21.在GDAL中，使用VRT格式的主要优势是：A.压缩比高B.无需写实际像素C.支持矢栅混合D.支持事务答案：B解析：VRT为XML元文件，可动态拼接，无需复制像素。22.下列关于空间SQL的索引使用，最可能触发全表扫描的是：A.WHEREST_Contains(geom,ST_Point(120,30))B.WHEREST_DWithin(geom,ST_Point(120,30)::geography,100)C.WHEREST_Intersects(geom,'POLYGON(...)')D.WHEREST_Area(geom)>1000答案：D解析：ST_Area计算值无索引，无法使用空间索引。23.在ArcGISPro中，使用“深度学习”工具需安装的框架是：A.PyTorchB.TensorFlowC.fast.aiD.ArcGISAPIforPython答案：B解析：Esri官方支持TensorFlow后端。24.下列关于WMTS的TileMatrix限制，正确的是：A.同一TileMatrixSet内所有层级像素尺寸相同B.顶层矩阵必须覆盖全球C.必须采用256×256块D.支持自定义比例分母答案：D解析：WMTS允许自定义TileMatrix，块大小与像素尺寸可变。25.在QGIS中，使用“几何检查”插件无法发现：A.重复节点B.自相交C.多部件空洞D.属性缺失答案：D解析：插件仅检查几何错误，属性缺失由字段计算器发现。26.下列关于空间数据版本控制（Geogig）的描述，正确的是：A.基于集中式锁机制B.使用类似Git的三树结构C.仅支持ShapefileD.无法合并冲突答案：B解析：Geogig采用Git语义，支持分支合并。27.在PostGIS中，执行ST_ClusterDBSCAN窗口函数时，参数eps=1000,minpoints=10，单位：A.总是米B.取决于几何SRID单位C.总是度D.可指定为参数答案：B解析：eps单位与几何坐标系一致，地理类型需用度或转geography。28.下列关于激光雷达LAS文件的说法，正确的是：A.LAS1.4支持波形数据B.LAZ为无损压缩C.点记录格式0包含RGBD.版本号存储于文件头第4字节答案：B解析：LAZ基于LASzip无损压缩，1.4支持波形，格式0无RGB。29.在GDAL中，使用gdalbuildvrtseparate的作用是：A.创建金字塔B.按波段分离堆叠C.创建虚拟堆叠多波段D.重投影答案：C解析：separate将多个单波段文件虚拟堆叠为多波段VRT。30.下列关于空间数据脱敏的最佳实践，错误的是：A.随机偏移50–500mB.使用K匿名C.删除所有坐标小数位D.聚合到网格答案：C解析：删除小数位仍可通过拓扑恢复，非安全脱敏。二、多项选择题（每题2分，共20分。每题有两个或两个以上正确答案，多选、少选、错选均不得分）31.下列哪些格式支持存储要素时间戳（datetime）字段：A.ShapefileB.GeoPackageC.GeoJSOND.FileGeodatabase答案：B、C、D解析：Shapefile的.dbf仅支持Date类型，不含时分秒。32.在PostGIS中，可用于计算地理（椭球）距离的函数包括：A.ST_DistanceB.ST_DistanceSphereC.ST_DistanceSpheroidD.ST_Length(geography)答案：B、C、D解析：ST_Distance默认笛卡尔，需geography或球体函数。33.下列哪些操作会改变GeoTIFF的像素值：A.gdalwarpt_srsB.gdal_translateoutsize50%50%C.gdal_edit.pya_srsD.gdalwarptr1010答案：A、B、D解析：a_srs仅改元数据，不采样像素。34.关于OGCAPIFeaturesPart2，支持的查询参数有：A.datetimeB.bboxC.limitD.cqltext答案：A、B、C解析：cqltext为Part3扩展。35.下列哪些软件可直接打开CloudOptimizedGeoTIFF无需下载整文件：A.QGISB.ArcGISProC.GoogleEarthProD.GDAL答案：A、B、D解析：GoogleEarthPro不支持HTTPRange读取COG。36.在三维Tiles中，哪些文件属于核心元数据：A.tileset.jsonB.subtree.jsonC.glTFD.batchTable.json答案：A、B解析：glTF为几何载体，batchTable为属性，非核心元数据。37.下列哪些方法可用于Shapefile字符编码转换：A.ogr2ogrlcoENCODING=UTF8B.iconv修改.dbfC.QGIS“另存为”指定编码D.ArcMap“导出”更改代码页答案：A、C、D解析：直接修改.dbf需同步.shp，iconv易破坏文件。38.下列哪些QGIS表达式可返回多边形边界框中心：A.centroid($geometry)B.point_on_surface($geometry)C.bounds($geometry)D.x(centroid($geometry))答案：A、B解析：bounds返回矩形几何，非中心点。39.下列哪些属于空间大数据分布式计算框架：A.GeoMesaB.SedonaC.RasterFramesD.GRASSGIS答案：A、B、C解析：GRASS为单机，无原生分布式。40.下列哪些做法可提升WMS服务出图速度：A.预生成瓦片B.启用MapCacheC.简化图层样式D.降低源数据分辨率答案：A、B、C、D解析：四项均可减少实时计算量。三、案例分析题（共30分）41.某市自然资源局计划将1980–2020年土地利用矢量数据（共300GBShapefile）迁移至开源空间数据库，要求支持历史版本、快速空间查询及多用户并发编辑。请回答：（1）选择何种数据库与存储方案？（3分）答案：PostGIS+GeoPackage混合方案；历史数据用GeoPackage分年归档，活跃数据用PostGIS分区表。解析：PostGIS支持并发与索引，GeoPackage便于离线归档与版本快照。（2）说明数据迁移的技术路线与关键工具。（5分）答案：①使用ogr2ogr批量转Shapefile至PostGIS，lcoPRECISION=no加速；②建立分区表按年份+行政区哈希分区；③创建GIST索引并执行VACUUMANALYZE；④使用pgAgent定时备份至GeoPackage；⑤通过Geogig建立分支，支持版本回退。解析：分区减少索引膨胀，Geogig提供Gitlike历史。（3）如何验证数据一致性？（2分）答案：采用PostGIS的ST_SymDifference与原Shapefile对比，面积差<0.001%视为一致；属性行数与sum(dtype)校验。解析：几何与属性双重校验，避免丢要素。42.某在线地图需发布全球10m分辨率DEM，数据量1.2TB，要求浏览器端秒级加载任意视角三维地形。现有单机GeoServer，带宽1Gbps。请设计技术方案：（1）选择何种切片与格式？（2分）答案：采用CesiumTerrain格式（quantizedmesh）+TMS瓦块结构，层级0–14。解析：quantizedmesh支持Cesium客户端高效LOD与法线。（2）说明后端处理流程。（4分）答案：①使用gdalbuildvrt构建全球VRT；②gdal_translate生成COG，块大小512×512，LZW压缩；③使用cesiumterrainbuilder切片为quantizedmesh；④上传至S3对象存储，启用CDN；⑤配置Nginx反向代理，支持HTTP/2与br压缩。解析：COG作为中间产物，保证重投影与压缩；CDN降低延迟。（3）如何评估服务性能？（2分）答案：使用JMeter模拟1000并发，统计P95<500ms；CesiumInspector查看Tile加载时序，空缓存下首屏<3s。解析：并发与首屏双重指标，覆盖网络与渲染瓶颈。四、实操编程题（共20分）43.使用Python完成以下任务：给定某市POI点数据（EPSG:4490），需提取出距离地铁线路（EPSG:4527）500m范围内的餐饮类POI，并计算最近地铁站的站名与距离（米），结果输出为GeoPackage，坐标系保持CGCS2000。请补全代码并给出关键行注释。（10分）```pythonimportgeopandasasgpdfromshapely.opsimportnearest_points读取数据poi=gpd.read_file('poi.geojson').to_crs(epsg=4527)转投影坐标subway=gpd.read_file('subway.shp').to_crs(epsg=4527)筛选餐饮food=poi[poi['type']=='餐饮']500m缓冲buf=subway.buffer(500)joined=gpd.sjoin(food,subway[['name','geometry']],how='inner',predicate='within')计算最近站与距离defget_nearest_station(row,subway):point=row.geometrydistances=subway.distance(point)nearest=distances.idxmin()returnsubway.loc[nearest,'name'],distances.min()joined[['nearest_sta','dist_m']]=joined.apply(lambdarow:get_nearest_station(row,subway),axis=1,result_type='expand')转回4490并输出joined.to_crs(epsg=4490).to_file('result.gpkg',driver='GPKG')```答案要点：①投影转换需在缓冲前完成；②sjoin后仍须计算到最近站，避免多站匹配；③结果转回4490保持权威坐标系。44.使用PostGISSQL完成：将某表parcel中所有面积大于1000m²的多边形按“东西”方向一分为二，若结果仍大于1000m²则继续递归分割，直至满足条件。要求返回原fid、分割层级、新几何。（10分）```sqlWITHRECURSIVEsplitAS(SELECTfid,