Name: 使用单点光曲线绘制的地球样外行星的表面图 (视频) | JoVE Journal
Uploaded: 2020-05-10T12:30:00.000+00:00
Duration: 6 min 48 s
Description: 学习利用单点光曲线绘制地球样系外行星表面图的方法，该技术通过地球作为代理模型重建二维表面地图，为评估系外行星宜居性提供关键空间分辨能力。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

本工作部分由美国国家航空航天局（NASA）合同支持的加州理工学院喷气推进实验室资助。YLY 感谢华盛顿大学虚拟行星实验室提供的支持。

1. 编程设置
<ol>
	<li>为附带代码设置编程环境。需要使用 Linux 操作系统的计算机，因为 healpy 软件包在 Windows 系统上不可用。该代码的计算开销不大，普通个人计算机即可运行本实验方案。</li>
	<li>按照说明（https://docs.anaconda.com/anaconda/install/linux/）在系统中安装 Python 3.7 版本的 Anaconda，然后在终端中使用以下命令设置编程环境： 
	$ conda create --name myenv python=3.7 
	$ conda activate myenv 
	$ conda install anaconda 
	$ conda install healpy 
	注意：这些步骤所需时间可能因硬件性能和网络速度而异，耗时可达数十分钟。前两行命令中的环境名称“myenv”可替换为任意其他字符串。</li>
</ol>
2. 从观测中获取多波长光变曲线和观测几何结构
<ol>
	<li>在观测几何中，为每个对应的时间帧包含次恒星点和次观测点的经度和纬度。 
	若要使用以下附带的代码，请确保这两个文件的格式与 LightCurve.csv 和 Geometry.csv 相同。</li>
	<li>运行 PlotTimeSeries.py 以可视化数据并检查其质量。将生成两张图：LightCurve.png 和 Geometry.png（补充图1-2）。本步骤及后续绘图代码中的参数在应用于不同观测数据时可能需要调整。 
	$ python PlotTimeSeries.py LightCurve 
	$ python PlotTimeSeries.py Geometry</li>
</ol>
3. 从光变曲线中提取表面信息
<ol>
	<li>对系外行星的时间分辨多波长反照率光变曲线进行中心化处理，并通过每个波长对应的标准化偏差进行归一化。这使得每个通道具有相等的重要性。 
	<img alt="标准分数公式，R't,k=(Rt,k-μk)/σk；数据归一化与分析方程。" src="/files/ftp_upload/60951/60951equ01.jpg" data-altupdatedat="1788250015000" data-alt="方程 1"> 
	其中 R’t,k 和 Rt,k 分别表示在第 t 个时间步和第 k 个波长处的缩放后反照率和观测反照率；μk 和 σk 分别表示第 k 个波长处反照率时间序列的均值和标准差。
	<ol>
		<li>运行 Normalize.py 对光变曲线 Rt,k 进行归一化处理。输出结果保存在 NormalizedLightCurve.csv 中。 
		$ python Normalize.py</li>
	</ol>
	</li>
	<li>运行 PlotTimeSeries.py 以可视化归一化后的光变曲线。将生成图像 NormalizedLightCurve.png（补充图 3）。 
	$ python PlotTimeSeries.py NormalizedLightCurve</li>
	<li>对缩放后的反照率光变曲线应用奇异值分解（SVD），以提取主要的主成分（PCs）及其对应的时间序列。 
	<img alt="矩阵分解方程 R'=UΣV^T，数学表达式，用于教学目的。" src="/files/ftp_upload/60951/60951equ02.jpg" data-altupdatedat="1788250015000" data-alt="方程 2"> 
	等式左侧，T 和 K 分别为总时间步数和观测波长数；R’ 为缩放后反照率观测值的矩阵，其 (t,k) 元素为 R’t,k。等式右侧，V 的列为 PCs，即定义 SVD 投影空间的正交单位向量；Σ 为对角矩阵，其 (k,k) 元素为沿由 V 第 k 列定义的第 k 轴方向上缩放光变曲线的标准差；U 的列为 V 中各 PC 对应的时间序列。
	<ol>
		<li>运行 SingularValueDecomposition.py 对 R’ 进行分解。结果中的 U、Σ 和 VT 分别保存在输出文件 U.csv、SingularValue.csv 和 V_T.csv 中。 
		$ python SingularValueDecomposition.py</li>
	</ol>
	</li>
	<li>使用 PlotTimeSeries.py 和 PlotSVD.py 可视化 SVD 结果。将生成三幅图像：U.png、Sigma.png 和 V_T.png（补充图 4–6）。 
	$ python PlotTimeSeries.py U 
	$ python PlotSVD.py</li>
	<li>分析各主成分（PCs）的贡献及其对应的时间序列，以确定包含地表信息的主成分。
	<ol>
		<li>比较 Σ 对角线上的奇异值。一个类似地球且部分被云覆盖的系外行星预期具有两个相近的主要奇异值。 
		注意：Σ 可能包含少于或多于两个主要奇异值，具体情况将在下文讨论。</li>
		<li>比较两个主要 PC 的时间序列模式。包含地表信息的 PC 通常具有比另一个更规则的形状。由于经度不对称性以及地表在连续两天中以微小变化重现，其对应的时间序列往往表现出近似恒定的每日变化。</li>
		<li>使用 Lomb-Scargle 周期图9,10 计算两个主要 PC 的周期性，以确认 PC 的选择。包含地表信息的 PC 在功率密度谱中倾向于在自转周期处出现更高的峰值。</li>
		<li>运行 Periodogram.py 获取每个 PC 时间序列的功率谱。功率谱保存在 Periodogram.csv 中。 
		$ python Periodogram.py</li>
		<li>运行 PlotPeriodogram.py 以可视化这些周期图并确认 PC 的选择。将生成图像 Periodogram.png（补充图 7）。当前绘图代码添加了虚线，表示年周期、半年周期、日周期和半日周期作为参考，当应用于其他观测时可能需要调整。 
		$ python PlotPeriodogram.py</li>
		<li>选择包含地表信息的主成分 vj 及其对应的时间序列 uj。 
		<img alt="向量提取的矩阵方程；数学公式；用于教学目的。" src="/files/ftp_upload/60951/60951equ03v2.jpg" data-altupdatedat="1788250015000" data-alt="方程 3"> 
		<img alt="矩阵方程 $u_j = U[:,j]$，数学公式，线性代数概念。" src="/files/ftp_upload/60951/60951equ04v2.jpg" data-altupdatedat="1788250015000" data-alt="方程 4"> 
		其中 V[:,j] 和 U[:,j] 分别为 V 和 U 的第 j 列；j 为步骤 3.3 推断出的包含地表信息的 PC 的索引。</li>
	</ol>
	</li>
</ol>
4. 构建行星表面地图
<ol>
	<li>使用分层等面积等纬度像素化（Hierarchical Equal Area iso-Latitude Pixelization, HEALPix）11 方法对反演图进行像素化。该方法将行星的球面划分为面积相等且分布均匀的像素。记第 p 个像素的未知值为 xp。
	<ol>
		<li>运行 HEALPixRandom.py 以可视化像素化方法。将生成一幅图像 HEALPixRandom.png（补充图 8）。可修改第 17 行的参数 Nside 以实现不同分辨率。此步骤可能需要数秒至数分钟，具体耗时取决于分辨率。 
		$ python HEALPixRandom.py</li>
	</ol>
	</li>
	<li>利用观测几何关系，计算在第 t 个时间步长时第 p 个像素的权重 wt,p。 
	<img alt="静态平衡公式，wt,p = {ct cos(αt,p) cos(βt,p), 0}，条件方程示意图。" src="/files/ftp_upload/60951/60951equ05.jpg" data-altupdatedat="1788250015000" data-alt="方程 5"> 
	其中 αt,p、βt,p 分别为第 t 个时间步长时第 p 个像素处的太阳天顶角和航天器天顶角；ct 是第 t 次观测的归一化因子，使得每个时间步长的总权重之和为 1。 
	注：本步骤假设几何信息已知，或可通过其他分析获得，相关内容将在下文讨论。
	<ol>
		<li>运行 ComputeWeight.py 以计算 wt,p。修改第 23 行的 Nside 值，可获得不同分辨率的反演图。由于数据量较大，输出结果保存为 W.npz 文件。 
		$ python ComputeWeight.py</li>
	</ol>
	</li>
	<li>使用 PlotWeight.py 可视化这些权重。将在文件夹 Weight 中生成多个图像，每个时间步长对应一幅图。合并这些图像可得到 补充视频 1，展示各像素权重随时间的变化情况。由于可视化数量庞大，此步骤可能需要数小时才能完成。 
	$ python PlotWeight.py</li>
	<li>结合几何信息与观测数据，构建线性回归问题。 
	<img alt="矩阵方程 W[T×P]x[P×1]=u[T×1]，数学公式，线性代数概念。" src="/files/ftp_upload/60951/60951equ06.jpg" data-altupdatedat="1788250015000" data-alt="方程 6"> 
	其中 P 为反演像素的总数；W 为权重矩阵，其第 (t,p) 个元素为 wt,p；x 为待求解量，其第 p 个元素为 xp。 
	采用 L-2 范数正则化求解该线性回归问题。 
	<img alt="用于实验分析中数据拟合的线性回归方程。" src="/files/ftp_upload/60951/60951equ07.jpg" data-altupdatedat="1788250015000" data-alt="方程 7"> 
	其中 I 为单位矩阵，λ 为正则化参数。 
	注：当 T~104 且 P~3*103 时，λ 取 10-3 是一个合理值。应通过比较正则化平方误差 e 中两项的大小来调整 λ，如下所示。 
	<img alt="含范数、λ 的优化方程；数学表达式。" src="/files/ftp_upload/60951/60951equ08.jpg" data-altupdatedat="1788250015000" data-alt="方程 8">
	<ol>
		<li>运行 LinearRegression.py 求解该线性回归问题。x 的结果保存在文件 PixelValue.csv 中。修改第 16 行的 λ 值可调整正则化强度。 
		$ python LinearRegression.py</li>
	</ol>
	</li>
	<li>根据 HEALPix 的映射规则，将 x 转换为二维表面图。
	<ol>
		<li>运行 PlotMap.py，使用不同正则化参数构建反演图。当前设置下将生成三幅图像：Map_-2.png、Map_-3.png 和 Map_-4.png（补充图 9）。像素索引与其在图上位置的对应关系详见 HEALPix 文档11。此步骤耗时约数十秒。 
		$ python PolotMap.py</li>
	</ol>
	</li>
</ol>
5. 估算反演图的不确定性
<ol>
	<li>在第 4.3 步中，将线性回归问题重写为以 x 的“真实值”z 和观测噪声 ε 表示的形式。<img alt="矩阵方程 $W_{[T \times P]}Z_{[P \times 1]} + \epsilon_{[T \times 1]} = u_{t[T \times 1]}$。" src="/files/ftp_upload/60951/60951equ09.jpg" data-altupdatedat="1788250015000" data-alt="方程 9">
	<ol>
		<li>假设 ε 服从高斯分布 N (0, σ2I[T*T])，并估计其协方差。当反演得到的图谱固定时，T-P 是观测值 uj 的自由度。 
		<img alt="静态平衡；方差公式 σ²=(uj-Wx)ᵀ(uj-Wx)/(T-P)；数学方程。" src="/files/ftp_upload/60951/60951equ10.jpg" data-altupdatedat="1788250015000" data-alt="方程 10"></li>
		<li>合并第 4.4 步和第 5.1 步中的方程，得到关于 x 的高斯向量。 
		<img alt="正则化方程，x = (WᵀW + λI)⁻¹(WᵀW)z + (WᵀW + λI)⁻¹Wᵀε，数学表达式。" src="/files/ftp_upload/60951/60951equ11.jpg" data-altupdatedat="1788250015000" data-alt="方程 11"></li>
		<li>计算 x 的期望值和协方差矩阵。 
		<img alt="正则化线性回归方程；统计分析方法；E[x] 和 Cov[x] 的公式。" src="/files/ftp_upload/60951/60951equ12.jpg" data-altupdatedat="1788250015000" data-alt="方程 12"></li>
		<li>将 x 中每个元素的不确定性表示为 Cov[x] 对角线上对应元素的平方根。 
		<img alt="静态平衡方程 ep=√Diag[Cov[x]]p 数学公式分析图示" src="/files/ftp_upload/60951/60951equ13.jpg" data-altupdatedat="1788250015000" data-alt="方程 13"> 
		其中 ep 是 xp 的不确定性；Diag[Cov[x]]p 是 Cov[x] 对角线上的第 p 个元素。</li>
		<li>运行 Covariance.py 计算 x 的协方差矩阵。由于数据量较大，结果将保存为 Covariance.npz。该步骤所需时间从数十秒到数分钟不等，具体取决于 W 的大小。 
		$ python Covariance.py</li>
	</ol>
	</li>
	<li>根据 HEALPix 的映射规则，将 ep 转换为反演得到的二维图谱。
	<ol>
		<li>运行 PlotCovariance.py 以可视化 Cov[x] 并将不确定性 ep 映射到反演图谱上。将生成两张图像：Covariance.png 和 Uncertainty.png（补充图 10-11）。 
		$ python PlotCovariance.py</li>
	</ol>
	</li>
</ol>

<ol><li>Schwieterman, E. W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Exoplanet+Biosignatures%3A+A+Review+of+Remotely+Detectable+Signs+of+Life%5BTitle%5D)+AND+%22Astrobiology%22%5BJournal%5D)">Exoplanet Biosignatures: A Review of Remotely Detectable Signs of Life.</a> Astrobiology. 18 (6), 663-708 (2018).</li>
<li>Campbell, B., Walker, G. A. H., Yang, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+Search+for+Substellar+Companions+to+Solar-type+Stars%5BTitle%5D)+AND+%22The+Astrophysical+Journal%22%5BJournal%5D)">A Search for Substellar Companions to Solar-type Stars.</a> The Astrophysical Journal. 331, 902(1988).</li>
<li>NASA. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=%22NASA+Exoplanet+Archive+%282019%29+Confirmed+Planets+Table%22">NASA Exoplanet Archive (2019) Confirmed Planets Table</a>. , (2019).</li>
<li>Gillon, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Seven+temperate+terrestrial+planets+around+the+nearby+ultracool+dwarf+star+TRAPPIST-1%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Seven temperate terrestrial planets around the nearby ultracool dwarf star TRAPPIST-1.</a> Nature. 542 (7642), 456-460 (2017).</li>
<li>Kawahara, H., Fujii, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Global+Mapping+of+Earth-like+Exoplanets+from+Scattered+Light+Curves%5BTitle%5D)+AND+%22The+Astrophysical+Journal%22%5BJournal%5D)">Global Mapping of Earth-like Exoplanets from Scattered Light Curves.</a> The Astrophysical Journal. 720 (2), 1333(2010).</li>
<li>Fujii, Y., Kawahara, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Mapping+Earth+Analogs+from+Photometric+Variability%3A+Spin-Orbit+Tomography+for+Planets+in+Inclined+Orbits%5BTitle%5D)+AND+%22The+Astrophysical+Journal%22%5BJournal%5D)">Mapping Earth Analogs from Photometric Variability: Spin-Orbit Tomography for Planets in Inclined Orbits.</a> The Astrophysical Journal. 755 (2), 101(2012).</li>
<li>Cowan, N. B., Fujii, Y. Mapping Exoplanets. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aCowan%2C+NB+%22Handbook+of+Exoplanets%22">Handbook of Exoplanets</a>. , Springer, Cham. (2018).</li>
<li>Farr, B., Farr, W. M., Cowan, N. B., Haggard, H. M., Robinson, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((exocartographer%3A+A+Bayesian+Framework+for+Mapping+Exoplanets+in+Reflected+Light%5BTitle%5D)+AND+%22The+Astronomical+Journal%22%5BJournal%5D)">exocartographer: A Bayesian Framework for Mapping Exoplanets in Reflected Light.</a> The Astronomical Journal. 156 (4), 146(2018).</li>
<li>Lomb, N. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Least-Squares+Frequency+Analysis+of+Unequally+Spaced+Data%5BTitle%5D)+AND+%22Astrophysics+and+Space+Science%22%5BJournal%5D)">Least-Squares Frequency Analysis of Unequally Spaced Data.</a> Astrophysics and Space Science. 39 (2), 447(1976).</li>
<li>Scargle, J. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Studies+in+astronomical+time+series+analysis.+II.+Statistical+aspects+of+spectral+analysis+of+unevenly+spaced+data%5BTitle%5D)+AND+%22The+Astrophysical+Journal%22%5BJournal%5D)">Studies in astronomical time series analysis. II. Statistical aspects of spectral analysis of unevenly spaced data.</a> The Astrophysical Journal. 263, 835(1982).</li>
<li>Górski, K. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((HEALPix%3A+A+Framework+for+High-Resolution+Discretization+and+Fast+Analysis+of+Data+Distributed+on+the+Sphere%5BTitle%5D)+AND+%22The+Astrophysical+Journal%22%5BJournal%5D)">HEALPix: A Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere.</a> The Astrophysical Journal. 622 (2), 759(2005).</li>
<li>Jiang, J. H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Using+Deep+Space+Climate+Observatory+Measurements+to+Study+the+Earth+as+an+Exoplanet%5BTitle%5D)+AND+%22The+Astronomical+Journal%22%5BJournal%5D)">Using Deep Space Climate Observatory Measurements to Study the Earth as an Exoplanet.</a> The Astronomical Journal. 156 (1), 26(2018).</li>
<li>Fan, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Earth+as+an+Exoplanet%3A+A+Two-dimensional+Alien+Map%5BTitle%5D)+AND+%22The+Astrophysical+Journal+Letters%22%5BJournal%5D)">Earth as an Exoplanet: A Two-dimensional Alien Map.</a> The Astrophysical Journal Letters. 882 (1), 1(2019).</li>
<li>Cowan, N. B., Strait, T. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Determining+Reflectance+Spectra+of+Surfaces+and+Clouds+on+Exoplanets%5BTitle%5D)+AND+%22The+Astrophysical+Journal+Letters%22%5BJournal%5D)">Determining Reflectance Spectra of Surfaces and Clouds on Exoplanets.</a> The Astrophysical Journal Letters. 765 (1), 17(2013).</li>
<li>Fujii, Y., Lustig-Yaeger, J., Cowan, N. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Rotational+Spectral+Unmixing+of+Exoplanets%3A+Degeneracies+between+Surface+Colors+and+Geography%5BTitle%5D)+AND+%22The+Astronomical+Journal%22%5BJournal%5D)">Rotational Spectral Unmixing of Exoplanets: Degeneracies between Surface Colors and Geography.</a> The Astronomical Journal. 154 (5), 189(2017).</li>
<li>Kawahara, H., Fujii, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Mapping+Clouds+and+Terrain+of+Earth-like+Planets+from+Photomertic+Variability%3A+Demonstration+with+Planets+in+Face-on+Orbits%5BTitle%5D)+AND+%22The+Astrophysical+Journal+Letters%22%5BJournal%5D)">Mapping Clouds and Terrain of Earth-like Planets from Photomertic Variability: Demonstration with Planets in Face-on Orbits.</a> The Astrophysical Journal Letters. 739 (2), 62(2011).</li>
<li>Kawahara, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Frequency+Modulation+of+Directly+Imaged+Exoplanets%3A+Geometric+Effect+as+a+Probe+of+Planetary+Obliquity%5BTitle%5D)+AND+%22The+Astrophysical+Journal%22%5BJournal%5D)">Frequency Modulation of Directly Imaged Exoplanets: Geometric Effect as a Probe of Planetary Obliquity.</a> The Astrophysical Journal. 822 (2), 112(2016).</li>
<li>Schwartz, J. C., Sekowski, C., Haggard, H. M., Pall ́e, E., Cowan, N. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Inferring+planetary+obliquity+using+rotational+and+orbital+photometry%5BTitle%5D)+AND+%22Monthly+Notices+of+the+Royal+Astronomical+Society%22%5BJournal%5D)">Inferring planetary obliquity using rotational and orbital photometry.</a> Monthly Notices of the Royal Astronomical Society. 457 (1), 926-938 (2016).</li>
</ol>

作者无任何利益冲突需要披露。

该方案的一个关键要求是从光变曲线中提取表面信息的可行性，而这取决于云层覆盖情况。在步骤3.5.1中，不同系外行星的主成分（PC）相对值可能有所不同。以地球为例，前两个主成分主导了光变曲线的变化，分别对应于与地表无关的云层和地表信息（Fan et al. 2019）13。它们的奇异值大小相近，因此可通过步骤3.5.2和3.5.3将地表信息分离出来。对于未来对系外行星的观测，在极端情况下，若行星完全被云层覆盖或完全无云，则在步骤3.3的奇异值分解（SVD）中仅会出现一个主导的主成分。此时必须进行光谱分析，以解释该主成分的物理含义，因为云层和地表的成分不同。如果主导主成分对应的是地表信息，则仍可继续执行步骤4和5；若其对应的是云层，则可得出结论：地表信息被云层遮挡，因此无法利用给定波长下的光变曲线提取地表信息，此时地表成像不可行。此外，也可能存在第三个甚至第四个具有可比性的主导主成分，其可能对应另一层云或大尺度水文过程，但只要能够提取地表信息，这些情况不会影响本方法后续步骤的有效性。
如 Cowan & Strait (2013)<sup ...

确定宜居世界是天体生物学的最终目标之一1。自首次探测到系外行星以来2，迄今为止已确认了超过4000颗系外行星3，其中包括若干类地行星（例如TRAPPIST-1e）4。这些行星的轨道和行星特性与地球相似，因此可能具备宜居性。在此背景下，基于有限的观测数据评估其宜居性至关重要。根据地球生命的认知，地质系统和气候系统对宜居性具有关键作用，因而可作为生物标志物。原则上，即使行星在空间上无法被分辨为超过单个像素点，这些系统的特征仍可能通过远距离观测获得。在此情况下，通过单像素光变曲线识别地质特征和气候系统，对于评估系外行星的宜居性至关重要。对这些系外行星进行表面测绘已成为紧迫任务。
尽管观测几何与光谱特征之间存在复杂关联，但系外行星表面的信息仍包含在其随时间变化的单点光变曲线中，这些光变曲线可以从远处获取，并通过足够多的观测数据推导得出。然而，由于云层的影响，对可能宜居的类地系外行星进行二维（2D）表面成像仍具挑战性。研究人员已基于模拟光变曲线和已知光谱开发并测试了多种重建二维表面图的方法5,6,7,8，但这些方法尚未应用于真实的观测数据。此外，在当前及近期将来的系外行星观测分析中，当行星表面成分尚未被较好约束时，采用特征光谱的假设可能会引发争议。
本文展示了一种针对类地系外行星的表面映射技术。我们采用奇异值分解（SVD）方法，无需假设任何特定光谱特征，即可评估并分离多波长光变曲线中来自不同来源的信息。结合观测几何构型，我们利用时间分辨但空间上混合的表面信息，实现了表面图的重建。为演示该方法，本文分析了深空气候观测站/地球多色成像相机（DSCOVR/EPIC；www.nesdis.noaa.gov/DSCOVR/spacecraft.html）获取的为期两年的地球多波长单点观测数据。由于目前可获得的系外行星观测数据尚不充分，我们以地球作为系外行星的代理样本，用于评估该方法的有效性。本文附带了示例代码，该代码基于 Python 3.7 环境开发，使用了 anaconda 和 healpy 软件包，但该方案的数学原理也可在其他编程环境中实现（例如 IDL 或 MATLAB）。

<table><tbody><tr><td>配备 anaconda 和 healpy 包的 Python 3.7</td><td></td><td></td><td>其他编程环境（例如 IDL 或 MATLAB）也可使用。</td></tr></tbody></table>

利用单点光变曲线对类地系外行星进行表面测绘

Surface Mapping of Earth-like Exoplanets using Single Point Light Curves

通过单点观测实现对系外行星特征的空间分辨，对于评估系外行星的潜在宜居性至关重要。本实验方案的最终目标是确定这些行星世界是否具有地质特征和/或气候系统。我们提出一种方法，从多波段单点光变曲线中提取信息并反演表面图。该方法利用奇异值分解（SVD）分离导致光变曲线变化的各个来源，并推断部分多云气候系统的存在。通过对SVD所得时间序列的分析，可在无需假设任何光谱特性的前提下，推断主成分（PCs）的物理归属。结合观测视角几何关系，若发现某一主成分包含表面信息，则有可能重建表面图。由于像素几何结构与光谱信息的卷积所导致的退化性问题会影响重建表面图的质量，因此需要引入正则化方法加以约束。为演示该方案，本文分析了地球的多波段光变曲线（作为类比系外行星的代理样本），并将结果与真实情况对比，以展示该方案的性能与局限性。本研究为未来推广至系外行星应用提供了基准参考。

我们利用地球的多波长单点光变曲线来演示本方案，并将结果与真实情况对比以评估地表反演的质量。所使用的观测数据来自位于地球与太阳之间第一拉格朗日点（L1）附近的 DSCOVR/EPIC 卫星，该卫星对地球被阳光照射的一面在十个波长下进行成像。本次演示使用了两年（2016 年和 2017 年）的观测数据，与 Jiang 等（2018）12 和 Fan 等（2019）13 所用数据相同，更多关于观测的细节可在上述文献中找到。2017 年 2 月 8 日世界时 9:27 的一次观测样本如 图 1 所示。地球图像被积分成单个点，以模拟“外星人”或遥远观测者所获得的光变曲线观测，这些观测者无法将地球的空间分辨率提高到超过一个像素。因此，生成了包含约 10,000 个时间步长的多波长单点系外行星光变曲线，作为本方案的输入数据。
在完成步骤3后，我们在多波长光变曲线中发现两个主要的主成分（PC），其中第二个主成分（PC2）包含地表信息。根据步骤3...

该方案从系外行星的光变曲线中提取信息，并构建其表面图。通过使用地球的光变曲线（作为类比系外行星）来演示这一方法。

从单点观测中实现对系外行星特征的空间分辨，对于评估系外行星的潜在宜居性至关重要。最终……

该方案可用于从单点观测中空间解析系外行星的特征，对于评估系外行星潜在的宜居性至关重要。该技术可用于重建类地系外行星的二维表面图，并且是首个利用地球作为代理通过真实观测进行测试的技术。该技术的数学原理简单明了，可轻松调整以适用于其他观测。使用者无需严格遵循编码脚本。该技术的可视化演示至关重要，因为一图胜千言。在配置好附加代码的编程环境后，输入命令以在系统中安装带有 Python 3.7 的 Anaconda。编程环境设置完成后，获取多波段光变曲线，查看观测数据的几何构型，并运行 plot time series.py 命令以可视化数据并检查其质量。然后输入命令以生成几何图形。为提取光变曲线表面信息，请运行 normalize.py 命令。输出结果将保存在 normalized light curve.csv 文件中。输入命令以可视化归一化的光变曲线，将生成一幅归一化光变曲线图。输入命令以分解归一化的光变曲线。生成的时间序列、奇异值和主成分将保存为相应的 csv 格式输出文件。使用命令可视化奇异值分解结果。将生成奇异值和主成分的图表。分析主成分的贡献度及其对应的时间序列，以确定包含表面信息的主成分，并比较奇异值矩阵对角线上的奇异值。类地部分多云系外行星预期具有两个相近的主要奇异值。为了确认主成分的选择，请输入命令以获取每个主成分时间序列的功率谱。功率谱将保存在 periodogram.csv 文件中。输入命令以可视化周期图，并确认主成分的选择。将生成一个周期图。当前的绘图代码会添加虚线，分别表示年周期、半年周期、日周期和半日周期，以供参考。请选择包含地表信息的主成分及其对应的时间序列。为了构建行星表面图，使用 HEALPix random 命令来可视化像素化方法。将生成一幅 HEALPix 随机图像。第 17 行中的参数 n_subside 可以根据不同的分辨率进行修改。为了确定每个像素的权重，请输入命令。由于文件较大，输出结果将保存为 w.npz。根据所获取地图的其他分辨率，将第 23 行的 n subside 值更改为适当的数值。使用 plot weight.py 命令可视化权重。将在权重文件夹中生成多个图像。合并这些图像可以展示每个像素的权重随时间的变化情况。请使用线性回归方法。使用 py 命令求解线性回归问题。像素值的结果将保存在 pixel value.csv 文件中。可根据需要在线16处修改Lambda的取值，以实现不同强度的正则化。然后运行plot_map.py命令，使用不同的正则化参数构建重构图谱。将生成三幅图。每个图中像素索引与其位置之间的关系在HEALPix文档中有详细说明。为计算每个像素的协方差矩阵，请运行协方差计算。py 命令。由于文件较大，结果将保存在 covariance.npz 中。为了可视化协方差矩阵并将不确定性映射到反演结果图中，请运行 plot covariance.py 命令。将生成一幅协方差图和一幅不确定性图。此处展示了2017年2月8日协调世界时927的地球多波长观测样本。此处可观察到多波长光变曲线的两个主要主成分的时间序列。第二主成分的时间序列表现出更规则的形态，其功率谱中每日变化近似恒定，且昼夜周期比第一主成分更显著。随后可以构建该代理系外行星的表面分布图，图中每个像素点的数值对应第二主成分的值。与地球的真实情况相比，在将地表信息与云层分离后，重建的分布图成功恢复了所有主要大陆。由于南半球海洋上空的云层覆盖较多，其重建结果较北半球更差。每个像素值的不确定性约为反演图中数值的10%，表明表面映射质量良好，结果为正向。将本方案应用于未来分析的关键要求是确认可以从光变曲线中提取表面信息。该技术在系外行星表面映射中具有基准作用，未来可通过其他去补偿方法及正则化方法结合新的观测数据进一步改进。

surface-mapping-earth-like-exoplanets-using-single-point-light

Research

科研

JoVE Journal

Engineering

工程

使用单点光曲线绘制的地球样外行星的表面图

来自天文物体的光必须穿过地球的湍流大气，才能被地面望远镜成像。为了实现接近最大理论角分辨率的直接成像，必须采用先进的技术，例如Robo-AO自适应光学系统所使用的技术。...

使用Robo-AO将可见宇宙聚焦清晰

在本视频方案中，我们使用经改造的双光子显微镜结合轨道追踪方法，对活细胞内荧光标记的溶酶体进行高速、三维的动态追踪。...

在改进的双光子显微镜中进行三维轨道追踪：用于细胞内囊泡追踪的应用

本文介绍了密集排列颗粒的离散随机介质中光的多重散射的数值与实验方法。这些方法被用于解释灶神星（4号小行星）和楚留莫夫－格拉希门克彗星（67P）的观测结果。...

行星风化层中光的散射与吸收

本方案描述了 ElectroMap 的设置与使用方法。ElectroMap 是一个基于 MATLAB 的开源软件平台，用于分析心脏光学标测数据。ElectroMap...

使用 ElectroMap 进行光学标测数据的高通量分析

冰冻圈中的碳通量目前几乎尚未得到评估，但在气候变化方面至关重要。本文展示了一种新型原型装置，该装置基于激光诱导荧光发射（L.I.F.E.）技术，在原位条件下提供高光谱和高空间分辨率的数据，用于捕捉表生冰川环境中的光合潜力。...

激光诱导荧光发射（L.I.F.E.）作为一种新型非侵入式工具用于冰冻圈生境中生物标志物的原位测量

本文介绍了一种用于测试月球表面大规模射电阵列成像能力的模拟框架。讨论了主要的噪声组分，并详细阐述了软件处理流程，包括如何针对新的科学应用对其进行定制化修改。...

模拟月球表面大型射电阵列成像

本方案涉及转染cAMP传感器和一种光遗传学蛋白bPAC-nLuc，以精确追踪其在细胞内的分布及对光刺激的响应。利用点扫描系统构建cAMP响应图谱的创新方法，有望推动光遗传学蛋白在不同研究领域的应用进展。...

使用光刺激网格绘制光遗传学蛋白的细胞分布图

本文介绍了一种利用布里渊显微光谱法在三维基质中对活体球状体进行力学表征的实验方案。该全光学方法可实现微米级分辨率的球状体成像，并提供定量的力学特性信息。该技术可用于力学表型分析，例如，用于确定三维微环境中肿瘤球状体的病理状态。...

利用布里渊光谱技术对类球体进行力学图谱分析

本方案涵盖一种四步大范围成像调查方法，用于提取珊瑚礁群落的结构复杂性、群落组成以及种群人口统计学指标。该方案的每一步均优先确保所采集图像的质量以及对原始图像的集成访问。...

利用大面积成像测量水下环境的结构、组成及变化

白光显微干涉测量法是一种光学、非接触式且快速的表面形貌测量方法。本文展示了该方法在机械磨损分析中的应用，可用于分析摩擦学测试样品上的磨损痕迹；同时在材料科学中用于测定离子束溅射或激光烧蚀的体积和深度。...

通过白光干涉术表征表面改性：在离子溅射、激光烧蚀和摩擦学实验中的应用

学习利用单点光曲线绘制地球样系外行星表面图的方法，该技术通过地球作为代理模型重建二维表面地图，为评估系外行星宜居性提供关键空间分辨能力。

利用单点光变曲线对类地系外行星进行表面测绘

探索更多视频