2025年地理信息科学考试试卷及答案

2025年地理信息科学考试试卷及答案一、单项选择题（每题1分，共20分。每题只有一个正确答案，请将正确选项字母填在括号内）1.在Web墨卡托投影（EPSG:3857）中，赤道处1个像素对应地面实际距离约为：A.0.1m  B.0.3m  C.0.6m  D.1.2m【答案】B【解析】Web墨卡托在第0级（整幅世界地图256×256像素）时，赤道周长40075017m，每像素≈40075017/256≈156543m；每放大一级分辨率减半，第19级分辨率≈156543/2^19≈0.3m。2.利用Sentinel2影像计算NDVI时，若B4中心波长为665nm，B8中心波长为842nm，则NDVI表达式为：A.(B8B4)/(B8+B4)  B.(B4B8)/(B4+B8)  C.(B8B2)/(B8+B2)  D.(B8B4)/B8【答案】A【解析】NDVI=(NIRRed)/(NIR+Red)，B8为近红外，B4为红光。3.在PostGIS中，执行ST_Area(geom)返回值的单位由：A.数据库编码决定  B.坐标参考单位决定  C.总是度²  D.总是m²【答案】B【解析】ST_Area直接计算几何坐标单位下的面积，若geom为地理坐标则单位为度²，投影坐标则为投影单位。4.下列关于Delaunay三角网的描述正确的是：A.任意四点共圆  B.最小角最大化  C.边长总和最小  D.所有三角形面积相等【答案】B【解析】Delaunay三角剖分满足“空圆”准则，同时最大化最小角，避免狭长三角形。5.在移动GIS中，影响GPS单点定位精度的主要误差源不包括：A.电离层延迟  B.多路径效应  C.卫星钟差  D.地图符号大小【答案】D【解析】地图符号大小属于可视化范畴，不影响测距观测值。6.利用无人机航测生成DSM时，以下哪种滤波算法最适合城市区域：A.渐进形态学滤波  B.基于坡度的滤波  C.基于高程阈值  D.中值滤波【答案】A【解析】渐进形态学可在保留建筑物边缘的同时逐步去除地面点，适合复杂城市表面。7.在QGIS中，将矢量图层坐标系从WGS84转换到CGCS2000，应使用：A.图层属性→样式  B.图层属性→源→设置CRS  C.处理工具箱→重投影  D.保存为→选择编码【答案】C【解析】重投影需使用“处理工具箱→矢量通用→重投影”，仅设置CRS不会改变几何坐标。8.下列哪种空间索引结构最适合高维点云范围查询：A.Rtree  B.Quadtree  C.KDtree  D.Btree【答案】C【解析】KDtree对高维k近邻查询效率高于Rtree，Rtree更适合二维外包框。9.在Python中使用rasterio读取GeoTIFF，获取仿射变换参数的接口是：A.dataset.crs  B.dataset.transform  C.dataset.nodata  D.dataset.dtypes【答案】B【解析】transform返回Affine对象，包含六参数仿射变换。10.根据《地理信息数据安全分类分级指南》，涉密测绘成果的核心识别要素是：A.比例尺大于1:5万  B.包含等高距小于5m  C.精度优于±1m  D.包含国家秘密坐标系参数【答案】D【解析】坐标系参数、控制点成果属于国家绝密级，直接决定数据安全等级。11.在SuperMap中，将线对象平滑为B样条曲线，需使用：A.拓扑处理  B.栅格分析  C.矢量编辑→光滑  D.空间查询【答案】C【解析】光滑操作可将折线转换为B样条或贝塞尔曲线。12.利用激光雷达回波次数信息提取电力线，主要依据：A.回波强度  B.多次回波间隔  C.首次回波高程突变  D.末次回波密度【答案】C【解析】电力线直径小，激光首次回波高程明显高于周围植被末次回波。13.在GDAL中，将栅格重采样方法设为“lanczos”对应的枚举值是：A.GRA_NearestNeighbour  B.GRA_Bilinear  C.GRA_Cubic  D.GRA_Lanczos【答案】D【解析】GDALWarpOptions的eResampleAlg字段直接支持GRA_Lanczos。14.下列关于H3空间索引的说法错误的是：A.由Uber开源  B.采用六边形网格  C.层级使用GeoHash编码  D.支持球面层次剖分【答案】C【解析】H3使用六边形+五边形剖分，编码为64位整数，非GeoHash。15.在遥感影像分类中，使用光谱角制图（SAM）衡量像元相似性，其取值范围：A.[0°,45°]  B.[0°,90°]  C.[0°,180°]  D.[0°,360°]【答案】B【解析】SAM计算两向量夹角，最小0°最大90°，夹角越小越相似。16.在PostgreSQL中，执行“SELECTST_3DDistance(a.geom,b.geom)FROMtba,tbb”时，需确保：A.安装postgis_sfcgal  B.安装postgis_topology  C.安装postgis_raster  D.安装postgis_tiger【答案】A【解析】ST_3DDistance需要SFCGAL库支持三维几何。17.利用InSAR技术监测地面沉降，若两次成像间隔ΔT=12天，相位变化Δφ=2π，则形变量为：A.λ/2  B.λ/4  C.λ  D.2λ【答案】A【解析】InSAR一个完整2π相位周期对应雷达波长λ的一半形变量。18.在Web地图切片中，GoogleXYZ规范的第0级瓦片数量为：A.1  B.2  C.4  D.8【答案】A【解析】第0级仅1张256×256像素瓦片覆盖全球。19.下列哪项不是OGC标准服务：A.WMS  B.WFS  C.WMTS  D.SHP【答案】D【解析】SHP是EsriShapefile格式，非服务接口。20.在深度学习中，用于遥感影像建筑物提取的UNet网络，其跳跃连接主要解决：A.梯度消失  B.下采样信息丢失  C.过拟合  D.计算量过大【答案】B【解析】跳跃连接将编码器高分辨率特征与解码器融合，恢复空间细节。二、多项选择题（每题2分，共20分。每题有两个或两个以上正确答案，多选少选均不得分）21.下列哪些属于空间数据质量元素：A.位置精度  B.属性精度  C.逻辑一致性  D.时间精度  E.地图配色【答案】ABCD【解析】ISO19113定义空间数据质量元素包括位置、属性、逻辑、时间、完整性和语义精度。22.在ArcGISPro中，以下哪些工具可用于矢量数据综合：A.简化建筑物  B.聚合面  C.平滑线  D.消除面  E.波段合成【答案】ABCD【解析】波段合成属于栅格处理，其余为综合工具。23.关于激光雷达点云分类，以下哪些特征可用于机器学习：A.高程  B.回波强度  C.全波形能量  D.点密度  E.扫描角【答案】ABCDE【解析】五类特征均能反映地物差异，提高分类精度。24.下列哪些投影属于等角投影：A.高斯克吕格  B.Web墨卡托  C.UTM  D.阿尔伯斯等积  E.Lambert等角圆锥【答案】ABCE【解析】阿尔伯斯为等积投影，其余保持角度不变。25.在Python中，使用geopandas.sjoin实现空间连接，可支持的谓词包括：A.intersects  B.within  C.contains  D.nearest  E.crosses【答案】ABCE【解析】nearest需使用sjoin_nearest，sjoin原生不支持。26.下列哪些现象会导致遥感影像出现几何畸变：A.地球曲率  B.地形起伏  C.传感器倾斜  D.大气散射  E.扫描镜速度变化【答案】ABCE【解析】大气散射影响辐射而非几何。27.在PostGIS中，以下哪些函数返回几何类型为线串：A.ST_MakeLine  B.ST_Intersection(两线相交)  C.ST_Boundary(面)  D.ST_Centroid  E.ST_Union(多点)【答案】ABC【解析】ST_Centroid返回点，ST_Union(多点)返回多点。28.下列哪些属于时空数据模型：A.时空立方体  B.事件驱动模型  C.三域模型  D.面向对象模型  E.场模型【答案】ABC【解析】面向对象与场模型描述静态空间，非时空专用。29.在GDAL中，构建VRT虚拟栅格可实现：A.多幅影像无缝镶嵌  B.波段选取  C.重投影  D.金字塔构建  E.矢量裁剪【答案】ABCD【解析】矢量裁剪需使用warp或mask。30.下列哪些技术可提高GPS高程精度：A.差分GPS  B.星基增强  C.精密单点定位  D.气压计辅助  E.增加观测卫星数【答案】ABCDE【解析】五项均可削弱或改正高程误差。三、判断题（每题1分，共10分。正确打“√”，错误打“×”）31.WGS84与CGCS2000在厘米级精度上可以视为一致。  【答案】×【解析】两框架历元差异导致厘米级偏差，需七参数转换。32.在QGIS中，字段计算器使用“$area”返回的单位始终为平方米。  【答案】×【解析】$area使用图层CRS单位，地理坐标下为度²。33.KDtree构建时间复杂度为O(nlogn)。  【答案】√【解析】标准KDtree递归划分，平均O(nlogn)。34.激光雷达波长越短，穿透植被能力越强。  【答案】×【解析】短波长（如532nm）易被叶片散射，近红外（1064nm）穿透更强。35.GeoPackage可同时存储矢量、栅格、属性表与元数据。  【答案】√【解析】GeoPackage基于SQLite，支持多表与多种数据类型。36.在遥感影像中，热红外波段可用于反演地表温度。  【答案】√【解析】热红外（10–12μm）对温度敏感，常用劈窗算法反演。37.空间自相关Moran’sI取值范围为[1,1]，0表示完全随机。  【答案】√【解析】Moran’sI定义满足该范围，0对应空间随机。38.使用无人机航测时，增加旁向重叠度可提高DEM精度但对DOM无影响。  【答案】×【解析】旁向重叠改善空三匹配，同样减少DOM接缝误差。39.在PostGIS中，执行ST_Transform(geom,4326)前必须确保geom已定义SRID。  【答案】√【解析】未定义SRID时ST_Transform无法知道源坐标系。40.卷积神经网络中，空洞卷积可扩大感受野而不增加参数。  【答案】√【解析】空洞卷积通过间隔采样实现，参数量与标准卷积相同。四、填空题（每空2分，共20分）41.在GDAL中，使用命令行将GeoTIFF转换为COG（云优化GeoTIFF）需添加参数__________。【答案】ofCOG或coTILED=YEScoCOMPRESS=LZWcoCOPY_SRC_OVERVIEWS=YES【解析】COG需内部瓦片、概览与压缩。42.高斯克吕格投影3°分带，中央经线117°E对应的带号是__________。【答案】39【解析】带号=中央经线/3，117/3=39。43.在Python中，geopandas读取Shapefile后，查看坐标参考系属性名为__________。【答案】.crs【解析】GeoDataFrame的crs属性存储pyproj.CRS对象。44.根据《测绘法》，基础测绘成果实行__________制度。【答案】汇交【解析】测绘成果需依法汇交至测绘地理信息主管部门。45.在激光雷达点云分类中，ISPRS标准将“桥梁”归为__________类别。【答案】人工地物【解析】ISPRS测试集将桥梁、建筑物等统一为“artifact”或“manmade”。46.利用Sentinel1SAR影像提取水体，常用__________极化组合增强水体与陆地对比度。【答案】VH/VV【解析】VH对水体表面粗糙度敏感，比值可抑制植被噪声。47.在PostGIS中，创建空间索引的SQL语句为CREATEINDEXidx_nameONtableUSING________(geom)。【答案】GIST【解析】PostGIS默认GiST索引支持二维几何。48.在Web地图中，瓦片行列号采用__________方向递增。【答案】谷歌XYZ（或ZXY）【解析】谷歌规范原点在西经180°、北纬85°，行号自上而下递增。49.在遥感影像中，蓝光波段中心波长约为__________nm。【答案】480【解析】Landsat8B2中心波长480nm，常用于水体穿透。50.在深度学习中，DiceLoss常用于解决__________问题。【答案】类别不平衡【解析】DiceLoss对前景背景像素数量不敏感，适合小目标提取。五、简答题（每题8分，共40分）51.简述利用InSAR技术监测城市地面沉降的主要流程，并指出两个关键误差源及改正方法。【答案】流程：1.选取合适SAR影像对（时间基线短、空间基线小）；2.影像精配准（亚像素级）；3.干涉图生成与去平地效应；4.相位滤波与解缠；5.地理编码与形变转换（φ→Δr）；6.时间序列分析（SBAS或PSInSAR）获取沉降速率。关键误差源：1.大气延迟——利用GNSS水汽数据或外部气象模型（如ERA5）进行大气改正；2.轨道误差——采用精密轨道（POD）或引入GCP做多项式拟合改正。52.说明如何利用多光谱影像与激光雷达点云进行城市建筑物提取，给出特征组合与分类器选择理由。【答案】步骤：1.影像预处理：辐射定标、大气校正、正射校正；2.点云滤波：渐进形态学去除地面，保留非地面点；3.特征级融合：将点云高程、强度、归一化高程（nDSM）与影像NDVI、纹理（GLCM）重采样至同一分辨率；4.特征组合：光谱+高度+强度+纹理+形状（面积、长宽比）；5.分类器：采用随机森林，原因：可处理高维特征、对噪声稳健、提供特征重要性；6.后处理：形态学闭运算去除小洞，规则化建筑物边缘。优势：高度特征有效区分建筑物与树木，光谱特征抑制裸地，组合精度可达95%以上。53.解释“空间数据引擎（SDE）”在大型企业GIS中的核心作用，并比较ArcSDE与PostGIS的架构差异。【答案】核心作用：1.统一存储：将多源矢量、栅格、元数据存入RDBMS，实现版本管理、事务、备份；2.并发控制：采用乐观/悲观锁、版本表（A、D表）支持多用户同时编辑；3.空间索引：利用Rtree或格网索引加速查询；4.SQL接口：提供ST_Geometry类型，使空间查询成为标准SQL。架构差异：ArcSDE为中间件模式，位于ArcGIS与Oracle/SQLServer之间，需安装SDE服务进程，proprietary；PostGIS为开源插件，直接嵌入PostgreSQL，无中间层，使用GiST索引，遵循OGC标准，更轻量、可扩展。54.给出利用Python实现“矢量面要素自动旋转最小外接矩形（MBR）”的代码思路，并说明如何计算长轴方向。【答案】思路：1.使用geopandas读取面要素；2.对每条geometry调用shapely.minimum_rotated_rectangle，返回Polygon；3.提取外接矩形四个顶点，按顺序排列；4.计算长轴向量：顶点0→顶点1，利用math.atan2(dy,dx)获取方位角θ（弧度）；5.将θ转为度，若θ<0则+180，得到0–180°范围长轴方向；6.将角度、面积、长宽比写入新字段，输出Shapefile。关键代码：```pythonimportgeopandasasgpd,mathgdf=gpd.read_file('input.shp')defget_angle(geom):rect=geom.minimum_rotated_rectanglecoords=list(rect.exterior.coords)dx=coords[1][0]coords[0][0]dy=coords[1][1]coords[0][1]angle=math.degrees(math.atan2(dy,dx))returnangleifangle>=0elseangle+180gdf['angle']=gdf.geometry.apply(get_angle)gdf.to_file('output.shp')```55.描述“数字孪生城市”中GIS与BIM融合的技术路径，并列举两项典型应用场景。【答案】技术路径：1.数据标准：采用IFC4GeoRef将BIM模型地理参考，导出CityGML或3DTiles；2.坐标转换：通过七参数将BIM局部坐标系转为CGCS2000；3.语义映射：建立IFC实体到CityGMLLOD3/LOD4的对应表，保留材质、属性；4.集成平台：基于Cesium或UnrealEngine加载3DTiles，叠加GIS服务（WMTS、WFS）；5.实时IoT：通过MQTT接入传感器，驱动BIM构件状态（如阀门开关）与GIS空间分析联动。应用场景：1.市政管网爆管模拟：GIS定位爆管点，BIM精确定位阀门，联动关闭影响范围；2.应急疏散：GIS计算建筑出入口容量，BIM提供室内路径，实时优化疏散路线。六、综合应用题（共40分）56.（20分）某市计划利用2025年6月获取的0.2m分辨率无人机影像与同期激光雷达点云，完成主城区1:500地形图更新。测区面积120km²，平均点云密度50pts/m²，影像倾斜角<5°，控制点采用CORS网络RTK测量，平面精度±2cm，高程±3cm。（1）设计航测技术方案，包括航线布设、像控点密度、点云分类指标；（8分）（2）给出DEM与DOM生产流程，指出关键质量检查项；（6分）（3）估算数据总量（原始影像+点云+成果），并说明存储与备份策略。（6分）【答案】（1）航线布设：航向重叠80%，旁向重叠70%，满足1:500规范；分区飞行，单航线长≤4km，转弯半径>500m，避开高楼电磁干扰；像控点密度：每km²15点，共1800点，采用“周边+中间”布点，拐角、高低处加密；点云分类指标：地面点高程中误差≤5cm，建筑物边界平面误差≤0.3m，植被剔除率≥98%。（2）生产流程：空三：使用AgisoftMetashape，导入POS、控制点，平差报告重投影误差<0.5像素；点云：分类→地面→建筑物→植被→其他，人工质检；DEM：地面点构TIN，0.5m格网，双线性内插，检查高程精度（RTK检核点高程中误差≤8cm）；DOM：真正射纠正，消除建筑物投影差，拼接线选阴影少区域，色彩均衡，分辨率0.2m；质量检查：空三报告、DEM高程检核、DOM拼接线目视、图幅接边差≤1mm（图上）。（3）数据量估算：原始影像：单张像元数8000×