rpb/rpc文件说明与matlab读取

已于 2024-01-12 14:16:27 修改
阅读量1.4k 收藏 21
点赞数 16
分类专栏： Earth Observation CV 文章标签： matlab 开发语言
于 2024-01-12 14:13:59 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/dou3516/article/details/135551098
版权

什么是rpb/rpc文件？

        rpb文件是用来存储用于遥感数据几何校正的RPC（Rational Polynomial Coefficients ）模型的文件。类似的还有RPC文件，rpb与rpc文件只是格式不同，但包含的信息一致。其用于从图像坐标转换到地理坐标，L1级别的遥感数据是没有经过地理编码及几何校正的，因此可以通过数据包内的RPC文件进行数据几何校正（地理编码）。

        RPC参数是有理函多项式函数模型的参数，可以将遥感影像的图像坐标（Line，Sample）即（行号，列号）与其对应的大地坐标（Lat，Long， Height）即（纬度，经度，高程）用比值多项式一一对应起来。RPC的概念得到了很大的关注是因为RPC不仅可以将传感器成像的物理模型隐藏起来，同时可以通过RPC模型校正出几何精度更高的卫星数据。不同数据厂商使用不同的方法生成RPC文件，但目前大多数都采取的非线性最小二乘法（Non linear Least square）。在成像期间，要拟合的数据由物理传感器模型生成。        

        更详细的了解参考：

        卫星RPC有理多项式模型原理与正反算实现（python）-CSDN博客
        看懂卫星遥感数据RPC文件-CSDN博客

rpb/rpc文件格式

        rpb/rpc文件通常以文本文件或xml文件格式存储，对于国产高分数据多采用文本文件格式。

        rpb文件***.rpb格式示例：

satId = "XXX";
bandId = "XXX";
SpecId = "XXX";
BEGIN_GROUP = IMAGE
	errBias =   1.0;
	errRand =    0.0;
	lineOffset = 	+8.787000000000000e+03
	sampOffset = 	+8.787000000000000e+03
	latOffset = 	+5.197788782488139e+01
	longOffset = 	+1.267276217581132e+02
	heightOffset = 	-3.665824527965528e+01
	lineScale = 	+8.787500000000000e+03
	sampScale = 	+8.787500000000000e+03
	latScale = 	+1.000000000000000e+01
	longScale = 	+1.000000000000000e+01
	heightScale = 	+5.036755360952578e+03
	lineNumCoef = (
		-2.489085287819830e-04,
		-1.114244712474867e+01,
		-6.168967195179680e+01,
		+7.117708104229436e-05,
		-1.395349820950307e-01,
		+8.696093392117951e-03,
		-9.575461257281814e-04,
		+7.734382012664702e-01,
		+2.574838134723002e-02,
		+5.070051498763076e-06,
		-5.809725628331042e-03,
		+5.140801824748853e-03,
		+4.289598645908305e-03,
		+2.294752920112858e-03,
		-9.015843957850683e-04,
		+1.035394535926532e-02,
		+1.232218463012144e-02,
		+3.230339260864579e-02,
		+1.047590959343009e-03,
		+2.012768440253724e-08);
	lineDenCoef = (
		+1.000000000000000e+00,
		+4.180010565216353e-04,
		+1.324655197657377e-02,
		-3.795304837771012e-05,
		+2.545365596961364e-01,
		-3.899350380745237e-04,
		+5.589868166444578e-05,
		+1.188681625956534e-03,
		+6.863241299655094e-01,
		-2.002690706130058e-04,
		-2.072288733039329e-04,
		+6.133069607160820e-03,
		-1.449422536681941e-03,
		-2.011485529779256e-06,
		+3.500482182174969e-03,
		+4.198143375537498e-04,
		-2.285845110305829e-06,
		-2.172369339810132e-04,
		-4.144345265331140e-04,
		-1.836223330071302e-09);
	sampNumCoef = (
		+3.116290493890969e-03,
		+3.762565560709448e+01,
		-1.780287250778219e+01,
		+4.667388215545534e-02,
		-1.418392610803636e+00,
		-4.579691619531623e-02,
		+3.030408501996984e-02,
		-5.911326661297898e-01,
		-2.781904253558995e+00,
		-9.782108549302553e-05,
		+1.239834648584903e-03,
		+1.999566996539736e-01,
		+6.939464036194087e-02,
		-4.737705433260331e-02,
		-2.858026398839704e-01,
		-9.898030356627813e-03,
		+2.240107912228240e-02,
		+2.526433759914681e-03,
		+4.945279318567586e-03,
		-5.870434133733310e-05);
	sampDenCoef = (
		+1.000000000000000e+00,
		+3.227373471194835e-01,
		+3.944306307994457e-02,
		-9.509654358749953e-03,
		+1.207778166635094e+01,
		-3.232572945854702e-02,
		+1.355313365456889e-02,
		-1.277420643579939e+01,
		-2.837944635154219e+00,
		-1.255397137925662e-03,
		+8.474582054034840e-02,
		+8.587278629665410e-03,
		+1.984850060313677e-02,
		-6.307604213724449e-04,
		+2.407506405453826e-02,
		-1.315252027876808e-02,
		+2.805758450086238e-05,
		-9.639541191512180e-02,
		-1.882857688509176e-02,
		+1.022672251679239e-05);
END_GROUP = IMAGE
END;

        rpc文件***.rpc格式示例：

LINE_OFF :+8.787000000000000e+03  pixels
SAMP_OFF : +8.787000000000000e+03  pixels
LAT_OFF : +5.197788782488139e+01   degrees
LONG_OFF:  +1.267276217581132e+02   degrees
HEIGHT_OFF: -3.665824527965528e+01   meters
LINE_SCALE: +8.787500000000000e+03  pixels
SAMP_SCALE:  +8.787500000000000e+03  pixels
LAT_SCALE: +1.000000000000000e+01   degrees
LONG_SCALE: +1.000000000000000e+01   degrees
HEIGHT_SCALE :+5.036755360952578e+03   meters
LINE_NUM_COEFF_1:-2.489085287819830e-04
LINE_NUM_COEFF_2:-1.114244712474867e+01
LINE_NUM_COEFF_3:-6.168967195179680e+01
LINE_NUM_COEFF_4:+7.117708104229436e-05
LINE_NUM_COEFF_5:-1.395349820950307e-01
LINE_NUM_COEFF_6:+8.696093392117951e-03
LINE_NUM_COEFF_7:-9.575461257281814e-04
LINE_NUM_COEFF_8:+7.734382012664702e-01
LINE_NUM_COEFF_9:+2.574838134723002e-02
LINE_NUM_COEFF_10:+5.070051498763076e-06
LINE_NUM_COEFF_11:-5.809725628331042e-03
LINE_NUM_COEFF_12:+5.140801824748853e-03
LINE_NUM_COEFF_13:+4.289598645908305e-03
LINE_NUM_COEFF_14:+2.294752920112858e-03
LINE_NUM_COEFF_15:-9.015843957850683e-04
LINE_NUM_COEFF_16:+1.035394535926532e-02
LINE_NUM_COEFF_17:+1.232218463012144e-02
LINE_NUM_COEFF_18:+3.230339260864579e-02
LINE_NUM_COEFF_19:+1.047590959343009e-03
LINE_NUM_COEFF_20:+2.012768440253724e-08
LINE_DEN_COEFF_1:+1.000000000000000e+00
LINE_DEN_COEFF_2:+4.180010565216353e-04
LINE_DEN_COEFF_3:+1.324655197657377e-02
LINE_DEN_COEFF_4:-3.795304837771012e-05
LINE_DEN_COEFF_5:+2.545365596961364e-01
LINE_DEN_COEFF_6:-3.899350380745237e-04
LINE_DEN_COEFF_7:+5.589868166444578e-05
LINE_DEN_COEFF_8:+1.188681625956534e-03
LINE_DEN_COEFF_9:+6.863241299655094e-01
LINE_DEN_COEFF_10:-2.002690706130058e-04
LINE_DEN_COEFF_11:-2.072288733039329e-04
LINE_DEN_COEFF_12:+6.133069607160820e-03
LINE_DEN_COEFF_13:-1.449422536681941e-03
LINE_DEN_COEFF_14:-2.011485529779256e-06
LINE_DEN_COEFF_15:+3.500482182174969e-03
LINE_DEN_COEFF_16:+4.198143375537498e-04
LINE_DEN_COEFF_17:-2.285845110305829e-06
LINE_DEN_COEFF_18:-2.172369339810132e-04
LINE_DEN_COEFF_19:-4.144345265331140e-04
LINE_DEN_COEFF_20:-1.836223330071302e-09
SAMP_NUM_COEFF_1:+3.116290493890969e-03
SAMP_NUM_COEFF_2:+3.762565560709448e+01
SAMP_NUM_COEFF_3:-1.780287250778219e+01
SAMP_NUM_COEFF_4:+4.667388215545534e-02
SAMP_NUM_COEFF_5:-1.418392610803636e+00
SAMP_NUM_COEFF_6:-4.579691619531623e-02
SAMP_NUM_COEFF_7:+3.030408501996984e-02
SAMP_NUM_COEFF_8:-5.911326661297898e-01
SAMP_NUM_COEFF_9:-2.781904253558995e+00
SAMP_NUM_COEFF_10:-9.782108549302553e-05
SAMP_NUM_COEFF_11:+1.239834648584903e-03
SAMP_NUM_COEFF_12:+1.999566996539736e-01
SAMP_NUM_COEFF_13:+6.939464036194087e-02
SAMP_NUM_COEFF_14:-4.737705433260331e-02
SAMP_NUM_COEFF_15:-2.858026398839704e-01
SAMP_NUM_COEFF_16:-9.898030356627813e-03
SAMP_NUM_COEFF_17:+2.240107912228240e-02
SAMP_NUM_COEFF_18:+2.526433759914681e-03
SAMP_NUM_COEFF_19:+4.945279318567586e-03
SAMP_NUM_COEFF_20:-5.870434133733310e-05
SAMP_DEN_COEFF_1:+1.000000000000000e+00
SAMP_DEN_COEFF_2:+3.227373471194835e-01
SAMP_DEN_COEFF_3:+3.944306307994457e-02
SAMP_DEN_COEFF_4:-9.509654358749953e-03
SAMP_DEN_COEFF_5:+1.207778166635094e+01
SAMP_DEN_COEFF_6:-3.232572945854702e-02
SAMP_DEN_COEFF_7:+1.355313365456889e-02
SAMP_DEN_COEFF_8:-1.277420643579939e+01
SAMP_DEN_COEFF_9:-2.837944635154219e+00
SAMP_DEN_COEFF_10:-1.255397137925662e-03
SAMP_DEN_COEFF_11:+8.474582054034840e-02
SAMP_DEN_COEFF_12:+8.587278629665410e-03
SAMP_DEN_COEFF_13:+1.984850060313677e-02
SAMP_DEN_COEFF_14:-6.307604213724449e-04
SAMP_DEN_COEFF_15:+2.407506405453826e-02
SAMP_DEN_COEFF_16:-1.315252027876808e-02
SAMP_DEN_COEFF_17:+2.805758450086238e-05
SAMP_DEN_COEFF_18:-9.639541191512180e-02
SAMP_DEN_COEFF_19:-1.882857688509176e-02
SAMP_DEN_COEFF_20:+1.022672251679239e-05

rpb文件读取（matlab）

        定义读取函数 rpb_reader.m：

function RPC=rpb_reader(f_in_name)
   %author  dou3516
   % 2024-1-13
   % 参考：https://blog.csdn.net/weixin_43955546/article/details/123702554
   % 读取遥感RPB文件，输入经纬度，获取经纬度值对象的像素坐标。
   %f_in_name： 输入RPB文件名称

    % 1 读取RPB文件内容
    fid1=fopen(f_in_name,'r');
    
    read_one_line = fgetl(fid1);
    read_one_line = fgetl(fid1);
    read_one_line = fgetl(fid1);
    read_one_line = fgetl(fid1);
    read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',});
    errBias=str2double(temp(2));
    
    read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',});
    errRand=str2double(temp(2));
    
    read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',});
    RPC.line_off =str2double(temp(2));
    
    read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',});
    RPC.samp_off=str2double(temp(2));
    
    read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',});
    RPC.lat_off =str2double(temp(2));
    
    read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',','});
    RPC.lon_off=str2double(temp(2));
    
    read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',','});
    RPC.height_off=str2double(temp(2));
    
    
    read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',','});
    RPC.line_scale=str2double(temp(2));
    
    read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',','});
    RPC.samp_scale=str2double(temp(2));
    
    read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',','});
    RPC.lat_scale =str2double(temp(2));
    
    read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',','});
    RPC.lon_scale =str2double(temp(2));
    
    read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',','});
    RPC.height_scale =str2double(temp(2));
    
    read_one_line = fgetl(fid1); 
    for i=1:20
         read_one_line = fgetl(fid1);
        temp=strsplit(read_one_line,{';','=',')',','});
        RPC.line_num_coeff(i)=str2double(temp(1));
    end
    read_one_line = fgetl(fid1); 
    for i=1:20
         read_one_line = fgetl(fid1);
        temp=strsplit(read_one_line,{';','=',')',','});
        RPC.line_den_coeff(i)=str2double(temp(1));
    end
     read_one_line = fgetl(fid1); 
    for i=1:20
         read_one_line = fgetl(fid1);
        temp=strsplit(read_one_line,{';','=',')',','});
        RPC.samp_num_coeff(i)=str2double(temp(1));
    end
     read_one_line = fgetl(fid1); 
    for i=1:20
         read_one_line = fgetl(fid1);
        temp=strsplit(read_one_line,{';','=',')',','});
        RPC.samp_den_coeff(i)=str2double(temp(1));
    end
    fclose(fid1);

        示例读取代码：

fRPB = '***.rpb';
RPC = rpb_reader(fRPB);

        读取到的RPC信息示例（注意不同于上述rpb文件的实际信息，仅做示例）如下：

参考：

遥感 RPC， RPB文件相关信息-CSDN博客

there2belief
关注
  • 16
    点赞
  • 踩
  • 21
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
高分影像批处理第三回——RPC文件与几何校正
weixin_43762038的博客
03-16 1万+
前言 前面两章部分完成了图像的基本组织和大气校正，可是对于遥感影像地理相关数据，精确的地理坐标也是非常关键的。L1级数据，可惜，并没有帮我们处理地理坐标系统。无论是打开ENVI还是用gdal工具查看图像元数据，其坐标系统都是错误的或者说是未读入的。因此本章节来解决这个问题。
matlab开发-颜色谱创建不连续绘图颜色
08-24
matlab开发-颜色谱创建不连续绘图颜色。为绘图创建n个不同的三元素RPB向量
影像RPC和GCP校正
03-09
如果影像存在rpb文件,使用此程序可以自行完成校正, 无校正文件需要改造程序构造相关输入即可, 基于GDAL开发,亦可用于 跨平台
遥感 RPC， RPB文件相关信息
花藻の博客
03-24 5090
RPC相关信息1. rpc和rpb的区别2. c++(visual_studio)安装gdal（1）安装配置：（2）一个版本不合适导致的报错。3. gdal的使用（1）（2）python+gdal rpc批量正射校正4.gdal rpc正反算（1）从行列号推算其经纬度。5.介绍rpc参数很详细的一篇文章
RPC校正
stone_tigerLI的博客
12-24 3898
RPC是什么？ RPC（rational polynomial coefficients），有理多项式系数。RPC是传感器严格几何模型的拟合形式。 是根据卫星平台载荷测量的平台运行轨迹参数、姿态参数、传感器安装参数及传感器内部几何参数构建的像-地关系几何模型[1]^{[1]}[1]。 RPC计算过程 P = (Latitude - LAT_OFF) / LAT_SCALE L = (Longitude - LONG_OFF) / LONG_SCALE H =
使用GDAL自带函数读取RPC文件的方法
maybcsdn的博客
03-21 6689
RPC为有理多项式系数（Rational Polynomial Coefficient），在遥感影像预处理中用于影像像素坐标（行列号）与对应同名点物方坐标（经度，纬度，高程）之间的相互计算。RPC文件，顾名思义，就是包含有理多项式系数的文本文件。RPC文件通常有两种格式：一种为图1所示，称为RPC格式的文件，另一种就是图2所示的，称为RPB格式的文件。 图1 RPC文件图2 RPB文件 ...
给geotiff影像创建一个RPB文件
花藻の博客
05-29 363
给有地理坐标系的影像，创建一个不依赖地形的rpb文件。
解密计算机网络中的RPB与Pruning技术：构建高效通信网络的利器
最新发布
m0_72410588的博客
08-05 684
RPB技术是一种广播传输方式，其核心思想是通过反向路径选择，将数据包从源主机发送到目标主机的所有节点上。具体而言，RPB技术采用了反向路径选择的策略，首先将数据包从源主机发送到相邻节点，然后通过反向路径选择方法，在每个节点处选择一个出边接口将数据包继续转发，直至到达目标主机。Pruning技术是一种网络优化方法，通过删除网络中不必要的连接和节点，减少网络的规模和复杂度。其核心思想是基于网络的结构和特性，去除冗余的连接和节点，以提升网络的性能和效率。
卫星影像的RPC参数
ch206265的博客
07-09 1万+
什么是RPC数据 RPC（Rational Polynomial Coefficient）有理多项式系数 共线条件方程作为摄影测量与遥感的核心,也是传感器模型的理论基础。通过共线方程建立严格的几何成像模型需要...
看懂卫星遥感数据RPC文件
BellaLYW的博客
10-12 5942
持续更新有关使用遥感图像RPC文件相关内容包括不限于原理、说明、使用、代码
卫片影像RPC参数
fb_help的专栏
01-18 2916
卫片影像RPC参数 什么是RPC模型 RPC（RFC）参数是使用有理函数模型来模拟地面物方到卫片像方的映射关系。 为什么使用RPC的模型 传统的中心投影遵循透视映射的几何模型。可以用内外方位元素来描述这种中心投影的成像模型，但卫片影像一般采用线阵CCD推扫的成像方式，或者任何其它的非中心投影的成像方法。不同于中心投影，CCD推扫的成像方式在一条线上是满足中心投影的，但线与线之间却不满足中心投影。对于这种成像方式，有严格的线阵几何模型来描述，也可以用通用的有理函数模型（RFM，RPC模型就是有理函数模型的系数
基于RPC方式的文件传输应用开发
06-25
设计并实现一个简单的 RPC分布应用。实现一个基于RPC 方式的文件传输应用，启动服务器和客户端验证结果。
rpb.rar_engineering_rpb
09-24
industrial engineering
BioSwob.Ecoaro.gaO1rpb
03-18
BioSwob.Ecoaro.gaO1rpb
colorSpectrum：创建 N 个不同的绘图颜色：创建 N 个不同的三元素 RPB 向量用于绘图-matlab开发
05-30
避免在 matlab 图中很难看到的颜色，例如黄色和青色。 三个简单的方向： 获取要绘制的 N 条曲线或线：数据。 类型：C = colorSpectrum(N); 情节：图，等一下，for i = 1:N,plot(Data(i,:),'Color',C(i,:)),end,...
国家2000坐标系CGCS2000（EPSG4490和EPSG4479）参数
热门推荐
豆芽菜
04-27 3万+
EPSG：4490 China Geodetic Coordinate System 2000 (CGCS2000)经纬度坐标系 Unit:degree (supplier to define representation) Geodetic CRS:China Geodetic Coordinate System 2000 Datum:China 2000 Ellipsoid:...
使用Google Earth Engine（一）：提取下载MODIS、Landsat点数据
豆芽菜
09-07 3万+
现在可以使用Google Earth Engine 快捷的提取MODIS、Landsat等产品的点（或区域）数据到自己的Google Drive下载。PS:GEE可以做全球尺度长时间序列的数据分析出图等，将来必是地球科学领域的大杀器。这里下载数据不是仅运行下面代码就行，还需要注意两点： 1）那么如何找到数据集呢？在搜索栏输入产品名，便可快速定位到产品，点击产品名右侧的import可...
Landsat 7 ETM+影像条带修复ENVI补丁
豆芽菜
05-10 2万+
本文介绍两个常用的Landsat 7ETM+影像条带修复ENVI补丁的下载和简单使用。ENVI补丁sav文件的使用方法是放置到程序安装目录的save_add文件夹下，比如在ENVI4.7下是D:\Program Files\ITT\IDL71\products\envi47\save_add。 1）landsat_gapfill.sav（比较好用） 下载地址：http://pan.baidu.
ENVI栅格裸数据生成shp和kml矢量文件
豆芽菜
06-19 2万+
ENVI栅格裸数据生成shp和kml矢量文件。生成kml文件需要用到ArcMap，但是ArcMap不能直接加载ENVI裸数据，需要在ENVI里将栅格数据转为shp或另存为ArcViewRaster。 前者生成kml文件的流程如下： 1)        ENVI打开栅格裸数据； 2)        选择Vector -> Raster to Vector，获得evf格式的矢量数据； 3)  
RPB文件和RPC文件的不同
03-27
RPB文件和RPC文件是两种不同的文件格式，主要区别如下： 1. RPB文件：RPB是Raster Processing Bitmap的缩写，是一种基于栅格数据的图像文件格式。RPB文件中存储的是像元值，每个像元代表图像中的一个点。RPB文件适用于存储具有一定分辨率的图像数据，如卫星图像、遥感图像等。 2. RPC文件：RPC是Rational Polynomial Coefficients的缩写，是一种基于多项式模型的地理参考文件格式。RPC文件中存储了多项式系数和影像的位置信息，可以将像素坐标转换为地理坐标。RPC文件适用于将图像数据与地理坐标进行关联，实现地理信息系统（GIS）中的空间分析。 综上所述，RPB文件和RPC文件是两种不同的文件格式，RPB文件主要用于存储图像数据，RPC文件主要用于存储地理参考信息。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值