海浪信息检索使用模拟紧凑Radarsat-2极化SAR

文摘

本文的主要目的是演示的功能紧凑极化合成孔径雷达(SAR) (CP)来获取海洋波场参数。Souyris北国”和重建的算法用于执行CP SAR伪quad-polarized数据海洋表面下的循环线性接收(CTLR)和传输π/ 4模式。结果表明,CP重建、北部的算法有更好的收敛能力比Souyris”。此外,重建精度的调查表明,CTLR模式优于π/ 4模式下,海洋表面的重建。因此,得出的结论是,北国CP CTLR模式数据的重建参数的算法适应检索海浪信息。海浪斜率光谱和其他主波参数也从重建CP数据计算并与原位测量国家数据浮标中心(NDBC)匹配浮标。CP与浮标SAR-based波场信息输出的比较也显示良好的协议控制波高的情况下,波方向,和波的周期,偏差为0.36米,17.96°,和0.88 s,分别。

1。介绍

海浪是最重要的海洋表面特性在海洋领域的物理和海岸工程。海洋表面参数,如波高、波方向和波长,对海冰监测很重要,船舶检测、海洋污染预报,和在其他领域1,2]。波的斜率光谱,它提供了关于二维方向分布的信息波的能量密度,海浪信息检索是重要的。近年来,遥感作为地球表面的有效工具监测逐渐在海洋信息检索中发挥着关键作用。所有的遥感传感器、SAR系统主动卫星传感器的分辨率足够高,提供的信息海洋波谱,所以这些系统起着关键的作用在海洋表面波的测量3- - - - - -5]。与原位浮标、特别行政区可以在各种天气条件或在远离海岸的地区和全球可以测量海浪不断6]。除了提高对海洋的性质的理解,SAR观测可以用来开发操作波模型的预测精度。因此,重要的是研究SAR海洋表面进行监测的能力,尤其是海浪检索。

近几十年来,卫星,如ERS-1 ERS-2,萨特,Radarsat-1/2, Sentinel-1 A / B,中国GF-3进行许多SAR系统用于监测海洋表面。几波检索算法也为海洋开发的信息检索。例如,k . Hasselmann和美国Hasselmann [7,8)提出了马克斯普朗克研究所(MPI)算法的基础上,封闭的非线性变换海洋光谱和SAR图像光谱之间的关系;波的输出模型(WAM)作为初始假想输入谱方法。Mastenbroek和de Valk9)开发了一种半参数检索算法(SPRA)基于Hasselmann非线性理论(7)通过将风力波和增加风浪资料由同步散射仪测量和膨胀提供信息是来自SAR信号的其余部分。王(10)提出了一种半经验的算法有客观估计的波高从萨特ASAR波模式图像没有任何先验知识。邵et al。11有效波高)开发一种半经验的函数(Hs)和平均波的周期(Tmw)检索从c波段VV极化SAR Sentinel-1。随着极化技术的发展,quad-polarized SAR海洋表面也被应用于监测。Zhang et al。12)提供一个协同方法检索两个海洋表面波和风从星载SAR quad-polarization字段,考虑海浪谱之间的非线性映射关系和SAR图像光谱,为了获取风场和波浪方向谱增效剂。舒勒和李13,14)提出了一个算法获取海浪谱还使用quad-polarized SAR图像和调查结果使用L -和p波段SAR数据。然而,这种方法只适用于L和P乐队。通过修改极化取向角,他等。15,16)开发了一种新的检索算法,可以应用于c波段quad-polarized SAR数据。他等的方法也被验证使用Radarsat-2 quad-polarized SAR数据Zhang et al。17,18]。

此外,尽管quad-polarized SAR可以提供丰富的极化信息开发海洋信息检索的过程中,有限的宽度也成为大规模监控quad-polarized SAR的缺点。克服的局限性的CP SAR,一种新型的SAR传感器叫做紧凑p (CP) SAR已经发明,它提供了一种新的方法来测量海浪和其潜在的用于海洋表面监测逐渐变得清晰(19]。与quad-polarized SAR相比,CP SAR在功率较低,但仍能够重建quad-polarized数据,从而克服缺点的狭窄范围和高功率quad-polarized SAR极化单所提供的信息不足或双极化SAR (20.,21]。CP SAR可以,事实上,被视为双极化之间的妥协和quad-polarized SAR这片宽度增加,同时减少了运行功率。迄今为止,印度RISAT-1推出2012年4月,日本ALOS-2 SAR卫星于2014年5月,进行CP SAR。在未来,计划Radarsat-2还将支持CP SAR的卫星星座。尽管CP SAR进行现有的SAR平台,仍然没有精确校准的数据,可用于研究,特别是在海浪的情况下测量。因此,重要研究潜在的CP SAR海浪信息的检索。在先前的研究中,Nunziata et al。22]研究了海面散射极化模型,利用混合极化(HP)模式来解释海洋表面极化参数。阴et al。23)提出了一种新的扩展X-Bragg模型并使用c波段quad-polarized SAR数据收购SIR-C / X-SAR和Radarsat-2惠普和效仿π/ 4 CP模式。Geldsetzer et al。24)提出,入射角之间的依赖,风速和风向的CP SAR参数提出了开放的水。选择20 CP参数c波段地球物理模型相关函数(CMOD)输出。从极化Radarsat-2 CP参数模拟数据和仿真数据可以在等待Radarsat星座任务(RCM)。李等人。25)提出一个新的经验重建算法修改N基于数据集的2000多Radarsat-2 quad-polarization图像和集中的浮标观测。这是证明,交叉极化促进风速检索,提高了结果的准确性,尤其是高风速。总之,CP SAR的能力重建伪quad-pol数据和海洋风检索在详细调查;然而,海洋表面重建的潜力尚未进行讨论,更不用说波信息检索。

本文使用两个quad-polarized SAR图像重建CP伪quad-polarized数据和演示的潜力CP SAR海浪信息检索。本文中使用的数据集所示部分2。节3,我们简要介绍了CP理论和海浪的检索方法。节4,我们调查两个CP算法用于海洋表面重建和评估海洋表面重建的最佳CP模式。FP和CP海浪检索结果与NDBC浮标数据。结论给出了部分5。

2。数据集

2.1。SAR数据

两个场景Radarsat-2 quad-polarized图像被用于这项研究。两个图像采集参数都是列在表中1以及研究区域的数据所示1(一)和1 (b)。HH -高压和VV-polarized单看起来复杂的图像(SLC) SAR图像也如图1(一)和1 (b)。

2.2。浮标数据

为了验证CP SAR波检索结果,NDBC浮标数据可以从国家海洋和大气管理局(NOAA)收集与SAR图像。主波参数提取的原始测量数据和表中列出2。在所有的波测量数据,占主导地位的波高(DWH),占主导地位的波的周期(劳务和退休金部)和主波方向(DWD)是用于验证检索结果。

由于有限的集合NDBC浮标,只有一个浮标与SAR图像场景2的配置。所以场景1执行CP的图像重建和海洋坡谱检索方法,和场景2添加验证方法。

3所示。理论和方法

3.1。CP重建

CP理论是一个新开发的方法,可以产生伪quad-polarized双极化SAR数据的数据。偏振信息可以重建在满足精度的降低传感器操作能力(26,27]。CP有三种主要模式:CTLR、双圆极化(DCP)和π/ 4模式。

3.1.1。CTLR模式

在CTLR模式下,微波在左旋或右旋圆偏振传输和接收在HH和VV极化。散射向量表示为

有关埃尔米特协方差矩阵 在哪里代表矩阵转置操作,代表共轭矩阵操作,〈〉代表平均操作,复数的实部,极化散射矩阵和p, q是水平或垂直极化通道,是虚数单位,是复数的虚部。也明显这个矩阵可以简化基于假设反射是对称的。

为了重建一个伪quad-polarized协方差矩阵的CP CTLR模式,也需要迭代解决。该协方差矩阵可以表示为

它是注意到在(4从CP)元素2×2的协方差矩阵。

3.1.2。DCP模式

各区总监模式中,右旋圆偏振的微波传输和接收在左旋或右旋圆偏振。散射向量表示为

有关埃尔米特协方差矩阵

应该注意的是,最后一项可以设置为零的假设的基础上反射对称,协方差矩阵可以简化

为了重建一个伪quad-polarized CP DCP的协方差矩阵模式,需要迭代解决。该协方差矩阵可以表示为(28]

它是注意到在(8从CP)元素2×2的协方差矩阵。

3.1.3。π/ 4模式

在π/ 4模式下,微波传输在线性极化45°倾角和接收在HH和VV极化。散射向量表示为

有关埃尔米特协方差矩阵

应该注意的是,最后一项可以设置为零的假设的基础上反射对称,协方差矩阵可以简化

为了重建伪quad-polarized CP的协方差矩阵π/ 4模式下,需要迭代解决。的伪quad-polarized协方差矩阵π/ 4模式可以表示为(28]

它是注意到在(12从CP)元素2×2的协方差矩阵。

还应该指出的是,散射向量和是相关的

的偏振信息CTLR和DCP模式将是相似的,因为上面的线性关系,这意味着这两个模式的CP重建结果将是相同的(29日,30.]。

3.2。波谱检索算法

线性极化后向散射截面,quad-polarized协方差矩阵用于海洋表面波极化签名,如下(31日]:

它是注意到在(14从FP)元素3×3的协方差矩阵。因为波斜率最大的贡献是获得只使用线性偏振,这里,我们将椭圆率设置为零,所以线性极化后向散射截面任意极化取向角可以派生(14,15]: 在哪里 , ,和是水平的线性极化后向散射横截面,十字架,分别和垂直偏振。水平和垂直偏振之间的相关性,可以直接源自于协方差矩阵。

在海浪的线性调制理论,海面高度可以视为一个正弦波叠加和雷达后向散射截面表示为(16,32,33] 在哪里海洋表面的高度,是波数,是波空间向量,是波角频率,本系列的复共轭。是SAR的调制传递函数,包括倾斜调制,流体力学调制、极化角调制,调制范围和方位方向,和速度聚束调制。在所有调制函数,只有倾斜和极化取向角调制依赖于雷达极化(16]。此外,波斜率诱导极化取向角调制可以忽略的水平或垂直极化(16]。由(18),对偏振不敏感的条款消失和方程可以简化为16] 在哪里和分别的复共轭系列;和是由(16] 在哪里雷达入射波数和吗是迎角。和是由(16)(16]。通过合并方程,简单的代数计算可以概括为(16] 在哪里 和 分别和方位波斜率范围。参数和给出了(8)[16]。用一个适当的选择 ,斜率和方位范围波斜率可以直接获得。波坡谱,另一波参数可以从波估计斜率。

4所示。结果与讨论

4.1。重建算法的比较

除了反射对称假设,copolarized和正交极化响应之间的关系是最主要的条件限制CP重建方程的解决方案(21,34]。

尽管Souyris et al。(21和Nord et al。34)算法用于重建伪quad-polarized数据,他们有不同的定义 。Souyris等人使用4的值 ,而北方等人计算值

因此,我们使用了两种重建算法来计算斜率之间的速度曲线 和 (图2)。重要的是要注意这一点在图2(一个)等于4在图2 (b) 是由 的价值,从FP数据计算。我们发现图的融合2(一个)远远不如图吗2 (b);相关系数, ,分别为0.53和0.93,。确定最优算法,我们还计算统计的nonconvergence像素在重建过程中。结果如表所示3。nonconvergent像素的数量发现当使用两个子图象的海洋表面是北方的算法比Souyris低得多的算法。北方的算法,因此,更适合重建伪quad-polarized海洋表面的数据。

4.2。评价CP重建精度

CP特区有三种主要模式:CTLR模式,DCP模式,π/ 4模式。每个CP模式都有不同的取向和椭圆率角,它描述波传播和接收状态。因为DCP和CTLR模式之间的线性关系,这两个模式有相似的极化信息;因此,我们只讨论π/ 4模式和CTLR模式。图3散点图显示详细的派生伪quad-pol原始quad-pol的结果符合。图中显示的值和FP和CP之间都集中在一对一的线,这意味着和 ,重建的伪quad-polarized数据很好。然而,重建显示了一个轻微的低估CTLR模式和的高估π/ 4模式,这意味着重建不一样好吗和 。的重建精度低的原因可能与疲软的体积在海洋表面散射。之前的记录,CP重建算法通常应用于田野和森林,在散射散射包含一个强大的总量组件。然而,交叉极化信息不是直接用于我们的论文,因此,重建的性能低劣不会影响检索结果。

更详细地检查重建性能,我们也选择7种极化参数,包括copolarization反向散射系数 , ,交叉极化后向散射系数 ,copolarization相关系数 ,和相干矩阵特征值 , ,和 ,对于每个参数,计算CTLR和重建精度π/ 4模式。重建精度 用于描述偏差的数量在CP重建。的重建精度 , , , , , ,和表中列出4。我们使用了描述“高”、“媒介”,和“低”来描述偏差,在“高”被定义为小于10%,“媒介”−10%与30%,和“低”被定义为超过30%。表4表明,CTLR和π/ 4模式,重建精度的参数 , , ,和与值高,CTLR模式优于的吗π/ 4模式。因此,我们决定使用CRTL模式为随后的海洋表面参数重建和海浪斜率光谱检索。

4.3。检索海浪斜率光谱

在本节中,我们描述北方的算法是用来重建伪quad-pol CTLR模式中的数据。检索到的海浪斜率基于重建的伪quad-polarized数据的光谱图所示4。与海浪斜率光谱来源于FP数据,结果CP CTLR模式这一研究获得的满足。同样清楚的是,波斜率光谱来源于FP相比提供了更多的细节,来自CP的数据。然而,这两个波光谱斜率总体上是相似的。

为了使波的定量评价检索结果,占主导地位的主波参数波高、波方向,和波的周期也从波斜率光谱中提取并与NDBC浮标数据。主波参数的表达式 在哪里代表了波的均方根斜率,代表了波长,和 ,分别代表的组件在方位方向和范围,是波数占主导地位,重力加速度,水的深度。

图5显示四子图象的位置执行检索结果统计,和180°模棱两可的波斜率频谱被移除根据浮标测量风向。结果见表5显示的平均错误占主导地位的波高、波方向,和波的周期来源于FP数据是0.14,17.12和0.61,而来自CP数据是0.36,17.96,和0.88,分别。值得注意的是,尽管平均错误主要的波高,波方向,和波的周期来源于FP数据远远大于CP,检索到的CP SAR对海浪也满意的测量结果。这也是事实使用CP SAR在这项研究中获得的检索精度有限的重建CP伪quad-pol蠕升尤其由于精度的损失 。在未来的研究中,检索精度获得使用CP SAR需要改善。

5。结论

在本文中,我们调查了CP海浪测量的能力和潜力。一个Radarsat-2 quad-polarized图像被用来重建CP伪quad-polarized数据和验证潜在的海洋坡谱检索。本研究的主要结论如下。

海洋表面波,北方比Souyris重建算法更有效。因为它的收敛能力和迭代时间,我们建议采用北方的CP伪quad-polarized数据的重建算法。

CTLR和π/ 4模式,重建精度的参数 , , ,和高而 , ,和这是低的。检索的海浪斜率光谱使用CP模式,CTLR模式优越的π/ 4模式,因此检索海浪斜率的光谱,我们建议使用CTLR模式。

CP的潜在海浪斜率光谱检索已经证明。主波参数的值也使用原位NDBC浮标数据进行验证。比较CP SAR-based波场信息与浮标输出显示良好的协议主要的波高、波方向,和波的周期,偏差为0.36米,17.96°,和0.88 s,分别。总的来说,检索结果的海浪斜率光谱和海浪参数使用CP SAR是令人满意的。此外,只考虑单一波SAR图像进行检索算法研究中,海浪信息检索的普遍性验证CP SAR仍然需要更多的结果支持。在未来的研究中,我们将收集更多quad-pol SAR匹配与浮标绘制波参数精确FP和CP SAR的区别。

的利益冲突

作者宣称他们没有关于这篇文章的出版的利益冲突。

确认

作者欣然承认提供的金融支持中国国家重点研发项目(2016 yfb0502504 2016 yfb0502500)和国家科学基金会的中国批准号。41301500,41431174,61471358,41401427。这篇文章在一定程度上是由中国科学院奖学金资助。作者还想感谢张教授从第一海洋学研究所,国家海洋局,支持研究CP理论。作者还想承认扁博士和刘博士从遥感与数字地球抛光原稿。

引用

1. c·j·陈,李解释道,胡,j .魏”提高产品开放海洋卫星数据与计划纠正剩余错误在遥感反射率,”海洋地球物理研究杂志》上,卷121,不。6,3866 - 3886年,2016页。视图:出版商的网站|谷歌学术搜索
2. j·陈,x,周b和d,“推导有色溶解有机物吸收系数与神经准解析算法从海洋颜色,”海洋地球物理研究杂志》上,卷122,不。11日,第8556 - 8543页,2017年。视图:出版商的网站|谷歌学术搜索
3. s . k .查图尔维迪·Shanmugam c·s·杨和Guven,“海浪参数的检测使用合成孔径雷达(SAR)数据,”期刊导航,卷66,不。2、283 - 293年,2013页。视图:出版商的网站|谷歌学术搜索
4. d·金,”风检索从x波段SAR图像海洋散射模型,用数值”韩国《遥感,25卷,不。3、243 - 253年,2009页。视图:谷歌学术搜索
5. m .杉:Shiroto k .大内,“估计海浪高度使用合成孔径雷达的极化率数据,”2011年IEEE国际地球科学和遥感研讨会(雪茄烟)温哥华,页2821 - 2824年,公元前,加拿大,2011年7月。视图:出版商的网站|谷歌学术搜索
6. p·w·瓦尚、h·e·Krogstad和j·s·帕特森“机载和星载合成孔径雷达观测海浪,“硕士,32卷,不。1,第112 - 83页,1994。视图:出版商的网站|谷歌学术搜索
7. k . Hasselmann和美国Hasselmann海浪谱的非线性映射成一个合成孔径雷达图像谱及其反演,”海洋地球物理研究杂志》上,卷96,不。C6, 10713 - 10729年,1991页。视图:出版商的网站|谷歌学术搜索
8. Hasselmann和k . Hasselmann”计算和参数化的非线性能量转移光谱台引力波。第一部分:新方法有效地计算精确的非线性转移积分,“《物理海洋学,15卷,不。11日,第1377 - 1369页,1985年。视图:出版商的网站|谷歌学术搜索
9. c . Mastenbroek和c . f . de Valk“半参数算法检索从合成孔径雷达海浪谱,”海洋地球物理研究杂志》上,卷105,不。C2, 3497 - 3516年,2000页。视图:出版商的网站|谷歌学术搜索
10. 朱j . h . Wang, j·杨和c·施“半经验的算法对SAR波高检索及其验证使用环境ASAR数据波模式”Oceanologica学报没有,卷。31日。3,59 - 66年,2012页。视图:出版商的网站|谷歌学术搜索
11. w·邵z张、李x和h·李,“海浪参数从Sentinel-1 SAR图像检索,“遥感,8卷,不。9日,2016年。视图:出版商的网站|谷歌学术搜索
12. 张,x, w . Perrie, y,他“协同测量从SAR海洋风和海浪,“海洋地球物理研究杂志》上,卷120,不。9日,第6184 - 6164页,2015年。视图:出版商的网站|谷歌学术搜索
13. d·l·舒勒和j·s·李,“微波技术提高定向海浪谱的测量,”国际遥感杂志》上,16卷,不。2、199 - 215年,1995页。视图:出版商的网站|谷歌学术搜索
14. d·l·舒勒j·s·李,d . Kasilingam Pottier大肠,“海洋表面测量斜坡和波光谱使用极化SAR图像数据”环境遥感,卷91,不。2、198 - 211年,2004页。视图:出版商的网站|谷歌学术搜索
15. y, h .沈,w . Perrie“极化SAR遥感海浪的。”大气和海洋技术杂志》上,23卷,不。12日,第1773 - 1768页,2006年。视图:出版商的网站|谷歌学术搜索
16. y, w . Perrie, t·谢,问:邹,“从线性偏振SAR海浪谱。”IEEE地球科学和遥感,42卷,不。11日,第2631 - 2623页,2004年。视图:出版商的网站|谷歌学术搜索
17. 张b、w . Perrie和y,他“验证RADARSAT-2完全偏振SAR测量海洋表面波,“海洋地球物理研究杂志》上,卷115,不。C6, 302 - 315年,2010页。视图:出版商的网站|谷歌学术搜索
18. y, b . Zhang, w . Perrie”验证RADARSAT-2极化SAR测量海浪”2009年IEEE国际地球科学和遥感研讨会第三,页3 - 168 - 171,开普敦,南非,2009年7月。视图:出版商的网站|谷歌学术搜索
19. f r . l . Paes Nunziata m .阶段”的能力hybrid-polarity特性观察金属海上目标,“IEEE海洋工程》杂志上第41卷。。2、346 - 361年,2016页。视图:出版商的网站|谷歌学术搜索
20. r·k·兰尼“紧凑的偏振测定角度看,”IEEE地球科学和遥感信卷。160年,12 - 18,2011页。视图:谷歌学术搜索
21. J.-C。Souyris, p .大修,r . Fjortoft s Mingot, js。李,“紧凑型偏振测定基于对称性的地球物理媒体:π/ 4模式”,IEEE地球科学和遥感,43卷,不。3、634 - 646年,2005页。视图:出版商的网站|谷歌学术搜索
22. f . Nunziata m .阶段,李x”海上浮油观察使用hybrid-polarity SAR架构。”IEEE海洋工程》杂志上,40卷,不。2、426 - 440年,2015页。视图:出版商的网站|谷歌学术搜索
23. j .阴j .杨z s周和j·歌,“扩展的布拉格散射模型船和油污观测方法使用偏振SAR紧凑,”IEEE选择杂志的主题应用地球观测和遥感,8卷,不。8,3760 - 3772年,2015页。视图:出版商的网站|谷歌学术搜索
24. t . Geldsetzer f•夏博诺、m . Arkett和t . Zagon“海洋风能研究使用模拟RCM compact-polarimetry特区,“加拿大遥感杂志》上第41卷。。5,418 - 430年,2015页。视图:出版商的网站|谷歌学术搜索
25. w·h·李,j . Wu Perrie, y,他“风速检索从hybrid-pol紧凑极化合成孔径雷达图像,”IEEE海洋工程》杂志上,没有。99年1 - 12,2017页。视图:出版商的网站|谷歌学术搜索
26. r . Shirvany m .夏伯特,j . y . Tourneret”船和漏油检测使用线性极化的程度和混合动力/ dual-pol SAR紧凑,”IEEE选择杂志的主题应用地球观测和遥感,5卷,不。3、885 - 892年,2012页。视图:出版商的网站|谷歌学术搜索
27. r . Touzi和f•夏博诺”,要求hybrid-compact SAR校准,”2014年IEEE地球科学和遥感研讨会魁北克市,页1109 - 1112年,QC,加拿大,2014年7月。视图:出版商的网站|谷歌学术搜索
28. s . Boularbah m . Ouarzeddine, a . Belhadj-Aissa”调查能力的紧凑偏振测定模式重建完整的偏振测定模式使用RADARSAT2数据,”先进的电磁学,1卷,不。1,p。19日,2012。视图:出版商的网站|谷歌学术搜索
29. 张x, j . Zhang m . Liu和j .孟”评估c波段紧凑的偏振测定SAR海冰分类,“Oceanologica学报,35卷,不。5,79 - 88年,2016页。视图:出版商的网站|谷歌学术搜索
30. p c . Dubois-Fernandez j . c . Souyris s Angelliaume f . Garestier,“小型偏振测定选择星载SAR在低频率,“IEEE地球科学和遥感,46卷,不。10日,3208 - 3222年,2008页。视图:出版商的网站|谷歌学术搜索
31. j . j . van Zyl h·a·Zebker, c . Elachi”成像雷达极化签名:理论和观察,”无线电科学,22卷,不。4、529 - 543年,1987页。视图:出版商的网站|谷歌学术搜索
32. js。李·d·l·舒勒,t·l·安斯沃思“偏振SAR数据地形方位补偿斜率变化,“IEEE地球科学和遥感,38卷,不。5,2153 - 2163年,2000页。视图:出版商的网站|谷歌学术搜索
33. js。李·d·l·舒勒,t·l·安斯沃思e . Krogager d . Kasilingam, W.-M。不过,“估计雷达极化取向变化引起的地形斜坡上,“IEEE地球科学和遥感,40卷,不。1,30-41,2002页。视图:出版商的网站|谷歌学术搜索
34. m . e .北部、t·l·安斯沃思j·s·李,和n·j·s·斯泰西,”比较紧凑的偏振合成孔径雷达模式”,IEEE地球科学和遥感卷,47号1,第188 - 174页,2009。视图:出版商的网站|谷歌学术搜索

版权

版权©2018王小陈et al。这是一个开放分布式下文章知识共享归属许可,它允许无限制的使用、分配和复制在任何媒介,提供最初的工作是正确引用。