GIS探索之路

图中的图表展示了标准正态分布曲线（红虚线）、使用Silverman经验法则计算带宽后的核密度估计曲线（蓝实线）以及数据点（红竖线）。从图中可以看出，在多峰分布的数据场景下，相比基于Silverman经验法则的核密度估计，使用ISJ算法得到的曲线能更好地拟合数据的多个峰值，体现了其在处理多峰数据上的优势。图中的图表展示了核函数（Kernel，虚线）、核密度估计（KDE，实线）和数据点（红点）之间的关系。图中的图表展示了数据点（红点）、不同带宽下的核函数（虚线）以及最终的核密度估计曲线（实线）。

使用QGIS Processing工具将采集的高程点转换为Mesh，格式为PLY。目前QGIS支持加载Mesh及3DTiles。

使用`CGAL`将地形曲面网格体栅格化。

以计算植被`NDVI`指数为例，介绍`QGIS`栅格的相关操作以及应用。

qgis矢量数据处理，字段编辑、字段计算器、几何形状有效性检查、Z值提取以及多几何炸开等等操作。

Scoops3D &GIS &DEM 使用DEM 对滑坡稳定性分析

主要使用材质世界偏移对静态网格的顶点位置进行修改，CesiumForUE 地形压平测试。

.NET 环境下非常简单的GIS地图开发库。

本教程将展示如何以编程方式从网络地图（通常称为瓦片地图）瓦片服务器下载地图图像，对其进行地理参考（设置坐标系）并将其保存为GeoTIFF。

使用python GDAL 对下载好的瓦片（一定处理后）进行拼接。

使用python GDAL进行简单的坐标系转换，暂时不考虑不同基准坐标系转换的精度问题。

通过CesiumJS 沙盒快速测试CesiumJS的一些功能，免去安装开发环境的困恼。

使用QGIS快速构建任意地区三维可视化地形效果图及其他。

使用QGIS加载常见在线地图服务及相应插件使用和导出感兴趣地区的数据。

PLY，多边形文件格式，表示存储描述为多边形集合的图形对象的 3D 文件格式。这种文件格式的目的是建立一种简单易用的文件类型，该文件类型足够通用，可用于各种模型。PLY 文件格式有 ASCII 和二进制格式，用于紧凑存储和快速保存和加载。PLY 格式的对象由一组顶点、面和其他元素以及可以附加到这些元素的颜色和法线方向等属性来描述。由 PLY 格式表示的对象可以是各种来源的结果，例如手动数字化对象、建模应用程序中的多边形对象、范围数据、立方体中的三角形、地形数据和辐射模型。

SLPK bin 文件与draco 解压缩SLPK bin 文件与draco 解压缩0.bin1.bindraco 解压示例draco 压缩示例参考SLPK bin 文件与draco 解压缩0.binSLPK 二进制0.bin 早前使用的几何文件存储格式，没有进行draco压缩处理。4个字节 表示顶点的数量顶点数据法线数据纹理数据颜色数据FeatureIdFaceRange解析代码： ifstream fin; fin.open("0.bin", ios::binary);

OSGB数据是一个四叉树结构，本文测试使用OSG加载OSGB倾斜摄影测量数据。

使用 `cesium terrain builder`工具制作3DTiles地形。

地理配准（Georeferencing） 是指将某个环境与物理空间中的位置相关联。该术语常用于地理信息系统（GIS）领域，以描述将物理地图或地图的栅格图像与空间位置相关联的过程。一旦对环境进行地理配准后，虚幻引擎就可以表示实际地理坐标（如纬度/经度或UTM坐标）。

UE4 集成 Proj4 进行 GIS 坐标转换操作。

CityJSON 数字孪生模型

3D 城市数据库是一个免费开放的 3D 地理数据库，借助于强大的标准空间关系数据库，可以对虚拟 3D 城市的空间数据模型进行存储、表示和管理。

GDAL 根据DOM生成DEM根据DOM影像生成指定高度的地形高度图，分辨率可指定。安装GDAL ToolsGDAL 版本3.2.0。 .\vcpkg.exe install gdal[tools]:x64-windows --recurseGDAL gdal_creategdal_create.exe 用于生成栅格。gdal_create -outsize 20 20 -a_srs EPSG:4326 -a_ullr 2 50 3 49 -burn 10 out.tif创建宽高为2

QuantizedMesh2STL Demo解析 cesium 地形 Quantized-Mesh数据，生成STL。请求地形瓦片http://example.com/tiles/z/x/y.terrain解析下载数据及解析1_0_0.terrain->1_0_0.stlbounds:(-180.0, -90.0, -90.0, 0.0)https://assets.cesium.com/1/1/0/0.terrain?v=1.2.0&extensions=octvertex

Atlas - A OpenSource 3D GIS Visualization ProjectAtlas 是一个三维地图可视化分析项目。Atlas 结合了 QGIS 和 谷歌地球 的主要特性，利用 osgEarth 来将二维地图和三维模型展示在同一三维空间内。程序使用 Qt 作为框架，以提供简洁、可定制的 UI，并且实现了一个高扩展性的插件系统。该项目的主要目的是为地理应用提供更好的可视化解决方案，并在此基础上支持复杂的空间分析、数据处理等，以方便工程应用和学术研究。该程序在 Windows 1

空间填充曲线从数学的角度上看，可以将空间填充曲线看成是一种把N维空间数据转换到1维连续空间上的映射函数。实际上，存储磁盘是一维的存储设备，而空间数据是多维数据，不存在天然的一维顺序。因此，为了使空间上邻近的元素映射也尽可能是一维直线上接近的点，提出了许多的映射方法。最常用的方法包括Z Ordering、Hilbert 曲线等等。皮亚诺曲线希尔伯特曲线z ordering curveS2参考3DTILES_bounding_volume_S2Debug S2 in CesiumJS

3D Tiles Next将地理空间数据流式传输到元宇宙世界Metaverse。Extensions高效地流式传输语义元数据通过空间索引对大规模的地理空间数据进行可视化模拟和分析统一的生态系统： glTF content 扩展属性数据的扩展空间索引的扩展KD-Tree, loose quadtree, and octree3DTILES_implicit_tilingS2 Subdivision集成gltf生态b3dmi3dmpnts参考I

RTree 3D 空间查询给定3D空间的包围盒，查询与包围盒特定空间关系的几何对象。本文使用的是LibSpatialIndex 开源库。自定义数据从 索引（int tileid） 到几何的映射#ifndef RTREE_BUILDER_3D_HPP#define RTREE_BUILDER_3D_HPP#include <iostream>#include <spatialindex/SpatialIndex.h>/*3D 空间包围盒*/struct mbb_

gdal处理landsat8合成真彩色影像tiflandsat8 影像的4，3，2波段分别对应于可见光的红光、绿光和蓝光波段，本文使用开源GIS 库jupyter python gdal 合成真彩色影像。下载landsat8数据在地理空间数据云上下载北京地区的landsat8影像数据codeimport osfrom osgeo import gdalimport numpy as npband1_fn='data/LC81230322021250LGN00/LC08_L1TP_123

OSGeolive 开源GIS 使用OSGeoLive是直接从可启动DVD或USB拇指驱动器运行的Linux发行版。也可以加载在VMware Workstation环境中运行的预制虚拟机磁盘文件(vmdk)。它允许您在不安装任何东西的情况下尝试各种开源地理空间软件。OSGeoLive PackagesOSGeoLive 14 软件清单VMware 安装OSGeolive vmdk自定义创建一个Linux虚拟机，在创建硬盘的时候选择现有虚拟硬盘即OSGeolive的vmdk文件。安装 VMw

mage3D韩国人开发的一款3DGIS 平台, 曾用名Gaia3D(Gaia 大地女神)，现取名mago3D。mago3D是下一代3D平台，可以将AEC（Architecture、Engineering、Construction）建筑、工程和建造领域的管理和可视化与传统的3D空间信息（3D GIS）相结合起来。[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-mnG3h8Su-1633780843516)(https://note.youdao.com/yws/res/16

Geoserver 发布 wmts瓦片服务准备数据全球的带投影栅格数据，可以利用SXEarth下载拼接而来，分辨率较低。安装部署Geoserver前提是安装Java环境。下载免安装的Geoserver压缩包解压即可。启动服务：注意：需要检查JAVA_HOME变量是否配置成功。配置工作区、数据源和图层登录Geoserver 服务器界面：http://localhost:8080/geoserver/web/新建一个工作区：global_map为工作区添加数据存储：新建图层：配

几种常见的地图瓦片服务按投影类型分为 WGS84 （EPSG：4326）和Web墨卡托（EPSG：3857）。WGS84 属于椭球体便于计算量测，常用于GPS或其他测图相关；Web墨卡托属于标准的球体，便于显示，常用于Web地图显示。按照服务类型分为：WMTS 和TMS。OpenStreetMapWeb墨卡托投影瓦片规则（0，0，0）https://a.tile.openstreetmap.org/0/0/0.png（1，0，0）https://a.tile.openstreetm

Unreal GISsupport gis (datum wgs84)support glob 3d earthmap tms wms(local vector maps raster)https://normalvector.com/tutorials/low_poly_world/building_the_globe/train hight (tms wms)gis basic tools (ppi los roller)51World creates digital twin

三维七参数转换模型三维仿射变换的一般形式[XYZ]=[ΔXΔYΔZ]+(1+m)[1θz−θy−θz1θxθy−θx1][X0Y0Z0]\begin{bmatrix} {X}\\ {Y}\\ {Z} \end{bmatrix}=\begin{bmatrix} {\Delta X}\\ {\Delta Y}\\ {\Delta Z} \end{bmatrix}+(1+m)\begin{bmatrix} {1}&{\theta_z}&{-\theta_y} \\ {-\the

二维四参数转换模型仿射变换的一般形式[x2y2]=[ΔxΔy]+(1+m)[cos⁡α−sin⁡αsin⁡αcos⁡α][x1y1]\begin{bmatrix}x_2\\y_2 \end{bmatrix}=\begin{bmatrix}\Delta x \\ \Delta y\end{bmatrix} + (1+m)\begin{bmatrix}\cos\alpha & -\sin\alpha \\ \sin\alpha & \cos\alpha\end{

投影带和中央经线计算/** * 投影带的带号计算：通过用户输入的经纬度及定义的投影带宽度，得到该经纬度下投影带的带号 * @param long 输入 用户输入的经度 * @param projectedBandWidth 输入 用户定义的投影带宽度 * @returns {number} 返回 经纬度下投影带的带号 */function bandNumberCalculation(long, projectedBandWidth) { var num; if (lon

var BJ54_JSON = { //3度+带号 "BJ54_zone_3_25": "2401", "BJ54_zone_3_26": "2402", "BJ54_zone_3_27": "2403", "BJ54_zone_3_28": "2404", "BJ54_zone_3_29": "2405", "BJ54_zone_3_30": "2406", "BJ54_zone_3_31": "2407", "BJ54_zone_3

proj4测绘坐标系转换相同基准面之间投影坐标和地理坐标的转换pj_transform projPJ pj_merc,pj_latlong; double x,y; //定义墨卡托投影坐标系 pj_merc=pj_init_plus("+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"); //定义地理坐标系 pj_latlong = pj_in

概况来说GIS坐标系统分为一下三类：1、平面直角坐标系（xy）（笛卡尔坐标系）：各类投影坐标（UTM，高斯投影等等）。2、球面坐标（或是地理坐标`Geographic`）（亦或是大地坐标`Geodetic`）（亦或是参心坐标）（亦或是经纬度BLH）：北京54，西安80等等都是属于这一类。3、空间直角坐标系（地心坐标系XYZ）（`Geocentric`）:WGS84，CGS2000等等。如上图所示，从上到下依次为两种不同基准下的平面坐标、地理坐标和地心坐标的相互转换方式。...