GIS
关注
分享
扫一扫
文章平均质量分 64
卡尔曼和玻尔兹曼谁曼
面朝大海，春暖花开
展开
-
GEE导出图像到本地结果全部为空
今天在使用Google Earth Engine处理数据进行导出为GeoTIFF到Google云盘的时候，发现下载下来以后的图像值全部为空（NAN）。我尝试将结果加载在GEE的Code Editor提供的在线地图上进行显示，发现结果可以正常显示，图像都是有值的。后来我对图像的数据类型进行修改，发现导出以后还是没值。再后来我尝试在导出函数中设置CRS参数，导出结果正确。...原创 2022-08-01 12:11:50 · 1537 阅读 · 0 评论 -
QGIS中WKT转为可视化涂层
文章目录QGIS中WKT转为可视化涂层常见的几种用于互操作的空间数据标准格式使用QGIS的Python接口将WKT转为可视化涂层版权声明：本文为博主原创文章，转载请注明原文出处！作者：阿振写作时间：2021-01-21 19:00:05QGIS中WKT转为可视化涂层常见的几种用于互操作的空间数据标准格式WKT（Well-Known Text）是开放地理空间联盟OGC（Open Geospatial Consortium ）制定的一种文本标记语言，用于表示矢量几何对象及其空间参照系统等。W原创 2021-01-21 19:46:47 · 1675 阅读 · 0 评论 -
GIS空间分析之Clip
文章目录开篇案例介绍代码展示方法总结版权声明：本文为博主原创文章，转载请注明原文出处！作者：阿振写作时间：2020-06-16 周一 早安开篇接着上篇《GIS空间分析之Buffer》之后，这篇来介绍矢量数据的裁剪（Clip）。裁剪其实算不上是一种空间分析功能，就是一种简单矢量数据处理操作。一般，我们会使用一个代表研究区域或者感兴趣的区域（AOI，Area of Interest）的多边形数据去裁剪一个比较大的全区域的数据，得到的结果就是感兴趣区域中包含的数据。矢量数据的裁剪在平时工作中用得不原创 2020-06-16 08:48:39 · 2416 阅读 · 2 评论 -
GIS空间分析之Buffer
文章目录开篇案例介绍代码展示方法总结版权声明：本文为博主原创文章，转载请注明原文出处！作者：阿振写作时间：2020-06-14 周末 夏夜开篇GIS空间分析是通过对GIS系统中的空间地物的空间位置以及分布形态等空间特性进行分析推理等得到额外有用信息的过程。GIS空间分析包含广泛的内容，是GIS系统的核心功能。从这篇博文开始，我们会简单介绍几种GIS系统中最常见最简单的空间分析功能，并使用Python的GDAL API进行实现。首先，我们来看一下开源GIS软件QGIS中提供的几种简单的空间分原创 2020-06-14 23:12:15 · 5913 阅读 · 0 评论 -
矢量数据空间查询
文章目录开篇案例一案例说明代码演示方法总结案例二案例说明代码演示方法总结版权声明：本文为博主原创文章，转载请注明原文出处！作者：阿振写作时间：2020-06-14 周天开篇在前面四篇博客中我们主要讲了对于空间矢量数据的属性数据的增删改查，在这篇博文中我们要讲解空间查询–GIS系统很重要的一项功能。空间查询就是根据地物的空间位置进行查询的一种数据检索方式。比如，我们要查询一条河流经的城市；一个公园内的所有路灯；离当前位置最近的公共卫生间等等都属于常用的空间查询。OGC简单要素规范定义了空间几何原创 2020-06-14 12:26:42 · 2429 阅读 · 0 评论 -
Shapefile属性操作之查
文章目录开篇案例介绍SQL查询方式代码展示方法总结遍历Feature要素方式代码展示方法总结版权声明：本文为博主原创文章，转载请注明原文出处！作者：阿振写作时间：2020-06-13 夜开篇在前面几篇博文中，我们分别介绍了矢量空间数据的属性数据的增（CREATE）删（DELETE）改（UPDATE）操作，这篇博文我们来聊聊属性数据的查询（Retrieve）操作。案例介绍我们还是使用之前的分省Shapefile数据，主要包含每个省的一些统计信息。下面以两个案例进行介绍：从给定数据中查询中原创 2020-06-13 23:40:59 · 1696 阅读 · 0 评论 -
Shapefile属性操作之改
文章目录开篇案例介绍代码展示方法总结版权声明：本文为博主原创文章，转载请注明原文出处！作者：阿振写作时间：2020-06-13开篇在上篇的《Shapfile属性操作之增》和《Shapefile属性操作之删》中我们分别介绍了对于空间矢量数据属性的增加（CREATE）和删除（DELETE）操作，这篇我们聊聊数据的更新操作（UPDATE）。案例介绍这里我们要处理的数据是中国地图分省的矢量Shapefile，是一个面状数据。该数据有一个NAME字段，给出了每个省的名称。这里我们要更新NAME字段，原创 2020-06-13 17:49:28 · 1927 阅读 · 0 评论 -
Shapefile属性操作之删
文章目录开篇案例介绍代码展示方法总结版权声明：本文为博主原创文章，转载请注明原文出处！作者：阿振写作时间：2020-06-13 又是一个周末开篇延续上篇的《Shapfile属性操作之增》，这篇我们来聊聊如何进行属性删除。使用的工具依旧是我们的GDAL库。案例介绍这里我们要处理的数据是中国地图分省的矢量Shapefile，是一个面状数据。在上篇中我们给该数据添加了一个属性字段Abbr用以表示省的简称。这篇我们再把该字段给删除掉。删除的方法包括两部：首先，从属性表中找到该字段，然后删除该字段原创 2020-06-13 17:34:57 · 1891 阅读 · 0 评论 -
Shapefile属性操作之增
文章目录开篇案例介绍代码展示方法总结版权声明：本文为博主原创文章，转载请注明原文出处！作者：阿振写作时间：2020-06-13 又是一个周末开篇《Python空间数据处理实战》系列的博文好久都没有更新了，今天乘周末有点时间，补了个觉，然后写几篇博文。关于Python空间数据处理，如果大家有什么想看到的内容，欢迎评论区留言，我会增加一些大家需要的内容！今天我打算用四篇博客简要介绍一下如何对空间矢量数据的属性数据进行操作。对于属性数据，我们可以简单将其看作一个二维表格，学过数据库的朋友，可能会原创 2020-06-13 17:21:42 · 1481 阅读 · 2 评论 -
QGIS制图中面积小的区域不显示注记
文章目录QGIS制图中面积太小的区域不显示注记版权声明：本文为博主原创文章，转载请注明原文出处！作者：阿振写作时间：2020-05-24 周天QGIS制图中面积太小的区域不显示注记在使用QGIS进行制图的过程中，对于面积太小的区域有可能存在注记显示不出来的情况。比如在中国地图中，香港和澳门区域面积较小，就存在显示不了注记的情况。如下图，就会发现澳门没有显示出来。[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-A25XDtjl-1590280272837)(/im原创 2020-05-24 08:31:49 · 3871 阅读 · 0 评论 -
QGIS面数据融合以后有小的线段或者洞存在如何解决
文章目录QGIS面数据融合以后有小的线段或者洞存在如何解决版权声明：本文为博主原创文章，转载请注明原文出处！作者：阿振写作时间：2020-05-24 周末QGIS面数据融合以后有小的线段或者洞存在如何解决在QGIS中我们使用Dissolve工具进行面数据的融合，如下图对中国分省的矢量面数据进行融合得到国界[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-clYEtjNC-1590280206719)(/images/QGIS/屏幕快照2020-05-24上午7.15原创 2020-05-24 08:31:11 · 4546 阅读 · 0 评论 -
SNAP Java API处理Sentinel-1数据
版权声明：本文为博主原创文章，转载请注明原文出处！作者：阿振写作时间：2020年4月14日 周二使用SNAP JAVA API处理Sentinel-1数据SNAP软件使用Java语言开发，提供了Python接口snappy，官方教程中也多以Python接口进行示范。但是我在使用Python接口过程中，发现并不是很好用，你必须要同时懂Java语言才能很好地使用Python接口，在IDEA...原创 2020-04-14 10:05:47 · 2852 阅读 · 1 评论 -
Sentinel数据处理工具包SNAP Python开发环境搭建
文章目录Sentinel数据处理工具包SNAP Python开发环境搭建安装问题缺包问题JDK版本问题环境变量问题版权声明：本文为博主原创文章，转载请注明原文出处！写作时间：2020年12月4日周末Sentinel数据处理工具包SNAP Python开发环境搭建这篇博文主要探索安装SNAP工具包并且使用Python接口进行开发过程中搭建开发环境所踩的坑。不得不说欧空局SANP官方提供的...原创 2020-04-12 12:23:47 · 2759 阅读 · 0 评论 -
ESA SNAP工具包Java接口的使用
版权声明：本文为博主原创文章，转载请注明原文出处！写作时间：2010年4月10日晚ESA SNAP工具包Java接口的使用欧盟的SNAP工具包提供了对Sentinel卫星数据的处理功能。SANP使用Java编写，UI界面使用了NetBeans框架。如果要进行批处理的话，还提供了基于图的Graph Processing Framework (GPF)处理模式。此外，SNAP还贴心的提...原创 2020-04-10 23:00:15 · 2366 阅读 · 2 评论 -
macOS下GDAL Java开发环境搭建
文章目录macOS下GDAL Java开发环境搭建GDAL源码编译安装Maven安装本地JAR使用Java版GDAL示例版权声明：本文为博主原创文章，转载请注明原文出处！写作时间：2020年4月9日下午macOS下GDAL Java开发环境搭建今天在笔记本（macOS平台）上折腾了一下如何搭建GDAL的Java开发环境。虽然GDAL的Python接口更好用，但是有时候需要协同使用一些J...原创 2020-04-09 17:50:57 · 2381 阅读 · 5 评论 -
手把手教你使用QGIS制作地图
手把手教你使用QGIS制作地图@[TOC](手把手教你使用QGIS制作地图)手把手教你使用QGIS制作地图加载矢量数据加载背景底图美化矢量数据添加晕线切换到排版视图添加经纬度格网添加其他修饰元素地图输出版权声明：本文为博主原创文章，转载请注明原文出处！写作时间：2020年7月4日晚手把手教你使用QGIS制作地图QGIS是一款开源免费的地理信息系统软件，虽然比不上商业的ArcGIS软件，但...原创 2020-04-07 23:41:55 · 42069 阅读 · 12 评论 -
Fiona简介及Shapefile数据读取
Fiona简介及Shapefile数据读取作者：阿振邮箱：tanzhenyugis@163.com博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375修改时间：2018-06-06声明：本文为博主原创文章，转载请注明原文出处Fiona简介用GDAL的Python绑定API书写程序有没有一种仍然在写...原创 2018-06-07 12:00:17 · 7839 阅读 · 1 评论 -
使用Rasterio做投影变换
使用Rasterio做投影变换作者：阿振邮箱：tanzhenyugis@163.com博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375修改时间：2018-06-11声明：本文为博主原创文章，转载请注明原文出处思路分析在之前GDAL系列文章中的《栅格数据投影转换》提到过，做投影转换最重要的是计算...原创 2018-06-12 11:36:57 · 3388 阅读 · 1 评论 -
使用Fiona创建Shapefile矢量数据
使用Fiona创建Shapefile矢量数据作者：阿振 邮箱：tanzhenyugis@163.com博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375修改时间：2018-06-10声明：本文为博主原创文章，转载请注明原文出处基本思路使用Fiona写入Shapefile数据，主要是构建一个Schem...原创 2018-06-12 10:56:31 · 2854 阅读 · 0 评论 -
矢量数据投影转换
矢量数据投影转换作者：阿振邮箱：tanzhenyugis@163.com博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375修改时间：2018-06-03声明：本文为博主原创文章，转载请注明原文出处案例说明接着上一篇博文中，我们得到了WGS84坐标系下的中国省区图，而我们一般中国地图中使用的是割圆...原创 2018-06-04 06:16:55 · 7432 阅读 · 1 评论 -
创建Shapefile文件并写入数据
打开Shapefile文件的正确方式作者：阿振邮箱：tanzhenyugis@163.com博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375修改时间：2018-06-02声明：本文为博主原创文章，转载请注明原文出处基本思路使用GDAL创建Shapefile数据的基本步骤如下：使用osge...原创 2018-06-03 11:02:34 · 10432 阅读 · 2 评论 -
读取HDF或者NetCDF格式的栅格数据
读取HDF或者NetCDF格式的栅格数据作者：阿振邮箱：tanzhenyugis@163.com博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375修改时间：2018-05-17声明：本文为博主原创文章，转载请注明原文出处HDF和NetCDF简介HDFHDF（Hierarchical Da...原创 2018-05-26 09:24:05 · 6721 阅读 · 7 评论 -
栅格数据创建与保存
栅格数据创建与保存作者：阿振邮箱：tanzhenyugis@163.com博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375修改时间：2018-05-24声明：本文为博主原创文章，转载请注明原文出处思路与方法使用Python进行栅格数据处理，很多时候，我们会将GDAL的Dataset对象转化为N...原创 2018-05-26 02:23:27 · 4460 阅读 · 3 评论 -
地图投影
地图投影作者：阿振邮箱：tanzhenyugis@163.com博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375修改时间：2018-04-29声明：本文为博主原创文章，转载请注明原文出处图片来源于网络，如有侵权请私信删除什么是地图投影我们的地球是圆的，而我们的纸张是平面。为了将地球绘...原创 2018-04-30 03:26:57 · 3091 阅读 · 0 评论 -
使用Rasterio读取栅格数据
使用Rasterio读取栅格数据作者：阿振 邮箱：tanzhenyugis@163.com 博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375 修改时间：2018-06-06 声明：本文为博主原创文章，转载请注明原文出处Rasterio简介有没有觉得用GDAL的Python绑定书写的代码很不Pyth...原创 2018-06-08 03:34:40 · 12643 阅读 · 2 评论 -
栅格数据投影转换
栅格数据投影转换作者：阿振邮箱：tanzhenyugis@163.com博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375修改时间：2018-06-01声明：本文为博主原创文章，转载请注明原文出处使用GDAL提供的命令行工具进行转换GDAL提供了gdalwarp命令可以方便地让我们进行影像拼接，...原创 2018-06-01 23:24:15 · 15936 阅读 · 5 评论 -
打开Shapefile文件的正确方式
打开Shapefile文件的正确方式作者：阿振邮箱：tanzhenyugis@163.com博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375修改时间：2018-05-25声明：本文为博主原创文章，转载请注明原文出处Shapefile文件简介Shapefile文件是美国ESRI公司发布的文件格式...原创 2018-06-02 01:53:57 · 14794 阅读 · 0 评论 -
使用Python的netCDF4和matplotlib.basemap包进行气象数据的可视化
首先，安装netCDF4和Basemap，Windows下和Linux下会稍微有些不一样，请自行百度。 netCDF4：https://github.com/Unidata/netcdf4-python Basemap：http://matplotlib.org/basemap/users/installing.html下面是一个示例代码：from netCDF4 import Dataset原创 2016-03-04 20:51:34 · 30962 阅读 · 15 评论 -
Morton码
版权声明：本文为博主原创文章，转载请注明原文出处！写作时间：2019-07-08 17:11:27Morton码的计算Morton码是对栅格格网进行编码的一种算法，在Google中搜索Morton，搜索结果第一位是Wikipedia的Z-order Curve，这是因为Morton码编码结果展现为一种Z形的填充曲线。下面简要说一下如何计算四进制和十进制的Morton码。四进制Morto...原创 2019-07-08 17:36:23 · 18086 阅读 · 0 评论 -
栅格数据裁剪
文章目录使用GDAL命令使用Python代码版权声明：本文为博主原创文章，转载请注明原文出处！写作时间：2019-03-22在进行遥感影像处理的时候，我们经常需要进行裁剪的工作，来看看如何使用GDAL工具进行这项操作吧！参考资料：GDAL: gdalwarpGDAL: gdal_translateGDAL/OGR Python API使用GDAL命令GDAL提供了两个命令可以...原创 2019-03-23 03:57:13 · 3862 阅读 · 0 评论 -
使用GDAL读取Sentinel数据
使用GDAL读取Sentinel数据GDAL 2.1已经原生支持对于Sentinel数据的读取，我这里使用Sentinel-2光学卫星数据给出使用GDAL工具对其进行读取的方法。这里我们要大概知道Sentinel数据的组织。下载下来的Sentinel数据是一个ZIP压缩包，里面包含了JPEG2000格式的影像数据以及一些XML格式的元数据文件。GDAL将Sentinel数据看做一个数据...原创 2018-08-29 22:50:17 · 4242 阅读 · 15 评论 -
使用Rasterio创建栅格数据
使用Rasterio创建栅格数据作者：阿振 邮箱：tanzhenyugis@163.com博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375修改时间：2018-06-09声明：本文为博主原创文章，转载请注明原文出处方法描述使用Rasterio创建并写入栅格数据比GDAL还简单一些，基本使用到两个函数...原创 2018-06-10 10:26:08 · 3661 阅读 · 2 评论 -
GDAL简介
GDAL简介作者：阿振邮箱：tanzhenyugis@163.com博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375修改时间：2018-05-13声明：本文为博主原创文章，转载请注明原文出处Geospatial Data Abstraction Library （GDAL）是使用C/C++语言编写的用...原创 2018-05-14 07:30:52 · 23572 阅读 · 3 评论 -
空间数据
空间数据作者：阿振 邮箱：tanzhenyugis@163.com 博客：https://blog.youkuaiyun.com/theonegis/article/details/80089375 修改时间：2018-05-06 声明： 本文为博主原创文章，转载请注明原文出处矢量数据（Vector）矢量数据基于对象模型（object-based）的空间数据描述模型。矢量数...原创 2018-05-07 11:45:21 · 2046 阅读 · 0 评论 -
GDAL获取栅格数据各个像素对应的经纬度（Python版）
这里要使用的一个比较重要的函数是GetGeoTransform函数，GDAL官网的API中对GetGeoTransform的解释如下： CPLErr GDALDataset::GetGeoTransform ( double * padfTransform ) virtual Fetch the affine transformation coefficients.Fetches原创 2016-03-04 21:11:16 · 26494 阅读 · 20 评论 -
Neo4j Spatial数据导入
首先，安装neo4j数据库。我的开发环境是Ubuntu，安装过程参考官网：Neo4j Debian Packages，安装后配置：Post-installation tasks。然后安装Neo4j Spatial的插件，参加Github介绍：neo4j-contrib/spatial。即把下载的文件解压到安装目录的plugins目录下$NEO4J_HOME/plugins。接下来，我们使用Java代原创 2016-10-01 14:38:02 · 3332 阅读 · 1 评论 -
MongoDB导入Shapefile数据
两种解决方案： 一、将整个shapefile转为GeoJSON然后直接导入mongoDB数据库中 首先，将shapefile数据转为WGS84地理坐标，然后使用GDAL的命令行工具ogr2ogr进行格式的转换，转换命令如下： ogr2ogr -f geoJSON continents.json continents.shp 删除生成JSON文件的前两行{ "type": "FeatureCo原创 2016-10-01 14:06:35 · 3518 阅读 · 0 评论 -
Java GDAL开发环境搭建
首先从官网下载编译好的二进制包 GISInternals Support Site。 解压以后的文件如下： 将该目录添加到自己的环境变量Path中，比如我的路径是：E:\C++\Library\gdal\release-1800-gdal-2-1-0-mapserver-7-0-1\bin在Eclipse中新建一个Java工程。将下载好的GDAL中bin\gdal\java下面的文件拷贝到工原创 2016-09-19 14:35:47 · 9462 阅读 · 9 评论 -
ShapeFile数据到mongodb的导入
开发环境为： 系统环境 Linux gisweb2 4.4.0-36-generic #55~14.04.1-Ubuntu SMP Fri Aug 12 11:49:30 UTC 2016 x86_64 x86_64 x86_64 GNU/Linux mongodb版本 当前最新版本3.2.9具体mongodb的安装参照官方文档，强烈建议参照官方安装文档。网上博客这种资料良莠不齐，而且新版本可能原创 2016-09-17 21:38:01 · 2834 阅读 · 3 评论 -
根据天数计算具体日期
问题描述：Landsat数据的命名规范中对于日期是：年份+该年的第几天 具体参见：Landsat File Naming Convention 那如果通过天数计算具体的日期呢？ 下面给出计算源码（Python版本）： Python3下运行通过def is_leap_year(year): # 判断闰年的方法是该年能被4整除且不能被100整除，或者是可以被400整除 if ((原创 2016-09-03 17:17:53 · 4085 阅读 · 1 评论