摘要: 1)栅格数读取 from osgeo import gdal import numpy as np #打开栅格数据集,只读 dataset=gdal.Open("E:/DATA/Image.TIF",gdal.GA_ReadOnly) #读取图像信息 #数据格式 print("Driver: {}/ 阅读全文
posted @ 2022-04-03 10:39 MisakaYier 阅读(187) 评论(0) 推荐(0) 编辑
摘要: 1)批量处理中,使用os获取批量文件名称 如文件夹E:DATA下有此6个文件夹 import os #根目录路径 rt_file_path="E:/DATA" #获取根目录路径下文件/文件夹名称 rt_file_list=os.listdir(rt_file_path) Out: ['2001', 阅读全文
posted @ 2022-04-03 09:29 MisakaYier 阅读(108) 评论(0) 推荐(0) 编辑
摘要: 使用gdal写入栅格之后,应使用 out_dataset=None 来关闭数据集(这也是官方帮助文档使用的方法) 如果使用 out_dataset.FlushCache( ) 可能会存在问题 阅读全文
posted @ 2022-04-02 09:48 MisakaYier 阅读(104) 评论(0) 推荐(0) 编辑
摘要: 问题起因:gdal GetProjection()返回空值 原因:GDAL找不到proj库 解决方法:设置环境变量 阅读全文
posted @ 2022-04-01 16:27 MisakaYier 阅读(529) 评论(0) 推荐(0) 编辑
摘要: 数组a 1)检查数组各元素是否为NA,NA返回TRUE,非NA返回FALSE is.nan(a) 2)检查数组个各元素是否为Inf,Inf返回TRUE,非Inf返回FALSE is.infinite(a) 3)确定NA或Inf元素在数组中的位置 which(is.na(a)==TRUE)) whic 阅读全文
posted @ 2022-03-24 16:07 MisakaYier 阅读(812) 评论(0) 推荐(0) 编辑
摘要: #通过向量建立时间序列格式,frequency是频率,比如季度为1年4季,月为1年12月,start为开始的时间 Timeseries<-ts(Vector,frequency=4,start=c(2012,1)) #展示时间序列 tsdisplay(Timeseries) 阅读全文
posted @ 2022-03-18 15:42 MisakaYier 阅读(686) 评论(0) 推荐(0) 编辑
摘要: 1)R语言读取csv,txt数据 #读取txt文件,header文件首行是否包含变量名,sep表示文件中的分隔符,col.names列名称 col_names<-c("col1","col2","col3") data_txt<-read.table("traindata.txt",header=T 阅读全文
posted @ 2022-03-18 15:37 MisakaYier 阅读(2478) 评论(0) 推荐(0) 编辑
摘要: 1)下载程序包 2)解压之后,使用控制台进入程序所在目录 输入python setup.py install 如果使用Anaconda,则使用Anaconda的控制台输出命令 阅读全文
posted @ 2022-03-18 14:34 MisakaYier 阅读(54) 评论(0) 推荐(0) 编辑
摘要: 用pandas读取,读取之后为DataFrame格式 import pandas as pd traindata=pd.read_csv("traindata.csv") 阅读全文
posted @ 2022-03-18 11:19 MisakaYier 阅读(41) 评论(0) 推荐(0) 编辑
摘要: 使用for循环向向量中循环添加计算之后的变量: #建立空向量vector vector<-c() m=1 #向vector中循环添加m的计算结果 for (i in 1:10) { m=m+5 vector<-c(vector,m) } #Out vector=6 11 16 21 26 31 36 阅读全文
posted @ 2022-03-18 11:00 MisakaYier 阅读(288) 评论(0) 推荐(0) 编辑
点击右上角即可分享
微信分享提示