getpatch
import time
import os
import math
import sys
import os,os.path,shutil
import numpy as np
import cv2 img_in_path = 'F:/project/Breast/InBreast/INBreast/removeother/'
# img_in_path = 'F:/project/Breast/InBreast/INBreast/imgtest/'
outpathimgyes = 'F:/project/Breast/InBreast/INBreast/outimgpatch/calcification/'
outpathimgno = 'F:/project/Breast/InBreast/INBreast/outimgpatch/no/'
txtPath = 'F:/project/Breast/InBreast/INBreast/AllTXTall/'
# txtPath = 'F:/project/Breast/InBreast/INBreast/ALLtest/'
outpathtxt = 'F:/project/Breast/InBreast/INBreast/outtxtpatch/'
txtType = 'txt'
txtLists = os.listdir(txtPath) #列出文件夹下所有的目录与文件 # Read all points into a list(before 21:30,715) for filename in txtLists:
pp = 0
ppp = 0
with open(txtPath + filename, 'r') as file:
print(filename)
lines = file.readlines()
dataset = [[] for i in range(len(lines))]
for i in range(len(dataset)):
dataset[i][:] = (item for item in lines[i].strip().split(',')) # 逐行读取数据
#print("dateset:", dataset)
# print(type(dataset[0][0]))
# print(dataset.__sizeof__())
u = np.array(dataset)
list = np.zeros((u.shape[0], 2))
for i in range(u.shape[0]):
# print(u[i,0][0])
findNumber = u[i,0].find(" ")
# print(findNumber)
list[i,0] = round(float(u[i, 0][0:findNumber]))
list[i,1] = round(float(u[i, 0][findNumber+1:])) # Read the same name image(before 21:30,715)
img = img_in_path + filename[:-4] + '.png' # out_path
img = cv2.imread(img)
black = np.zeros((img.shape[0], img.shape[1]))
for i in range(img.shape[0]):
for j in range(img.shape[1]):
if img[i,j][0] == 0 and img[i,j][1] == 0 and img[i,j][2] == 0:
black[i][j] = 1
print(img.shape[0],img.shape[1])
# Cutting distance has been determined, Handling cutting details(before 22.30, 715)
row = img.shape[0]/56;
col = img.shape[1]/56;
print(row,col) # 待考虑问题:黑色像素
for i in range(int(row-1)):
# print(i)
for j in range(int(col-1)):
tt = 0
pb = np.sum(black[i * 56:i * 56 + 112,j * 56:j * 56 + 112]) / 12544
# print("top_left:", i*56, j*56)
# print("bottom_right:",i*56+112,j*56+112)
for k in range(u.shape[0]):
if j*56<list[k,0]<j*56+112 and i*56<list[k,1]<i*56+112:
tt = tt + 1
patch_label = np.zeros((tt, 4))
tt1 = 0
for k in range(u.shape[0]):
if j*56<list[k,0]<j*56+112 and i*56<list[k,1]<i*56+112:
patch_label[tt1,0] = list[k,0] - j*56
patch_label[tt1,1] = list[k,1] - i*56
patch_label[tt1,2] = list[k,0]
patch_label[tt1,3] = list[k,1]
tt1 = tt1 + 1
if tt>0 and pb == 0:
patch = img[i * 56:i * 56 + 112,j * 56:j * 56 + 112]
pp = pp + 1
# print(patch.shape)
outputimg = outpathimgyes + '/' + filename[:-4] + '/' + filename[:-4] + '_'+ str(pp) + '.png'
cv2.imwrite(outputimg, patch)
outputtxt1 = outpathtxt + '/' + filename[:-4] + '/' + filename[:-4] + '_' + str(pp) + '.txt'
np.savetxt(outputtxt1, patch_label, fmt='%f')
if tt==0 and pb == 0:
patch = img[i * 56:i * 56 + 112,j * 56:j * 56 + 112]
ppp = ppp + 1
# print(patch.shape)
outputimg1 = outpathimgno + '/' + filename[:-4] + '/' + filename[:-4] + '_' + str(ppp) + '.png'
cv2.imwrite(outputimg1, patch) print(img.shape[0], img.shape[1])
print(pp,ppp) for i in range(int(row-1)):
tt = 0
pb = np.sum(black[i * 56:i * 56 + 112, img.shape[1] - 113:img.shape[1]-1]) / 12544
# print("top_right:",i*56,img.shape[1]-1)
# print("bottom_left:", i * 56+112, img.shape[1] - 113)
for k in range(u.shape[0]):
if img.shape[1] - 113 < list[k, 0] < img.shape[1]-1 and i * 56 < list[k, 1] < i * 56 + 112:
tt = tt + 1
patch_label = np.zeros((tt, 4))
tt1 = 0
# print(tt)
for k in range(u.shape[0]):
if img.shape[1] - 113 < list[k, 0] < img.shape[1]-1 and i * 56 < list[k, 1] < i * 56 + 112:
patch_label[tt1, 0] = list[k, 0] - j * 56
patch_label[tt1, 1] = list[k, 1] - i * 56
patch_label[tt1, 2] = list[k, 0]
patch_label[tt1, 3] = list[k, 1]
tt1 = tt1 + 1
if tt > 0 and pb == 0.1:
patch = img[i * 56:i * 56 + 112, img.shape[1] - 113:img.shape[1]-1]
pp = pp + 1
# print(patch.shape)
outputimg = outpathimgyes + '/' + filename[:-4] + '/' + filename[:-4] + '_' + str(pp) + '.png'
cv2.imwrite(outputimg, patch)
outputtxt1 = outpathtxt + '/' + filename[:-4] + '/' + filename[:-4] + '_' + str(pp) + '.txt'
np.savetxt(outputtxt1, patch_label, fmt='%f')
if tt == 0 and pb == 0:
patch = img[i * 56:i * 56 + 112, img.shape[1] - 113:img.shape[1]-1]
ppp = ppp + 1
# print(patch.shape)
outputimg1 = outpathimgno + '/' + filename[:-4] + '/' + filename[:-4] + '_' + str(ppp) + '.png'
cv2.imwrite(outputimg1, patch)
print(img.shape[0], img.shape[1])
#
#
#
for j in range(int(col-1)):
# print("bottom_left:", img.shape[0]-1, j*56)
# print("top_right:", img.shape[0] - 113, j * 56+112, )
tt = 0
pb = np.sum(black[img.shape[0]-113:img.shape[0] - 1, j * 56:j * 56+112]) / 12544
# print("top_right:",i*56,img.shape[1]-1)
# print("bottom_left:", i * 56+112, img.shape[1] - 113)
for k in range(u.shape[0]):
if j * 56 < list[k, 0] < j * 56+112 and img.shape[0]-113 < list[k, 1] < img.shape[0]-1:
tt = tt + 1
patch_label = np.zeros((tt, 4))
tt1 = 0
# print(tt)
for k in range(u.shape[0]):
if j * 56 < list[k, 0] < j * 56+112 and img.shape[0]-113 < list[k, 1] < img.shape[0]-1:
patch_label[tt1, 0] = list[k, 0] - j * 56
patch_label[tt1, 1] = list[k, 1] - i * 56
patch_label[tt1, 2] = list[k, 0]
patch_label[tt1, 3] = list[k, 1]
tt1 = tt1 + 1
if tt > 0 and pb == 0:
patch = img[img.shape[0]-113:img.shape[0] - 1, j * 56:j * 56+112]
pp = pp + 1
# print(patch.shape)
outputimg = outpathimgyes + '/' + filename[:-4] + '/' + filename[:-4] + '_' + str(pp) + '.png'
cv2.imwrite(outputimg, patch)
outputtxt1 = outpathtxt + '/' + filename[:-4] + '/' + filename[:-4] + '_' + str(pp) + '.txt'
np.savetxt(outputtxt1, patch_label, fmt='%f')
if tt == 0 and pb == 0:
patch = img[img.shape[0]-113:img.shape[0] - 1, j * 56:j * 56+112]
ppp = ppp + 1
# print(patch.shape)
outputimg1 = outpathimgno + '/' + filename[:-4] + '/' + filename[:-4] + '_' + str(ppp) + '.png'
cv2.imwrite(outputimg1, patch)
# print(img.shape[0], img.shape[1]) # A = np.array([[1, 2], [3, 4], [5, 6]])
# print(np.sum(A))
print(pp)
getpatch的更多相关文章
- Android系统Recovery工作原理之使用update.zip升级过程分析(一)
通过分析update.zip包在具体Android系统升级的过程,来理解Android系统中Recovery模式服务的工作原理.我们先从update.zip包的制作开始,然后是Android系统的启动 ...
- Amazon API Gateway Importer整合过程小结
(1)需要将swagger json转换成amazon api gateway 所需要的格式(根据Method Request中 Request PathsURL Query String Param ...
- 深入了解Libgdx中间Skin分类
文不是直接翻译.. . 本文在Libgdx的官方wiki的基础上,加上一些自己的理解.所以,难免会一些甚至是非常多的理解非常片面的东西.写的不好,还请见谅.... 事实上 事实上.在LibGDX的官方 ...
- recovery 差分升级包制作超时
我们在对android系统升级的时候,可以减少升级包的大小,只升级差异部分,也就是差分包升级,相关的描述可以参考:http://blog.csdn.net/csdn66_2016/article/de ...
- recovery 差分升级包制作超时【转】
本文转载自:https://blog.csdn.net/csdn66_2016/article/details/73800349 我们在对android系统升级的时候,可以减少升级包的大小,只升级差异 ...
- centos 6.4安装杀毒软件clamAV 0.98[转]
原文出处: http://dnuser.blog.51cto.com/4863891/1303829 1.查看系统版本 [root@local]# lsb_release -a LSB Versi ...
- Android系统Recovery工作原理之使用update.zip升级过程分析(一)---update.zip包的制作【转】
本文转载自:http://blog.csdn.net/mu0206mu/article/details/7399822 这篇及以后的篇幅将通过分析update.zip包在具体Android系统升级的过 ...
随机推荐
- 201871010111-刘佳华《面向对象程序设计(java)》第十四周学习总结
201871010111-刘佳华<面向对象程序设计(java)>第十四周学习总结 实验十二 Swing图形界面组件(一) 实验时间 2019-11-29 第一部分:基础知识总结 1.设计 ...
- 面向对象程序设计(JAVA) 第13周学习指导及要求
2019面向对象程序设计(Java)第13周学习指导及要求 (2019.11.19-2019.11.25) 学习目标 (1) 掌握事件处理的基本原理,理解其用途: (2) 掌握AWT事件模型的工作 ...
- HDL的三种描述方式
结构化描述 结构化描述方式是最原始的描述方式,是抽象级别最低的描述方式,但同时也是最接近于实际的硬件结构的描述方式.结构化的描述方式,思路就像在面包板上搭建数字电路一样,唯一的不同点就是我们通过HDL ...
- Java程序猿怎么才能月薪过万?
每一个略微有点长进的人,都应该把作业里的前三名作为自己斗争的政策和对手.你离成为冠军Java程序员还有多远,看完这篇你就知道了. 软件工程师的作业生涯里,知识有一个三年的半衰期.这意味着三年后,你所具 ...
- 最近公共祖先(LCA)基础模板(倍增法)
之前在澡堂学过这么个东西,听课时理解非常透彻,然后做题时是这种状态: 因为并没有切板子题,最近切掉以后看同桌,他默默地说了一句话: 我是什么时候A的来着... 我当时就心态爆炸... 现在来进行简单整 ...
- 如何将Excel表批量赋值到ArcGIS属性表
情景再现 现需要将Excel表信息批量赋值(不是挂接)到Shp文件的属性表,两张表的字段.记录数一模一样,至于为什么会出现这样的问题,咱也不敢问,只有想个法子把它搞定! 原始的Excel信息表共57列 ...
- [学习笔记] Manacher与PAM
\(1\) Manacher 挺短,背是挺好背的 Manacher用于求回文串长度.思想大概就是: 1.加入字符集之外的识别字符(比如#)分隔开原来相邻的字母,这样所有的回文串都变成了以某个字符为中心 ...
- vb.net 带有一个参数的线程
For Each _row As DataGridViewRow In datagridview.Rows 'searchRecords_refreshRow(_row) ' 上面以前,直接运行函数, ...
- css的伪元素 ::after ::before 和 图标字体的使用
浅谈css的伪元素::after和::before css中的::after和::before已经被大量地使用在我们日常开发中了,使用他们可以使我们的文档结构更加简洁.但是很多人对::after和 ...
- js常用但是不容易记住的代码
<!-- iframe 自适应高度度 --><iframe src="__CONTROLLER__/showlist" frameborder="0&q ...