USGS-ISIS处理LROC-NAC(EDR)数据，添加地理投影，转换为.tif格式

已于 2023-04-26 09:11:37 修改
阅读量1.7k 收藏 21
点赞数 4
文章标签： linux ubuntu python
于 2023-04-26 01:05:17 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_31721595/article/details/130365130
版权

目录

1.LROC-NAC数据下载

2.利用ISIS处理LROC-NAC(EDR)数据

2.1 将NAC-EDR转换为ISIS的cube格式

2.2 附加相机位置信息

2.3 辐射校正

2.4 回波校正

2.5 光度校正（可选）

2.6 创建一个地图文件mosaic.map

2.7 地图投影

2.8 镶嵌

2.9 配准（可选）

3.利用GDAL转换为地理参考格式（.tif）

参考资料

1.LROC-NAC数据下载

数据介绍：Lunar Reconnaissance Orbiter Camera (LROC) Narrow Angle Cameras (NACs)——月球轨道窄角相机，是月球勘测轨道器(LRO)上的一对高分辨率推扫式(pushbroom)相机。两台相机（NAC-Left，NAC-Right）略有重叠地平行扫描月面，并同时获取图像。典型的LROC NAC图像的像素分辨率在0.5-2.0 m/px之间，部分图像的空间分辨率可达0.25-0.3 m/px，是目前分辨率最高的月球影像数据集。一个典型的NAC成对图像大约有10,000个像素宽（东西向）和52,000个像素长（南北向），对应的正常覆盖范围大小为5x26公里到20x100公里左右。NAC数据集由NASA's Planetary Data System（PDS）提供，可供公众下载。

数据下载网址：

注意：获取的数据为.IMG格式，为了读取和提取图像的地理信息，需要将其转换为(.geotiff or .gtif)格式，该数据一般包括：EDR/CDR/RDR，我们下载EDR数据，CDR数据的处理比EDR更复杂。

  1. EDR：未经过处理的原始NAC数据；
  2. CDR：经过光度校正但无地理投影的NAC数据；
  3. RDR：经过几何校正和地理投影的NAC数据

打开第一个网址后如图所示，可根据NAC的（轨道号、相机角度、入射角和发射角、分辨率和图像ID）查询数据，Product ID（eg.M1172873803LE，M1172873803RE）下载对应的EDR数据。

2.利用ISIS处理LROC-NAC(EDR)数据

我们使用美国地质调查局United States Geological Survey (USGS)开发的成像仪和光谱仪集成软件Integrated Software for Imagers and Spectrometers (ISIS)处理数据。

ISIS的下载安装参考：深空探测处理软件USGS-ISIS的下载安装_qq_31721595的博客-CSDN博客

2.1 将NAC-EDR转换为ISIS的cube格式

（1）首先找到你所下载的数据文件的位置，然后打开终端，利用conda进入ISIS，如图：

（2）将M1172873803LE.IMG，M1172873803RE.IMG转换为ISIS的cube格式，命令如下：

lronac2isis from=M1172873803LE.IMG to=M1172873803L.raw.cub

2.2 附加相机位置信息

ISIS使用Navigation and Ancillary Information Facility SPICE system在观测过程中确定相机位置，以便将图像正确投影到地面上。该步只要在光度校正前执行即可，开始时执行此步骤便于立即查询图像的几何信息，或在进行额外处理前确定是否有所需要的kernels可用。命令如下：

spiceinit from=M1172873803L.raw.cub spksmithed=true web=true
spiceinit from=M1172873803R.raw.cub spksmithed=true web=true

对M1172873803R做同样处理，如图：

2.3 辐射校正

lronaccal from=M1172873803L.raw.cub to=M1172873803L.cal.cub
lronaccal from=M1172873803R.raw.cub to=M1172873803R.cal.cub

2.4 回波校正

lronacecho from=M1172873803L.cal.cub to=M1172873803L.echo.cub
lronacecho from=M1172873803R.cal.cub to=M1172873803R.echo.cub

注：回波校正后会生成伪影，利用trim去除。

如果对图像求和(在EDR标签中，查找“CROSSTRACK_SUMMING = 2”，或者检查图像宽度是否为2532像素)，将这些数字除以2:修剪量应该是23和13 ：

NAC-L：

trim from=M1172873803L.echo.cub to=M1172873803L.tr.cub left=46 right=26

NAC-R：

trim from=M1172873803R.echo.cub to=M1172873803R.tr.cub left=26 right=46

2.5 光度校正（可选）

这一步不是必须的

photomet from=M1172873803L.tr.cub to=M1172873803L.pho.cub frompvl=basicpho.pvl
photomet from=M1172873803R.tr.cub to=M1172873803R.pho.cub frompvl=basicpho.pvl

不过，若要进行光度校正首先需要新建一个basicpho.pvl文件，该文件内容为：

Object = NormalizationModel
	Group = Algorithm
 		Name = Mixed
 		Incmat=87.0
 		Albedo=0.05
 		Incref=0.0
 		Thresh=10E30
 	EndGroup
EndObject
Object = PhotometricModel
 	Group = Algorithm
 		Name = HapkeHen
 		Wh=0.278906221
 		Hh=0.0438
 		B0=1.0
 		Theta=10.0
 		HG1=-0.229315607
 		HG2=0.0
 	EndGroup
EndObject

2.6 创建一个地图文件mosaic.map

首先将所有想要组合成一个mosaic的图像的路径放入自建的文本文件"mosaic_level_1.lis"中，如图：

然后，利用mosrange命令 生成地图文件。此外，如果知道感兴趣区域的确切边界或中心点，还可用maptemplate，但这里用的mosrange命令，使用的是等矩形投影equirectangular，若纬度在70°以上或靠近极点的区域，利用polar stereographic。

注：因为所有发布的LRO产品使用的都是0-360的经度值，而ESRI程序（如ArcMap）经度值在-180-180，所以若要在ArcMap中使用地图投影产品，还需在后面添加命令，londom=180。

mosrange fromlist=mosaic_level_1.lis to=mosaic.map projection=equirectangular londom=180

2.7 地图投影

cam2map from=M1172873803L.pho.cub to=M1172873803L.map.cub map=mosaic.map warpalgorithm=forwardpatch patchsize=50
cam2map from=M1172873803R.pho.cub to=M1172873803R.map.cub map=mosaic.map warpalgorithm=forwardpatch patchsize=50

2.8 镶嵌

ls M1172873803?.map.cub > M1172873803.lis

automos fromlist=M1172873803.lis mos=M1172873803.mosaic.cub

运行完第一句命令后，生成M1172873803.lis文件，如图：

然后运行第二条命令，但是报错了，目前还不清楚是什么原因，如图： 

但我们可以将得到的 M1172873803L.map.cub 和 M1172873803R.map.cub 放在ArcMap中进行镶嵌，见第3节。

2.9 配准（可选）

有时一对NAC图像可能存在明显的偏移，需要使用corge来对齐两个图像。

map2map from=M1172873803R.map.cub to=M1172873803R.Lbounds.cub map=M1172873803L.map.cub matchmap=true

下一步，首先需要创建nac_pair_align.def文件，内容如下：

Object = AutoRegistration
 	Group = Algorithm
		Name = MaximumCorrelation
 		Tolerance = 0.7
 		SubpixelAccuracy = True
 	End_Group
 	
	Group = PatternChip
 		Samples = 20
 		Lines = 20
 	End_Group

# If you determine that the vertical or horizontal offset is >50 pixels, adjust the SearchChip 
# size to be at least `20+offset*2`, where 20 is the PatternChip size.
 
	Group = SearchChip
 		Samples = 120
 		Lines = 120
 	End_Group
 
	Group = SurfaceModel
 		DistanceTolerance = 1.5
 		WindowSize = 5
 	End_Group
End_Object
End

然后，运行如下命令，会得到 SampleAverage 和 LineAverage，用于下一步translate命令 ： 

coreg from=M1172873803R.Lbounds.cub match=M1172873803L.map.cub deffile=nac_pair_align.def

利用translate命令纠正偏移， 将$sampleaverage 和 $lineaverage换为上图corge输出的值，-0.58和-1.13.

translate from=M1172873803R.map.cub to=M1172873803R.map.shifted.cub strans=$sampleaverage ltrans=$lineaverage

然后再对M1172873803L.map.cub、M1172873803R.map.shifted.cub运行2.8中的镶嵌命令，但依然报错，所以我直接跳过2.8、2.9的命令，利用ArcMap拼接影像。

 最后需要的就是这两个数据文件：

3.利用GDAL转换为地理参考格式（.tif）

因为未在虚拟机中安装GDAL，我在本机电脑的Pycharm中利用Python进行转换，前提是也需要安装gdal包，或从osgeo中导入gdal包。

from osgeo import gdal

cubdata = r"D:\Moon_Data\processedData\M1096293859L.map.cub"   # 经上面步骤得到的.cub文件路径
tifdata = r"D:\Moon_Data\processedData\M1096293859L.tif"       # .tif文件路径

gdal.Translate(tifdata, cubdata)

最终得到下图的.tif文件，并在ArcMap中打开它们，可以看到这些.tif图像带有投影信息了： 

  

参考资料

为ISIS处理LROC-NAC的官方资料，上述过程均参考这些资料，详细步骤请参考如下资料！

  1. 链接：https://pan.baidu.com/s/1EwGj3xdbeES6-9LQVmh-2g 
    提取码：6e8i
  2. 链接：https://pan.baidu.com/s/1cyIBdrzbh4wBuKh15gHP9A 
    提取码：o7fz
qq_31721595
关注
  • 4
    点赞
  • 踩
  • 21
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
LROC2FGearthview：用于FlightGear地球视图的月亮纹理转换器
02-16
FlightGear Earthview的月球纹理生成器 此代码使您可以从梦幻般的公共领域LROC图像创建Moon纹理，RGB和。 LRO（C）代表NASA / GSFC /亚利桑那州立大学的月球侦察轨道卫星（相机）。 请查看网页或页面以获取更多详细信息，并直接访问数据。 跑步 您需要有效安装imagemagick（API6> =），可从（ ）获得的normalmap二进制文件和的有效安装。 然后启动bash脚本： $ ./moontextures moon heights 8k 该脚本改编自Chris Blues的脚本，该脚本生成了地球纹理，有关更多信息，请 学分 Chris_Blues（ ） LROC图片，NASA / GSFC /亚利桑那州（ ） 执照
isis安装指南（Liunx版本）+ ASP
qqwowo99的博客
01-14 4034
ISIS（Integrated Software for Imagers and Spectrometers ）是一个数字图像处理软件包，用于处理当前和过去 NASA 和国际行星任务收集的图像。 例如，从PDS下载的LRO月球影像数据源文件为IMG格式，需要使用ISIS进行预处理，并生成遥感图像常用TIF格式 Ames Stereo Pipeline（ASP）主要功能是从深空探测轨道器立体像对自动生成三维地形（DEM）等，“ISIS3+ASP” 这一流程已经成为深空探测领域构建立体像对的基本流程 环境要求：
CV LRO mission （LROC introduction）
dbzr15157的博客
12-02 209
美国观测轨道计划，卫星发射与2009年6月28日，是美国重返月球计划的一部分。空间器携带7中传感器：Lunar reconnaissance orbiter camera(LROC)Lunar Orbiter Laser altimeter(LOLA)Lunar Exploration Neutron Detector(LEND)Diviner Lunar Radiometer Ex...
【infiniband】SDR、DDR、QDR、FDR、EDR、HDR、NDR
数字人生
09-12 8270
- SDR（Single Data Rate） - DDR（Double Data Rate） - QDR（Quad Data Rate） - FDR（Fourteen Data Rate） - EDR（Enhanced Data Rate）
深空探测处理软件USGS-ISIS的下载安装
qq_31721595的博客
04-17 1342
深空探测处理软件USGS-ISIS下载安装教程
【DEM格式转换】转换国标DEM数据格式为通用格式
小黑先森的博客
05-10 2万+
使用ArcGIS工具转换国标.dem文件为通用.tif格式 *.dem是一种比较常见的DEM数据格式，其有两种文件组织方式，即NSDTF-DEM和USGS-DEM 1。 NSDTF-DEM NSDTF-DEM是一种明码的中国国家标准空间数据的交换格式，遵从国家标准GB/T 17798-2007《地理空间数据交换格式》中的格网数据组织规范。 USGS-DEM USGS-DEM是一种由美国...
Sentinel-2 哨兵二号数据（Level-1C）下载及预处理教程
热门推荐
C_hqpxd的博客
05-15 3万+
欧空局（ESA）仅发布了哨兵2号（S2）的L1C级多光谱数据（MSI），L1C级数据是经过几何精校正的正射影像，并没有进行辐射定标和大气校正。同时，ESA还对S2 L2A级数据进行了定义，L2A级数据主要包含经过辐射定标和大气校正的大气底层反射率数据（Bottom-of-Atmosphere corrected reflectance），但这个L2A数据需要用户根据需求自行生产，为此，ESA发布了专门生产L2A级数据的插件Sen2cor。 ...
usgs.rar_USGS_usgs-dem
09-19
【标题】"usgs.rar_USGS_usgs-dem" 指的是一个与美国地质调查局（United States Geological...通过解压并编译这个项目，用户不仅可以了解USGS的DEM数据格式，还能学习到如何使用C++处理GIS数据，是一个很好的实践案例。
dataRetrieval：此R包旨在直接从Web服务获取USGS或EPA水质样本数据，流数据和元数据。 参见：http：usgs-r.github.iodataRetrieval
02-05
4. **保存数据**：若需长期存储，可以使用R的`write.csv`或`write.table`函数将数据保存为CSV或其他格式的文件。 5. **可视化**：R中丰富的可视化包，如`ggplot2`，可以帮助我们将数据以图表形式展示出来，便于理解...
aringler-usgs-mlabBob.zip
08-08
"aringler-usgs-mlabBob.zip" 是一个压缩包文件，从标签 "matlab" 我们可以推断，这个文件很可能包含了与 MATLAB 相关的代码、数据或者项目。MATLAB 是一种专用于数值计算、数据分析和算法开发的高级编程环境。它...
美国USGS的数据说明---DEM数据格式读取等
05-23
美国USGS的数据说明---DEM数据格式读取等
Quake-Report:Android应用程序可从USGS网站获取实时数据
05-12
API通常提供JSON或XML格式的数据，应用解析这些数据并将其转化为易于处理的格式。 2. JSON解析：JSON（JavaScript Object Notation）是一种轻量级的数据交换格式，便于人阅读和编写，同时也方便机器解析和生成。在...
wdfn-blog:USGS国家水数据博客
04-28
USGS国家水数据博客 提交准则 适用/对更多受众有用。 内容提供了其他地方没有的其他信息。 与USGS员工或USGS数据用户相关 展示与USGS WMA或USGS WMA产品组合相关的数据库，应用程序，软件或通用做法。 提交更新 ...
NAC配置与管理——1
kaoa000的专栏——非淡泊无以明志，非宁静无以致远
05-29 9335
在计算机网络安全管理中，用户的网络接入控制（NAC）是必须充分考虑的一件大事，因为现在的网络安全隐患主要不是来自外网，而是内网。在用户接入控制方面，最直接、最有效的方法就是基于接入设备接口的各种用户认证方法，如802.1x认证、MAC认证和Portal认证。它们针对不同的用户需求和实际的网络环境提供的几种基于端口的实用接入控制方案。相对前面的AAA方案来说，此处基于接口的接入控制方法更直接，直接在...
ROC FROC SROC曲线
weixin_41043240的博客
05-10 1万+
1. ROC曲线 ROC（Receiver Operating Characteristic Curves）曲线是显示Classification模型真正率和假正率之间折中的一种图形化方法。 解读ROC图的一些概念定义:： 真正（True Positive , TP）被模型预测为正的正样本； 假负（False Negative , FN）被模型预测为负的正样本； 假正（False Pos...
ROC曲线、FROC曲线、SROC曲线和LROC曲线
Acmer_future_victor的博客
03-05 4539
该文章为转载文章，原文链接为https://blog.csdn.net/weixin_41043240/article/details/80265577 1. ROC曲线 ROC（Receiver Operating Characteristic Curves）曲线是显示Classification模型真正率和假正率之间折中的一种图形化方法。 解读ROC图的一些概念定义:：...
USGS官方的各种卫星数据产品的详细说明（很多英文缩写都可在这里找到详细信息，如ETM+和TM的意思，也可以下载）
fenghuayoushi的专栏
11-10 1万+
如题了：USGS官方的各种卫星数据产品的详细说明（很多英文缩写都可在这里找到详细信息，如ETM+和TM的意思，也可以下载） USGS是什么？ （United States Geological Survey，简称：USGS），美国地质调查局。很厉害的一个官方网站。想想，Landsat卫星数据的官方网站就在USGS旗下，就知道有多牛B了。美国政府内政部辖下的科学机构，是内政部里唯一一个纯
【持续集成_05课_Linux部署SonarQube及结合开发项目部署】
最新发布
weixin_42333261的博客
07-12 281
前置条件：sonarQube不能使用root账号进行启动，所以需要创建普通用户及。3）CMD上传qube文件-不能传到home路径下哦。2）添加用户、组名、密码。4）检查并解压上传的包。
USGS-ISIS的下载安装
10-07
USGS-ISIS是美国地质调查局（USGS）开发的一个用于处理行星地质图像和数据的软件包。要下载和安装USGS-ISIS，您可以按照以下步骤进行操作： 1. 访问USGS-ISIS的官方网站：https://isis.astrogeology.usgs.gov/7.0.0/Application/index.html。 2. 在网站上找到下载链接，并选择适合您操作系统的版本进行下载。 3. 下载完成后，按照下载页面上的说明，执行安装程序。 4. 在安装过程中，您可能需要选择安装路径和其他选项。请根据个人需求进行选择。 5. 完成安装后，您需要配置ISISROOT环境变量。打开编辑器，找到配置文件，并在文件末尾添加以下内容：`ISISROOT=/home/v/ISIS3/isis`（这里的isis路径根据实际情况进行修改）。 6. 保存文件并退出编辑器。 7. 打开终端或命令提示符，输入以下命令以应用配置变更：`source ~/.bashrc`（根据您的操作系统和Shell可能会有所不同）。 8. 接下来，您可以通过命令行输入`isis3Startup.sh`来启动USGS-ISIS。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值