给高分系列卫星影像添加地理坐标
问题陈述:两张高分影像。影像A:具有地理坐标信息XML;影像B:不具备地理坐标信息。将影像A的地理信息赋给影像B。
应用的场景:图像mask标签无地理坐标,有原始图像,对mask赋值。
addpos-v1.py
import sys
from threading import currentThread
import numpy as np
# from Py6S import *
from osgeo import gdal,gdal_array
def parse_rpc_file(rpc_file):
# rpc_file:.rpc文件的绝对路径
# rpc_dict:符号RPC域下的16个关键字的字典
# 参考网址:http://geotiff.maptools.org/rpc_prop.html;
# https://www.osgeo.cn/gdal/development/rfc/rfc22_rpc.html
rpc_dict = {
}
with open(rpc_file) as f:
text = f.read()
# .rpc文件中的RPC关键词
words = ['errBias', 'errRand', 'lineOffset', 'sampOffset', 'latOffset',
'longOffset', 'heightOffset', 'lineScale', 'sampScale', 'latScale',
'longScale', 'heightScale', 'lineNumCoef', 'lineDenCoef','sampNumCoef', 'sampDenCoef',]
# GDAL库对应的RPC关键词
keys = ['ERR_BIAS', 'ERR_RAND', 'LINE_OFF', 'SAMP_OFF', 'LAT_OFF', 'LONG_OFF',
'HEIGHT_OFF', 'LINE_SCALE', 'SAMP_SCALE', 'LAT_SCALE',
'LONG_SCALE', 'HEIGHT_SCALE', 'LINE_NUM_COEFF', 'LINE_DEN_COEFF',
'SAMP_NUM_COEFF', 'SAMP_DEN_COEFF']
for old, new in zip(words, keys):
text = text