1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
// 图像分割 16-bit/float 变体
// watershed / kmeans / SLIC 支持 Image16 / ImageF
// ========== 辅助函数 ==========
///|
fn read_u16_le_seg(data : Bytes, offset : Int) -> Int {
data[offset].to_int() + data[offset + 1].to_int() * 256
}
///|
fn read_float_le_seg(data : Bytes, offset : Int) -> Float {
let b0 = data[offset].to_int()
let b1 = data[offset + 1].to_int()
let b2 = data[offset + 2].to_int()
let b3 = data[offset + 3].to_int()
let bits = b0 | (b1 << 8) | (b2 << 16) | (b3 << 24)
Float::reinterpret_from_int(bits)
}
///|
fn clamp_int_seg(v : Int, lo : Int, hi : Int) -> Int {
if v < lo {
lo
} else if v > hi {
hi
} else {
v
}
}
///|
fn img16_to_img8(img : @types.Image16) -> @types.Image {
let w = img.width
let h = img.height
let ch = img.channels
let data8 : Array[Byte] = Array::make(w * h * ch, b'\x00')
for i in 0..<(w * h * ch) {
let v16 = read_u16_le_seg(img.data, i * 2)
data8[i] = (v16 / 257).to_byte()
}
@types.Image::{
width: w,
height: h,
channels: ch,
data: Bytes::from_array(data8),
}
}
///|
fn imgf_to_img8(img : @types.ImageF) -> @types.Image {
let w = img.width
let h = img.height
let ch = img.channels
let data8 : Array[Byte] = Array::make(w * h * ch, b'\x00')
for i in 0..<(w * h * ch) {
let vf = read_float_le_seg(img.data, i * 4)
let v8 = clamp_int_seg((vf * 255.0F).to_int(), 0, 255)
data8[i] = v8.to_byte()
}
@types.Image::{
width: w,
height: h,
channels: ch,
data: Bytes::from_array(data8),
}
}
// ========== Watershed 16-bit ==========
///|
/// 16-bit Watershed 分割。
pub fn watershed_16(
img : @types.Image16,
markers : Array[Int],
) -> Array[Int] raise @types.LoadError {
let img8 = img16_to_img8(img)
watershed(img8, markers)
}
///|
/// 16-bit 自动 Watershed 分割。
pub fn watershed_auto_16(
img : @types.Image16,
) -> (Array[Int], Int) raise @types.LoadError {
let img8 = img16_to_img8(img)
watershed_auto(img8)
}
// ========== Watershed float ==========
///|
/// float Watershed 分割。
pub fn watershed_f(
img : @types.ImageF,
markers : Array[Int],
) -> Array[Int] raise @types.LoadError {
let img8 = imgf_to_img8(img)
watershed(img8, markers)
}
///|
/// float 自动 Watershed 分割。
pub fn watershed_auto_f(
img : @types.ImageF,
) -> (Array[Int], Int) raise @types.LoadError {
let img8 = imgf_to_img8(img)
watershed_auto(img8)
}
// ========== K-means 16-bit ==========
///|
/// 16-bit K-means 图像分割。
pub fn kmeans_segment_16(
img : @types.Image16,
k : Int,
max_iters? : Int = 20,
) -> (SegmentLabelImage, Array[SegmentRegion]) raise @types.LoadError {
let img8 = img16_to_img8(img)
kmeans_segment(img8, k, max_iters~)
}
// ========== K-means float ==========
///|
/// float K-means 图像分割。
pub fn kmeans_segment_f(
img : @types.ImageF,
k : Int,
max_iters? : Int = 20,
) -> (SegmentLabelImage, Array[SegmentRegion]) raise @types.LoadError {
let img8 = imgf_to_img8(img)
kmeans_segment(img8, k, max_iters~)
}
// ========== SLIC 16-bit ==========
///|
/// 16-bit SLIC 超像素分割。
pub fn slic_16(
img : @types.Image16,
k : Int,
m : Int,
max_iters : Int,
) -> SuperpixelResult raise @types.LoadError {
let img8 = img16_to_img8(img)
slic(img8, k, m, max_iters)
}
// ========== SLIC float ==========
///|
/// float SLIC 超像素分割。
pub fn slic_f(
img : @types.ImageF,
k : Int,
m : Int,
max_iters : Int,
) -> SuperpixelResult raise @types.LoadError {
let img8 = imgf_to_img8(img)
slic(img8, k, m, max_iters)
}