GIS小虫的博客

将本地数据同步到内网服务器，使用的postgis，表含空间字段。

对于ArcGIS server的mapserver，可以方便的获取图层字段信息，如下http://cache1.arcgisonline.cn/arcgis/rest/services/SimpleFeature/ChinaBoundaryLine/MapServer/0?f=pjsongeoserver的要素图层，可以通过WFS服务的DescribeFeatureType请求获取字段信息（范围，坐标参考等信息可以通过wms的GetCapabilities）例子图层信息如下请求地址http

需要GDAL源码OpenCLVisual StudioSQLite（可选）PROJ（可选）1、查看自己电脑的显卡配置支不支持OpenCL方法很多，可以自行百度。我是下载TechPowerUp GPU-Z（点击链接进入下载页）查看OpenCL项目前有没打勾2、安装OpenCL可以使用Intel，Nvidia官网都能下，也能用CUDA里带的。我用的Intel家的Intel® SDK for OpenCL，选择符合你本地系统的SDK下载（需要自行注册账号），安装。安装过程没啥需要注意的，

问题源于使用4528投影按照官网示例定义view时地图一片空白，于是开始了对Proj4js和投影转换的抽丝剥茧，这里记录下了解的知识和最后的结论。官网示例和Proj4js的使用OL示例见wms-image-custom-proj这里简单介绍下proj4，我认知的proj4是基于一些约定参数组成的用于定义坐标参考的一组参数，和WKT相比更简约，省略了WKT里的完整参数定义。比如4528的wktPROJCRS["CGCS2000 / 3-degree Gauss-Kruger zone 40", .

我是传了个样式SLD上来，预览图例的时候发现中文都变成了小方块经过一番研究，发现这个图例是每次现请求geoserver现生成的，可以通过设置请求参数里的字体为中文字体来生成中文图例如http://localhost:8080/geoserver/wms?REQUEST=GetLegendGraphic&VERSION=1.0.0&FORMAT=image/png&WIDTH=20&HEIGHT=20&LAYER=basic-platform:crops_typ

geoserver可以通过GetCapabilities接口获取wms图层的原始信息，比如图层列表，图层的范围等。但是随着发布的图层数量增加，这个返回的结果会很大，这里记了一些优化方法以减小其体积。设置空间参考列表不设置的话，你会看到这么一大堆无用信息设置方法如下这里直接填写需要的空间参考的SRID需要的话也可以同时勾选Output bounding box for every supported CRS，选择输出所需坐标系的图层边界。效果如下十分方便加载图层后设置extent或者缩放

根据openlayers官方的popup实例改的，使用Vue2+TS，功能比较简单OLPopop.vue<template> <div ref="container" class="ol-popup"> <a href="#" ref="closer" class="ol-popup-closer"></a> <slot></slot> </div></template>impo

python GDAL生成COG

OpenLayers WebGL瓦片图层简介和样式

COG的简介，生成，使用

还没有解决，先记一下，以后再研究用GDAL生成了个cog，生成中使用band.ComputeStatistics计算了波段信息，现在用gdal.Info可以获取到波段的最大最小值信息但是geotif.js获取不到元信息。geotifjs源码里解析主函数应该是这个我想要的那个tag，GDAL_METADATA，对照如下结合查到的资料，应该是保存了需要的元信息的，这些信息作为IFD保存在了tif文件里，geotifjs在parseFileDirectoryAt函数里读取IFD然而读不到。所

JSTS是一个符合OGC标准的，包含空间拓扑功能的JavaScript类库。和turf类似，区别在于，JSTS是通过原始JTS Java 源代码翻译来的，保留了JTSAPI；而turf来自MapBox。个人感觉JSTS的文档不如turf明了(一些网址还404了)，但是可以直接适用OpenLayers（以下简称OL）的几何对象，对JTS API熟悉的话在OL工程里可能会更方便。JSTS的参考文档比较分散，GitHub库里有比较全面的简介，这个地址里有一些简单的示例，JTS API直接用的JavaAPI，.

python安装GDAL推荐使用conda

用到的插件是ol-wind，作者除了Openlayers风场插件外还有其他平台的，详见wind-layer获取数据git里有如何获取数据的介绍，这里记一下我的实践过程。数据来源是nomads，见下图。git上curl命令里的下载地址现在好像不好使了，所以我用的是GFS 1.00 Degree，点击grib filter可以进入后续的筛选。一路往后点，到这个页面上面应该是一些参数啥的，总之就照着git上curl命令里的参数选或者修改curl命令里的数据源，把filter_gfs.pl改.

leaflet加载瓦片图层时，若缩放到没有瓦片的层级或者超出默认的minZoom、maxZoom，会一片空白。比如图层服务最高只有12级，地图缩放到13级就没了。若想再目标层级显示最近层级的瓦片的放大或缩小，比如缩放到13级就把12级的瓦片放大一点，可以设置maxNativeZoom/minNativeZoom// 超出NativeZoom的不再请求，而是放大原有的// 超出Zoom的不显示const tile = L.tileLayer(url, { minZoom: 0, minNative

在http://download.gisinternals.com/sdk.php下载所需版本的压缩包将下好的压缩包里release-1928-x64/bin/目录下的东西解压到jdk/bin目录下（也有说放jre/bin目录下或者java工程根目录下的，还没尝试）将release-1928-x64/bin/gdal/java/下的gdalalljni.dll或者还有其他dll拷贝到java工程根目录，gdal.jar拷贝到java工程里并用maven引用<dependenc..

前言OpenLayers最近版本（6.11.0）上看到了使用WebGLTile图层加载GeoTIFF的示例，功能强大，不仅可以在前端直接显示tif影像，还可以做分波段彩色合成，对比度拉伸等色彩上的调整。简单试了下，发现数据源上存在一定的限制WebGLTile的source属性接受DataTileSource和TileImage两种source类型，结合示例里的数据源发现WebGLTile图层只接受单张tif或者XYZ切片格式数据源，目前还不支持WMS地图服务，虽然GeoServer的WMS服务支持输出

如题

背景使用openLayers 6.5.0，html2canvas 1.3.3，使用geoserver发布的地图服务。问题使用html2canvas打印时，地图一片空白，控制台报错’Unable to clone canvas as it is tainted’解决1、加载地图服务，配置source时添加配置：crossOrigin: 'anonymous'2、同时geoserver或者放geoserver的tomcat需要配置跨域，否则服务会报跨域的错坑不要在html2canvas里配置

发布连接动态空间的服务可选数据库，矢量数据(shp)文件夹，栅格数据(tif)文件夹。下面示例的是PostgreSQL数据库。这里的workSpaceId是JS调用的关键属性。前端代码API版本4.16data Source官方配置文档参考Data Source object// 返回的是一个引用的mapImageLayerreturn arcgisFuns.mapImageLayer({ id: 'mapImageLayer_test', url, sublay

之前添加外部的图片作为symbol时图片路径一直有问题。有两种方式，一种是图片放到nginx或tomcat里，在symbol的url里使用绝对路径，这种适合图片可能动态增加或修改的情况。另一种如下 private addPoint(angle: number) { this.view.graphics.removeAll(); const symbol = { type: 'picture-marker', url: require('../assets/im

4X去掉了之前的ArcGISDynamicMapServiceLayer，替换的MapImageLayer更能更强大一些，但是MapImageLayer没有图形的鼠标交互，对需要添加图形交互的MapServer图层需要用FeatureLayer加到map对象上。 // 包含图层url等信息的对象数组 layerList.reverse().forEach((ele: any) => { let layer; if (ele.type === 'feature

根据geotools自己封装的工具类，主要有shp文件的内容读取，转geojson，要素的增删改，新shp生成需要的依赖这里主要用的是geotools的依赖，版本是<geotools.version>23.2</geotools.version>还用到了fastjson，版本1.2.57<dependency> <groupId>com.alibaba</groupId> <artifactId>fastjson

什么是GeoToolsGeoTools是一个Java类库，它是构建在OGC标准之上的，提供了很多标准类和方法来处理空间数据，是开源空间数据处理的主要工具。什么是JTSJTS是加拿大的 Vivid Solutions 做的一套开放源码的 Java API。它提供了一套空间数据操作的核心算法,为在兼容OGC标准的空间对象模型中进行基础的几何操作提供2D空间谓词API。（目前已经切换为org.locationtech开源组织了）。JTS解决了空间对象之间拓扑关系的判定和计算，并提供很多有用的算法来解决对象

合并两个shp文件用到的工具UnionFeatureCollection需要的依赖 <dependency> <groupId>org.geotools</groupId> <artifactId>gt-process</artifactId> <version>${geotools.version}</version>

这一系列翻译自openlayers官网的WorkShop。OL官网提供了多个系列教程供开发者学习参考，其中QuickStart是面向初学者的hello world，Tutorials提供了构建OL应用的一些基础知识，WorkShop（本系列）详细介绍了一些入门向的高阶应用，最后是APIDocs，适合开发时查阅接口。教程中需要下载的资源可以在WorkShop原网站获得链接。在本模块中，我们将创建一个显示用户 GPS 位置和航向的移动地图。此示例的目的是展示如何将 OpenLayers 与浏览器 API.

1.数据准备需要准备发布的栅格数据（tif格式，不支持img），以及包含所有栅格数据边界的shp文件。2.服务发布添加数据存储，选择imageMosaic工作区和名称看着填写，url点击浏览，选择上一步包含了影像和shp的文件夹，确定->保存新建图层，选择刚才建立的数据源。图层设置里，保证坐标参考系统该有的值都有（一般没问题，都会自动从数据源获取），覆盖参数看情况设置，全部保持默认也行。点击保存，就可以预览了。完成后数据所在的目录会新增以数据源为名的...

这一系列翻译自openlayers官网的WorkShop。OL官网提供了多个系列教程供开发者学习参考，其中QuickStart是面向初学者的hello world，Tutorials提供了构建OL应用的一些基础知识，WorkShop（本系列）详细介绍了一些入门向的高阶应用，最后是APIDocs，适合开发时查阅接口。教程中需要下载的资源可以在WorkShop原网站获得链接。

本文整理了地图坐标和投影的基础知识，一方面是自己备忘，同时面向非gis/测绘行业的普通开发人员做个基础知识铺垫，对常见问题有个了解。1、不同坐标系的存在，投影的意义。2、

栅格数据集的适量信息记录在数据库中 数据名_cat表的shape字段里。行政区划建了个要素类(FeatureClass)。注意和栅格数据集坐标系需要一致。插一下矢量84转2000参考（这个文章）重点是需要先转成ITRF2000再在ArcCatalog里重新定义坐标系再回ArcMap里转。行政区划要素类这里有名称和代码，当然还有几何属性shape下面上SQL这里先把cat表里的shape和行政区划的shape使用sde.st_intersects做相交，获取行政区划表里相交为true的县级行政