利用ctx图像制作马赛克vrt的基本cli
circ的Python项目详细描述
自2018年以来,一直在制作不错的ctx马赛克
circ通过读取ctx图像元数据来选择 与用户提供的边界框相交,然后执行一些基本过滤和 随机图像还原以生成覆盖bbox的合理最小图像列表 低发射角。
使用的图像来自asu mars space flight facility viewer.mars.asu.edu。
安装
目前,您必须已安装gdal和wget,并在您的路径中可用。 然后下载或克隆repo,然后运行python setup.py install
用法
Usage: circ make-vrt MINX MINY MAXX MAXY [NAME][EM_TOL][DRY_RUN] circ make-vrt --minx MINX --miny MINY --maxx MAXX --maxy MAXY [--name NAME][--em-tol EM_TOL][--dry-run DRY_RUN]指定最小和最大经度(最小值)和最小纬度和最大纬度。 默认情况下,允许的最大发射角为5,但这可以降低到较低的值。 失去了很多空间覆盖。尝试使用–dry run标志和不同的发射角值 观察bbox的图像计数是如何变化的。
因此,要在盖尔陨石坑周围制作CTX马赛克,安装后只需运行:
circ make-vrt 136.0 -7.0 139.5 -3.5 --name gale --em_tol 1.0
这将创建一个名为gale的文件夹,其中将下载一组ctx图像。发射角公差也降低到1.0(小于或等于) 这里传递的bbox比实际需要的要大一点,因此将下载一些额外的图像。 如果运行此示例,则值得为每个坐标将其收缩半到四分之三度。
最后创建的vrt图像可以直接用于qgis等应用程序,或者使用其他gdal命令行工具(如so:
gdal_translate -co COMPRESS=JPEG -co TILED=YES -co COPY_SRC_OVERVIEWS=YES gale.vrt gale_jpeg.tif gdaladdo -r average --config COMPRESS_OVERVIEW JPEG --config JPEG_QUALITY_OVERVIEW 85 gale_jpeg.tif 248
要查看我们实际节省的空间,这里是gale文件夹占用的总空间:
du -sh ./gale
5.0G ./gale
这是最终jpeg压缩tif版本的空间大小:
du -sh ./gale_jpeg.tif 817M ./gale
此时,您可以删除vrt和tiff文件夹以节省磁盘空间,并使用其他gdal命令将马赛克剪裁到边框。