04 2022 档案
摘要：//设置要下载的范围 var roi_China=ee.Geometry.Rectangle(73,3,135,52); Map.centerObject(roi_China, 4); Map.addLayer(roi_China, {}, 'China'); //获取矢量数据集，使用filter来
阅读全文
摘要：//设置要下载影像的区域，生成roi var roi_China=ee.Geometry.Rectangle(73,3,135,52); //将生成的下载区域roi在地图上显示 Map.centerObject(roi_China, 4); Map.addLayer(roi_China, {}, '
阅读全文
摘要：1）GetTransform()返回6个参数，分别为： (左上角x坐标[0]， 水平分辨率[1]，旋转参数[2]， 左上角y坐标[3]，旋转参数[4]，竖直分辨率[5]) 遥感影像如果没有发生旋转，即上北下南，则旋转参数[2]与旋转参数[4]均为0 返回的参数对应的是影像当前坐标系下的坐标值，比如W
阅读全文
摘要：不同矩阵a,b对应位置元素运算 1）合并矩阵：merge=np.array([a,b]) 1）求平均np.mean(merge,axis=0) 2）求标准差np.std(merge,axis=0) 3）求差值，直接a-b 4）求分数，直接a/b 关于axis运算方向的问题，在看了网上的一些解析与实验
阅读全文
摘要：#导入支持库 import sys #读取数组所占内存，单位为MB sys.getsizeof(array)/1024/1024
阅读全文
摘要：获取数据类型位深,由于获取可能出现不同大小写的表述方法，而python数据格式为小写，所以使用lower()统一成小写 DataType=gdal.GetDataTypeName(dataset.GetRasterBand(1).DataType).lower()
阅读全文
摘要：1）栅格数读取 from osgeo import gdal import numpy as np #打开栅格数据集，只读 dataset=gdal.Open("E:/DATA/Image.TIF",gdal.GA_ReadOnly) #读取图像信息 #数据格式 print("Driver: {}/
阅读全文
摘要：1）批量处理中，使用os获取批量文件名称 如文件夹E:DATA下有此6个文件夹 import os #根目录路径 rt_file_path="E:/DATA" #获取根目录路径下文件/文件夹名称 rt_file_list=os.listdir(rt_file_path) Out： ['2001',
阅读全文
摘要：使用gdal写入栅格之后，应使用 out_dataset=None 来关闭数据集（这也是官方帮助文档使用的方法） 如果使用 out_dataset.FlushCache( ) 可能会存在问题
阅读全文
摘要：问题起因：gdal GetProjection()返回空值 原因：GDAL找不到proj库 解决方法：设置环境变量
阅读全文
