高分3号介绍及PIE使用


edited by nrzheng,2022.2.12

参考链接:123456

GF3介绍

高分三号卫星是我国首颗分辨率达到1米的雷达卫星,高分三号卫星数据是C频段多极化合成孔径雷达数据(SAR数据),发射于2016年8月,卫星具有多个拍摄模式。

1. 分辨率

高分三号卫星提供的SAR图像可以提供1米至500米分辨率、10公里至650公里幅宽的微波遥感数据,用于服务于海洋环境监测与权益维护、灾害监测与评估、水利设施监测与水资源评价管理、气象研究及其他多个领域。

2. C波段的好处

根据功能和使命的不同,地球微波遥感探测卫星可以划分为L、 S、C、X等频段,频率依次由低到高。频率越低,穿透力越强,越高则穿透力越差。打个比方,假如用这四种频段看一棵树,X频段只能看到树梢,L、 S频段只能看到树根,C频段则介于两者之间,既能看到地表,也能看到树干。因此,在对海洋、气象、减灾等环境和目标进行探测时,C频段的优势更为突出。

(补充)

在雷达行业,雷达工作频率划分成为若干波段,由低到高的顺序依次是高频(HF),甚高频(VHF)、超高频(UHF,也称作P)、L波段、S波段、C波段、X波段,Ku波段、K波段、Ka波段、U波段、V波段和W波段。
在这里插入图片描述

标志字符全程/名称来源频率范围中心频率中心频率对应波长
HF高频:high frequency3~30MHz————
VHF甚高频:very high frequency30~300MHz————
UHF/P特高频:ultra high frequency300~1000MHz————
L长波:long wave1~2GHz1.5GHz20cm
S短波:short wave2~4GHz3GHz10cm
CS和X的折中:compromise4~8GHz6GHz5cm
X用于火控雷达,X来自于瞄准镜的“准心”8~12GHz10GHz3cm
Kukurz-under12~18GHz15GHz2cm
Kkurz,德语“短”(容易被水蒸气吸收,不能再雨雾天气使用,所以有Ka、Ku)18~27GHz————
Kakurz-above27~40GHz35GHz8mm
U40~60GHz50GHz6mm
V60~80GHz70GHz4mm
W80~110GHz100GHz3mm

各波段雷达功能:

  • 波长短,分辨率高,穿透性差,容易被吸收

  • 波长长,分辨率低,穿透性强

  • 搜索雷达,工作频段(VHF、UHF、L)

  • 搜索跟踪雷达,工作波段(L、S、C)

  • 火控、成像雷达,工作波段(C、X、Ku)

  • 弹载雷达,工作频段(X、Ku、K、Ka、V、U、W)

3. 极化方式

高分三号多极化包含单极化、双极化及全极化

4. 成像模式

共包括6种条带成像模式、1种滑块聚束成像模式、3种扫描成像模式、1种波成像模式、2种扩展入射角成像模式。

高分三号卫星是迄今为止成像模式(12种)最多的SAR卫星,它不仅涵盖了传统的条带、扫描成像模式,而且可在聚束、条带、扫描、波浪、全球观测、高低入射角等多种成像模式下实现自由切换,既可以探地,又可以观海,达到“一星多用”的效果。

(我们数据的成像模式是UFS/SL/QPS1)

序号成像模式入射角 (°)视数A*E分辨率 (m)成像带宽 (km)极化方式
中文名英文名标称方位向距离向标称范围
1滑块聚束Spotslight (SL)20~501*111.0~1.50.9~2.510*1010*10可选单极化
2超精细条带Ultra-fine stripmap (UFS)20~501*1332.5~53030可选单极化
3精细条带1Fine stripmap (FSI)19~501*1554~65050可选双极化
4精细条带2Wide fine stripmap (FSII)19~501*210108~1210095~110可选双极化
5标准条带Standard stripmap (SS)17~503*2252515~3013095~150可选双极化
6窄幅扫描Narrow ScanSAR (NSC)17~501*65050~6030~60300300可选双极化
7宽幅扫描Wide ScanSAR (WSC)17~501*810010050~110500500可选双极化
8全球观测模式Global observation (GLO)17~532* (2~4)500500350~700650650可选双极化
9全极化条带1Quad-pol stripmap (QPSI)20~411*1886~93020~35全极化
10全极化条带2Wide quad-pol stripmap (QPSII)20~383*2252515~304035~50全极化
11波成像模式Wave (WAV)20~411*210108~125*55*5全极化
12扩展入射角低入射角Extended incidence angle (EXT)10~203*2252515~30130120~150可选双极化
高入射角50~6020~308070~90可选双极化

5. 产品分级

(不确定,我们的数据是L1A级)

产品级别产品形式定义
L1A复数据产品(SLC)根据卫星参数,进行成像处理、相对辐射校正后获得的斜距复数据产品,提供斜地转换系数;复数据产品保留幅度、相位、极化信息。
L1B单视图像产品(SLP)根据卫星参数,进行成像处理、相对辐射校正的图像数据斜距产品。
多视图像产品(MLP)根据卫星参数,进行成像处理、多视处理、相对辐射校正、拼接后获得的图像数据产品。
L2系统几何校正产品(SGC)根据卫星下传的姿轨数据,进行几何定位、地图投影、重采样后获得的系统级几何校正产品。

6. 单幅产品组成

序号名称文件格式文件命名简单介绍说明
1SAR1级影像数据Tiffsar.tif位深度16,用gdal读没有地理信息影像数据、RPC参数文件、影像元信息文件、浏览图、和拇指图以景号为唯一标识。详细介绍见下面的“各个文件介绍”
2RPC参数文件.RPB/.RPCsar.rpc保存的坐标信息,应该就是用来对图像进行地理编码的
3影像元信息文件XMLsar.meta.xml存放平台、传感器、GPS、影像、采集过程等信息
4入射角文件XMLsar.incidence.xml存放拍摄时的入射角
5浏览图文件JPEGsar.jpg
6入射角文件JPEGsar.shumb.jpg

各个文件介绍:

  • rpc文件:用于几何校正的RPC模型。简单来说,RPC模型就是有理函数纠正函数,是将地面点大地坐标与其对应的像点坐标用比值多项式关联起来,这就像数字摄影测量学上在外场用单反相机拍张照片,并求出其内外方位元素、已知对应像点坐标的大地坐标值(一般三对以上)将相片的所有像点坐标转换为大地坐标的求解过程。举个简单例子,用ENVI打开无法支持对应影像RPC文件的数据(比如GF-2)就会显示无参考坐标系统,即没有对应的坐标信息,软件需要自动链接RPC文件才得以显示影像的坐标信息;或者进行RPC校正后的影像就带有坐标信息了,之后影像需要根据问题本身、精度要求、用途等考虑是否需要进行GCPs控制点校正以及是否需要生产数字正射影像产品(DOM)。

7. GF3数据应用

行业应用领域监测内容
海洋海岛海岸带动态监测海岸线变迁监测
海岸带滨海湿地监测
海岛监测
围填海监测
浮筏养殖区监测
海洋环境灾害监测海面风场
海浪
北极海冰分布
海洋权益维护信息服务船舶态势监测
船舶识别
海上油气平台监测
减灾专题应用水体提取
建筑物群倒塌信息提取
道路损毁提取
水利专题应用地表水体监测
地表土壤含水量监测
洪涝灾害监测
干旱灾害监测
水利工程监测
气象降水量监测

GF3数据使用——PIE地理编码

PIE官网(软件下载,注册码都有)(下载的里面有用户手册,还是好用的)

1. 最原始数据展示

原始的sar.tiff(展示的是HH极化)
在这里插入图片描述

  • 位深度为16(window显示)
  • matlab读出来是(w, h, 2)的 int16的图像,无法show,单独一个通道show出来没东西。
  • python不能用PIL读出来。
  • python能用gdal读出来,但是 没有地理信息,可以show单独一个通道
  • 二者读出来数值一样,都是有正有负有零的矩阵

2. 数据导入

在这里插入图片描述
导入成功:
在这里插入图片描述

这个时候输出路径底下的文件就会有:
在这里插入图片描述

  • 用matlab读不出来
  • 用python的PIL读不出来
  • 只能用gdal读出来,但是 没有地理信息,没法show,因为是复数
  • 是(h, w)的数据,数据类型是complex64

3. 复数数据转换

选择转换的数据,选上一步的输出
在这里插入图片描述
在这里插入图片描述
然后就会自己加载进来了,这时候:
在这里插入图片描述
输出文件夹变成:
在这里插入图片描述

  • matlab读完是(h, w)的 single数据(数值都很小)
  • python用PIL可以读
  • python用gdal可以读, 没有地理信息。读出来是跟matlab一样的,python的show出来很黑,因为数值小
  • matlab的show如下跟window的有些许不一样,因为是single,各自的处理不一样:
    在这里插入图片描述

4. 地理编码

选择上一步输出的复数据影像,全极化就四张都选
在这里插入图片描述
在这里插入图片描述
地理编码结束后如下:
在这里插入图片描述- 地理编码之后图片变大了??h和w都变大了好多

  • matlab读出来是(h*, w*)single 数据,show出来如下,数值很小
  • python的PIL和gdal都可以读, 有地理信息
    在这里插入图片描述

对应电子地图下载

使用bigemap

  • 5
    点赞
  • 52
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 6
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

damonzheng46

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值