YuWanwanyayaya-CSDN博客

转载 R语言十倍交叉验证

theta.fit <- function(x,y){lsfit(x,y)} #最小二乘拟合结果。cat(“Original R-square =”,r,“\n”) #将多个参数连接起来。#对数据states所有预测变量进行回归，然后用shrinkage（）函数做10重验证。library(bootstrap) #crossval函数实现k重交叉验证。

2023-04-18 22:55:30 555

BandMath工具能够方便的执行图像中各个波段的加、减、乘、除、三角函数、指数、对数等数学函数计算，也可以使用IDL编写的函数。BandMath工具使用函数都是基于IDL的数据组运算符。IDL的数组运算符使用方便且功能强大。它们可以对图像中的每一个像元进行单独检验和处理，而且避免了FOR循环的使用（不允许在波段运算中使用）。数组运算符包含关系运算符（LT、LE、EQ、NE、GE、GT）、Boolean运算符（AND、OR、NOT、XOR）和最小值、最大值运算符（<、>）。这些特殊的运算符对图

2021-04-19 09:57:37 6708

原创 ArcGIS与EXCEL结合，对度分秒进行转换

一 Arcgis可以提取矢量面的中心点坐标打开ArcMap软件，添加面矢量数据，先为面统计面积，转换成点后该字段会保留，点文件不再能统计面积在工具栏中选择“ArcToolbox”按钮，打开ArcToolbox工具箱，在其中选择“数据管理工具”——“要素”——“要素转点”，选择输入要素，设置输出要素“内部”为可选选项，勾选时转点一定在面内，不勾选则可能出现在面外转换结果如图所示二统计点的经纬度信息3. 为点图层添加X 和Y两个字段，分别计算经纬度4. 右击X-计算几何-点的

2021-04-14 15:45:45 6311

转载 ArcGIS属性表批量替换

转载def func(name): ##附表00湿地 if name==('S#-2#=t'.decode('utf-8')): a='S2t' return a.decode('utf-8') else: return name1、简述实验课上，老师让我们依据出台的《国土空间调查、规划、用途管制用地用海分类说明》更改原有三调数据分类名中的属性字段值，如上图所示，原始的三调数据中的分类名：冰川及永久积累，依据新出台的文件《国土空间调查、规划、用途管制用地用海分类说明》，要求把其

2021-04-14 15:04:03 6099

转载 Python批量计算NDVI

PYTHON批量计算NDVI做了少量修改，剔除了异常值，执行代价时需要更换影像对应波段及文件存储位置import osimport numpy as npfrom osgeo import gdalimport globlist_tif = glob.glob("F://2020TreeClassification//NDVIPY//MASKGF//*.tif")out_path = 'F://2020TreeClassification//NDVIPY//NDVIGF'for tif

2020-12-12 17:05:08 2641 5

原创 Python批量计算Pearson相关系数和RMSE

import mathimport xlrdfrom sklearn.metrics import mean_squared_error # 用于评估模型workbook = xlrd.open_workbook('E:\Pycharm\PycharmLearning\RMSE_Test\\RMSE_20190415.xlsx', 'r')sheet1 = workbook.sheet_by_index(0)GF_B1 = sheet1.col_values(5)[1:] # 第6列为变换

2020-10-21 10:48:11 2632 1

YuWanwanyayaya的博客