@Cy86
@jdhao
@chashao
@paopjian
老哥们, 辛苦啊...
我现在是直接用真实的患者 CT 片子写..
不能外传....

另外, 我看有人已经上 pydicom 了...
另外我请教一下, 我现在处理的 dicom 格式不是通用的 dicom 格式...
是排版好之后的 dcm 文件.
很多信息不能直接获取..最重要的 pydicom.pixel_array 和 pixelData 是没有的..
我是通过 pydicom
dcm = pydicom.read_file(filePath)
imageInfo = dcm.ReferencedImageBoxSequence[0].BasicGrayscaleImageSequence[0]
这种方式获取的 pixelData 的 bytes
有人能科普一下这个问题吗?

@wcsjtu 我先尝试 numpy 的方案....再看预编译或者 jit 类的...
谢谢~

各位的代码, 够我消化一夜了....
太感谢了...太感谢了!
@hsfzxjy
@yoohwzy
@stein42

@dlsflh
我是直接写个 tensor 扔给 torch 跑，我知道它会自动用 SIMD 或者 CUDA (如果启动 gpu 的话)

@lmshl 这大哥的做法解释了你的疑问..

@stein42 卧槽卧槽, 果然大神多....
我需要消化一下代码...
这之前都没用过 numpy...感谢感谢..

另外请教一下,
buffer = bytes([208, 4]) * (width * height)
文件内容不只是 [208, 4] 还有其他的灰度值和通道值, 这一步不理解..

@threebr 我的理解是, python 要通过对应的接口和 numpy 通信, 而 numpy 本身对于下标赋值这种操作的过程是很复杂的.
这基本就是昨天看到的一篇文章的原话. 作者还列出了对应的处理过程, 因为是 C, 直接放弃阅读了..

@hsfzxjy 是吗? 一直很怕上 C 相关的东西...哈哈
nodejs 可以直接读数据库获取路径, 去转换完了写回对应字段...
解决问题来说, 应该这个方法是最快的..
但是不服气啊...60 秒和 2.3 秒的差距...

@yoohwzy
@lookStupiToForce 感谢回复..先尝试 np.stack 的用法...
再看 google 结果...
对..如果不用转换, 直接读取生成, 是正确的思路..只是当时没找到方法, 就先硬上了..
先解决,再优化~

@lizytalk 试过.. 仅限于 @jit 装饰. 一堆警告和报错.
主要是不支持一些方法. 包括不限于 numpy 的...
网上翻的教程, 所有的数组转 uint8 类型, 都是 numpy 直接干的...
其他我就不会了...有其他方法我也想试试.... 不用 numpy 转换, 或许可以直接上 pypy

@MoYi123
我读出来的 bytes 数据处理成 uint8 后. 是这样的形式:
[208, 4, 208, 4, 208, 4...196, 8]
不懂图像处理的知识, 我的理解就是, 一个灰度, 一个 Alpha 通道值(透明值?).
第一次处理数组是要把上面的数组, 改成如下形式: <- 暂且叫 生成数组阶段.
就是再赋值两次灰度..写入数组.构建成如下形式:

[208, 208, 208, 4, 208, 208, 208, 4.......]

再下一步就是把 208 这个灰度值变成 rgb 的形式.
然后再通过一个循环变成 pillow 支持的格式, 如下: <- 暂且叫 数组转换阶段吧.
[
[ [r, g, b, a], [r, g, b, a], [r, g, b, a], ],
[ [r, g, b, a], [r, g, b, a], [r, g, b, a], ],
[ [r, g, b, a], [r, g, b, a], [r, g, b, a], ],
]

我补充了点文字..你看一下能理解吗?

@dlsflh 感谢回复..
我猜..他们只是读...并不涉及到数组大范围的修改, 或者大量的下标赋值操作吧?所以并没有那么慢..
哦对, 另外可能还直接调用了 GPU 做计算...而我目前不能这么做..服务器的显卡很垃圾.


@hsfzxjy cython....是不是不如直接 nodejs 搞一个了...研究以下打包发布, python 直接拉起一个 node 的进程后台转换..这个东西有开发周期的约束, 快到了..

所有的电子设备, 我都是分期上的...
尤其是 apple 产品...24 期挺好的..

@lostberryzz 干眼症可能看绿色比较舒服吧? 或者散光? 或者近视? 反正我都有
.羡慕吗.
我干眼症严重的时候, 就是墨绿色的背景色... 亮色都太刺眼, 感觉看不清楚....

@nielinjie 受教了...我又学会了一个专业名词..
可以去忽悠客户了

你昨晚买域名的时候是不是上市计划都做好了?哈哈哈

@ysc3839 噢...受教了

uuid 是不是有个限定命名空间的? python 是 uuid5
能够确保在命名空间里是唯一的.

@fzcf 嗯...挺厉害的.
国内东莞长安那一带的指纹考勤厂, 也很多朝鲜的工程师..

分享一些经历, 希望对你有帮助...
朋友高中时候做了 "磨平角膜" 的那个手术, 治疗近视...专业名称是叫 飞秒 吗? 还是飞秒是更先进的了?
大夫嘱咐要卧床休息 30 天, 不能用眼, 比如看电视, 手机, 玩电脑, 看书 都不行..
他忍了 2 周, 然后觉得没事了...开始看电视, 跑网吧和我们一起玩.. 那个时候 iphone3 都没发布.. 所以手机没啥玩的..
最强智能机可能是索爱 K750c 或者 moto E680i ??? 有手机估计更完蛋了.
然后玩了不到一周就说眼睛不舒服, 看不清东西..
现在戴的眼镜和瓶酒瓶底差不多..

后来看到过一些相关资料, 说现在有些手术不需要长时间休息?? 不知道和你了解的手术是否是一个类型的..
那个年代我记得这种手术刚出来没多久... 大约是 2004 年的样子吧?

然后我老婆当年也做了类似的手术, 但是她就很老实的休息了 20 多天, 恢复的比较好...
他俩做的时间相差不了 2 年... 老婆现在视力没问题, 没有什么不良反应.

以上, 两条例子供你参考...希望有帮助

另外, 平时可以吃一些维生素 A..
我之前有段时间眼睛突然不舒服, 长时间看代码就感觉很空洞, 看不清楚, 眼睛干涩. 查了一下大概是干眼症吧? 我还有光诱型偏头痛, 那段时间动不动就犯...
后来陆续吃了 30 斤胡萝卜, 没事就坐着啃, 补充维生素 A, 叶黄素, 注意休息...
现在应该是好过来了, 大概上海疫情封城期间的事...5 月吧? 8 月就基本好了..
现在没事就出去转转, 抬头看天空和远处...