saltgiser的博客

本文介绍了一种Python解析HGT高程数据文件的方法。HGT是美国地质调查局提供的二进制高程数据格式，但缺乏坐标信息无法直接使用。文章详细解析了HGT文件结构，包括命名规则（如N36E117表示36-37°N,117-118°E）、数据存储方式（2字节有符号整数）和多种分辨率（1角秒、3角秒等）。重点提供了一个完整的Python函数，通过文件名解析地理范围、判断分辨率、读取二进制数据并转换为numpy数组，处理无效值后生成可直接用于分析的高程网格数据。该函数支持不同分辨率，最终输出包含高程网格、经纬度范围

本文探讨了Python实现选址方案时距离分析的高效优化方法。传统矢量计算面临效率低下和结果不稳定的问题，而通过借鉴ArcMap的思路，采用栅格化方法可显著提升性能。具体步骤包括数据预处理、线要素合并与栅格化、基于栅格的距离变换等。优化后，计算效率大幅提升，同时保持合理精度，500米分辨率下平均偏差小于250米，完全满足宏观选址需求。这种栅格化方法为大规模线要素距离计算提供了有效解决方案。

（分辨率此处有个坑，大致面积在350平方公里左右可以下载约10m分辨率的数据，下载后实测大概市15m左右分辨率，如果范围过大，则没办法下载10m左右的数据）背景 ：最近有需求需要哨兵2号真彩色图片，传统方式是下载各波段原始数据自己进行合成，但新的官网可以下载已有的合成过后的数据产品，因此为了偷懒，直接下载已经合成好的True Color Image。根据以上信息进行python编码，我的geojson数据都是在一个文件夹下，命名有规则，然后可以进行用户认证信息的输入，下载图片的日期则是手动写死。

背景：拥有一组包含“开始时间”和“结束时间”的数据，可用于计算每条记录的持续时间。同时，每条记录还包含一个“面积”字段。为了更直观地展示每条数据的持续时间与面积，希望通过图形化方式进行可视化表达。，并在每条线上标注对应的序号以便识别。每条记录绘制为一条表示持续时间的横向线段，解决思路：设计一张时间轴图，