saltgiser-CSDN博客

原创 HGT 文件解析：从二进制数据到高程网格的实现方案

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

2025-09-25 11:24:27 655

原创 OpenLayers 绘图面积校验：为什么 drawend 里删不掉图形？

本文探讨了OpenLayers地图开发中实现图形面积校验功能时遇到的事件选择问题。作者最初在drawend事件中校验面积并删除超限图形，但因图形尚未存入数据源而报错。后改为在addfeature事件中处理，此时图形已存入数据源，操作得以正常执行。文章总结出"时机原则"：绘制动作相关逻辑用drawend事件，图形数据操作必须用addfeature事件。类似网购流程，只有商品入库后才能进行后续操作。此外还提到可通过实时计算或虚拟图形等替代方案解决该问题。

2025-08-21 20:06:55 497

原创基于栅格的线要素距离计算：从矢量困境到高效优化（Python实现）

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

2025-08-20 17:28:06 616

原创 WEBGIS-WKT 几何合并工具

在地理信息系统（GIS）应用开发中，经常需要处理各种空间几何数据。WKT（Well-Known Text）作为一种开放的文本标记语言，被广泛用于表示矢量几何对象，如点、线、面以及更复杂的 MultiPolygon 等类型。然而，现有工具要么功能单一，要么操作复杂，缺乏一个轻量、直观且支持多种几何类型的合并工具。为此，我们开发了这款基于 Web 技术的 WKT 几何合并工具，旨在为开发者和 GIS 从业者提供便捷的几何合并解决方案。

2025-08-01 16:42:30 888

原创 python批量下载哨兵2号TrueColor图片

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

2025-07-29 11:08:36 503

原创 python使用matplotlib库生成持续时间图表

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

2025-05-30 16:59:03 638

原创 openlayers 实现多个面组成一个多面并获取范围

本文介绍了一个用于获取多个面要素组成范围的函数getPolygonsExtent。该函数能够处理Polygon和MultiPolygon类型的要素，通过收集所有面的坐标，构建一个新的MultiPolygon对象，并返回其地理范围。适用于需要计算多个面要素整体范围的地图缩放场景，解决了从geojson获取多个面要素时的范围计算问题。

2025-05-28 10:55:09 351

原创 ArcMap中对矢量数据某个字段进行归一化

摘要：在网格统计中，使用自带字段计算器进行归一化处理时，结果往往不够精确（如最大值归一化值仅为0.99379）。为解决此问题，本文采用arcpy代码实现归一化计算，通过公式(x-min)/(max-min)在ArcMap Python工具中运行，获得了准确的归一化结果。该方法有效提升了数据处理的精度和可靠性。

2025-05-27 15:45:05 720

原创 Arcmap 模型构建器(Model Builder)中带有后缀的文件重命名方式

在输出结果的命名时，栅格数据名称（如“xxx.tif”）未去除“.tif”后缀，致使输出文件名为“xxx.tif.dbf，出现了命名错误。第一种方法：使用了参考文献1中的方案，使用模型构建器中的 Calculate Value工具进行文件命名操作，测试能够跑通，但需要注意的拖入设置好 Calculate Value工具参数时，前面栅格迭代数据要重新设置下，因为之前已经跑过一次了。在ArcMap中，模型构建器（Model Builder）是一种强大的工具，可以自动化复杂的地理数据处理流程，并支持批量操作。

2025-04-28 11:53:57 2540 1