中误差是衡量观测精度的一种数字标准,亦称“标准差”或“均方根差”。在相同观测条件下的一组真误差平方平均值的平方根。

因真误差不易求得,所以通常用最小二乘法求得的观测值改正数来代替真误差。它是观测值与真值偏差的平方和观测次数n比值的平方根。

中误差的大小反映了该组观测值精度的高低,因此,通常称中误差为观测值的中误差。

**1.问题:**用中误差检验影像精度

2.已知:

3.关键思路:

(1)影像转换位数组

(2)excel数据转换成二维数组

4.过程:

import math#计算两个数组的均方根均方根差
def rmse(records_real, records_predict):if len(records_real) == len(records_predict):return math.sqrt(((records_real - records_predict) ** 2).mean())else:return None
import numpy as np
import xlrd#excel表格数据转换为二维数组
def excel(input_excel):table = xlrd.open_workbook(input_excel).sheets()[0]  # 获取第一个sheet表row = table.nrows  # 行数col = table.ncols  # 列数data = np.zeros((row, col))  # 生成一个nrows行ncols列,且元素均为0的初始矩阵for x in range(col):cols = np.matrix(table.col_values(x))  # 转换为矩阵data[:, x] = cols  # 按列把数据存进矩阵中return data

import sys
from PIL import Imageif __name__ == '__main__':reload(sys)sys.setdefaultencoding("utf-8")inputfile = "I:\\PycharmProject\\ceshi.tif"input_excel = "I:\\PycharmProject\\ceshi.xls"##1.输入单波段的待检验影像,并把影像灰度值转换为数组arrarr = np.array(Image.open(inputfile))##2.读取观测值表为数组brrbrr =excel(input_excel)##3.利用arr与brr计算均方根差rmse_val = rmse(arr, brr)print(rmse_val)return rmse_val

参照文章:python例用中误差检验精度

python利用中误差检验影像精度相关推荐

  1. 利用Python爬取全国250m精度的人口数据、房价数据等数据 | CSDN博文精选

    作者 | 中原百科来源 | CSDN博客 (一) 我的第一篇博客写的就是爬取人口数据基于腾讯位置大数据平台的全球移动定位数据获取(Python爬取),精度是1000m,后来有朋友和我说有个网站开放过2 ...

  2. 利用Python爬取全国250m精度的人口数据

    此次以GeoQ(智图)为基础,利用Python爬取全国250m精度的人口数据(GeoQ)这个网站开放过250m精度的人口分布数据,而且人口分布有年龄分段等属性.先得注册登录到达创建地图的界面. 看人口 ...

  3. 利用Python爬取全国250m精度的人口数据(GeoQ)、房价数据和公交站(线路)等数据

    前言 文的文字及图片来源于网络,仅供学习.交流使用,不具有任何商业用途,版权归原作者所有,如有问题请及时联系我们以作处理. 作者: 中原百科 GIS大师兄 PS:如有需要Python学习资料的小伙伴可 ...

  4. 利用python进行T检验

    引入所需的包 from scipy import stats import numpy as np 注:ttest_1samp, ttest_ind, ttest_rel均进行双侧检验 H0:μ=μ0 ...

  5. python构建t检验(Student’s t-test)

    python构建t检验(Student's t-test) t检验,亦称student t检验(Student's t test),主要用于样本含量较小(例如n < 30),总体标准差σ未知的正 ...

  6. python利用什么写模板_Python利用逻辑回归分类实现模板

    Logistic Regression Classifier逻辑回归主要思想就是用最大似然概率方法构建出方程,为最大化方程,利用牛顿梯度上升求解方程参数. 优点:计算代价不高,易于理解和实现. 缺点: ...

  7. python从random生成列表_详解Python利用random生成一个列表内的随机数

    详解Python利用random生成一个列表内的随机数 首先,需要导入random模块: import random 随机取1-33之间的1个随机数,可能重复: random.choice(range ...

  8. python/gdal处理遥感影像(读取、投影转换、裁剪、建立图像金字塔等)

    python/gdal处理遥感影像(读取.投影转换.裁剪.建立图像金字塔等) gdal库简单介绍 python使用gdal 一.安装python环境 二.安装gdal库 三.使用gdal处理遥感影像 ...

  9. python进行 t 检验

    t检验主要是针对正态总体均值的假设检验,即检验样本的均值与某个值的差异,或者两个样本的均值是否有差异等. 其不需要事先知道总体的方差,并且在少量样本情况下也可以进行检验. python进行t检验使用s ...

最新文章

  1. Quartz学习笔记
  2. 简单实现AJAX: ASP.NET2.0 中回调的实现及常见问题的解决
  3. 给NavigationCtrl 增强动画.
  4. python中斐波那契数列_斐波那契数列–在Python,JavaScript,C ++,Java和Swift中进行了解释...
  5. Liunx之nginx代理
  6. @scheduled注解配置时间_Spring Boot中使用@Scheduled创建定时任务
  7. 曾创下IRR62%的超高战绩的VC离开中国:一切因人而来因人而去
  8. 海词词典android v3.1.2新版发布 英语学习必备工具,海词词典Android V3.1.2新版发布 英语学习必备工具...
  9. linux打开txt文件命令_linux系统文件及常用命令
  10. ad16怎么画弧线_ad15软件中如何画圆弧 多点画圆弧软件
  11. 风场可视化与原理剖析
  12. 微信小程序api基础应用
  13. 关于Android如何集成QQ登录及分享
  14. 保护接地 功能接地 等电位接地
  15. 怎么把word转换成excel表格的格式
  16. Ubuntu零基础教学-GParted磁盘分区工具使用|超级详细,手把手教学
  17. Integer源码详解
  18. mall商城微服务遇见的坑
  19. 西浦计算机系怎么样,西交利物浦大学
  20. python怎么运用函数_如何在python中使用step函数

热门文章

  1. Netherlands central buffer of UPS european delivery center
  2. 《建筑工程定额与预算》
  3. 图像评价指标——PSNR、SSIM、FID、KID
  4. Win10下用SCP命令免密码上传、下载阿里云服务器上的文件
  5. WSL挂载移动硬盘U盘
  6. DellR720服务器上安装EXSI6.5.0全教程
  7. 如何规避苹果审查,实现app store上的app版本强制更新
  8. lisp注册注册机源码_[求助]请教在lisp中加入注册码
  9. ln x的matlab表示,matlab中ln怎么表示
  10. AutoCad 修改多段线