
GIS
cwr888
GIS老油条,有任何GIS问题很乐意一起交流(v:junco--),有项目更好。
展开
-
地图可视化上盖蒙版效果案例记录
地图可视化项目蒙版效果实现案例记录原创 2025-03-24 13:54:16 · 164 阅读 · 0 评论 -
RedHat8安装postgresql15和 postgis3.4.4记录及遇到的问题总结
RedHat8安装postgresql15和 postgis3.4.4错误: 问题: 冲突的请求 - 没有东西可提供 hdf5(postgis34_15-3.4.4-1PGDG.rhel8.x86_64 需要) - 没有东西可提供 xerces-c(postgis34_15-3.4.4-1PGDG.rhel8.x86_64 需要)(尝试添加 '--skip-broken' 来跳过无法安装的软件包 或 '--nobest' 来不只使用软件包的最佳候选)原创 2025-02-14 09:29:19 · 972 阅读 · 0 评论 -
通过leaflet的L.geoJSON添加图层,没图形显示也没报错
通过leaflet的L.polygon加载面数据正常显示,但通过L.geoJSON添加面图层数据却不显示也没报错。原创 2022-08-18 18:05:06 · 1579 阅读 · 0 评论 -
leaflet加载geoJSON图层数据时报Invalid LatLng object
leaflet加载geoJSON图层数据时报:Uncaught Error:Invalid LatLng object错误原创 2022-08-16 21:07:55 · 2938 阅读 · 0 评论 -
python cartopy.io.shapereader Incomplete shapefile definition in
问题:python 通过cartopy.io.shapereader读取shp文件,出现Incomplete shapefile definition in 'c:\aaa.dbf'ditubianjie = ccrs.PlateCarree() # 定义投影转换fname = 'aaa.dbf'features = shpreader.Reader(file).geometries()# 绘制地图ax.add_geometries(features , ditubianjie, facec原创 2022-05-17 19:20:17 · 3824 阅读 · 5 评论 -
Leaflet加turf生成等值线图(或色斑图)并单击显示值
Leaflet加turf生成色斑图并单击显示值范围或精确值参考代码如下:index.html<!DOCTYPE html><html lang="en" ><head> <meta charset="UTF-8"> <title>CodePen - turfjs-interpolate-isoband</title> <link rel='stylesheet' hr.......原创 2022-05-11 19:46:40 · 5332 阅读 · 4 评论 -
python读取nc文件数据生成专业气象图
python读取nc文件数据生成专业气象图如色斑图、等值线、风向标,如下图所示:代码如下:# -*- coding: utf-8 -*-"""Created on Tue Apr 12 10:58:15 2022@author: junco"""import numpy as npimport matplotlib.pyplot as pltimport matplotlib as mplimport scipy.ndimage as ndimageimport xarra.原创 2022-04-25 18:55:16 · 7056 阅读 · 0 评论 -
leaflet加载leaflet-canvas-marker图层时出现偏移问题解决
问题背景: 用leaflet加载高德地图,对地图进行操作(如:放大、缩小或移动)后,再动态加载leaflet-canvas-marker图层,图层第一次(初始)叠加在地图上有偏移,对地图进行操作(如:放大、缩小或移动)后,偏移消失。[特别注意要先对地图进行操作,再加载leaflet-canvas-marker图层,否则问题不能重现]。如下图: 第一次(初始)叠加位置: 对地图进行操作后位置: 初步判定为leaflet-canvas-marker的一个Bug。问题解决:动态加...原创 2022-04-25 18:29:24 · 5077 阅读 · 0 评论 -
气象NC扇形经纬网格转换成前端要求的等经纬网格
气象NC扇形经纬网格转换成前端leaflet-vector-scalar.js要求的等经纬网格背景:最近从气象局拿到文件格式为NC的气象文件(包括温度、湿度、风、气压、雨量等),需要读取其中的温度数据并在前端展示。用专业软件打开一看,图是扇形的,经纬度间隔也不相等。问题:需要读文件格式为NC的气象文件中的温度数据,并转换为leaflet-vector-scalar.js要求的等经纬网格JSON数据。解决:通过JAVA调用netcdf包读取NC文件,将NC扇形经纬网格转换成等经纬网格,转.原创 2022-03-08 18:09:58 · 1478 阅读 · 0 评论 -
将ArcGIS Server的EsriJSON转化为GDB文件,支持带弧的面和线
将ArcGIS Server的EsriJSON转化为GDB文件,支持带弧的面和线# -*- coding: utf-8 -*-# ---------------------------------------------------------------------------# JSONToGDB.py# Created on: 2022-02-12 17:00:00.00000# # Description:# ------------------------------------原创 2022-02-13 11:08:31 · 606 阅读 · 0 评论 -
ArcGIS实现按比例尺进行地图全图导出一张完整的大比例地图
ArcGIS实现按比例尺进行地图全图导出一张完整的大比例地图1、启动 ArcMap2、打开一个mxd地图文档3、切换到布局视图框,视图-布局视图,或点击左下角的布局视图4、右键"视图"-"数据框属性"-"数据框"-"范围",选择要输出的图片固定比例5、点击"文件"-"页面和打印设置"。主要是选择"根据页面大小的变化按比例缩放地图元素"的复选框,然后指定页面宽、高(特别注意宽高可以手动指定,不是只能选标准A4等纸张大小)。6、点击"文件"-"导出地图...原创 2022-01-19 18:26:20 · 26442 阅读 · 0 评论 -
关于postgresql数据库pgrouting计算最短路径突然变慢问题排查记录
问题:用postgresql数据库pgrouting插件开发的最短路径规划服务,两个月来一直正常运行,今天突然变得很慢了起来。问题排查:1、首先路网数据没有修改,排除数据问题;2、查看各结点服务器服务正常,排除服务问题;3、查看数据库运行状态,也正常,排降数据库程序问题;4、查看系统,CPU资源正常,硬盘空间正常,但发现异常内存快满了,如下图所示:尝试解决:杀掉占内存的无用进程或增加内存,本次解决方案将占内存大的mysql数据库占用的内存调小。结果:最短路径规划服务速度恢复正常。结论:pos原创 2021-11-12 17:48:59 · 909 阅读 · 0 评论 -
Openlayers显示-360到360的网格图层
Openlayers显示-360到360的网格图层,效果如下图:实现代码如下:<!DOCTYPE html><html> <head> <meta charset="UTF-8"> <title>Openlayers显示-360到360的网格图层</title> <link rel="stylesheet" href="https://cdn.jsdelivr.net/gh/openlayers/ope.原创 2021-11-04 17:53:27 · 528 阅读 · 0 评论 -
Openlayer地图折线上点的信息错位标注方法
Openlayer地图折线上点的信息错位标注方法,尽量不压到线,如下图所示:解决思路如下:一、起点处理起点与第二点标牌方向相反二、中间点处理以中间点为坐标系原点,求中间点连接的两条边线在第几象限,并求出每天条边线与y轴正半轴之间的夹角。1、两个边不在一条直线上,标牌标在中间点连接的边线的两个夹角的角度最大的一边(1)两个边在同一象限case 同象限1、1 标牌标在象限3case 同象限2、2 标牌标在象限4case 同象限3、3 标牌标在象限1case 同象限...原创 2021-10-22 17:16:45 · 556 阅读 · 0 评论 -
Openlayers前端复用Turf.js生成平头Buffer
Openlayers前端复用Turf.js生成平头Buffer,实现如下效果:代码如下:<!doctype html><html lang="en"> <head> <link rel="stylesheet" href="https://cdn.jsdelivr.net/gh/openlayers/openlayers.github.io@master/en/v6.2.1/css/ol.css" type="text/cs.原创 2021-10-19 16:21:55 · 1221 阅读 · 0 评论 -
Geoserver清除“资源缓存“和重新载入“配置和目录“报错Error creating bean with name ‘gwcFacade‘ defined
Geoserver清除"资源缓存"和重新载入"配置和目录"后报错如下Error creating bean with name 'gwcFacade' defined in URL [jar:file:/geoserver219/webapps/geoserver/WEB-INF/lib/gs-gwc-2.19.0.jar!/applicationContext.xml]: Cannot resolve reference to bean 'DiskQuotaMonitor' while s.原创 2021-09-27 14:59:09 · 1840 阅读 · 0 评论 -
Geoserver 编辑SQL视图时由参数名引起的错误
Geoserver 编辑SQL视图时由参数名引起的错误问题如下:ava.lang.ClassCastException: java.lang.Double cannot be cast to java.lang.String java.lang.Double cannot be cast to java.lang.String编辑SQL视图中给定的错误参数名angle如下:解决办法,修改参数名为inAngle后错误问题解决,修改参数名后编辑SQL视图如下:总结,可能angle为geoserve.原创 2021-08-12 17:53:20 · 600 阅读 · 0 评论 -
Geoserver通过SLD配置字段值在面边线内侧显示
Geoserver通过SLD配置字段值在面边线内侧显示,如下图所示:参考配置SLD代码如下: <sld:Rule> <sld:Name>aaa</sld:Name> <sld:Title>aaa</sld:Title> <sld:MinScaleDenominator>100000</sld:MinScaleDenominator> .原创 2021-07-30 11:40:53 · 648 阅读 · 0 评论 -
Geoserver通过SLD配置道路线样式,随着放大加粗
一、需求,Geoserver通过SLD配置道路线样式,并随着放大加粗,随着缩小变细效果如下图:二、解决,通过uom 属性指定测量单位实现SE 1.1规范中给出的测量单位,如下:http://www.opengeospatial.org/se/units/metrehttp://www.opengeospatial.org/se/units/foothttp://www.opengeospatial.org/se/units/pixel以下SLD显示 uom 用于指定 L转载 2021-07-23 16:30:56 · 1227 阅读 · 0 评论 -
Geoserver通过SLD配置地图服务
配置如下线样式:代码如下: <sld:PolygonSymbolizer> <Geometry> <ogc:Function name="boundary"> <ogc:PropertyName>the_geom</ogc:PropertyName> </ogc:Funct原创 2021-07-23 16:02:27 · 442 阅读 · 0 评论 -
Geoserver点图层设置SVG图标显示不全问题记录
一、问题。Geoserver点图层设置SVG图标时,有些图标在地图上显示时会显示不全。如下图所示:SVG图标如下:<?xml version="1.0" encoding="utf-8"?><svg version="1.1" id="layer_1" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink" x="0px" y="0px" viewBox="0 0 ...原创 2021-07-14 09:31:01 · 1977 阅读 · 0 评论 -
Arcgis出图时提示未能加载文件或程序集stdole, Version=7.0.3300.0错误
Arcgis出图时提示未能加载文件或程序集stdole, Version=7.0.3300.0错误Arcgis开发的net程序,部署到新安装服务器时,调用出图功能时提示错误:未能加载文件或程序集“stdole, Version=7.0.3300.0, Culture=neutral, PublicKeyToken=b03f5f7f11d50a3a”或它的某一个依赖项。系统找不到指定的文件。解决方案:1、服务器安装windows office软件2、服务器安装Arcgis Desktop软件...原创 2021-03-24 18:43:16 · 1881 阅读 · 0 评论 -
有些CAD通过Arcgis程序读取后,发现面积不对
有些CAD通过Arcgis程序读取后,发现面积不对,后来发现一种情况是CAD绘图时将多个相连的实体拼成了多个镂空图形,导致Arcgis程序(或其它图形程序)读取后构图不对导致面积不对。如下图所示。1、CAD图形将多个相连的实体拼成了多个镂空图形,如下图所示:解决方法参考:1、通过ArcMap将CAD转成SHAPE。2、通过ArcMap将相连的实体合并,如下图所示:3、再通过ArcMap将SHAPE转成CAD,Arcgis程序(或其它图形程序)读取构图将正确同时面积也会正确。...原创 2021-03-09 20:05:31 · 8962 阅读 · 0 评论 -
ArcGIS编辑数据,处了调用AE、AO接口外,还可以通过SQL实现
ArcGIS编辑关系数据库SDE中数据时,处了调用通用的AE、AO接口外,还可以通过SQL语句实现,参考SQL语句如下:1、查询语句select d.objectid, d.NAME,sde.st_astext (d.shape) as coords from SDE.DISTRICT d2、创建语句INSERT INTO SDE.POLYGON (OBJECTID, shape) VALUES (642, sde.st_geometry('polygon((487353.04230000 24原创 2021-02-22 16:16:26 · 230 阅读 · 0 评论 -
开源PostgreSQL、POSTGIS、GeoServer环境搭建及数据管理与发布笔记
开源PostgreSQL、POSTGIS、GeoServer环境搭建及数据管理与发布笔记一、PostgreSQL下载、安装。进入官方网站根据需要选择版本。PostgreSQL官方下载地址:https://www.enterprisedb.com/downloads/postgres-postgresql-downloadsPostgreSQL安装步骤参考:https://blog.youkuaiyun.com/antma/article/details/83579920二、POSTGIS下载、安装。由于用P原创 2021-02-03 19:09:57 · 657 阅读 · 0 评论 -
ARCGIS环空洞和岛图形坐标顺序问题保存或几何修复后情况记录
ARCGIS环空洞图形中,不管空洞和岛坐标顺序是顺时针还是逆时针,保存或者几何修复后,空洞坐标顺序都会调整为逆时针,岛坐标顺序都会调整为顺时针。即如下四个坐标串:1、环 顺时针-逆时针-顺时针-逆时针[[[[105844.32,41880.61], [105902.08,41857.77], [105928.56,41915.48], [105869.59,41941.94]],[[105874.32,41921.74], [105909.97,41913.04], [105898.84,41877.原创 2020-12-09 15:03:01 · 973 阅读 · 0 评论 -
Arcmap数据框的空间参考改变导致坐标值改变问题
Arcmap数据框的空间参考改变导致坐标值改变问题1、Arcmap中MXD设置坐标系为2435时,编辑数据时查看数据值如下:2、将Arcmap中MXD设置坐标系为4547时,编辑数据时查看数据值如下:因此相同的矢量数据,设置不同的坐标系,编辑时显示的值也不一样。因此做数据编辑和数据校验时一点要确保arcmap数据框的空间参考和图层的空间参考要一致。...原创 2020-11-24 19:20:18 · 1534 阅读 · 0 评论 -
经纬度坐标与高斯坐标的转换JAVA代码(无带号)
经纬度坐标与高斯坐标的转换JAVA代码(无带号),以中央子午线114.0为例。 // 由高斯投影坐标反算成经纬度 public static double[] GaussToBL(double X, double Y) { double[] output = new double[2]; double longitude1,latitude1, longitude0, X0,Y0, xval,yval; //NN曲率半径,测量学里面用N表示 //M为子午线弧长,测量学里用原创 2020-08-27 20:15:00 · 1693 阅读 · 7 评论