遥感
文章平均质量分 58
星梦云庄
A doctor without technology
展开
-
SAR 学习资料
声明:以下内容引用,记录一下,以免找不到。转载 2023-05-17 09:11:03 · 287 阅读 · 0 评论 -
ArcPy批量合并矢量shape文件
Python+ArcPy批量合并矢量shape文件原创 2023-02-21 22:50:55 · 1441 阅读 · 0 评论 -
matlab同时画n多条曲线,设置颜色及图例
matlab同时画n多条曲线,设置颜色及图例转载 2023-02-13 21:30:07 · 3962 阅读 · 0 评论 -
使用ArcPy对栅格影像和矢量文件进行批量投影
在ArcGIS中,可以利用ArcToolbox中的数据管理下的投影工具对栅格影像数据和矢量文件进行投影转换,也可以进行批量投影操作。本代码和ArcGIS实现的功能一样,都是对栅格影像和矢量文件进行批量投影,所调用的工具是Project_management 和ProjectRaster_management 函数。代码如下:#!/usr/bin/python# -*- coding: utf-8 -*-import osimport arcpyarcpy.CheckOutExten..原创 2022-01-25 12:52:06 · 4804 阅读 · 7 评论 -
最新数据处理 之 栅格数据批量归一化
网上下载的数据集,有时候需要进行归一化处理,在 ArcGIS 中逐一处理比较费事,可以考虑用 ArcPy 进行批量归一化。代码如下:# -*- coding: utf-8 -*-import osimport arcpyfrom arcpy import envfrom arcpy.sa import *arcpy.CheckOutExtension("Spatial")env.workspace = "D:/Download/"save_path = "D:/Save/"for i原创 2021-12-28 14:16:32 · 3307 阅读 · 0 评论 -
最新数据处理 之 最新DSM(AW3D30)数据批量拼接--文末附数据获取方式
本文主要针对JAXA发布的最新版本3.1&3.2(AW3D30)数据进行数据拼接处理。JAXA利用ALOS的PRISM立体模式制作了全球的DSM数据,并免费分布了30米版本产品。详细的数据介绍和下载:AW3D30与DEM数据一样,数据是按照格网大小进行分发的,如 N028E090.zip,就是28N, 90E区域的DSM数据。数据下载完成后,解压出来,里面包括了 "_DSM.tif","_MSK.tif","_STK.tif" 三种tif格式的数据,这里 DSM 数据为例介绍批量拼接。原创 2021-12-23 13:43:03 · 3817 阅读 · 39 评论 -
最新数据处理 之 1:100万与1:25万全国基础地理数据库批量合并
背景数据来源:全国地理信息资源目录服务系统链接:全国地理信息资源目录服务系统这个网站是由自然资源部主管, 国家基础地理信息中心进行运行维护的,因此数据是比较权威的。目前,网站上提供了以下三类数据: 30米全球地表覆盖数据,数据整体现势性达到2010年 1:100万全国基础地理数据库,数据整体现势性达到2015年 1:25万全国基础地理数据库,数据整体现势性达到2015年 个人可以进行注册,然后查询和下载,下载的时候,需要提供单位和下载的目的。提交后,很快就可以进行下原创 2021-12-21 17:52:55 · 2682 阅读 · 10 评论 -
深度学习、遥感领域常用精度评价指标盘点
最近几年,以深度学习技术为代表的新型前沿技术在图像处理领域,取得了巨大的成功,从图像处理到场景理解,从图像分类到图像实时分割,深度学习无处不在。在遥感领域,深度学习技术更是火热。本篇博文就对常用的深度学习精度评价指标进行总结梳理,旨在为模型算法的精度评价指标选取提供一个参考。1、精确率、召回率和F1精确率(Precision),又可以称为正确率,主要是指识别出来的目标有多少是识别准确的,衡量的是识别结果的正确率;召回率(Recall),又可以成为查全率,主要是指所有正确的目标有多少被识别出来原创 2020-11-14 15:40:15 · 9060 阅读 · 0 评论 -
使用Python批量下载哨兵1号数据
写在前面好久没更新了,不是忘了,是懒!正好最近想下载哨兵1号数据来反演一下地表形变,于是乎,就撰写此文记录一下坎坷的路程。首先在Open Access Hub (copernicus.eu)网站上搜寻了一番,选了一些数据(此处略去)。选的数据很多,要是一个一个点击下载,实在太过于麻烦,直觉告诉我肯定有自动批量下载的方式,于是bing了一下,果然有。大概浏览了一下,无非就那么几种,要么是通过https://search.asf.alaska.edu/这个网站上自动生成的python...原创 2021-09-01 16:12:20 · 1892 阅读 · 6 评论 -
我国高分系列卫星遥感影像介绍
继上一篇介绍《遥感图像处理》的文章之后,本篇文章对我国的高分系列卫星遥感影像进行简单的整理。高分系列卫星是在高分专项的支持下,也就是高分辨率对地观测系统重大专项,由国防科技工业局牵头,组织实施建设的一系列高分辨率对地观测卫星,肩负着我国民用高分辨率遥感数据实现国产化的使命。高分专项,于2010年批准启动实施,到2020年,已经覆盖了从全色、多光谱到高光谱,从光学到雷达,从太阳同步轨道到地球同步轨道等多种类型,构成了一个具有高空间分辨率、高时间分辨率和高光谱分辨率能力的对地观测系统。高分一号原创 2021-02-03 21:40:07 · 10221 阅读 · 7 评论 -
遥感图像处理
2021年第一篇文章基本概念遥感(Remote Sensing):遥远的感知,是一种远距离不直接接触物体而获取其信息的探测技术,主要是以电磁波为媒介,包括从紫外--可见光--红外--微波的范围。DN值(Digital Number ):遥感影像像元亮度值,记录地物的灰度值。地表反射率(Surface Albedo):地面反射辐射量与入射辐射量之比,表征地面对太阳辐射的吸收和反射能力。反射率越大,地面吸收太阳辐射越少;反射率越小,地面吸收太阳辐射越多。表观反射率(Apparent Refl原创 2021-02-01 14:52:16 · 10369 阅读 · 4 评论 -
如何利用Python对DEM数据进行显示(对之前代码的更正),附带利用Python计算坡度
1、更正---图像显示代码在之前的某一篇文章中,《如何利用Python对遥感影像进行显示》,阐述了如何利用Python对遥感影像进行显示,最近利用代码对DEM数据进行显示时,输入DEM数据,显示的结果却是全黑的,仔细看了一下之前的代码,发现了一个小的BUG,如下:这句有点问题,当 (im_max-im_min) 值比较大,大于255时,255/(im_max-im_min) 结果就被python默认为int型,即0。这样,后面乘以任何数字都为0,导致整个显示图像都为0。因此,该部分代码应更正为:原创 2020-12-26 17:12:35 · 3063 阅读 · 20 评论 -
如何利用ArcPy批量裁剪栅格(包括外边框和实际边界范围)
最近2周,由于项目上事情很多,自己又偶感风寒,博客没有及时更新。今天稍微好一些,更新安排上。今天,遇到一个问题:有很多范围大小不一样的矢量文件,如何根据这些矢量的范围裁剪对应的栅格影像?这个问题,利用ArcGIS软件的Toolbox可以很好的解决,但批量的话,操作起来有点费力,所以这里写了一个批量裁剪的代码。用到的库主要是ArcPy,但需要在ArcGIS自带的IDLE中运行,也可以在ArcMap界面中运行。批量裁剪代码如下:import sys,osimport arcpy..原创 2020-12-07 16:55:43 · 2073 阅读 · 4 评论 -
利用Pyinstaller打包代码时出现小错误的解决方案
最近几天,由于项目需求,需要对代码进行打包封装,但打包完后,运行时,出现了一些小错误,这里进行简单的介绍一下。在之前的博文中《如何在Windows10环境下,对Tensorflow代码进行打包封装(利用PyQt5制作界面)》已经介绍了如何利用PyQt制作界面,如何利用Pyinstaller对代码进行封装打包,这里不作介绍。这篇文章主要对打包过程中存在的部分比较隐秘的小问题进行汇总。问题一:GDAL命令行下运行成功,但打包后运行出错之前的博文已经介绍了一些利用GDAL库进行读取影像、创建矢量、矢栅原创 2020-11-24 11:47:19 · 1818 阅读 · 4 评论 -
如何使用Python中的GDAL库循环写入矢量数据
该文章为《如何使用Python中的GDAL库创建矢量文件》后续在《如何使用Python中的GDAL库创建矢量文件》文章中,介绍了利用GDAL库创建点、线、面矢量文件,文章中只介绍了单个数据的创建情况,那么面对有很多个数据时,该怎么创建呢?问题来源:现有一些点坐标,目的是将其构建成一个点矢量文件,直接上代码。单个点数据情况:from osgeo import ogr,osr Longitude = 116.4 Latitude = 39.9 def create_point():原创 2020-11-22 17:40:31 · 900 阅读 · 3 评论 -
如何利用Python对遥感影像进行显示
遥感影像一般具有多个波段,比较常见的影像一般是4波段多光谱影像,比如高分一号、高分二号、资源三号等。这些影像数据一般体量较大,有的几百兆,有的多达几十G,格式一般是16位无符号整型,一般看图软件无法打开显示,需要ArcGIS、ENVI等专业的软件进行查看,有时候很不方便。这篇博客就简单的介绍一下,如何利用Python对遥感影像进行显示,需要用到的库为GDAL和Opencv。正文Python中,一般的绘图库都无法处理或显示遥感影像,例如Matplotlib、Opencv、Scipy等,很多连遥原创 2020-11-20 21:39:31 · 9039 阅读 · 25 评论 -
如何使用Python中的GDAL库进行矢量数据与栅格数据相互转换
今天继续更新GDAL相关内容,主要是进行矢量数据与栅格数据的相互转换。矢栅转换在ArcGIS软件中很好实现,有专门的工具可以使用,这里主要讲述如何利用Python中的GDAL库完成矢栅转换,直接上代码。矢量转栅格from osgeo import ogr,gdaldef shape_to_raster(shapefile, rasterfile, savefile): data = gdal.Open(rasterfile, gdal.GA_ReadOnly) x_res =原创 2020-11-19 12:44:39 · 2543 阅读 · 12 评论 -
如何使用Python中的GDAL库创建矢量文件
继《如何使用Python中的GDAL库对遥感影像进行读取和存储》之后,本文再简单地介绍一下,如何使用Python中的GDAL库创建矢量文件。这里,矢量文件指的是点、线、面文件,shape格式,可以使用 ArcGIS 等软件读取编辑。需要用到的Python库依然为GDAL,安装完GDAL后,就可以直接导入使用了。废话不多说,直接上代码。创建点文件from osgeo import ogr,osrLongitude = 116.4 Latitude = 39.9def creat原创 2020-11-18 13:33:06 · 3618 阅读 · 3 评论 -
如何使用Python中的GDAL库对遥感影像进行读取和存储
与自然影像稍有不同,遥感影像一般有多个波段(波段数大于4),数据格式多为tiff,比如国产高分一号/二号遥感影像,多光谱影像有4个波段,分别为蓝(1)、绿(2)、红(3)以及近红外(4)波段,真彩色是321显示,假彩色是432显示。利用遥感或GIS处理软件,如ENVI、ERDAS、ArcGIS等,可以很方便快捷地实现对图像进行读取、处理等操作。然而,很多时候我们需要对大量的图像数据进行批量处理操作,这时候就需要我们自己编写代码来实现。这篇博客就简单地介绍,如何在Python中实现遥感影像的读取和存储原创 2020-11-17 14:03:29 · 14501 阅读 · 32 评论