数字旗手

电气化、自动化、数字化、智能化、智慧化

0%

ImagePy解析：8 -- 由新建图像谈起（引出IPy、ImagePlus、Canvas、ImageManager）

Posted on Edited on In Disqus:

%%%%% 更新日志 %%%%%
2019-10-11更新：增加ImageManager管理器解析
%%%%%%%%%%%%%%%%%%%%

参考文献：
创建图像
ImagePy开发文档 —— 常用汇总
ImagePy开发文档 —— 图像封装类

新建图像的插件

最简单的新建图像的插件代码如下：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
from imagepy.core.engine import Free
from imagepy import IPy
import numpy as np

class Plugin(Free):
    title = 'New Image Demo'
    para = {'name':'new image','w':300, 'h':300}
    view = [(str, 'name', 'name',''),
            (int, 'w', (1,2048), 0,  'width', 'pix'),
            (int, 'h', (1,2048), 0,  'height', 'pix')]

    def run(self, para = None):
        imgs = [np.zeros((para['h'], para['w']), dtype=np.uint8)]
        IPy.show_img(imgs, para['name'])

该插件的写法可以参考之前的解析文章，在run()方法中接收para的宽度w和高度h的数值，然后通过numpy的zeros函数生成了一个全是0的array数组，注意imgs变量又对该array数组加入了一个方括号，即将其转化为python的列表（其实从后面可以看出，该列表的长度就是整个图像栈有多少个slices），所以imgs的具体数值就是：
1
2
3
4
5
6
7
imgs =  [array([[0, 0, 0, ..., 0, 0, 0],
       [0, 0, 0, ..., 0, 0, 0],
       [0, 0, 0, ..., 0, 0, 0],
       ...,
       [0, 0, 0, ..., 0, 0, 0],
       [0, 0, 0, ..., 0, 0, 0],
       [0, 0, 0, ..., 0, 0, 0]], dtype=uint8)]

然后将imgs和用户输入的图像名称传入IPy的show_img()方法。到这里就可以看出，ImagePy是将传统的numpy数组（一般的图像库都是将图像读作numpy数组）进行了进一步的封装，这个封装类就是ImagePlus，后面再具体分析ImagePlus。

IPy常用工具集

上面将imgs传入了IPy的show_img()方法，实际IPy是一些常用功能的汇总，比如展示一幅图像，获取当前图像窗口等功能，其实这些功能是分布在各个功能组建中的，IPy只是重新引用和整理，并且IPy处于最顶层目录，随处都可以方便的引用这样以来，提高了开发效率。（语出上面的参看文献）
具体看show_img()方法：

1
2
3
4
5
6
7
8
9
10
11
12
13
def showimg(imgs, title):
    print('show img')
    from .core import ImagePlus
    ips = ImagePlus(imgs, title)
    showips(ips)

pub.subscribe(showimg, 'showimg')
def show_img(imgs, title):
    if uimode()=='no':
        from .core import manager, ImagePlus
        from .ui.canvasframe import VirturlCanvas
        frame = VirturlCanvas(ImagePlus(imgs, title))
    else:wx.CallAfter(pub.sendMessage, 'showimg', imgs=imgs, title=title)

可以看出，该方法又通过Publisher/Subscriber机制调用了showimg()方法，这里就引出了ImagePlus封装类：
1
ips = ImagePlus(imgs, title)

ImagePlus接收imgs和title，组装成了ImagePlus类型的ips。

ImagePlus图像封装类

ImagePy 采用 Numpy 数组作为图像数据结构，Numpy 非常强大，但是要让它能够支持我们的全部特性（索引色支持，撤销，选区支持，多通道支持，图像栈支持），还需要一些其他的数据给予必要的辅助，因而我们有了ImagePlus类（语出上面的参考文献）。
下面看看ImagePlus收到imgs后干了什么：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
def set_imgs(self, imgs):
    self.is3d = not isinstance(imgs, list)
    self.scrchanged = True
    self.snap = None
    self.imgs = imgs
    self.size = self.imgs[0].shape[:2]
    self.height, self.width = self.size
    self.imgtype = get_img_type(self.imgs)
    if self.imgs[0].ndim==2: self.channels = 1
    else: self.channels = self.imgs[0].shape[2]
    self.dtype = self.imgs[0].dtype
    if self.dtype == np.uint8: self.range = (0, 255)
    else: self.range = self.get_updown('all', 'one')
    if self.dtype == np.uint8:
        self.chan_range = [(0, 255)] * self.channels
    else: self.chan_range = self.get_updown('all', 'all', step=512)
    self.chan = (0, [0,1,2])[self.channels==3]

针对于上述图像，经过初始化后，其对应的ImagePlus的各个属性值为：
1
2
3
4
5
6
7
8
9
self.is3d =  False
self.size =  (300, 300)
self.type =  8-bit
self.imgs[0].shape =  (300, 300)
self.channels =  1
self.dtype =  uint8
self.chan =  0
self.chan_range =  [(0, 255)]
self.range =  (0, 255)

如果是使用ImagePy源码里的new_plg，即：
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
import wx,os
from imagepy.ui.canvasframe import CanvasFrame
import numpy as np
from imagepy import IPy

from imagepy.core.engine import Free

class Plugin(Free):
    title = 'New'
    para = {'name':'Undefined','width':300, 'height':300, 'type':'8-bit','slice':1}
    view = [(str, 'name', 'name', ''),
            (int, 'width',  (1,10240), 0,  'width', 'pix'),
            (int, 'height', (1,10240), 0,  'height', 'pix'),
            (list, 'type', ['8-bit','RGB'], str, 'Type', ''),
            (int, 'slice',  (1,2048), 0,  'slice', '')]

    #process
    def run(self, para = None):
        w, h = para['width'], para['height']
        channels = (1,3)[para['type']=='RGB']
        slices = para['slice']
        shape = (h,w,channels) if channels!=1 else (h,w)
        imgs = [np.zeros(shape, dtype=np.uint8) for i in range(slices)]
        IPy.show_img(imgs, para['name'])

且此时设定Type为RGB、slice为2，那么此时图像的各个ImagePlus属性值为：
1
2
3
4
5
6
7
8
9
self.is3d =  False
self.size =  (300, 300)
self.type =  rgb
self.imgs[0].shape =  (300, 300, 3)
self.channels =  3
self.dtype =  uint8
self.chan =  [0, 1, 2]
self.chan_range =  [(0, 255), (0, 255), (0, 255)]
self.range =  (0, 255)

这里插一句，上面的插件的run()函数中判断channels的值的语句为：
1
channels = (1,3)[para['type']=='RGB']

这是Python的一种特殊的条件表达式，称为元组条件表达式，即(X, Y)[C]，如果C为False，则返回X，如果C为True，则返回Y，这是因为在Python中，True等于1，而False等于0，这就相当于在元组中使用0和1来选取数据。
元组条件表达式在这里是没有问题的，但下面的两篇文章有个提醒，在使用这样的条件表达式时，需要注意如下问题，即元组是提前构建数据，则X和Y会提前运行，然后再用True(1)/False(0)来索引到数据。如果X和Y是比较耗资源的运算，则会开销比较大，另外如果有异常运算，也会直接报错。此时建议使用X if C else Y的条件表达式形式。
Python中的三元操作
三元运算符

回到两张图像的对比，可以看出对于一张8-bit的只有一个slice的图像与一张RGB的有两个slices的图像来说，两者的图像类型type、图像矩阵的shape、通道标识chan、通道范围chan_range都是不同的。
这里再插一句，对于ImagePlus类的range和img属性，ImagePy采取了既非普通属性、又非方法调用的方式来设定和读取，即property装饰器的方法，这样的好处有两点：（1）相比于普通的属性，这样可以进行参数检查，保证安全性和合理性；（2）相比于set()和get()这样的配对函数，这样的方式又更加简洁。其实原理就是property装饰器将一个getter方法变成了属性，同时生成了一个setter装饰器，负责将setter方法变成属性赋值。对比img和range这两个属性，可以发现，range有property和相应的setter装饰器，而img只有property装饰器，这表明range是可读可写，而img是只读方式。这里的知识点见廖大爷的如下教程：
使用@property

另外，这个img属性的具体定义：

1
2
@property
def img(self):return self.imgs[self.cur]

可以看出，img具体调用的是哪个图像还与当前的游标有关。

Canvas画布

继续回到IPy，可以看出构造了ImagePlus类型的ips后，就继续调用IPy的showips()函数：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
showips(ips)

def showips(ips):
    if uimode()=='ipy':
        from .ui.canvasframe import CanvasPanel
        canvasp = CanvasPanel(curapp.canvasnb)
        canvasp.set_ips(ips)
        curapp.canvasnb.add_page( canvasp, ips)
        #canvasp.canvas.initBuffer()
        curapp.auimgr.Update()

    elif uimode()=='ij':
        from .ui.canvasframe import CanvasFrame
        frame = CanvasFrame(curapp)
        frame.set_ips(ips)
        frame.Show()

根据当前的UI界面是ImagePy风格还是ImageJ风格，来以不同的方式显示图像。其实ImageJ风格的CanvasFrame也是调用了ImagePy自定义的CanvasPanel，它是一张张图像分成不同的窗口显示，而ImagePy风格的显示方式是以标签页的形式聚合显示不同图像，即传入CanvasPanel的是curapp.canvasnb，查看这个canvasnb即可知道，它是CanvasNoteBook类型的数据，其派生自wx.lib.agw.aui.AuiNotebook。
CanvasPanel又调用了同级目录canvas下的Canvas类，前面已经提到，该类已经从ImagePy中解耦，可以作为独立组件运行。
Canvas类的运行过程乍一看不好理解，其实它是利用了wxPython的SizeEvent：
1
2
3
4
5
    self.bindEvents()
def bindEvents(self):
    for event, handler in [ \
            (wx.EVT_SIZE, self.on_size)]:
        self.Bind(event, handler)

即当窗口绘制时，就会调用on_size()函数，然后就会执行一系列的绘图等操作。
ImagePy的Canvas是个定制的wxPython Panel，它显示图像的基本原理是利用wxPython的GDI接口，具体是wxPython的BufferedDC，以及其DrawBitmap()函数。

当创建Canvas对象后，接着就是将图像传入其中，即：

1
2
3
4
5
# CanvasPanel设置ips
canvasp.set_ips(ips)

# CanvasFrame设置ips，实际也是调用的CanvasPanel
frame.set_ips(ips)

然后就是将这些Panel显示出来：
1
2
3
4
5
6
# imagepy风格的Panel显示
curapp.canvasnb.add_page( canvasp, ips)
curapp.auimgr.Update()

# imagej风格的Panel显示
frame.Show()

以imagej风格为例（因为它相对简单，直接Show()出来，没有标签页那些东西需要分析），当Show()方法（这是wxPython的每个frame显示时都要调用的方法）被调用时，就会产生“激活”事件（wx.ActivateEvent），即：
1
2
3
4
5
self.Bind(wx.EVT_ACTIVATE, self.on_valid)
def on_valid(self, event):
    if event.GetActive():
        ImageManager.add(self.canvaspanel.ips)
        WindowsManager.add(self.canvaspanel)

一旦被激活，可以看出，此时接着就会将当前的ips和当前的panel加入到ImageManager管理器和WindowsManager管理器中。
下面就以ImageManager为例详细分析一下管理器的运行路径。

ImageManager管理器

管理器是 ImagePy 里很重要的一个概念，作为插件系统，各个部件都是松散耦合。而管理器的作用就是全局协调。（语出上面的参考文献）
先来看一下ImageManager的全貌：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
class ImageManager:
    imgs = []
    @classmethod
    def add(cls, ips):
        …       
    @classmethod
    def remove(cls, ips):
        …           
    @classmethod
    def get(cls, title=None):
        …         
    @classmethod
    def get_titles(cls):
        …       
    @classmethod
    def name(cls, name):
        …

首先注意到这个类与之前类的定义的一个很大不同就是每个方法前面都加了classmethod修饰符，这样的好处是调用这些方法时不需要实例化，不需要 self 参数，但第一个参数需要是表示自身类的 cls 参数，可以来调用类的属性，类的方法，实例化对象等。参考资料有：
Python classmethod 修饰符
飘逸的python - @staticmethod和@classmethod的作用与区别
正确理解Python中的 @staticmethod@classmethod方法

然后具体看一下add()方法：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
@classmethod
def add(cls, ips):
    print("ips in add() = ", ips)
    cls.remove(ips)
    callback = lambda a: cls.remove(a())
    def callback(a):
        cls.remove(a())
    print('~~~ ~~~ image add!')
    cls.imgs.insert(0, weakref.ref(ips, callback))
    print("~~~ ~~~ ImageManager = ", cls.imgs)

@classmethod
def remove(cls, ips):
    for i in cls.imgs:
        if i() == ips:
            print("ips in remove() = ", ips)
            print("~~~ ~~~ !!! remove imgs !!!")
            cls.imgs.remove(i)

add()首先就调用了该类的remove()方法，在这个方法里会判断该图像ips是否在ImageManager的imgs列表中，如果有则删除。
插一句，这个remove方法不只是在add()中调用，在关闭某张图像时，也会调用，因为：
1
2
3
4
5
6
7
self.Bind(wx.EVT_CLOSE, self.on_close)
def on_close(self, event):
    WindowsManager.remove(self.canvaspanel)
    ImageManager.remove(self.canvaspanel.ips)
    self.canvaspanel.ips = None
    self.canvaspanel.back = None
    event.Skip()

然后定义了一个lambda表达式，虽然它是匿名函数，但是这里还是给它取了个名字callback，那其实这里有个问题需要后续研究一下，即这个lambda表达式与下面接着的显式定义的callback()函数有何不同。。。这个地方没明白，留作后续详细解读。。
插一句lambda表达式的知识点，见：
什么时候使用Lambda函数？

然后就是将当前图像ips插入到imgs列表中，这里用到了weakref弱引用，参考资料见：
weakref —- 弱引用
python中的weakref(弱引用)/)
3.9. weakref — 实现对象的弱引用
引用上面资料中对weakref的用途：

python使用weakref模块可以实现对象的弱引用，python中的垃圾回收机制是通过对象的引用计数来触发的，当对象间有相互引用时，垃圾回收会失败，这时可以使用weakref模块实现对象的弱引用，给对象设置弱引用并不会增加对象的引用计数，所以不会影响内存的回收，同时也正是因为弱引用不会增加对象的 引用计数，弱引用并不能保证原对象存在，弱引用的初衷是为了在进行大对象处理时实现缓存机制，比如从而节省资源占用，比如当你需要处理大量的二进制图片文件对象时，为了方便读写，我们需要进行一个图片名称到图片对象的映射，如果使用字典进行存储对象和名称的映射的存储，那么在程序运行期间，所有的图片对象都不会被销毁，因为此时对象已经作为dict对象中的key或者value，这就造成了过多的资源消耗，若此时使用weakref模块的WeakKeyDictionary和WeakValueDictionary类对对象进行弱引用，则可以减少资源占用，方便进行内存回收。

回调函数callback的用途：在引用 「死亡」 并且没有其他引用指向原始对象的时候将引用对象作为参数传递给回调函数。这个功能的一种应用方法是从缓存中删除弱引用对象。

有一点需要额外注意的是，使用weakref.ref()获得对象引用后，下面调用时需要加上括号才能获得原对象，如果是使用weakref.proxy()获得引用，则调用时不需要加上括号。

下面做一个测试，看看上面的代码中的print信息都有啥。
具体操作是这样的：先新建一个图像，名为Undefined，然后再新建一个图像，名为Undefined-1，然后关闭Undefined-1。
操作下来，print的信息如下：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
ips in add() =  <imagepy.core.wraper.imageplus.ImagePlus object at 0x000002914889F5C0>
~~~ ~~~ image add!
~~~ ~~~ ImageManager =  [<weakref at 0x0000029148895548; to 'ImagePlus' at 0x000002914889F5C0>]

ips in add() =  <imagepy.core.wraper.imageplus.ImagePlus object at 0x0000029148883198>
~~~ ~~~ image add!
~~~ ~~~ ImageManager =  [<weakref at 0x00000291447469A8; to 'ImagePlus' at 0x0000029148883198>, <weakref at 0x0000029148895548; to 'ImagePlus' at 0x000002914889F5C0>]
 
ips in remove() =  <imagepy.core.wraper.imageplus.ImagePlus object at 0x0000029148883198>
~~~ ~~~ !!! remove imgs !!!
ips in add() =  <imagepy.core.wraper.imageplus.ImagePlus object at 0x000002914889F5C0>
ips in remove() =  <imagepy.core.wraper.imageplus.ImagePlus object at 0x000002914889F5C0>
~~~ ~~~ !!! remove imgs !!!
~~~ ~~~ image add!
~~~ ~~~ ImageManager =  [<weakref at 0x0000029148895548; to 'ImagePlus' at 0x000002914889F5C0>]

从打印信息就可以一目了然地知道运行路径，即先add一张图像，然后再add另一张图像，这样ImageManager中就有两张图像引用，然后再remove第二张图像，这里注意，其实remove掉第二张照片后，第一张照片又会被add一下，然后再被remove，然后再insert到列表中，就像最后几行所示。原因就是因为第二张图像被删掉后，第一张图像的Canvas组件又被activated激活了，然后就会再调用add()方法。
其实不光是删掉图像后会激活另一张图像，也可以通过点击另一张图像来激活它，就会发现上面的过程又会重复一遍，即先删除图像，再添加它，这样带来的效果就是当前被激活的图像始终在ImageManager的imgs列表的第一个。
设计巧妙啊！