- 博客(29)
- 资源 (2)
- 收藏
- 关注
RSS订阅原创 电脑端OneNote2019（win10系统）登录微软账号登陆错误报错误码0x80860010错误-已解决
之前登录都没有问题，今天登录点击同步就无法登录了，手机端也不显示更新的笔记。出现以下这样的错误，尝试了两个方法，法一不行，法二解决了。1.法一：重置oneNote此方法非常简单，两张截图：开始菜单——OneNote——应用设置。重置这个方法对我不管用。2.法二：右击：我的电脑，找到网络和Internet。删除删除高级选项——重置——确定win+r 快捷键调出运行窗口，输入services.msc，查看 微软服务是否正常启动。 如...
2022-02-14 16:07:31
4104
1
原创 numpy.ndarray size changed, may indicate binary incompatibil
在笔记本电脑运行程序时报错：numpy.ndarray size changed, may indicate binary incompatibility.
2021-05-31 16:58:13
3342
原创 python地理数据可视化
参考1.知乎：https://www.zhihu.com/question/33783546/answer/7759464012.GitHub上代码下载：https://github.com/nancygaooo/Visualization-Practice
2021-05-26 17:58:42
2294
原创 Basemap安装
1.下载编译好的安装包Basemap无法使用pip install basemap进行安装，可以下载编译好的安装包，网址：https://www.lfd.uci.edu/~gohlke/pythonlibs/。在basemap安装之前需要先安装pyproj（我的电脑之前已经安装过了）。我的python版本是3.9的，所以选择 basemap-1.2.2-cp39-cp39-win_amd64.whl 。2.安装在cmd中把路径定位到上述安装包所在的文件夹，然后输入：pip in...
2021-05-26 10:15:11
2185
原创 .ipynb文件的使用问题
下载了一些文件，后缀名为.ipynb，没见过这种格式的，学习一下。一、关于.ipynb文件文件全称为ipython notebook，是由jupyter插件生成的文件。 可以在线编辑在线运行，在浏览器中打开二、.ipynb文件的使用，需要在Windows下安装.ipynb文件1.pip 安装（1）打开cmd输入：pip install jupyter安装的东西挺多的：（2）pip list查看2.使用命令进入命令行输入：jupyter notebook
2021-05-25 10:55:57
1406
原创 pip的安装与卸载
一、安装1.去官网（https://pypi.org/project/pip/#files）下载安装包2.解压到某一文件夹，如： E:\1_RuanJianBao\PythonRelevant\pip-21.1.2.tar\dist\pip-21.1.2\pip-21.1.23.在cmd命令行中输入（1）E：（2）在 E:\1_RuanJianBao\PythonRelevant\pip-21.1.2.tar\dist\pip-21.1.2\pip-21.1.2 下输...
2021-05-25 10:36:35
1077
原创 python对landsat8数据进行辐射校正
Landsat8 TI级数据已经过几何校正和地形校正，还需要做辐射校正和大气校正。1.辐射定标（1）代码倾向于TIF格式的影像，所以output format选择的是“GTIFF”格式。from osgeo import gdalfrom osgeo import gdal_arrayimport numpy as npfrom show import TwoPercentLinearfrom matplotlib import pyplot as pltimport cv2 a
2021-05-25 09:27:32
3348
14
原创 AtmosphericCorrection大气校正_landsat8
有关大气校正的代码在github上下了一些：https://github.com/zhaoguanhua/atmosphericcorrection/这里有多种数据（高分、landsat8、哨兵等）的大气校正代码。本文只是测试了landsat8：1.代码这是主文件（加了一些函数、方法的注解），还有部分文件就不一一显示了，在（https://github.com/zhaoguanhua/atmosphericcorrection/）直接下载即可。# glob模块参考： https://blog.
2021-04-20 10:10:21
3055
8
原创 Py6S模块安装过程（Win10环境，64位）
想了解一下用python做大气校正的过程，查了一下资料，需要安装Py6S模块。下面是自己的一个安装过程，仅供参考。一、前期准备工作1.安装MinGW软件（1）下载安装mingw-get-setup软件参考：https://blog.csdn.net/pdcxs007/article/details/8582559双击出现：（2）安装mingw软件参考：https://www.cnblogs.com/LIJIH/p/12533926.html（3）配置环境变
2021-04-15 13:56:22
1926
8
原创 python查看我国1990-2015年间的温度变化情况
偶然间看到一个博文（https://blog.csdn.net/kimol_justdo/article/details/113774790），正好手里有1990-2015年的年均气温TIF图，也想尝试一下。下图是所有文件：1.基本情况介绍下图是“中国1980年以来逐年年平均气温空间插值数据集”（详见资源环境科学与数据中心：http://www.resdc.cn/data.aspx?DATAID=228） ...
2021-04-12 16:07:25
437
原创 DEM渲染
一、Python实现1.二、ArcGIS方式实现1.山东地区的DEM图在ArcGIS的ArcToolbox中操作如下：按照顺序ArcToolbox -> 3DAnalyst Tools - > Raster Surface -> Hillshade，打开生成山体阴影的对话框，输入数据源选择DEM图层，其他参数默认，如图：结果如下图所示：参考：1.麻辣GIS：https://malagis.com/a...
2021-04-02 10:11:07
1919
3
原创 Python计算坡度坡向并输出二维、三维图
在上次代码的基础上做了一点儿修改，将定义的函数单独放在一个模块里面，主函数去单独调用该模块。DEMslopeAspect模块from osgeo import gdal,ogr,osrimport numpy as npimport mathimport datetime# Python matplotlib模块代码示例 https://vimsky.com/examples/detail/python-module-matplotlib.html# Axes3D是mpl_toolki
2021-03-29 15:45:35
2145
原创 PythonGDAL给图像设置投影坐标系
网站上下载的部分遥感数据只有地理坐标，没有投影坐标，本文章尝试给“只有地理坐标的图片”添加投影坐标。原始数据：山东地区ASTER GDEM30m数据，该数据只有WGS84地理坐标目标数据：给该数据添加UTM-50N1.确定你的目标投影坐标信息，就可以通过下述代码实现。（1）首先你得知道想添加的目标投影信息，可以通过网上查询的方式。此网站查的比较费劲，我反正没在这找到：https://developers.arcgis.com/javascript/3/jshelp/pcs.htm，该
2021-03-22 18:44:39
7469
9
原创 gdal.Buildvrt和gdal.Warp工具实现TIF影像拼接、矢量裁剪
看到GDAL库有很多工具可以直接使用，避免写长串的代码。这里算是记录一下自己的学习过程吧~比如gdalwarp，李老师的博客中写到的：（https://liminlu.blog.csdn.net/article/details/8531979）以上是用命令行的方式实现，自己敲了几行代码在pycharm中实现（毕竟不太习惯使用命令行的方式）：from osgeo import gdal,gdalnumericfrom osgeo.gdalconst import *# 查看DEM图的属
2021-03-15 17:28:33
5440
5
原创 Python批处理MODIS数据并计算NDVI
1.登录NASA官网下载MOD13Q1数据，红框标出来的是筛选条件，我的筛选条件列出符合要求的文件如下：NASA官网：https://ladsweb.modaps.eosdis.nasa.gov/MODIS数据的介绍：https://www.cnblogs.com/cuteshongshong/articles/3622855.html2.利用指定的MRT工具对MODIS数据进行批处理MRT下载地址：链接：https://pan.baidu.com/s/1aqD4UAhPQAWq83zq
2021-03-04 15:20:18
3388
3
原创 Python实现栅格数据拼接、裁剪等操作
想看看山东的高程如何。从地理空间数据云（http://www.gscloud.cn/）上下载了GDEMV2—30米分辨率的数据。
2021-03-02 08:33:25
2781
5
原创 利用Python中的GDAL和OGR模块实现shapefile对栅格DEM数据的裁剪
之前已经把31景DEM数据拼接完成，这篇文章，利用shapefile文件实现对山东DEM数据的裁剪。参考：https://blog.csdn.net/GISShiXiSheng/article/details/72843070?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522161423692216780271562057%2522%252C%2522scm%2522%253A%252220140713.130102334.pc%255Fb
2021-02-26 17:44:30
1952
1
原创 python地理数据处理相关的操作
ogrinfo：打印矢量图层的信息ogr2ogr：矢量数据格式转换gdalinfo：打印栅格图层的信息gdal_translate：栅格数据格式转换gdaladdo：对栅格数据集建立pyramidsgdalwrap：栅格数据集的投影变换gdal_merge：栅格图像的拼接...
2021-02-25 20:52:53
1540
原创 Python读取指定文件夹下指定类型数据的文件名并保存到TXT文件中
一、读取路径到指定TXT文件有31个DEM文件夹，每个文件夹下有3个文件：ASTGTM2_N34E115_dem.tif，ASTGTM2_N34E115_num.tif，README.tif。想读取每个DEM文件夹下的ASTGTM2_N34E115_dem.tif文件。代码如下：import os# 查询符合条件的文件f = open(r'D:\ProfessionalProfile\VectorRelevant\GDEMV2_30m_shandong\pathDEM.txt','w
2021-02-24 10:52:22
1141
原创 TypeError: cannot unpack non-iterable int object查找指定文件夹下指定文件类型的数量
想统计一下指定文件夹下，指定文件类型的数量，出现了这个错误：TypeError: cannot unpack non-iterable int object（见下图）找了一下是因为第14行出错，将 j,k=0 改为j,k=0，0就可以了。解释：报错的原因是函数返回值的数量不一致，查看函数返回值数量和调用函数时接收返回值的数量是不是一致，修改一致即可。import os# 查询符合条件的文件j,k=0def file(root, ext): for i in os.l...
2021-02-24 09:55:26
5191
原创 Python地理数据处理相关模块的安装方法
主要模块有下面几个：1.GDAL模块在此网站下载（https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal）与Python版本相对应的版本，我的Python版本是3.9.0，下载的文件是 GDAL-3.2.1-cp39-cp39m-win_amd64.whl ，39代表的python版本是3.9.x。然后，利用：pip installGDAL-3.2.1-cp39-cp39m-win_amd64.whl 命令安装。可以通过 “pip list”...
2021-02-23 17:46:26
524
原创 点线面矢量数据创建
from osgeo import ogrimport matplotlib.pyplot as pltfrom ospybook.vectorplotter import VectorPlotter# 一、构建几何类型:点point = ogr.Geometry(ogr.wkbPoint)p1=point.AddPoint(3, 3) # 创建点x, y = point.GetX(), point.GetY() # 画图vp = VectorPlotter(True) ...
2021-02-23 10:54:44
604
原创 ospybook安装包及安装过程
Python学习的过程中用到了ospybook库，安装了一下。ospybook-1.0 的下载地址：https://www.manning.com/books/geoprocessing-with-python?a_bid=c3bae5be&a_aid=geopy，百度云地址：链接：https://pan.baidu.com/s/1SPxndlTgeSARgBKl1bbbZw 提取码：1slx下载后，安装：1.打开cmd，定位setup.py目录，打开命令行运行 ：...
2021-02-22 20:23:02
1145
原创 求助：矢量数据FID字段无法读取问题
今天尝试读取矢量数据，遇到了一些小问题，有的自己没能解决，记录下来。如有大神偶尔看到，帮助一下小白，不胜感激。一直不太理解“驱动”，找了一些资料算是了解一些，见文章：https://blog.csdn.net/shaxiaozilove/article/details/79313515以下代码也参照了好多大神的资料：参考：https://blog.csdn.net/zsc201825/article/details/91909356参考：https://www.cnblogs.com/l..
2021-02-22 15:31:08
1819
5
原创 将一副遥感影像裁剪为多幅图像
参考：https://blog.csdn.net/zsc201825/article/details/89359995以Landsat8某景影像为例：原图（底图）是假彩色的，上面单波段图像为裁剪的9景影像。代码如下：from osgeo import gdalimport math# 读取要裁剪的原始遥感影像in_ds = gdal.Open("D:\ProfessionalProfile\DataRelevant\L134036_20170808.tif")print("op
2021-02-20 11:29:40
810
1
原创 【已经解决】TypeError: read_img() missing 1 required positional argument: ‘filename‘
前两把读取遥感数据、存储遥感数据函数写成了一个类；然后再去调用这个类，出现了这个问题（TypeError: read_img() missing 1 required positional argument: 'filename'）：后来查询了一下，类需要实例化，完整代码如下：from classIMAGE import IMAGEimport numpy as npimage = IMAGE() #实例化类的对象# proj, geotrans, data1, row, co
2021-02-07 15:38:21
4328
3
原创 利用Python把遥感影像的某几个波段合成
刚刚开始想着手学习Python，决定从最基础的读写遥感影像开始。学习借鉴了网上很多前辈们的经验，自己出现了一些小问题写在这里，算是记录一下。这是USGS上下载的一景Landsat8影像，地点在青海。用ENVI截取了扎陵湖的一部分，真彩色合成图如下：想利用Python把nir,red,green三个波段合成假彩色，代码如下：import osimport numpy as npfrom osgeo import gdalclass IMAGE: # 读图像文件
2021-02-05 10:57:03
4899
支持向量机软件包
2018-02-05
支持向量机
2018-02-05
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人