![](https://img-blog.csdnimg.cn/20201014180756923.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
Python地理数据处理
文章平均质量分 74
见贤思齐547
这个作者很懒,什么都没留下…
展开
-
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 · 3492 阅读 · 15 评论 -
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 · 3255 阅读 · 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 · 475 阅读 · 0 评论 -
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 · 1993 阅读 · 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 · 2217 阅读 · 0 评论 -
Python计算DEM(tif格式)坡度和坡向
slopeaspect将背景值设为-300slopeArcgis原创 2021-03-26 16:35:02 · 7288 阅读 · 10 评论 -
PythonGDAL给图像设置投影坐标系
网站上下载的部分遥感数据只有地理坐标,没有投影坐标,本文章尝试给“只有地理坐标的图片”添加投影坐标。原始数据:山东地区ASTER GDEM30m数据,该数据只有WGS84地理坐标目标数据:给该数据添加UTM-50N1.确定你的目标投影坐标信息,就可以通过下述代码实现。(1)首先你得知道想添加的目标投影信息,可以通过网上查询的方式。此网站查的比较费劲,我反正没在这找到:https://developers.arcgis.com/javascript/3/jshelp/pcs.htm,该原创 2021-03-22 18:44:39 · 7630 阅读 · 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 · 5807 阅读 · 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 · 3583 阅读 · 3 评论 -
Python实现栅格数据拼接、裁剪等操作
想看看山东的高程如何。从地理空间数据云(http://www.gscloud.cn/)上下载了GDEMV2—30米分辨率的数据。原创 2021-03-02 08:33:25 · 2872 阅读 · 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 · 2013 阅读 · 1 评论 -
python地理数据处理相关的操作
ogrinfo:打印矢量图层的信息ogr2ogr:矢量数据格式转换gdalinfo:打印栅格图层的信息gdal_translate:栅格数据格式转换gdaladdo:对栅格数据集建立pyramidsgdalwrap:栅格数据集的投影变换gdal_merge:栅格图像的拼接...原创 2021-02-25 20:52:53 · 1594 阅读 · 0 评论 -
利用Python把遥感影像的某几个波段合成
刚刚开始想着手学习Python,决定从最基础的读写遥感影像开始。学习借鉴了网上很多前辈们的经验,自己出现了一些小问题写在这里,算是记录一下。这是USGS上下载的一景Landsat8影像,地点在青海。用ENVI截取了扎陵湖的一部分,真彩色合成图如下:想利用Python把nir,red,green三个波段合成假彩色,代码如下:import osimport numpy as npfrom osgeo import gdalclass IMAGE: # 读图像文件原创 2021-02-05 10:57:03 · 5196 阅读 · 0 评论 -
点线面矢量数据创建
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 · 643 阅读 · 0 评论 -
将一副遥感影像裁剪为多幅图像
参考: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 · 836 阅读 · 1 评论 -
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 · 5321 阅读 · 0 评论 -
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 · 556 阅读 · 0 评论