遥感
文章平均质量分 69
独孤尚亮dugushangliang
古人之观于天地、山川、草木、虫鱼、鸟兽,往往有得,以其求思之深而无不在也。夫夷以近,则游者众;险以远,则至者少。而世之奇伟、瑰怪,非常之观,常在于险远,而人之所罕至焉,故非有志者不能至也。有志矣,不随以止也,然力不足者,亦不能至也。有志与力,而又不随以怠,至于幽暗昏惑而无物以相之,亦不能至也。然力足以至焉,于人为可讥,而在己为有悔;尽吾志也而不能至者,可以无悔矣,其孰能讥之乎?
展开
-
PIE-Engine使用哨兵1影像提取洪水范围
1 背景洪涝灾害是突发事件,具有持续时间短、危害性大等特征。基于卫星遥感技术,人们可实现快速、准确的洪涝淹没范围提取,直观地显示其空间分布及其动态变化和发展规律,从而使得其在洪涝灾害监测中发挥着越来越重要的作用。光学传感器,如MODIS、Landsat TM等均被用于提取洪水淹没范围。但是洪灾发生时,天气往往多云多雨,光学传感器波段不能穿透云层,因此无法获取有效地面信息。星载合成孔径雷达(SAR)具有全天候、全天时的工作能力和宽覆盖的优势,不受云、雨、雾的影响,可在夜间工作,尤其适合洪涝淹没范围监测。原创 2022-01-26 18:36:18 · 4627 阅读 · 2 评论 -
栅格的二元回归
有时候我们会用到残差趋势法,例如以植被覆盖度为因变量 、以气温和降水为自变量,逐像元建立二元线性回归模型 ,逐像元得到回归方程的系数;其次,利用气温和降水数据以及回归模型的系数,建立模型模拟得到气候影响下的植被覆盖度的预测值;最后,基于遥感影像获得的植被覆盖度观测值与基于回归模型模拟得到气候影响下的预测值做差值计算,得到的结果即为植被覆盖度残差,表示了人类活动对植被覆盖的影响。今天分享一下栅格的二元回归系数计算方法。1 二元回归系数计算二元回归系数系数的计算公式网上书上有计算公式,这里不再赘述。这里介绍原创 2021-12-30 14:46:55 · 3269 阅读 · 3 评论 -
重标极差分析 Hurst指数计算
在时间序列分析中,有时候会用到hurst指数,今天分享Hurst指数的计算方法。1 介绍本节介绍出自《地理数学方法:基础和应用》一书【第 21 章 时间序列的 R/S 分析】 R/S 分析是一种基于长程相关思想的时间序列分析方法。 这种方法由 H. E. Hurst 于 1965 年最先提出,后来伴随着非线性理论的发展而成长起来。 Hurst 原本是剑桥大学物理学博士,对埃及尼罗河( Nile)进行了长达 60 年的观测,记录了尼罗河水位原创 2021-12-30 14:46:11 · 7024 阅读 · 3 评论 -
计算两个栅格的相关系数
有时候我们需要计算两个栅格的相关系数,判断相关性,例如计算NDVI和降水的相关系数,NDVI和温度的相关系数。今天分享一下计算两个栅格相关系数的计算方法。1 相关系数计算相关系数的计算公式网上书上有计算公式,这里不再赘述。这里介绍一下Python的numpy库计算相关系数,使用np.corrcoef()函数,示例如下。import numpy as npx1 = np.array([9.6,17.1,64.8,40.9,136.3,182.5,78.3,3.7,26.,0.4])x2 = np.a原创 2021-12-30 14:45:13 · 8898 阅读 · 6 评论 -
esri土地覆盖全国拼接
1 数据介绍Esri_Land_Cover_2020_10m数据来源是Esri公司,基于10m哨兵影像数据,使用深度学习方法制作做的全球土地覆盖数据。该数据集一共分类十类,分别如下所示:更多介绍参考【10米ESRI土地覆盖数据介绍与下载教程】。2 数据处理首先下载全球的esri土地覆盖文件,在此基础上操作。2.1 找出中国范围的分块2.1.1 生成影像边界框根据影像范围生成影像的边界框shp,便于筛选。python代码如下。 这里使用pyshp创建矢量,pyshp需要是1.2.3版本。原创 2021-12-23 18:43:19 · 1316 阅读 · 0 评论 -
中国分省30米DEM(NASA 2020版)
“中国分省NASA DEM 30米分辨率2020版”01 数据介绍NASA DEM 发布于 2020年2月18日,是美国国家航空航天局(NASA)全新的全球30米分辨率DEM数据(NASA DEM)。是当前最新的一套数据,具有质量好、覆盖范围广等优点。官网地址:https://earthdata.nasa.gov/esds/competitive-programs/measures/nasadem更多介绍可以阅读:数据分享︱全国30 m分辨率DEM数据(ASTER GD...原创 2021-12-05 21:15:29 · 3698 阅读 · 0 评论 -
中国范围NPP-VIIRS逐年夜间灯光数据(2013-2020年)
“本文分享中国范围的NPP VIIRS Annual VNL V2夜间灯光影像”夜间灯光数据最近几年得到了广泛的应用。NPP VIIRS数据是常用的数据。NPP VIIRS数据官网下载地址https://eogdata.mines.edu/products/vnl/。由于Annual VNL V2逐年的数据比较大,下载比较费时,因此我下载了2013-2020年的数据,共有三百多G,并裁剪出了中国区域,方便国内用户获取。01 数据介绍进入官网后找到...原创 2021-12-05 21:12:09 · 7643 阅读 · 6 评论 -
ESA WorldCover和Globeland30土地覆盖数据的比较
1 数据介绍 本文将从几点角度比较一下ESA WorldCover和Globeland30两种土地覆盖产品。 ESA WorldCover是欧空局2021年最新发布的根据哨兵数据处理的全球10米分辨率土地覆盖数据,官方精度74%。分为林地(10),灌木(20),草地(30),耕地(40),建筑(50),裸地/稀疏植被区(60),雪和冰(70),开阔水域(80),草本湿地(90),红树林(95),苔藓(100)。官网地址:https://esa-w原创 2021-11-29 09:27:32 · 4128 阅读 · 0 评论 -
ubuntu20.04安装qgis
一开始使用官网的安装教程https://www.qgis.org/en/site/forusers/alldownloads.html#debian-ubuntu安装失败,报错鉴于仓库 'https://qgis.org/ubuntu focal InRelease' 不支持 'i386' 体系结构,跳不知怎么解决,后来使用了Ubuntugis安装就好了,按照focal过滤,Ubuntugis不但能安装qgis,也能安装其他的gis库例如gdal等,参见https://launchpad.ne..原创 2021-08-22 00:29:26 · 2600 阅读 · 1 评论 -
基于高分辨率影像的潍坊大棚遥感提取
1 背景 朋友,当你坐在经过潍坊的火车上的时候,是否被窗外的风景震惊过?那一望无际的大棚,像一片海洋,又像一片草原,连接天际,一眼望不到边。在卫星地图上,白茫茫一片,蔚为壮观。本文基于0.6米的高分辨率遥感影像提取潍坊的大棚分布。2 数据使用古戈影像的18级数据,0.6米分辨率,共三百多G。3 方法方法和思路都比较简单。使用深度学习的方法,首先选典型样本,其次训练模型,最后应用模型。对于效果不好的区域,再选择一部分样本对模型进行微调。样本的选取是十分重要...原创 2021-08-07 11:54:48 · 2347 阅读 · 2 评论 -
gdal怎样修改图像的投影和地理信息
最近有个需求,给图片定义投影、添加地理坐标信息,今天终于找到了,一起看一看。那就是学习一下gdal_edit.py这个文件geotransform和projection可以是从其他dataset读到。from osgeo import gdalds = gdal.OpenEx(tifpath, gdal.OF_RASTER | gdal.OF_UPDATE)ds.SetGeoTransform(trans)ds.SetProjection(proj)ds = None如果修改noda.原创 2021-06-18 16:01:02 · 913 阅读 · 1 评论 -
NPP VIIRS夜间灯光数据新下载地址
noaa的npp viirs夜间灯光下载原的地址为https://www.ngdc.noaa.gov/eog/viirs/download_dnb_composites.html。看官网发现数据只更新到2019年4月份。仔细观察上面黄底色的字,These products, based on VIIRS observations, are now publicly available thro...原创 2020-03-18 19:11:03 · 13463 阅读 · 14 评论 -
Google Earth Engine导出大影像时下载失败的一种解决方法
Google Earth Engine导出大影像时下载失败最近用gee导出了一些影像到drive云盘中,虽然不是很大,只有两百M左右,但是使用backup and sync同步软件来下载的时候却一直失败,不停的重试。同步软件不支持断点续传,一旦因为某原因中断,又要从头开始,感到十分的无奈经过多次教训,发现下载50M以内的小文件是没有问题的,因此考虑导出数据的时候使用分块来解决,最后下载这些...原创 2020-02-29 18:16:48 · 8386 阅读 · 9 评论 -
arcpy批量将tif影像添加到ArcGIS中同时指定RGB顺序
将一个文件夹里的所有的tif添加到arcgis中指定波段顺序,例如假彩色,真彩色,arcgis默认会使用前3个波段作为RGB显示。运行之前最好先建好金字塔,否则添加会很慢指定拉伸方式没有查到怎么实现# -*- coding:utf-8 -*-'''将一个文件夹里的所有的tif添加到arcgis中指定波段顺序,例如假彩色,真彩色'''import osimport sysr...原创 2020-02-06 11:01:57 · 3661 阅读 · 1 评论 -
Python将多个单波段tif合并为一个影像
将同一文件夹里的多个tif合并为一个tif或dat或其他的影像文件输入影像的行列数投影应保持一致dat格式可以给波段命名,tif格式看不到波段名称# -*- coding:utf-8 -*-import numpy as npimport gdalimport osimport sysreload(sys)sys.setdefaultencoding('utf8')#修改...原创 2020-02-06 09:26:45 · 6845 阅读 · 16 评论 -
arcpy将多个dbf合并到excel中
使用arcpy将多个dbf文件根据提供的字段提取相应的值并保存到txt中,然后再保存到excel,没有某个字段的dbf将设置默认值nodata# -*- coding:utf-8 -*-import arcpyimport osimport sysimport copy# in order to support Chinese languagereload(sys)sys.se...原创 2020-02-03 22:04:27 · 2622 阅读 · 3 评论 -
postgis使用wkt或wkb创建几何
使用wkt创建 ST_GeomFromText使用wkb创建 ST_GeomFromWKB导出为wkt ST_AsText导出为wkb ST_AsEWKBpostgres=# select ST_AsEWKB(ST_GeomFromWKB ('\x0101000020e6100000dfc0e44691145d40e355d636c5fb4340', 4326 ));警告: OGC...原创 2020-01-08 17:59:48 · 2661 阅读 · 0 评论 -
wkt转shp / wkt转shapefile / shp转wkt
使用ogr从wkt字符串创建shapefile输入为一个文本文件,每一行格式为a;b;c;wkt\n,分号分割,abc分别为三个字段的值,wkt为字符串。一行的示例dfas;defg;hjkpiuu;POLYGON((128.99 46.99略去一些))如果多个字段仅需略作修改即可将文本文件的每一行的wkt转为面,并将前3个值作为字段添加到shp中import osrimport ...原创 2020-01-07 17:00:29 · 5413 阅读 · 1 评论 -
Python arcpy检查矢量自相交
arcpy检查面矢量自相交或异相交的问题。基本思路如下:Python脚本如下,需在ArcGIS 10以上版本中运行 ,最后生成__WT.shp的矢量,即为问题矢量:import arcpy#打开覆盖写入arcpy.env.overwriteOutput=TrueA="E:\\zxj\\hh.shp"fold="E:\\zxj\\"index=A.rfind("\\")...原创 2018-07-23 17:36:02 · 4956 阅读 · 0 评论 -
利用影像的元数据文件生成落图的程序(xml/tfw)
有时候需要获取遥感影像的范围,即落图文件。利用高分影像的元数据文件生成则十分简便快捷。根据实际需求,编写了生成落图的程序,主要针对高分的影像。不同的传感器的数据存储的形式也不一定一样,针对相同的存储形式写了相应的函数。基于Python实现。考虑到不同的影像有不同的投影,最后每个xml文件对应一个shp。1 导入需要的库import osimport osgeo.ogr as ogri...原创 2018-07-18 21:42:40 · 2966 阅读 · 2 评论 -
全球30m土地覆盖产品(2010/2015)
最近苦苦寻找30m分辨率土地覆盖/利用而不可得,经袁道友提醒,找到了这个网站,免费的全球30m土地覆盖产品含2010/2015两期,向提供数据的科研工作者们表示敬意。而全国的一年的土地覆盖在地理国情监测云平台上卖24万元。http://data.ess.tsinghua.edu.cn/以2015年为例 ,部分文件如下:下载了80E50N一期,与同年的Landsat8数据比较如下...原创 2018-07-18 21:05:41 · 13906 阅读 · 20 评论 -
Python+gdal读写遥感影像
Python 和gdal读遥感影像,进行相应的波段运算后再写出。读为numpy数组,写出为envi格式的,只能单个波段读,略有缺陷。(部分代码摘自网络,向他们表示敬意)import numpy as npimport gdalimport osdef image(path): #读为一个numpy数组 dataset = gdal.Open(path) b...原创 2018-07-18 14:28:12 · 11745 阅读 · 9 评论 -
Python的Mann-Kendall非参数检验和计算Hurst指数
Mann-Kendall 检验法简称为 M-K 法, 是一种非参数统计检验方法, 可适用于不具有正态分布特征变量的趋势分析[38]。假定X1,X2,...Xn为时间序列变量[1],n为时间序列的长度,M-K 法定义统计量S为其中式中,xj、xk 分别为第j、k年对应的观测值,且j< k。标准化的检验统计量Z为当n≥10时,统计量S近似服从正态分布,在不考虑...原创 2018-07-20 17:14:40 · 14127 阅读 · 4 评论 -
Python熵权法确定权重
熵权法赋权是一种客观赋权方法, 在一些评价中, 通过对熵的计算确定权重, 就是根据各项评价指标值的差异程度, 确定各评价指标的权重。详细介绍及计算公式可参考文献[1]。主要步骤包括(1)原始数据矩阵进行标准化(2)定义熵(3)定义熵权。具体步骤也可参考https://blog.csdn.net/wangh0802/article/details/53981356。这里不再赘述。这里通过Pyth...原创 2018-07-26 20:34:17 · 17093 阅读 · 11 评论 -
Python gdal/ogr cookbook
http://pcjericks.github.io/py-gdalogr-cookbook/利用Python读写矢量、栅格数据的教程网站。原创 2018-07-30 15:58:37 · 923 阅读 · 0 评论 -
Python tkinter读查看遥感影像投影
Python tkinter读查看遥感影像投影结果如图 点击右上角选择,选择一个影像: 再次点选择按钮,选一个其他文件: 选择一个没有投影的:选择一个没有投影的影像: 退出。代码如下:# -*- coding: utf-8 -*-"""Created on Sun Aug 5 19:14:52 2018@author: Administra...原创 2018-08-05 19:59:07 · 1073 阅读 · 0 评论 -
Python opencv对多个相同大小图片裁剪出同一区域
arcgis同一个模板出的图,对每个裁剪出同一片区域,使得裁剪后的大小仍是相同的。word 、powerpoint似乎不太好实现。用Python和opencv很容易实现。import cv2img=cv2.imread(u'lun/aaa.jpg')#裁剪区域为行数的范围,列数的范围,波段数的范围img1=img[1900:3414,266:2364,:]cv2.imwrite(u'...原创 2018-08-21 15:27:57 · 2509 阅读 · 0 评论 -
Arcpy自动修补/填补polygon面的空洞/孔洞/洞/空白
使用arcpy填补孔洞,思路是读取所有坐标但排除空洞的坐标再创建矢量,测试使用三个图斑,一个没有洞,一个有一个洞,一个有三个洞,测试数据如下:程序如下:使用时需修改自己的输入和输出路径,不过属性都没了,如需以前的属性,可以再通过相交、连接等工具补充。import arcpydata=[]#infc为输入面矢量infc = "kong"outfc = "D:/wcg/test/...原创 2018-12-23 14:18:58 · 4632 阅读 · 0 评论 -
arcpy计算当前地图文档的每个数据框每个图层的面积
计算当前地图文档的每个数据框每个图层的面积,基于数据源的坐标系import arcpyimport osimport sysreload(sys)sys.setdefaultencoding('utf-8')#设置输出文件outputfile="D:\\aaa\\output.txt"#获取当前地图文档下的数据框mxd = arcpy.mapping.MapDocument(...原创 2019-01-07 21:42:58 · 1564 阅读 · 1 评论 -
Python做函数拟合
读一幅影像指定位置的一些点与另外一幅图相同位置的点做拟合,包括一次函数、二次函数、指数函数、幂函数、对数函数,计算R方,并绘制散点图,在图上显示保存为图片。# -*- coding: utf-8 -*-"""Created on Mon Jul 23 14:40:18 2018@author: Administrator"""import gdal import osimp...原创 2018-07-23 15:12:11 · 10196 阅读 · 0 评论