CN106959474B - 一种水云云粒子谱分布的空间反演分辨率优化方法 - Google Patents
一种水云云粒子谱分布的空间反演分辨率优化方法 Download PDFInfo
- Publication number
- CN106959474B CN106959474B CN201710161061.XA CN201710161061A CN106959474B CN 106959474 B CN106959474 B CN 106959474B CN 201710161061 A CN201710161061 A CN 201710161061A CN 106959474 B CN106959474 B CN 106959474B
- Authority
- CN
- China
- Prior art keywords
- cloud
- observation
- similarity
- resolution ratio
- cloud particle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000002245 particle Substances 0.000 title claims abstract description 106
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000003595 spectral effect Effects 0.000 title claims abstract description 28
- 238000005457 optimization Methods 0.000 title claims abstract description 27
- 230000008859 change Effects 0.000 claims abstract description 18
- 239000000284 extract Substances 0.000 claims abstract description 8
- 238000009434 installation Methods 0.000 claims abstract description 7
- 230000010287 polarization Effects 0.000 description 19
- 238000010586 diagram Methods 0.000 description 15
- 238000003860 storage Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000009826 distribution Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 229920006343 melt-processible rubber Polymers 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 4
- 239000000428 dust Substances 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000009828 non-uniform distribution Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000005427 atmospheric aerosol Substances 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V8/00—Prospecting or detecting by optical means
- G01V8/10—Detecting, e.g. by using light barriers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J4/00—Measuring polarisation of light
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本申请公开了一种水云云粒子谱分布的空间反演分辨率优化方法，解决反演精度差、计算效率低的问题。包含：计算偏振反射率查找表，包含偏振反射率的理论向量；从卫星历史观测数据计算像元的散射角和对应的偏振反射率观测值；设置空间反演分辨率，抽取散射角在145°‑165°范围内偏振反射率值构成观测向量；针对不同的云粒子有效半径和云粒子有效方差离散值，对理论向量进行曲线拟合得到多个理论曲线；对观测向量进行曲线拟合得到观测曲线；将理论曲线和观测曲线进行比较，得到与观测曲线最相似的理论曲线及相似度。改变空间反演分辨率，重复计算所述相似度，取所述相似度超过阈值时的最小反演单元所对应的空间反演分辨率，作为最佳空间反演分辨率。
Description
技术领域
本申请涉及遥感领域，尤其涉及水云云粒子谱分布反演的方法。
背景技术
地球表面约有20％-30％被水云覆盖，通过卫星数据反演得到的水云光学特性和微物理特性参数广泛应用于气候变化、遥感监测、天气预报、模式开发等领域。水云云粒子谱分布包含两个重要参数：云滴有效粒子半径和有效方差。
云粒子谱分布的反演，多角度偏振方法可以同时获得云粒子有效半径和有效方差信息。GF系列卫星搭载的大气气溶胶多角度偏振探测仪(DPC)能够提供多角度的偏振观测数据。
实际水云在水平方向上的非均匀分布，多角度偏振观测中，有效观测数据在有限的水平观测面积内，不一定完全落在偏振角度范围(145°-165°)，为了达到反演偏振角度范围要求，往往会以增加观测面积来增加观测数据，然而水云水平方向是非均匀的，观测面积的增加，会使云粒子有效半径、有效方差的反演精度降低。
发明内容
本发明提出一种水云云粒子谱分布的空间反演分辨率优化方法，解决反演精度差、计算效率低的问题。
本申请实施例提供的一种水云云粒子谱分布的空间反演分辨率优化方法，包括以下步骤：
步骤1、根据Mie散射理论计算偏振反射率随散射角、云粒子有效半径、云粒子有效方差变化的数据，组成查找表；所述查找表包含云粒子有效半径离散值和云粒子有效方差离散值条件下，偏振反射率随散射角变化的理论向量。
步骤2、从卫星历史观测数据中获取水云的多角度偏振观测信息，计算各个像元的散射角和对应的偏振反射率观测值。
步骤3、设置空间反演分辨率，在对应的反演单元内抽取散射角在145°-165°范围内任意像元的偏振反射率观测值，构成偏振反射率随散射角变化的观测向量；针对不同的云粒子有效半径离散值和云粒子有效方差离散值，对所述理论向量进行曲线拟合，得到多个理论曲线；对所述观测向量进行曲线拟合，得到观测曲线；将所述多个理论曲线和所述观测曲线进行比较，得到与所述观测曲线最相似的理论曲线及相似度。
步骤4、改变空间反演分辨率，重复上一步骤，更新所述相似度，得到所述相似度随空间反演分辨率变化的规律。取所述相似度超过阈值时的最小反演单元所对应的空间反演分辨率，作为最佳空间反演分辨率。
优选地，在本发明所述水云云粒子谱分布的空间反演分辨率优化方法的实施例中，所述空间反演分辨率所对应的反演单元，随空间反演分辨率改变而改变，依次为2×2像元、3×3像元、…、N×N像元，其中N为小于25的正整数。
优选地，本申请另一实施例提供的水云云粒子谱分布的空间反演分辨率优化方法，所述观测向量的元素数量不少于20。进一步优选地，所述观测向量相邻元素对应的散射角差异小于5°。
可选择地，所述水云云粒子谱分布的空间反演分辨率优化方法的实施例，还包含以下步骤：在任一所述空间反演分辨率条件下，改变所述观测向量中的元素数量，重复计算所述相似度，得到所述相似度随元素数量的变化规律，取所述相似度超过阈值时的最少元素数量。
可选择地，所述水云云粒子谱分布的空间反演分辨率优化方法的实施例，还包含以下步骤：改变所述观测向量中的任意相邻元素之间散射角差异，重复计算所述相似度，得到所述相似度随散射角差异的变化规律，取所述相似度超过阈值时的最大散射角差异。
在本申请所涉及的任意一项水云云粒子谱分布的空间反演分辨率优化方法的实施例中，所述云粒子有效半径离散值从5μm到20μm，间隔为0.5μm；所述云粒子有效方差离散值为0.01，0.02，0.05和0.1。
优选地，在本发明任意一项所述水云云粒子谱分布的空间反演分辨率优化方法的实施例中，所述相似度为相关系数；所述相似度超过阈值，是指所述相关系数大于预设的相关系数判决门限。作为另一可选择的方案，所述相似度为均方根误差；所述相似度超过阈值，是指所述均方根误差小于预设的均方根误差判决门限。
本申请实施例采用的上述至少一个技术方案能够达到以下有益效果：本发明针对现有被动卫星反演云粒子谱分布中，没有充分考虑水云在水平方向上非均匀性的缺陷，考虑了减少偏振角度上的观测数据，提高空间反演分辨率，减小了反演单元的面积，提高计算效率；降低水平非均匀性对反演结果精度的影响，提高云粒子谱分布反演精度。
附图说明
此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：
图1为在145°-165°范围内水云粒子偏振反射率随散射角、水云粒子有效半径、水云粒子有效方差变化的规律；
图2为水云粒子谱反演时的空间观测坐标、空间反演分辨率和反演单元示意图，其中，图2a为水云粒子谱反演时的空间观测坐标系示意图，图2b为空间反演分辨率和反演单元示意图；
图3为偏振反射率随散射角、云粒子有效半径、云粒子有效方差变化的数据查找表结构示意图；
图4为偏振反射率理论向量、观测向量、理论曲线、观测曲线关系的示意图；
图5为本发明方法确定小元素数量的实施例流程图；
图6为本发明方法确定最大散射角差异的实施例流程图；
图7为本发明方法确定最优空间反演分辨率的实施例流程图；
图8为本发明方法在最小元素数量和最大散射角差异条件下确定最优空间反演分辨率的实施例流程图；
图9为相似度随反演单元面积变化的规律、阈值示意图，其中，图9a相似度为相关系数，图9b相似度为均方根误差。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚，下面将结合本申请具体实施例及相应的附图对本申请技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
以下结合附图，详细说明本申请各实施例提供的技术方案。
图1为在145°-165°范围内水云粒子偏振反射率随散射角、水云粒子有效半径、水云粒子有效方差变化的规律。图中纵轴的数据为偏振反射相函数值，与偏振反射率是正相关的物理量。在145°-165°散射角范围内，偏振反射率(或偏振反射相函数)与偏振散射角的关系曲线称为“剩余虹”，剩余虹对水云云滴的有效粒子半径和有效方差敏感。例如图1a中，剩余虹波峰所在的散射角随着云滴有效粒子半径的减少而增加；图1b中，剩余虹的峰值随有效方差的减小而增大。由于随着有效半径和有效方差变化剩余虹的理论形状发生改变，因而可以通过实际观测的剩余虹反演水云云粒子谱(即推测实际水云粒子有效半径和有效方差)，水云粒子的偏振反射率就需要在观测散射角度范围内有足够的观测数据，用来拟合曲线。
图2为水云粒子谱反演时的空间观测坐标、空间反演分辨率和反演单元示意图；对水云粒子进行观测时，要建立空间观测坐标系，包含云平面、太阳方向、卫星观测方向(云粒子散射方向)、散射角等。图2a为水云粒子谱反演时的空间观测坐标系示意图；包含卫星观测天顶角θ’，卫星方位角φ’，太阳天顶角θ，太阳方位角φ，偏振散射角Θ。其中偏振散射角是太阳光入射方向和散射方向之间的夹角。利用当期卫星数据中的卫星方位角和太阳方位角计算象元每个角度观测的相对方位角；利用当期卫星数据中的卫星天顶角，太阳天顶角数据，和以上步骤中得到的相对方位角数据，计算象元每个角度观测的散射角Θ。
图2b为空间反演分辨率和反演单元示意图。在反演单元内，包含多个观测象元，假设在一个反演单元内的水云粒子分布相同，即各个像元中水云粒子有效半径相同、有效方差相同；采集反演单元内每个像元的历史卫星观测数据，包含卫星观测天顶角θ’，卫星方位角φ’，太阳天顶角θ，太阳方位角φ，多角度偏振反射率云、云识别和云相态数据组建样本集，对样本集中观测的散射角进行提取，获取每个象元落在145°至165°散射角范围内的观测值。然而，实际水云在水平方向上非均匀分布，在实际条件下进行多角度偏振观测，在有限的水平观测面积内及期望的散射角度范围内不一定存在有效数据。为了达到反演偏振角度范围要求，可以增加反演单元的面积来增加观测数据，然而假设水云粒子分布相同与实际水云水平非均匀的矛盾加大，观测面积的增加，也可能会使云粒子有效半径、有效方差的反演精度降低。
反演单元面积和反演分辨率是相对应的概念，面积越大分辨率越低。例如图2b中表示出云平面内4×4像元的反演分辨率及对应的反演单元面积、8×8像元的反演分辨率及对应的反演单元面积。假定反演单元面积中水云粒子有效半径相同、有效方差相同，通过比较观测的偏振反射率和理论的偏振反射率曲线，确定反演单元面积水云粒子有效半径、有效方差。如果每个像元的大小为6km×7km，8×8像元的反演单元面积大小则为48km×56km。4×4像元的反演单元面积大小则为24km×28km。在云平面内观测范围中，包含大量反演单元，在多个反演单元上获得的水云粒子参数值的集合，即云平面内观测范围中水云粒子分布状况。
图3为偏振反射率随散射角、云粒子有效半径、云粒子有效方差变化的数据查找表结构示意图。例如，在结构同中的横向，所述云粒子有效半径离散值从5μm到20μm，间隔为0.5μm；所述云粒子有效方差离散值为0.01，0.02，0.05和0.1。在结构图中的纵向，表示散射角离散值，从145°-165°间隔为0.5°。结构图中每一个最小方格，代表一个元素，包含一个偏振反射率值。
图4为偏振反射率理论向量、观测向量、理论曲线、观测曲线关系的示意图。图中，xi＝RM(Θi)为确定的云粒子有效半径reff＝reff0、云粒子有效方差veff＝veff0条件下，偏振反射率随离散的散射角Θi变化的理论值；yj＝RO(Θj)为在观测的云粒子有效半径、云粒子有效方差条件下，偏振反射率随离散的散射角Θj变化的理论值(其中，i，j为自然数)。
X＝{xi}为偏振反射率理论向量；Y＝{yj}为偏振反射率理论向量；
MPRs为X的离散值xi进行曲线拟合后得到的理论值曲线；
OPRs为Y的离散值yj进行曲线拟合后得到的观测值曲线；
当MPRs和OPRs最相似时，认为所述观测的云粒子有效半径≈reff0、所述观测的云粒子有效方差≈veff0。
所述MPRs和OPRs相似的程度，用衡量相似程度的参数(简称相似度)来表示。
具体地，例如以相关系数作为相似度，有T1＝correlate(OPRs,MPRs)，其中T1为相关系数；相关系数约大(趋近于1)表示两条曲线越相似。
或者，例如以均方根误差作为相似度，有其中T2为均方根误差；均方根误差越小(趋近于0)表示两条曲线越相似。SA表示散射角，N表示MPRs和OPRs在计算T2时的点数，每个点对应一个散射角。
图5为本发明方法确定小元素数量的实施例流程图。本实施例提供的一种水云云粒子谱分布的空间反演分辨率优化方法，包括以下步骤：
步骤101、根据Mie散射理论计算偏振反射率随散射角、云粒子有效半径、云粒子有效方差变化的数据，组成查找表；所述查找表包含云粒子有效半径离散值和云粒子有效方差离散值条件下，偏振反射率随散射角变化的理论向量。
步骤102、从卫星历史观测数据中获取水云的多角度偏振观测信息，计算各个像元的散射角和对应的偏振反射率观测值。
步骤103、设置空间反演分辨率，在对应的反演单元内抽取散射角在145°-165°范围内任意像元的偏振反射率观测值，构成偏振反射率随散射角变化的观测向量；针对不同的云粒子有效半径离散值和云粒子有效方差离散值，对所述理论向量进行曲线拟合，得到多个理论曲线；对所述观测向量进行曲线拟合，得到观测曲线；将所述多个理论曲线和所述观测曲线进行比较，得到与所述观测曲线最相似的理论曲线及相似度。
步骤104、在任一所述空间反演分辨率条件下，改变所述观测向量中的元素数量，重复计算所述相似度，得到所述相似度随元素数量的变化规律，取所述相似度超过阈值时的最少元素数量。
图6为本发明方法确定最大散射角差异的实施例流程图。本实施例提供的一种水云云粒子谱分布的空间反演分辨率优化方法，包括以下步骤：
步骤201、根据Mie散射理论计算偏振反射率随散射角、云粒子有效半径、云粒子有效方差变化的数据，组成查找表；所述查找表包含云粒子有效半径离散值和云粒子有效方差离散值条件下，偏振反射率随散射角变化的理论向量。
步骤202、从卫星历史观测数据中获取水云的多角度偏振观测信息，计算各个像元的散射角和对应的偏振反射率观测值。
步骤203、设置空间反演分辨率，在对应的反演单元内抽取散射角在145°-165°范围内任意像元的偏振反射率观测值，构成偏振反射率随散射角变化的观测向量；针对不同的云粒子有效半径离散值和云粒子有效方差离散值，对所述理论向量进行曲线拟合，得到多个理论曲线；对所述观测向量进行曲线拟合，得到观测曲线；将所述多个理论曲线和所述观测曲线进行比较，得到与所述观测曲线最相似的理论曲线及相似度。
步骤204、改变所述观测向量中的任意相邻元素之间散射角差异，重复计算所述相似度，得到所述相似度随散射角差异的变化规律，取所述相似度超过阈值时的最大散射角差异。
需要说明的是，按照图5-6实施例，优化的结果表明，在145°-165°偏振角度范围内，观测数值减少到20个，偏振角度间隔小于5°时，可以不影响水云粒子谱分布的提取精度，同时有效提高观测空间反演分辨率，减小水云在水平面上非均匀对反演结果精度的影响。
图7为本发明方法确定最优空间反演分辨率的实施例流程图。本实施例提供的一种水云云粒子谱分布的空间反演分辨率优化方法，包括以下步骤：
步骤301、根据Mie散射理论计算偏振反射率随散射角、云粒子有效半径、云粒子有效方差变化的数据，组成查找表；所述查找表包含云粒子有效半径离散值和云粒子有效方差离散值条件下，偏振反射率随散射角变化的理论向量。
步骤302、从卫星历史观测数据中获取水云的多角度偏振观测信息，计算各个像元的散射角和对应的偏振反射率观测值。
步骤303、设置空间反演分辨率，在对应的反演单元内抽取散射角在145°-165°范围内任意像元的偏振反射率观测值，构成偏振反射率随散射角变化的观测向量；针对不同的云粒子有效半径离散值和云粒子有效方差离散值，对所述理论向量进行曲线拟合，得到多个理论曲线；对所述观测向量进行曲线拟合，得到观测曲线；将所述多个理论曲线和所述观测曲线进行比较，得到与所述观测曲线最相似的理论曲线及相似度。
步骤304、改变空间反演分辨率，重复上一步骤，更新所述相似度，得到所述相似度随空间反演分辨率变化的规律。例如，在本发明所述水云云粒子谱分布的空间反演分辨率优化方法的实施例中，所述空间反演分辨率所对应的反演单元，随空间反演分辨率改变而改变，依次为2×2像元、3×3像元、…、直至N×N像元(N为正整数，一般地，N在小于25的范围内选取)。取所述相似度超过阈值时的最小反演单元所对应的空间反演分辨率，作为最佳空间反演分辨率。
作为进一步优化的实施例，还包含以下步骤
改变所述观测向量中的任意相邻元素之间散射角差异，重复计算所述相似度，得到所述相似度随散射角差异的变化规律，取所述相似度超过阈值时的最大散射角差异。
改变所述观测向量中的元素数量，重复计算所述相似度，得到所述相似度随元素数量的变化规律，取所述相似度超过阈值时的最少元素数量。
图8为本发明方法在最小元素数量和最大散射角差异条件下确定最优空间反演分辨率的实施例流程图
步骤401、根据Mie散射理论计算偏振反射率随散射角、云粒子有效半径、云粒子有效方差变化的数据，组成查找表；所述查找表包含云粒子有效半径离散值和云粒子有效方差离散值条件下，偏振反射率随散射角变化的理论向量。
步骤402、从卫星历史观测数据中获取水云的多角度偏振观测信息，计算各个像元的散射角和对应的偏振反射率观测值。
步骤403、设置空间反演分辨率，在对应的反演单元内抽取散射角在145°-165°范围内任意像元的偏振反射率观测值，构成偏振反射率随散射角变化的观测向量；
步骤404、判断：所述观测向量的元素数量不少于20、所述观测向量相邻元素对应的散射角差异小于5°。如果是，则对所述观测向量进行曲线拟合，得到观测曲线；
步骤405、针对不同的云粒子有效半径离散值和云粒子有效方差离散值，对所述理论向量进行曲线拟合，得到多个理论曲线；将所述多个理论曲线和所述观测曲线进行比较，得到与所述观测曲线最相似的理论曲线及相似度。
步骤406、改变空间反演分辨率，重复步骤403-405，更新所述相似度，得到所述相似度随空间反演分辨率变化的规律。取所述相似度超过阈值时的最小反演单元所对应的空间反演分辨率，作为最佳空间反演分辨率。
图9为相似度随反演单元面积变化的规律、阈值示意图。图9a相似度为相关系数；图9b相似度为均方根误差。优选地，在本发明任意一项所述水云云粒子谱分布的空间反演分辨率优化方法的实施例中，所述相似度为相关系数；所述相似度超过阈值，是指所述相关系数T1大于预设的相关系数判决门限P1。作为另一可选择的方案，所述相似度为均方根误差；所述相似度超过阈值，是指所述均方根误差T2小于预设的均方根误差判决门限P2。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (10)
1.一种水云云粒子谱分布的空间反演分辨率优化方法，其特征在于，包括以下步骤：
根据Mie散射理论计算偏振反射率随散射角、云粒子有效半径、云粒子有效方差变化的数据，组成查找表；所述查找表包含云粒子有效半径离散值和云粒子有效方差离散值条件下，偏振反射率随散射角变化的理论向量；
从卫星历史观测数据中获取水云的多角度偏振观测信息，计算各个像元的散射角和对应的偏振反射率观测值；
设置空间反演分辨率，在对应的反演单元内抽取散射角在145°-165°范围内任意像元的偏振反射率观测值，构成偏振反射率随散射角变化的观测向量；针对不同的云粒子有效半径离散值和云粒子有效方差离散值，对所述理论向量进行曲线拟合，得到多个理论曲线；对所述观测向量进行曲线拟合，得到观测曲线；将所述多个理论曲线和所述观测曲线进行比较，得到与所述观测曲线最相似的理论曲线及相似度；
改变空间反演分辨率，重复上一步骤，更新所述相似度，得到所述相似度随空间反演分辨率变化的规律；
取所述相似度超过阈值时的最小反演单元所对应的空间反演分辨率，作为最佳空间反演分辨率。
2.如权利要求1所述水云云粒子谱分布的空间反演分辨率优化方法，其特征在于，所述观测向量的元素数量不少于20。
3.如权利要求1所述水云云粒子谱分布的空间反演分辨率优化方法，其特征在于，所述观测向量相邻元素对应的散射角差异小于5°。
4.如权利要求1所述水云云粒子谱分布的空间反演分辨率优化方法，其特征在于，还包含以下步骤：
在任一所述空间反演分辨率条件下，改变所述观测向量中的元素数量，重复计算所述相似度，得到所述相似度随元素数量的变化规律，取所述相似度超过阈值时的最少元素数量。
5.如权利要求4所述水云云粒子谱分布的空间反演分辨率优化方法，其特征在于，还包含以下步骤；
改变所述观测向量中的任意相邻元素之间散射角差异，重复计算所述相似度，得到所述相似度随散射角差异的变化规律，取所述相似度超过阈值时的最大散射角差异。
6.如权利要求1-5任意一项所述水云云粒子谱分布的空间反演分辨率优化方法，其特征在于，所述云粒子有效半径离散值从5μm到20μm。
7.如权利要求1-5任意一项所述水云云粒子谱分布的空间反演分辨率优化方法，其特征在于，所述云粒子有效方差离散值为0.01，0.02，0.05和0.1。
8.如权利要求1-5任意一项所述水云云粒子谱分布的空间反演分辨率优化方法，其特征在于，所述相似度为相关系数；所述相似度超过阈值，是指所述相关系数大于预设的相关系数判决门限。
9.如权利要求1-5任意一项所述水云云粒子谱分布的空间反演分辨率优化方法，其特征在于，所述相似度为均方根误差；所述相似度超过阈值，是指所述均方根误差小于预设的均方根误差判决门限。
10.如权利要求1-5任意一项所述水云云粒子谱分布的空间反演分辨率优化方法，其特征在于，所述空间反演分辨率所对应的反演单元，随空间反演分辨率改变而改变，依次为2×2像元、3×3像元、…、直至N×N像元，其中N为小于25的正整数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710161061.XA CN106959474B (zh) | 2017-03-17 | 2017-03-17 | 一种水云云粒子谱分布的空间反演分辨率优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710161061.XA CN106959474B (zh) | 2017-03-17 | 2017-03-17 | 一种水云云粒子谱分布的空间反演分辨率优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106959474A CN106959474A (zh) | 2017-07-18 |
| CN106959474B true CN106959474B (zh) | 2018-10-30 |
Family
ID=59470370
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710161061.XA Active CN106959474B (zh) | 2017-03-17 | 2017-03-17 | 一种水云云粒子谱分布的空间反演分辨率优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106959474B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110411918B (zh) * | 2019-08-02 | 2020-09-08 | 中国科学院遥感与数字地球研究所 | 一种基于卫星偏振技术的pm2.5浓度遥感估算方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102298150B (zh) * | 2011-05-23 | 2012-11-14 | 北京师范大学 | 全球陆表宽波段发射率反演方法及系统 |
| CN103940515B (zh) * | 2013-11-14 | 2016-06-08 | 北京航空航天大学 | 一种波浪水面天空漫反射光偏振场模拟方法 |
| US9964482B2 (en) * | 2015-01-26 | 2018-05-08 | The United States Of America As Represented By The Secretary Of The Army | Detecting clouds using polarized sunlight |
| CN105954812B (zh) * | 2016-04-29 | 2018-09-07 | 中国科学院遥感与数字地球研究所 | 一种确定云相态的方法和设备 |
-
2017
- 2017-03-17 CN CN201710161061.XA patent/CN106959474B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106959474A (zh) | 2017-07-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9208612B2 (en) | Systems and methods that generate height map models for efficient three dimensional reconstruction from depth information | |
| Korolev | Reconstruction of the sizes of spherical particles from their shadow images. Part I: Theoretical considerations | |
| CN107895348B (zh) | 一种在散射环境中非均匀光场下的偏振图像复原方法 | |
| US8971640B1 (en) | Image alignment | |
| Hetem Jr et al. | Fractal 3-D simulations of molecular clouds | |
| CN109581383A (zh) | 各向异性海面雷达后向散射模拟方法及系统 | |
| CN109614459A (zh) | 应用于二维激光的地图构建回环检测方法及设备 | |
| CN113012063B (zh) | 一种动态点云修复方法、装置及计算机设备 | |
| CN103562966A (zh) | 景深栅格化 | |
| US9875575B2 (en) | Smoothing 3D models of objects to mitigate artifacts | |
| WO2015094956A1 (en) | Specularity determination from images | |
| CN106959474B (zh) | 一种水云云粒子谱分布的空间反演分辨率优化方法 | |
| CN105737806A (zh) | 基于星载微波辐射计反演海面盐度的方法和装置 | |
| CN105954812B (zh) | 一种确定云相态的方法和设备 | |
| CN110910435B (zh) | 建筑物点云提取方法、装置、计算机设备及可读存储介质 | |
| CN105701284A (zh) | 时变多尺度电大区域海面电磁散射矢量场的并行计算方法 | |
| CN108364349A (zh) | 墙面扫描建模方法、装置、终端及计算机可读存储介质 | |
| Huang et al. | Comparing TID simulations using 3-D ray tracing and mirror reflection | |
| CN109557588B (zh) | 一种煤矿井下二维矿震波速反演降维方法 | |
| CN106680157B (zh) | 一种确定大颗粒水云云粒子谱分布的方法 | |
| CN104574428B (zh) | 一种sar图像入射角估算方法 | |
| US8374391B1 (en) | Model of visibility of targets | |
| CN113724229B (zh) | 高程差的确定方法、装置及电子设备 | |
| Zhang et al. | A review for retrieving wind fields by spaceborne synthetic aperture Radar | |
| US20130287303A1 (en) | Display system with edge map conversion mechanism and method of operation thereof |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |