技术标签: python 图像16位转8位 机器学习 深度学习 人工智能 opencv
不同于一般的8位数据,数据库中的遥感图像都是16位的无符号整形数据,而且大部分集中在[0,500]的范围内,对于8位的RGB图像来说根本无法显示正确的图像。因此,必须对灰度级进行压缩,从65536个灰度级压缩到256个灰度级范围内。
但是,仅仅使用线性压缩是无法解决问题的,因为前面也说过,数据主要集中在[0,500]范围内,即它只占用了前面的几百个灰度级,后面的根本没用到。为此,我们必须丢掉出现概率小的灰度级,保留出现概率大的灰度级范围。
****1.读16位图像数据
def read_mat(input_file):
"""
读取影像
Inputs:
input_file:16位图像数据的路径
Outputs:
array_data:保存图像数据的numpy.array数组
rows:高度
cols:宽度
bands:深度
"""
result_data = scio.loadmat(input_file)
dataset= result_data['imgPS']
cols = dataset.shape[0]
rows = dataset.shape[1]
bands = dataset.shape[2]
array_data = np.zeros((rows, cols,bands))
for i in range(bands):
band = dataset[:,:,i]
array_data[:, :,i] = band
return array_data, rows, cols, bands
****2.直方图统计
def cumulativehistogram(array_data, band_min, band_max):
"""
累计直方图统计
Inputs:
array_data:16位图像数据
band_min:最小像素值
band_max:最大像素值
Outputs:
cutmax:累积分布函数(CDF)==0.98对应的灰度级
cutmin:累积分布函数(CDF)==0.02对应的灰度级
"""
# 逐波段统计最值
rows=array_data.shape[0]
cols=array_data.shape[1]
gray_level = int(band_max - band_min + 1)
gray_array = np.zeros(gray_level)
for row in range(rows):
for col in range(cols):
gray_array[int(array_data[row, col] - band_min)] += 1
count_percent2 = rows*cols * 0.02
count_percent98 = rows*cols* 0.98
cutmax = 0
cutmin = 0
for i in range(1, gray_level):
gray_array[i] += gray_array[i - 1]
if (gray_array[i] >= count_percent2 and gray_array[i - 1] <= count_percent2):
cutmin = i + band_min
if (gray_array[i] >= count_percent98 and gray_array[i - 1] <= count_percent98):
cutmax = i + band_min
return cutmin, cutmax
****3.灰度压缩
def compress(origin_16, output_8):
"""
Input:
origin_16:16位图像路径
Output:
output:8位图像路径
"""
array_data, rows, cols, bands = read_mat(origin_16)
compress_data = np.zeros(( rows, cols,bands))
for i in range(bands):
band_max = np.max(array_data[:, :,i])
band_min = np.min(array_data[:, :,i])
cutmin, cutmax = cumulativehistogram(array_data[:, :,i], band_min, band_max)
compress_scale = (cutmax - cutmin) / 255
for j in range(rows):
for k in range(cols):
if (array_data[j, k,i] < cutmin):
array_data[j, k,i] = cutmin
if (array_data[j, k,i] > cutmax):
array_data[j, k,i] = cutmax
compress_data[j, k,i] = (array_data[j, k,i] - cutmin) / compress_scale
write_mat(output_8, compress_data)
****5.生成8位图像
def write_mat(output_file,array_data):
output_file_order=["B","G","R","NIR"]
rows=array_data.shape[0]
cols=array_data.shape[1]
bands=array_data.shape[2]
for band in range(bands):
cv2.imwrite(output_file+output_file_order[band]+".tif", array_data[:,:,band])
文章浏览阅读3.2k次。本文研究全球与中国市场分布式光纤传感器的发展现状及未来发展趋势,分别从生产和消费的角度分析分布式光纤传感器的主要生产地区、主要消费地区以及主要的生产商。重点分析全球与中国市场的主要厂商产品特点、产品规格、不同规格产品的价格、产量、产值及全球和中国市场主要生产商的市场份额。主要生产商包括:FISO TechnologiesBrugg KabelSensor HighwayOmnisensAFL GlobalQinetiQ GroupLockheed MartinOSENSA Innovati_预计2026年中国分布式传感器市场规模有多大
文章浏览阅读1.1k次,点赞2次,收藏12次。常用组合逻辑电路结构——为IC设计的延时估计铺垫学习目的:估计模块间的delay,确保写的代码的timing 综合能给到多少HZ,以满足需求!_基4布斯算法代码
文章浏览阅读3.3k次,点赞3次,收藏5次。OpenAI Manager助手(基于SpringBoot和Vue)_chatgpt网页版
文章浏览阅读2.2k次。USACO自1992年举办,到目前为止已经举办了27届,目的是为了帮助美国信息学国家队选拔IOI的队员,目前逐渐发展为全球热门的线上赛事,成为美国大学申请条件下,含金量相当高的官方竞赛。USACO的比赛成绩可以助力计算机专业留学,越来越多的学生进入了康奈尔,麻省理工,普林斯顿,哈佛和耶鲁等大学,这些同学的共同点是他们都参加了美国计算机科学竞赛(USACO),并且取得过非常好的成绩。适合参赛人群USACO适合国内在读学生有意向申请美国大学的或者想锻炼自己编程能力的同学,高三学生也可以参加12月的第_usaco可以多次提交吗
文章浏览阅读394次。1.1 存储程序1.2 创建存储过程1.3 创建自定义函数1.3.1 示例1.4 自定义函数和存储过程的区别1.5 变量的使用1.6 定义条件和处理程序1.6.1 定义条件1.6.1.1 示例1.6.2 定义处理程序1.6.2.1 示例1.7 光标的使用1.7.1 声明光标1.7.2 打开光标1.7.3 使用光标1.7.4 关闭光标1.8 流程控制的使用1.8.1 IF语句1.8.2 CASE语句1.8.3 LOOP语句1.8.4 LEAVE语句1.8.5 ITERATE语句1.8.6 REPEAT语句。_mysql自定义函数和存储过程
文章浏览阅读188次。半导体二极管——集成电路最小组成单元。_本征半导体电流为0
文章浏览阅读2.8k次,点赞3次,收藏18次。游戏水面特效实现方式太多。咱们这边介绍的是一最简单的UV动画(无顶点位移),整个mesh由4个顶点构成。实现了水面效果(左图),不动代码稍微修改下参数和贴图可以实现岩浆效果(右图)。有要思路是1,uv按时间去做正弦波移动2,在1的基础上加个凹凸图混合uv3,在1、2的基础上加个水流方向4,加上对雾效的支持,如没必要请自行删除雾效代码(把包含fog的几行代码删除)S..._unity 岩浆shader
文章浏览阅读5k次。广义线性模型是线性模型的扩展,它通过连接函数建立响应变量的数学期望值与线性组合的预测变量之间的关系。广义线性模型拟合的形式为:其中g(μY)是条件均值的函数(称为连接函数)。另外,你可放松Y为正态分布的假设,改为Y 服从指数分布族中的一种分布即可。设定好连接函数和概率分布后,便可以通过最大似然估计的多次迭代推导出各参数值。在大部分情况下,线性模型就可以通过一系列连续型或类别型预测变量来预测正态分布的响应变量的工作。但是,有时候我们要进行非正态因变量的分析,例如:(1)类别型.._广义线性回归模型
文章浏览阅读69次。环境保护、 保护地球、 校园环保、垃圾分类、绿色家园、等网站的设计与制作。 总结了一些学生网页制作的经验:一般的网页需要融入以下知识点:div+css布局、浮动、定位、高级css、表格、表单及验证、js轮播图、音频 视频 Flash的应用、ul li、下拉导航栏、鼠标划过效果等知识点,网页的风格主题也很全面:如爱好、风景、校园、美食、动漫、游戏、咖啡、音乐、家乡、电影、名人、商城以及个人主页等主题,学生、新手可参考下方页面的布局和设计和HTML源码(有用点赞△) 一套A+的网_垃圾分类网页设计目标怎么写
文章浏览阅读614次,点赞7次,收藏11次。之前找到一个修改 exe 中 DLL地址 的方法, 不太好使,虽然能正确启动, 但无法改变 exe 的工作目录,这就影响了.Net 中很多获取 exe 执行目录来拼接的地址 ( 相对路径 ),比如 wwwroot 和 代码中相对目录还有一些复制到目录的普通文件 等等,它们的地址都会指向原来 exe 的目录, 而不是自定义的 “lib” 目录,根本原因就是没有修改 exe 的工作目录这次来搞一个启动程序,把 .net 的所有东西都放在一个文件夹,在文件夹同级的目录制作一个 exe._.net dll 全局目录
文章浏览阅读1.5k次。本文为转载,原博客地址:http://blog.csdn.net/hujingshuang/article/details/46910259简介 BRIEF是2010年的一篇名为《BRIEF:Binary Robust Independent Elementary Features》的文章中提出,BRIEF是对已检测到的特征点进行描述,它是一种二进制编码的描述子,摈弃了利用区域灰度..._breif description calculation 特征点
文章浏览阅读4.1k次,点赞21次,收藏79次。本文是《基于SpringBoot的房屋租赁管理系统》的配套原创说明文档,可以给应届毕业生提供格式撰写参考,也可以给开发类似系统的朋友们提供功能业务设计思路。_基于spring boot的房屋租赁系统论文