通过控制点生成RPC文件

1. 控制点生成RPC的软件

函数来源： https://www.orfeo-toolbox.org/CookBook-6.4/Applications/app_GenerateRPCSensorModel.html

需要下载LINUX版本的软件，才能使用这个功能

2. 一些便捷使用的函数

（1）无坐标的影像和参考的有地理坐标的影像生成控制点文件

这两张影像通过ENVI的map-regietration 手动选择一些同名点，然后对这个同名点进行整理

fclose all
clear all
clc


filename = 'G:\test\selected_43_points.pts';
[x1,y1,x2,y2]=textread(filename,'%f%f%f%f', 'headerlines ' , 5 );
x1=x1-6573;

landsat_name='H:\points_test\ls_landsat.tif';

[landsat_data,R_landsat_data]=geotiffread(landsat_name);


[lat,lon]=pix2latlon(R_landsat_data,y2,x2);


fid = fopen('G:\selected_43_points_new.txt', 'w');
for i=1:size(lat,1)
    fprintf(fid, '%s\n', [num2str(x1(i),'%.6f'),32,num2str(y1(i),'%.6f'),32,num2str(lon(i),'%.6f'),32,num2str(lat(i),'%.6f')]);
end

fclose(fid);

（2）生成的rpc文件，进行格式转换

% 读取生成的rpc文件，读取需要的信息
fclose all
clear all
clc



f_in_name='H:\ZY3_data\corona_test\output1.geom';
fid1=fopen(f_in_name,'r');

for i=1:37
        read_one_line = fgetl(fid1);  
end


read_one_line = fgetl(fid1);
temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
heightOffset=char(temp(2));


read_one_line = fgetl(fid1);
temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
heightScale =char(temp(2));

read_one_line = fgetl(fid1);

read_one_line = fgetl(fid1);
temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
latOffset =char(temp(2));

read_one_line = fgetl(fid1);
temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
latScale  =char(temp(2));

for i=1:20
     read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
    lineDenCoef{i}=temp(3);
end    

for i=1:20
     read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
    lineNumCoef{i}=temp(3);
end    

read_one_line = fgetl(fid1);
temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
lineOffset =char(temp(2));    


read_one_line = fgetl(fid1);
temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
linescale=char(temp(2));


read_one_line = fgetl(fid1);
read_one_line = fgetl(fid1);



read_one_line = fgetl(fid1);
temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
longOffset=char(temp(2));

read_one_line = fgetl(fid1);
temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
longScale =char(temp(2));

for i=1:15
read_one_line = fgetl(fid1);

end

for i=1:20
     read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
    sampDenCoef{i}=temp(4);
end    

 for i=1:20
     read_one_line = fgetl(fid1);
    temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
    sampNumCoef{i}=temp(4);
end   



read_one_line = fgetl(fid1);
temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
sampOffset=char(temp(3));


read_one_line = fgetl(fid1);
temp=strsplit(read_one_line,{';','=',')',',','p','d','m',':'});
samplescale=char(temp(3));

fclose(fid1);








% 2 将rpc文件写为rpb文件
f_out_name='H:\ZY3_data\corona_test\test.rpb';
fid2=fopen(f_out_name,'w');
fprintf(fid2,'satId = "XXX";');
fprintf(fid2,'\n');
fprintf(fid2,'bandId = "XXX";');
fprintf(fid2,'\n');
fprintf(fid2,'SpecId = "XXX";');
fprintf(fid2,'\n');
fprintf(fid2,'BEGIN_GROUP = IMAGE');
fprintf(fid2,'\n');
fprintf(fid2,'	errBias =   1.0;');
fprintf(fid2,'\n');
fprintf(fid2,'	errRand =    0.0;');
fprintf(fid2,'\n');

fprintf(fid2,['	lineOffset = ',strtrim(char(lineOffset))]);
fprintf(fid2,'\n');
fprintf(fid2,['	sampOffset = ',strtrim(char(sampOffset ))]);
fprintf(fid2,'\n');
fprintf(fid2,['	latOffset =  ',strtrim(char(latOffset))]);
fprintf(fid2,'\n');
fprintf(fid2,['	longOffset =  ',strtrim(char(longOffset ))]);
fprintf(fid2,'\n');
fprintf(fid2,['	heightOffset =  ',strtrim(char(heightOffset ))]);
fprintf(fid2,'\n');
fprintf(fid2,['	lineScale =  ',strtrim(char(linescale))]);
fprintf(fid2,'\n');
fprintf(fid2,['	sampScale =  ',strtrim(char(samplescale))]);
fprintf(fid2,'\n');
fprintf(fid2,['	latScale = ',strtrim(char(latScale))]);
fprintf(fid2,'\n');
fprintf(fid2,['	longScale =  ',strtrim(char(longScale ))]);
fprintf(fid2,'\n');
fprintf(fid2,['	heightScale = ',strtrim(char(heightScale ))]);
fprintf(fid2,'\n');

fprintf(fid2,'	lineNumCoef = (');
fprintf(fid2,'\n');
for i=1:19
    fprintf(fid2,['			',strtrim(char(lineNumCoef{i})),',']);
    fprintf(fid2,'\n');
end
fprintf(fid2,['			',strtrim(char(lineNumCoef{20})),')',';']);
fprintf(fid2,'\n');

fprintf(fid2,'	lineDenCoef = (');
fprintf(fid2,'\n');
for i=1:19
    fprintf(fid2,['			',strtrim(char(lineDenCoef{i})),',']);
    fprintf(fid2,'\n');
end
fprintf(fid2,['			',strtrim(char(lineDenCoef{20})),')',';']);
fprintf(fid2,'\n');

fprintf(fid2,'	sampNumCoef = (');
fprintf(fid2,'\n');
for i=1:19
    fprintf(fid2,['			',strtrim(char(sampNumCoef{i})),',']);
    fprintf(fid2,'\n');
end
fprintf(fid2,['			',strtrim(char(sampNumCoef{20})),')',';']);
fprintf(fid2,'\n');


fprintf(fid2,'	sampDenCoef = (');
fprintf(fid2,'\n');
for i=1:19
    fprintf(fid2,['			',strtrim(char(sampDenCoef{i})),',']);
    fprintf(fid2,'\n');
end
fprintf(fid2,['			',strtrim(char(sampDenCoef{20})),')',';']);
fprintf(fid2,'\n');

fprintf(fid2,'END_GROUP = IMAGE');
fprintf(fid2,'\n');
fprintf(fid2,'END;');
fclose(fid2);