兰小莫-CSDN博客

原创 axios封装geoserver wfs服务请求中cql_filter的intersect方法

1.cql支持属性查询或者空间关系查询过滤，具体可以参考cql reference。https://docs.geoserver.org/latest/en/user/filter/ecql_reference.html#filter-ecql-reference2.wfs是发布在地图服务器上矢量数据的要素服务类型，可以对wfs进行查询、新增、更新、删除等操作。首先可以去geoserver官网文档查看具体的来源https://docs.geoser...

2021-08-07 20:48:22 1030

原创 Java读写excel,excel转成json写入磁盘文件

pom读写excel主要的dependency：  <dependency> <groupId>org.apache.poi</groupId> <artifactId>poi</artifactId> <version>3.16</version&...

2020-11-01 21:58:46 412 2

原创 selenium自动化学习笔记（Java版）

背景：selenium是python爬虫，同时也是自动化测试使用的一个模块，当然不止有python语言版，本文就是Java版，因为python相关的博客很多了，而我个人更加倾向使用Java到实际的工作项目中。在几年前，尤其是没有做过web相关开发之前，曾经使用过不少python进行过爬虫练习，但都是零零星星地跟着其它博客的例子进行学习，改造成为自己的需求。不过，在经历了webgis开发，部分Java开发甚至安卓移动端自动化（用的是Auto.js，后面的文章有时间会写相关的内容）后，便有了想往P...

2020-10-18 18:18:11 486

原创 cesium实现绘制线、面贴地效果

此篇在vue里进行cesium整合开发，需要预备vue的开发基础。代码上本人也阅读了官网的例子，借鉴一部分，https://sandcastle.cesium.com/?src=Drawing%20on%20Terrain.html主要思路如下：1）鼠标左键单击事件发生，用两个数组来存点坐标2）鼠标移动事件，利用两个数组长度比较，左键单击事件发生时，两个长度一样，move...

2020-03-09 21:32:01 9209 6

原创 openlayers3在Vue中鼠标右键取消绘制方法

// 本文是介绍右键取消绘制的一种解决方式，以下代码是从其他项目代码中抽取关键部分仅供大致参考，并不是完整的实例代码，请读者注意。//初始化mapinitMap(){ let self=this; ... self.map=addmap(self.center,self.resolutions,self.resolution,s...

2019-10-14 12:56:08 1972

原创 openlayers3之jsts空间分析完整例子

接着上一篇，把上一篇相关的完整代码写出来。<!DOCTYPE html><html> <head> <title>jsts sample</title> <meta charset="utf-8"> <link rel="stylesheet" href="ol3/ol.css" ty...

2019-08-12 11:58:11 726

原创 openlayers3之jsts笔记

jsts是jts的js库的实现，官方找到的帮助都是jts Javadoc帮助文档，其实完全可以参照jts的帮助。一般地，我们只需要在前端读取格式转成geometry即可，那么很多来自于jts的API都是一样的。附上jts的链接http://locationtech.github.io/jts/javadoc/。下面举一个jsts与ol geometry互转以及使用jts api的例子，如此链接下...

2019-08-09 10:22:34 908

原创记geotools java开发坐标转换之一坑

最近接触到了一个关于geotools后端坐标转换的工作，手上想自己动手尝试一下，看着官方文档，一路照做下来，无论是建立Java project还是maven project，一样的代码还是一样的导包还是有问题，起先我以为一定要在maven中写代码，开始是maven的环境配置上出了一点问题，后来环境问题搞了很久才解决，到了后面还是报坐标转换先关的包的异常，于是再仔细地看看官方给的文档，...

2019-07-12 22:52:38 2967

原创线状图斑的平均宽计算

本来我是有可能完全不做三调的相关工作，三调主要还是抠图与地物类型判读占最大比重工作的。扫尾的时候，上面扔给了我一项任务，就是让我赶紧弄出个计算线性图斑宽度的工具。我此时的想法就是提取出中心线，等间距取点作近邻分析，求得平均值*2作为该线状地物的平均宽。与此同时，正好在用FME，发现有一个线性面中心线提取的转换器——centerlinereplacer。不妨试试效...

2018-10-28 18:23:27 5744 10

转载 Java+appium+Android模拟器自动化

基础开发环境Java JDK1.7/1.8，Java开发工具idea、eclipse等安装略 Android studio安装，Android SDK（Android 9.0），模拟器创建虚拟设备1.Android studio下载地址：https://developer.android.google.cn/studio/2.安装完成后，在SDK manager下载对应的SDK模拟器创建一个Android 9的手机点击AVD managerCreate virtual...

2020-11-14 15:14:20 269

原创 Linux笔记文档1

linux使用笔记1解压文件命令:.gz解压1：gunzip FileName.gz解压2：gzip -d FileName.gz压缩：gzip FileName.tar.gz 和 .tgz解压：tar zxvf FileName.tar.gz压缩：tar zcvf FileName.tar.gz DirName———————————————.bz2解压1：bzip2 -d...

2019-07-12 15:22:32 111

原创数据清洗笔记

数据采集:互联网爬虫，外业采集，已有的第三方库数据、部分标准化结构数据，统计资料……数据清洗:RDBMS数据清洗清洗的主要内容：缺失:根据其它字段生成或表连接生成或根据数学模型计算重复:去重，把冗余的数据清理脏数据:文本里有不符合要求的字符需要剔除，如特殊符号(标点符号)，空格，全角半角，乱码等数据挖掘:从已有数据，用数学模型，计算出的结果，常见有：机器学习深度学习 NLP等方法...

2019-04-19 22:49:31 393

原创 kettle环境部署

1．最新下载链接：https://sourceforge.net/projects/pentaho/files/Data%20Integration/7.1/pdi-ce-7.1.0.0-12.zip/download2.没有安装Java jdk的要安装，要求jdk在1.8以上。此步可百度Java jdk下载及安装教程3.下载完成解压，将data-integration文件夹拷到Java ...

2018-12-02 18:48:08 638

原创余弦相似文本算法——python实现

在地址匹配过程中，由于不同的数据源，地理坐标存在不少差异，但我们可以通过文本相似度去比较。越相似越说明两个地址是一样的。文本相似度算法很多，作为刚入门的一员，选择从余弦先入手，余弦通过向量公式计算应该是我们中学时代最熟的一个公式之一，选择它入门相对好理解些。两个文本cos计算从我的算法步骤说主要有4个：1、对两个文本进行分词，分词我用jieba分词的。2、获取两...

2018-10-30 22:46:12 1854 2

原创添加XY离散点数据

arcgis 添加数据里有个添加XY数据可以将txt、CSV、xlx按一定格式存储的坐标数据创建图层事件，对应的工具里有创建 XY 事件图层工具下面将介绍一种比较高级的一种方式，用ArcPy完成。数据是用txt存储的，格式为“X,Y”。这是一种插入游标的方法，向featureclass里写入几何。"shape"是几何字段，用来存储几何的。首先要创建一个要素类，再用游标插入完成。impo

2017-08-06 16:38:43 1983

原创要素内部几何空间查询

featureclass里几何字段（shape）存储有很多的几何，假如我们要的是某一个几何部分作空间查询，现有的工具只有图层查询图层的，没有几何查询图层的。ArcPy为我们提供searchcursor的函数，通过SQL条件语句查找出某部分的几何，再把这个几何用于空间查询是可行的，空间查询用SelectLayerByLocation函数。以下是用矢量格网中的

2017-08-06 16:27:26 1104

原创 ArcGIS 10.2 SDE连接postgresql配置

Arcgis 10.2 Sde连接postgresql：1、安装系统对应的postgre客户端，64位安装64,32位装32 ，10.2只能安装9.1和9.2的2、将arcgis安装位置下的desktop中 DatabaseSupport 中如安装9.1找9.1 window64中st_geometry.dll复制到postgresql安装目录下的lib下3、下载32位postg

2017-07-14 17:42:27 9557 2

原创 C#实现批量高斯投影正算、反算

批量计算有利于提高工作/学习效率，本文以EPSG提供《Coordinate Conversions and Transformations including Formulas》的高斯投影正算、反算算法写成c#代码为例。。 // 高斯投影反算，将高斯坐标反算出经纬度坐标class Gausstojw { private double atanh(d...

2017-07-14 13:21:52 6756

原创利用Python folium库调用OpenStreetMap制作GDP热力图

利用Python的开源模块Folium调用openstrertmap生成html，利用Python爬取的互联网地图的经纬度坐标与统计局公布的年度GDP数据结合，构造data（x,y,gdp）http://blog.csdn.net/qq_14906811/article/details/73717447import osimport foliumimport numpy as npda

2017-07-09 21:34:41 12433 7

原创 ArcPy批量掩膜裁剪栅格/图像

我在实习的时候，根据某些实际的工作需求，从批量裁剪矢量图层想到也想搞一个栅格/影像批量裁剪，在ArcGIS中掩膜裁剪栅格/图像是非常好用的，掩膜裁剪工具很好用，既能矢量裁剪，也能栅格/图像裁剪。实际工作中，可能会面临大批量数据处理，ArcGIS现有的工具箱不能满足要求的情况下，不得不想想构建模型迭代器，或者Python脚本用以处理批量处理的问题了。假如是栅格/影像标准图幅...

2017-07-07 20:54:21 4235

原创批量经纬度转成空间直角坐标系

直接上代码。。import numpy#从键盘自定义输入扁率，长半轴，椭球高（区域中心高）等#f=1/(298.257)f=float(input("f:"))#a=6378140.0a=float(input("a:"))L=0.0e=2*f-f*fB=0.0H=float(input("H:"))f1=open(r'C:\Users\suns\Desktop\POI

2017-06-25 18:53:56 4064 2

原创 Python爬取全国主要城市经纬度坐标

1、首先注册百度或高德地图的开发密钥。2、在Python搞爬虫，你可能需要这几个库，json，urllib，直接pip就行。import jsonfrom urllib.request import urlopen, quoteurl = 'http://api.map.baidu.com/geocoder/v2/'output = 'json'ak = 'UMvGLi

2017-06-25 14:57:56 13337 2

翻译 Rasterio 安装

想来了解更多的安装方式去安装文档阅读 https://pypi.python.org/pypi/rasterio/这里仅介绍以whl方式安装一、whl包准备：gdal、rasteriohttp://www.lfd.uci.edu/~gohlke/pythonlibs二、安装在安装rasterio之前，确保你已经安装了Python，最好是Python3.5

2017-06-04 17:01:39 6706 2

原创基于Rasterio与matplotlib的DEM三维展示

直接上代码了>>> from mpl_toolkits.mplot3d import Axes3D>>> import matplotlib.pyplot as plt>>> import numpy as np>>> import rasterio as rio>>> f=rio.open(r'C:\Users\suns\Desktop\Ex1\dem')>>> nrows

2017-06-04 16:55:41 1078

原创 geopandas简单绘图及横坐标调整

from shapely.geometry import * from geopandas import * import numpy as np import matplotlib.pyplot as plt xmin, xmax, ymin, ymax = 900000, 1080000, 120000, 280000 xc = (xmax - xmin) * np.random.

2017-05-29 19:02:38 8939

原创基于Python 的简单栅格图像边界提取

在GIS中，栅格属性里有关于栅格自身的信息，背景（nodata value）对于识别一张图像的边界像元尤为重要，我们目的只要把每行每列中的第一次出现不是nodata的像元和最后一次出现nodata的前一个像元就可以了。对于栅格，可以用ArcPy中的RasterToNumpyArray函数将将栅格转成numpy数组，然后就可以按照所想读取出每行列中首尾像元。以下是部分代码提取边界像元的核心算法，其实

2017-05-21 11:56:57 6313

翻译 Anaconda 3 4.2.0上安装geopandas

一、首先安装Anaconda3 4.2.0下载链接：https://repo.continuum.io/archive/.winzip/二、安装Geopandas安装教程请参照官方文档：http://geopandas.org/install.html按照官方安装文档的说明，geopandas库依赖：numpy、pandas、fiona、shapely、pyproj等库，在安装g

2017-05-19 18:01:48 12767 8

原创 Python3.4.2 win32 安装 gdal

一、软件准备：1、Python3.4.2 win32下载https://www.python.org/ftp/python/3.4.2/2、gdal-111-1600-core . msi 下载http://download.gisinternals.com/sdk/downloads/release-1600-gdal-1-11-mapserver-6-4/

2017-05-14 14:56:11 2813