地理国情普查培训dom生产说明

地理国情普查

正射影像生产技术规定说明国家基础信息中心2013.12.6主要内容需求资料情况试验结论主要内容指标项目对DOM需求地理国情普查是地理国情监测的基础，普查数据的精度对今后开展地理国情监测也十分重要。DOM用于地理国情普查底图，同时用于更新国家DOM产品数据库。几何精度和影像质量对地理国情普查成果的质量至关重要。技术规定参考标准兼顾DOM产品的标准规范（CH/T9009.3-2010）参考了五万建库和更新工程技术规定，但原来分辨率稍低、精度单一、数据量小。调研了土地二调底图生产技术指标（严格按规范，多种精度）。第一次高分覆盖，精度从严控制，保证普查精度，也便于以后利用。既有整景纠正成果，又有分幅的产品数据，整景纠正用于普查解译，分幅的按产品要求制作影像数据源分析多源影像数据

以优于1米分辨率影像为主（大部分0.5米）；较高的定轨定姿参数；新型传感器影像较多（WV、Pleiades等）

控制资料分析（1）1:50000更新区域外业像控点、航空影像空三加密成果、1:1万量测控制点西部测图区域

外业像控点、航空影像和SPOT5HRS立体加密成果五万更新/西部测图五万航片生产区域控制资料分析（2）1:1万地形图目前已覆盖166400幅，约占全国46%（全国约373749幅）,DRG本身以及覆盖区域的像控成果。1:1万地形图覆盖区域技术试验概况中心和三个直属局（陕西、龙江、四川）参加。试验数据包括：WorldView1、WorldView2、Quickbird、Geoeye、Pleiades、资源三号试验区域

考虑不同地形类别，同时有高精度检测数据

陕西：西安14景wv，2景QB，1景ZY3

龙江：齐齐哈尔9景wv，1景ZY3，呼玛县16景wv，不同侧视角（最大33.5）

四川：北川安县3景ZY3、2景WV、简阳1景Geoeye、攀枝花西部、德阳、洪雅县10景WV主要实验内容高分辨率卫星影像在不同地形类别、不同像控条件下，正射纠正精度情况DEM精度对正射纠正精度的影响（1万、5万）弱连接(交会)区域网平差的可行性及意义已有成果用于像控的利用效率现有商用软件的支持程度实验内容及结论（1）优于1米高分卫片不同控制点正射精度

平地丘陵地：1-2个控制点可以达到1:1万DOM精度，3-4个点最好，再多意义不大。ZY34个以上控制点也可以达到1:1万精度。

山地高山地：

1-2个控制点可以达到1:2.5万DOM精度。ZY34个以上控制点也可以达到1:2.5万精度。但侧视角较大（大于25度）时精度较差，需要4个点以上才能达到2.5万要求。

实验内容及结论（2）DEM分辨率和精度对高分卫片正射精度影像

1万DEM比5万DEM有明显精度提高DEM精度对侧视角较大影像精度影响较大弱连接条件下区域网平差工程可实施性可以明显减少控制点数量，但对DEM精度依赖较大使用前需要评估DEM精度情况侧视角大的需要单独纠正实验内容及结论（3）新数据源和新技术商用软件的支持程度

Erdas、PCI、PixelGrid、DPGrid、像素工厂等软件均支持新型卫片,ErdasLPS、PixelGrid、像素工厂支持区域网平差大数据量对现有生产流程的影响

工作量比以前多倍增加，最好有并行处理软硬件系统。已有成果用于像控选点的可行性

航片区域可用性大，卫片区域像控点少，总体上只有40%的点可以利用。和已有DOM配准：只适合优于1米DOM（备选方式）技术规定总体精度要求总体上控制在1:2.5万以上精度1:1万覆盖区域和航空摄影区域原则上按1:1万精度要求生产。其他区域原则上按1:2.5万精度要求生产。几何精度指标少数困难地区适当放宽，另行报批。接边差从严要求（分辨率高）不同影像数据源、不同控制数据源以及不同生产批次之间的正射影像接边限差放宽至倍地形类别1:10000成图精度影像平面中误差1:25000成图精度影像平面中误差平地5.0米12.5米丘陵地5.0米12.5米山地7.5米18.75米高山地7.5米18.75米地形类别1:10000成图精度影像图幅接边限差1:25000成图精度影像图幅接边限差平地5.0米12.5米丘陵地5.0米12.5米山地7.5米18.75米高山地7.5米18.75米影像源选择如果生产区域内同时存在多源影像数据，正射影像生产优先选用空间分辨率和光谱分辨率更高、时相更靠近生长季、现势性更新的影像，如果同时具有相当的卫星影像和航摄影像，优先采用卫星影像。影像源选择（西部地区）优先利用高分辨率卫星影像进行生产，下列情况可选择资源三号影像，或报项目部协调满足要求的高分辨率卫星影像。高分辨率卫星影像云量超过20%，影响地物判读。山地、高山地影像侧视角大于20度，平地、丘陵地影像侧视角大于25度，正射纠正精度不能满足“技术规定”中3.8节要求。高分辨率卫星影像的时相在冬季，影响地理国情普查地物解译。像控和DEM控制资料现有像控资料和空三加密成果外业像控特殊地区可以用现有DOM配准、周边影像做控制等方案，航片采用联合平差等方式。卫片按航片加密点要求，适当调整。DEM

优先选用1:1万DEM，不低于1:5万DEM加强资料收集收集作业区域各种像控资料，除高精度像控和加密成果外，包括国界界碑、CORS站点、原西部测图工程区域的HRG整景正射影像资料等，作为正射影像纠正的补充像控数据源。外业调绘核查过程中采集外业控制点西部地区要求每1:5万幅图不少于2个。其他地区可根据需要采集。西部地区平面精度按照表8中1:25000精度执行，高程精度按照表9中1:25000精度执行。控制点精度地形类别1:100001:25000平面误差高程误差平面误差高程误差平地1.00.12.50.4丘陵地1.00.252.50.5山地1.00.52.50.6高山地1.00.52.51.2控制点精度地形类别1:100001:25000平面误差高程误差平面误差高程误差平地2.51.05.01.0丘陵地2.51.05.01.5山地3.752.07.52.0高山地3.753.07.53.5内外业一体化生产流程从2次外业变为1次外业像控之前可以同步开展内业作业定向流程航片采用空三加密方式进行定向，高程精度可适当放宽；保证平面精度条件下，像控布设可放宽要求可采用新老资料联合平差，但老资料精度需满足要求卫片RPC模型多光谱和全色波段套合误差在一个像素内；允许弱连接区域网平差同轨同时相的，可采用先拼接再纠正方式；弱连接(交会)区域网平差精度景号

地形类别控制点数

检测点误差类型201103260354122011102904020120120529033918201204250256322012052903391720120425025629高山高山山地山地平地平地6个4个7个4个19个10个无控最大16.17016.87314.36817.97011.63110.778最小4.9574.7935.0476.1824.0816.364中误差11.25912.8839.04011.1527.2167.9636点最大13.83113.5857.96710.8914.3344.502最小2.2512.9753.7842.8102.4002.620中误差8.8677.2864.5835.8362.7093.23314点最大13.05913.6375.0207.6445.2124.008最小2.0303.5031.1881.2341.0572.080中误差8.0137.5523.5714.0772.9053.17926点最大13.24313.4054.5614.9773.21713.855最小1.1793.8772.3411.0641.5701.932中误差7.9747.5002.8772.9552.0042.967弱连接(交会)区域网平差精度控制点数量类型点数最大误差中误差4个控制点定向点42.91.95检查点45.93.768个控制点定向点82.91.81检查点0控制点点数中误差（米）最大误差（米）备注41.5582.64621.9002.940（对角分布）22.1593.014（一侧分布）12.4343.97104.1056.098明显减少像控点

提高接边精度在高山区1:5万DEM情况下，精度不稳定地形类别：平地丘陵地不同像控、地形、DEM纠正精度地形类别景号影像类型侧视角（度）检查点中误差（5万25米间距）（米）中误差（5万10米间距）（米）无控1点4点无控1点4点平地025629QB9.3297.76473.2223.068.4763.3222.795033917WV118.5357.0063.8943.20210.5923.9332.992山地033918WV118.3209.3584.6424.15610.8664.8024.006025632QB9.5299.1384.1324.0719.8694.3043.94高山地035233WV16.2108.4217.9263.1378.5595.5252.871035232WV16.5307.0825.8464.7676.9255.0094.388035230WV17.1336.6373.9253.6116.4113.7323.071035330WV111.4328.9136.2195.0987.9886.1364.832035329WV113.9279.8036.7985.1359.4996.3304.772035647WV123.43613.9069.2697.92613.7549.4107.644035648WV123.92511.4888.3475.93011.8148.9285.037035941WV123.93212.7187.9266.37913.7857.5666.280035413WV126.41410.997.8175.76311.6997.1085.644035412WV126.93112.1028.0745.10512.8328.6515.264035411WV127.53012.4019.8467.83712.6049.6097.621040201WV229.13028.40818.78911.13128.27414.73310.587部分实验精度地类景号类型侧视角检查点中误差（5万5米间距）（米）中误差（1万5米间距）（米）无控

1点

4点无控

1点

4点平地025629QB9.3298.4883.3122.9688.5883.2982.645033917WV118.53510.6813.3953.13411.0172.8862.206山地033918WV118.32010.8064.1414.43110.4853.8823.097025632QB9.5299.7074.1623.9319.6683.6363.128高山地035233WV16.2108.8373.7583.058035232WV16.5307.0144.9604.805035230WV17.1336.2333.6103.478035330WV111.4327.3855.4194.7586.7463.2771.951035329WV113.9279.4255.3604.8168.7463.9202.725035647WV123.43613.5599.3716.53010.6144.5173.956035648WV123.92512.1278.8705.4819.8243.7153.055035941WV123.93213.6028.2956.559035413WV126.41412.2047.8805.58112.8444.9643.550035412WV126.93113.0428.8285.66110.3695.0463.635035411WV127.53012.6688.2227.66111.5723.4182.159040201WV229.13027.50414.83710.14425.48211.9564.203DEM对大侧视角影像纠正影响侧视角29.1度地形类别

高山地1:5万DEM1:1万DEM控制点数控制点中误差检查点数检查点中误差控制点数控制点中误差检查点数检查点中误差11911.541192.9621.371810.7121.37182.6640.62169.7940.62161.4360.78148.5060.78141.44侧视角11.4度地形类别

高山地1:5万DEM1:1万DEM控制点数控制点中误差检查点数检查点中误差控制点数控制点中误差检查点数检查点中误差1265.591262.0520.74255.5620.74252.1640.84234.6640.84231.1360.76214.6960.76211.10误差单位：米DEM对大侧视角影像纠正影响侧视角26.9度地形类别

高山地1:5万DEM1:1万DEM控制点数控制点中误差检查点数检查点中误差控制点数控制点中误差检查点数检查点中误差1187.951182.9820.54177.8420.54172.6740.45156.1440.45152.1360.41136.5460.41132.24DEM精度大侧视角影像纠正精度影响明显（采用高精度DEM（1:1万）或者需要4个以上控制点）侧视角大于25度约占总景数15%-20%。误差单位：米卫片像控指标卫星影像类型纠正模型控制点个数侧视角<=25度侧视角>25度WorldView-1/2RPC模型≥2≥4GeoEye-1RPC模型≥2≥4QUICKBIRDRPC模型≥2≥4IKONOSRPC模型≥2≥4PleiadesRPC模型≥2≥4ZY-3RPC模型≥3≥6TH-1RPC模型≥6≥10法国SPOT5/6严格轨道模型≥3≥6印度P5RPC模型≥6≥10日本ALOSRPC模型≥6≥10每个加密分区不超过30景卫星影像，控制点数量不少于卫星影像景数的一半，最少不低于6个控制。相邻加密分区接边区域应该分布不少于2个共用控制点，并利用共用控制点进行接边检查，网间公共点平面较差不超过P15表10规定的指标要求。多光谱和全色之间纠正后配准误差在1个像素之内（全色）整景和分幅分别要求分幅：参考产品规范，Tif+TFW格式，大小不超过2G，匀光，统一三种分辨率，色彩模式等；统一外扩100个像素（和CH/T9009.3-2010不一样）分景：保留原影像指标（分辨率、光谱等），用于地理国情普查解译，tif或img格式均可；坐标均无带号，文件不设投影信息分幅要求优于（包含）1米，应按照1:25000分幅裁切，则一个1:50000图幅内应存在4个数据文件、4个TFW影像信息文件、1个元数据文件和1个投影信息文件。分辨率为2米，应按照1:50000分幅裁切，则一个1:50000图幅内应存在1个数据文件、1个TFW影像信息文件和1个元数据文件。投影信息文件（*.Xml）<?xmlversion="1.0"?><PAMDataset><SRS>PROJCS["CGCS_2000_GK_CM_105E",GEOGCS["GCS_CGCS_2000",DATUM["D_CGCS_2000",SPHEROID["CGCS_2000",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Gauss_Kruger"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",105.0],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]</SRS></PAMDataset>影像数据源统一规定便于生产统一，把适合本项目的数据源具体规定，包括时相，以及分辨率、色彩和时相选取原则；如果增加再补充规定不同影像规定了不同的像控要求影像分辨率和像素位分幅数据

地面分辨率采用0.5米、1.0米和2.0米三种规格。色彩模式为黑白（8位）、RGB彩色（24位）二种（P3表1和P4表2对卫星影像和航空影像进行统一要求）；分景数据

分辨率和像素位原则上和原始影像地面分辨率保持一致，为了统一对分辨率进行归纳（P16表12和P17表13）。分幅DOM分辨率对照表传感器类型原始地面分辨率正射影像分辨率WorldView-1/20.5米0.5米GeoEye-10.5米0.5米QUICKBIRD0.61米0.5米IKONOS1.0米1.0米Pleiades0.7米0.5米ZY-32.1米2.0米TH-12.0米2.0米SPOT5（HRG）2.5米2.0米SPOT6（1.5米）1.5米1.0米印度P52.5米2.0米日本ALOS2.5米2.0米按产品要求，统一三个规格（0.5米、1米和2米）个别利用已有DOM，分辨率较高（0.2米），可保留原分辨率，按1:1万分幅存储整景DOM分辨率对照表传感器类型原始地面分辨率正射影像分辨率WorldView-1/20.5米0.5米GeoEye-10.5米0.5米QUICKBIRD0.61米0.6米IKONOS1.0米1.0米Pleiades0.7米0.7米ZY-32.1米2.0米TH-12.0米2.0米SPOT5（HRG）2.5米2.5米SPOT6（1.5米）1.5米1.5米印度P52.5米2.5米日本ALOS2.5米2.5米只是对分辨率取整，便于统一。尽量保持原分辨率。同一图幅包含不同数据源包含的影像源分辨率高分辨率影像所占比例分幅成果分辨率0.5米1.0米超过30%0.5米低于30%1.0米1.0米2.0米超过30%1.0米低于30%2.0米0.5米2.0米超过90%0.5米70%-90%1.0米低于30%2.0米不同分辨率组合：按照1:25000标准分幅（7.5´×5´）为基本存储单元。如果图幅内多光谱影像所占比例超过图幅的10%，则整个图幅的按RGB彩色（24位）模式存储，如果所占比例低于图幅的10%，则按灰度（8位）模式存储。全色和多光谱组合：按照1:25000标准分幅（7.5´×5´）为基本存储单元。其他指标融合一般只对同一卫星遥感影像的多光谱数据和全色波段数据进行融合，建议用Pansharp。

多光谱所有波段均融合，分幅的选三个波段，接近真彩色。影像增强

匀光、拉伸等，只针对分幅，整景不增强处理。其他数学基础、命名、裁切影像扭曲变形的处理整景正射影像不需对扭曲变形的影像进行处理。分幅正射影像制作时尽量避免使用扭曲变形的影像，当影像扭曲变形影响地物的判读和采集时，需对该部分影像重新进行处理消除变形。由于地形陡峭引起的正射影像拉伸（拉花）不需处理。元数据参考CH/T1007-2001规范和1:5万更新规定，对一些过时的或不必要的进行精简DOM基本情况。数据源情况生产过程信息数据分发信息每个图幅一个文件，保存为XML格式（网上提供了录入工具程序）控制点影像采集与整理控制点影像数据包括控制点点位索引影像和控制点信息文件。索引影像为对应的卫星影像全色波段遥感影像（纠正前原始影像）上和原始航空影像上采集控制点信息文件为记录控制点坐标等信息的描述性文件包括像控点和加密点像控点格式[Controlpointinfo]Geographic:Lat=31°59′58.044″//地理坐标，纬度Geographic:Lon=113°21′35.551″//地理坐标，经度Geodetic:X=3544228.510（注）//大地投影坐标，南北方向Geodetic:Y=723021.310（注）//大地投影坐标，东西方向，Geodetic:Z=1120.310（注）//高程PointType=ControlPoint（TiePoint）//像控点（加密点）ImgResolution=0.50//影像地面分辨率，单位：米