close

Вход

Забыли?

вход по аккаунту

?

О восстановлении функции яркости изображения методом полной вариации.

код для вставкиСкачать
Известия Тульского государственного университета
Естественные науки. 2015. Вып. 4. С. 19–30
Математика
УДК 517.97, 519.63
О восстановлении функции яркости
изображения методом полной вариации
Н. Х. Т. Данг, С. Д. Двоенко
Аннотация. Рассмотрена задача восстановления функции яркости исходного незашумленного изображения на основе метода полной
вариации. В настоящее время для создания изображений широко
используется различное современное цифровое оборудование. В работе цифрового оборудования возможно появление различных оптических дефектов. В первую очередь, качество изображений зависит
от оптических датчиков, где из-за технологических ограничений на
изображениях появляется шум. Шум снижает качество изображений
и качество результата их обработки. Для эффективного устранения
шумов разных типов разработаны специальные методы. Для электронных микроскопических изображений характерными считаются
гауссовский и пуассоновский шумы, а также их линейная комбинация. Решение задачи восстановления функции яркости позволяет эффективно устранить шумы данных типов в реальных изображениях.
Ключевые слова: задача восстановления функции яркости, полная вариация, модель ROF, гауссовский шум, пуассоновский шум,
обработка изображений, электронные микроскопические изображения, уравнение Эйлера – Лагранжа.
Введение
В настоящее время для создания изображений широко используется различное современное цифровое оборудование. В работе цифрового оборудования возможно появление различных оптических дефектов. В первую
очередь, качество изображений зависит от оптических датчиков, где из-за
технологических ограничений на изображениях появляется шум. Шум снижает качество изображений и качество результата их обработки. Поэтому
проблема устранения шума на изображениях является актуальной.
В настоящее время разработано много методов устранения шумов для
случаев, когда тип шума известен. Например, метод полной вариации [1-13]
является известным и эффективным методом.
Впервые, по-видимому, концепция полной вариации была применена для
устранения шума в работах Рудина [1, 2, 12, 13]. Он предложил использовать полную вариацию для решения задачи обработки изображений. Но его
20
Н. Х. Т. Данг, С. Д. Двоенко
модель ROF предназначена только для устранения гауссовского шума [12,
13].
Другим популярным шумом является пуассоновский шум. Например,
такой шум появляется в рентгеновских снимках. Модель ROF не сможет
эффективно устранить такой шум. Поэтому в работе Ли [14] была построена
другая модель, известная как модифицированная модель ROF.
Оба типа шумов (гауссовский и пуассоновский) распространены, но их
комбинация также важна [15]. Она появляется на биомедицинских изображениях, например, на изображениях электронной микроскопии [16, 17].
В данной статье предлагается совместно применить модель ROF и модифицированную модель ROF для устранения комбинации шумов.
В экспериментах были использованы реальные изображения с натуральным неизвестным шумом. Предполагалось, что это смесь гауссовского и пуассоновского шумов. Качество обработки сравнивалось с другими методами
устранения шума: модель ROF, модифицированная модель ROF, медианный
фильтр [18], фильтр Винера [19], регуляризация Бельтрами (Beltrami) [20].
Для сравнения качества изображений до и после их обработки был использован известный критерий BRISQUE (Blind/Referenceless Image Spatial QUality
Evaluator) [21].
1. Модель устранения шума
Пусть в пространстве R2 задана ограниченная область Ω ⊂ R2 . Назовем
функции u(x, y) ∈ R2 и v(x, y) ∈ R2 , соответственно идеальным (без шума) и
реальным (зашумлённым) изображениями, где (x, y) ∈ Ω.
Если функция u гладкая, то ее полная вариация имеет вид
Z
VT [u] = |∇u|dxdy,
Ω
где ∇u = (ux , uy ) — градиент, ux = ∂u/∂x, uy = ∂u/∂y, |∇u| =
q
u2x + u2y . В
этой работе мы считаем, что полная вариация функции u ограничена VT [u] <
< ∞.
Согласно работам [1, 2, 12, 13], гладкость изображений можно характеризовать их полной вариацией. Полная вариация зашумлённых изображений
всегда больше полной вариации соответствующих гладких изображений.
Для решения задачи VT [u] → min необходимо ввести ограничение на
вариацию гауссовского шума:
Z
(v − u)2 dxdy = const.
Ω
О восстановлении функции яркости изображения методом полной вариации
21
В этих условиях модель ROF для устранения гауссовского шума на изображении имеет вид [12]


Z
Z
λ
(v − u)2 dxdy 
u∗ = arg min  |∇u|dxdy +
2
u
Ω
Ω
как решение задачи безусловной оптимизации, где λ > 0 – множитель
Лагранжа.
Для устранения пуассоновского шума на основе модели ROF была предложена другая модель [14]. Такая модель получается при решении задачи
VT [u] → min с ограничением:
Z
Z
ln(p(v|u))dxdy = (u − v ln(u))dxdy = const
Ω
Ω
как решение задачи безусловной оптимизации:


Z
Z
u∗ = arg min  |∇u|dxdy + β (u − v ln(u))dxdy  ,
u
Ω
Ω
где β > 0 – коэффициент регуляризации. Такая модель известна как модифицированная модель ROF для устранения пуассоновского шума.
Чтобы построить модель устранения смешанного шума, будем также
решать задачу устранения шума, основанную на свойстве гладкости полной
вариации VT [u] → min.
Определим ограничения в новой задаче. Предполагается, что при заданном изображении вариация шума постоянна (пуассоновский шум не изменяется, а гауссовский шум зависит только от вариации шума):
Z
ln(p(v|u))dxdy = const,
(1)
Ω
где p(v|u) – условная вероятность наблюдения реального изображения v при
заданном идеальном изображении u.
Рассмотрим гауссовский шум. Его плотность распределения с дисперсией
σ 2 определяется как:
¶
µ
√
(v − u)2
/(σ
2π).
p1 (v|u) = exp −
2σ 2
Плотность пуассоновского шума
p2 (v|u) = exp(−u) uv /v!.
Обратим внимание, что значения функций яркости изображения u и
v – это целые числа (например, для восьмибитового изображения интервал
яркости определяется значениями от 0 до 255).
22
Н. Х. Т. Данг, С. Д. Двоенко
Для устранения комбинации гауссовского и пуассоновского шумов рассмотрим следующую линейную комбинацию:
ln(p(v|u)) = λ1 ln(p1 (v|u)) + λ2 ln(p2 (v|u)),
где λ1 > 0, λ2 > 0, λ1 + λ2 = 1.
Согласно (1), получим задачу устранения шума с ограничениями:
R
 ∗
 u = arg min |∇u|dxdy,
Ω
¢
R ¡ λ1 u
2 + λ (u − v ln(u)) dxdy = κ ,
(v
−
u)

2
2
2σ
Ω
где κ — постоянное значение.
Сведем эту задачу к задаче безусловной оптимизации с использованием
функционала Лагранжа


Z
Z
Z
λ
1
L(u, τ )= |∇u|dxdy + τ  2 (v − u)2 dxdy+ λ2 (u − v ln(u))dxdy − κ ,
2σ
Ω
Ω
Ω
чтобы найти решение:
(u∗ , τ ∗ ) = arg min L(u, τ ),
u, τ
(2)
где τ > 0 — множитель Лагранжа.
В данной модели, если λ 1 = 0 и λ2 = 1, то при β = τ λ2 = τ будет получена
модифицированная модель ROF для устранения пуассоновского шума. Если
λ2 = 0 и λ 1 = 1, то при λ = τ λ 1 /(2σ 2 ) = τ /(2σ 2 ) будет получена модель ROF
для устранения гауссовского шума. Если λ 1 > 0, λ2 > 0, то будет получена
модель для устранения смеси гауссовского и пуассоновского шумов.
2. Дискретная модель устранения шума
Для решения задачи (2) можно применить метод множителей Лагранжа
[22, 23, 24]. Его применение в данной задаче оказывается не совсем удобным,
поэтому в данной работе использовано уравнение Эйлера – Лагранжа [24].
Пусть функция f (x, y) определена в ограниченной области Ω ⊂ R2 и
непрерывно дифференцируема до второго порядка по x и y при (x, y) ∈ Ω.
Пусть F (x, y, f, fx , fy ) – выпуклый функционал, где fx = ∂f /∂x, fy = ∂f /∂y.
Решение задачи оптимизации
Z
F (x, y, f, fx , fy )dxdy → min
Ω
удовлетворяет уравнению Эйлера – Лагранжа:
Ff (x, y, f, fx , fy ) −
∂
∂
Ffx (x, y, f, fx , fy ) −
Ff (x, y, f, fx , fy ) = 0,
∂x
∂y y
О восстановлении функции яркости изображения методом полной вариации
23
где Ff = ∂F/∂f , Ffx = ∂F/∂fx , Ffy = ∂F/∂fy .
Тогда решение задачи (2) удовлетворяет следующему уравнению Эйлера – Лагранжа:




∂ 
λ1
v
u
 − µ ∂  q uy
 = 0, (3)
q x
− 2 (v − u) + λ2 (1 − ) − µ
σ
u
∂x
∂y
2
2
2
2
u +u
u +u
x
y
x
y
где µ = 1/τ .
Представим уравнение (3) в следующем виде:
uxx u2y − 2ux uy uxy + u2x uyy
v
λ1
= 0,
(v
−
u)
−
λ
(1
−
)
+
µ
2
σ2
u
(u2x + u2y )3/2
где
uxx
∂2u
=
,
∂x2
uyy
∂2u
=
,
∂y 2
uxy
∂
=
∂x
µ
∂u
∂y
(4)
¶
= uyx .
Для получения дискретной модели (4) добавим искусственный параметр
времени u = u(x, y, t). Уравнение (4) соответствует следующему уравнению
диффузии:
ut =
uxx u2y − 2ux uy uxy + u2x uyy
v
∂u
λ1
.
= 2 (v − u) − λ2 (1 − ) + µ
∂t
σ
u
(u2x + u2y )3/2
(5)
Пусть размер изображения N1 × N2 . Тогда дискретная форма уравнения
(5) имеет вид
!
µ
v
λ
i
j
1
k
(vi j − ukij ) − λ2 (1 − k ) + µφkij ,
(6)
uk+1
i j = ui j + ξ
σ2
ui j
φkij =
∇xx (ukij )(∇y (ukij ))2
((∇x (ukij ))2 + (∇y (ukij ))2 )3/2
+
+
−2∇x (ukij )∇y (ukij )∇xy (ukij ) + (∇x (ukij ))2 ∇yy (ukij )
∇x (ukij ) =
((∇x (ukij ))2 + (∇y (ukij ))2 )3/2
uki+1,j − uki−1,j
,
2∆x
∇y (ukij ) =
uki,j+1 − uki,j−1
,
2∆y
∇xx (ukij ) =
uki+1,j − 2ukij + uki−1,j
,
(∆x)2
∇yy (ukij ) =
uki,j+1 − 2ukij + uki,j−1
,
(∆y)2
,
24
Н. Х. Т. Данг, С. Д. Двоенко
∇xy (ukij ) =
uki+1,j+1 − uki+1,j−1 − uki−1,j+1 + uki−1,j−1
,
4∆x∆y
uk0j = uk1j ; ukN1 +1,j = ukN1 ,j ; uki0 = uki1 ; uki,N2 +1 = uki,N2 ;
i = 1, . . . , N1 ; j = 1, . . . , N2 ;
k = 0, 1, . . . , K; ∆x = ∆y = 1; 0 < ξ < 1, где K — достаточно большое число,
K = 500.
3. Параметры модели устранения шума
Процедура (6) может быть использована для устранения шума на изображениях, если значения параметров λ1 , λ2 , µ, σ заданы. Часто на практике
эти параметры неизвестны, и их нужно оценить. Тогда параметры λ1 , λ2 , µ
в процессе (6) нужно представить как λk1 , λk2 , µk на каждой итерации k.
В новой процедуре такие параметры будут вычисляться на каждом шаге
итерации.
3.1. Оптимальные параметры λ1 и λ2
Пусть u, τ являются решением задачи (2). Тогда мы получим условие
∂L(u, τ )/∂u = 0. Данное условие позволяет вычислить оптимальные параметры линейной комбинации шумов λ 1 , λ2 :
R
(1 − uv )dxdy
Ω
R
, λ2 = 1 − λ 1 .
λ1 = 1 R
(v
−
u)dxdy
+ (1 − uv )dxdy
2
σ
Ω
Ω
Дискретная форма для вычисления параметров имеет вид
N1 P
N2
P
λk1
=
(1 −
i=1 j=1
N1 P
N2
P
(
i=1 j=1
vij −ukij
σ2
vij
ukij
)
+1−
,
vij
ukij
λk2 = 1 − λk1 ,
)
где k = 0, 1, . . . , K.
3.2. Оптимальный параметр µ
Для поиска оптимального параметра сглаживания µ умножим (3) на (v −
− u) и проинтегрируем по частям по области Ω. В итоге получим выражение
для оптимального параметра µ:
R
2
)dxdy
(− λσ21 (v − u)2 − λ2 (v−u)
u
Ω
.
µ= R q
ux vx +uy vy
)dxdy
( u2x + u2y − √
2
2
Ω
ux +uy
О восстановлении функции яркости изображения методом полной вариации
25
Его дискретная форма имеет вид
N1 P
N2
P
µk =
λk
(− σ12 (vij − ukij )2 − λk2
i=1 j=1
N1 P
N2
P
i=1 j=1
(vij −ukij )2
)
ukij
,
k
ηij
где
k
ηij
=
q
∇x (ukij )∇x (vij ) + ∇y (ukij )∇y (vij )
q
(∇x (ukij ))2 + (∇y (ukij ))2 −
,
(∇x (ukij ))2 + (∇y (ukij ))2
∇x (ukij ) =
uki+1,j − uki−1,j
,
2∆x
∇y (ukij ) =
uki,j+1 − uki,j−1
,
2∆y
k
∇x (vij
)=
k
k
vi+1,j
− vi−1,j
v,
2∆x
k
∇y (vij
)=
k
k
vi,j+1
− vi,j−1
,
2∆y
uk0j = uk1j ; ukN1 +1,j = ukN1 ,j ; uki0 = uki1 , uki,N2 +1 = uki,N2 ;
v0j = v1j ; vN1 +1,j = vN1 j ; vi0 = vi1 ; vi,N2 +1 = vi,N2 ;
i = 1, ..., N1 , j = 1, ..., N2 ; k = 0, 1, ..., K; ∆x = ∆y = 1.
3.3. Оптимальный параметр σ
Для вычисления параметра σ здесь использован метод Иммеркера [25]:
p
N1 X
N2
X
π/2
|uij ∗ Λ|,
(7)
σ=
6(N1 − 2)(N2 − 2)
i=1 j=1
Ã
1 −2 1
−2 4 −2
1 −2 1
конволюции (свертки):
где Λ =
!
– маска изображения. Оператор * – это оператор
uij ∗ Λ = ui−1,j−1 Λ33 + ui,j−1 Λ32 + ui+1,j−1 Λ31 + ui−1,j Λ23 +
+uij Λ22 + ui+1,j Λ21 + ui−1,j+1 Λ13 + ui,j+1 Λ12 + ui+1,j+1 Λ11 ,
где i = 1, . . . , N1 ; j = 1, . . . , N2 ; uij = 0, если i = 0, или j = 0, или i = N1 + 1,
или j = N2 + 1. Параметр σ вычисляется сразу же на первой итерации.
4. Поиск начального решения
Очевидно, что в локальной итерационной процедуре (6) результат обработки в общем случае зависит от начальных значений параметров λ01 , λ02 , µ0 .
26
Н. Х. Т. Данг, С. Д. Двоенко
Если сначала задать параметры λ01 , λ02 , µ0 , то их неудачные значения
определят не очень хорошие оценки uij , а через них также не очень хорошие
оценки параметров распределений.
Случайный выбор параметров λ01 , λ02 , µ0 также неприемлем, так как
фактически вносит дополнительный шум в изображение.
Очевидно, что начальные значения параметров λ01 , λ02 , µ0 должны быть,
по-возможности, достаточно близки к тем значениям, которые будут найдены. Поэтому параметры λ01 , λ02 , µ0 оцениваются в данной работе как средние
по соседним пикселам изображения с использованием метода Иммеркера.
5. Оценка качества изображений
Оценка качества реальных зашумленных изображений, для которых
неизвестно, как они выглядят без шума (отсутствуют т.н. опорные изображения) – это трудная теоретическая задача. В настоящий момент разработано
несколько методов для решения такой задачи. В данной статье используется критерий BRISQUE, предложенный Анишом (Anish) в работе [21]. Чем
больше значение критерия QBRISQUE, тем лучше качество изображения.
Оценка качества изображения на основе метода BRISQUE заключается
в выполнении следующих шагов:
1) поиск специальных коэффициентов, представляющих т.н. натуральную статическую сцену (НСС) для зашумленного изображения;
2) построение специального вектора, компоненты которого получены из
коэффициентов НСС;
3) применение методов регрессии для оценки значения критерия качества
изображения.
6. Эксперимент
Предложенная модель была протестирована на реальных изображениях
с неизвестным естественным шумом. Например, здесь показан результат
обработки изображения «aphid» (рис. 1а) из базы данных Университета
Альберты (Канада) [26]. Данное изображение получено на электронном микроскопе. На рис. 1б показан увеличенный фрагмент.
Результат устранения шума показан на рис. 1в. На рис. 1д также показан
увеличенный фрагмент обработанного изображения.
При обработке предполагалось, что шум на изображении представляет
собой смесь гауссовского и пуассоновского шумов в неизвестной пропорции.
Поэтому в ходе обработки были автоматически определены следующие значения параметров: λ1 = 0.4042, λ2 = 0.5958, µ = 0.5238, σ = 8.3884. Значение
QBRISQUE показано в табл. 1.
О восстановлении функции яркости изображения методом полной вариации
27
Рис. 1. Устранение шума на реальном изображении: а) – реальное
изображение; б) – увеличенный фрагмент реального изображения;
в) – изображение после устранения шума; д)– увеличенный фрагмент
изображения после устранения шума
Таблица 1
Сравнение качества методов устранения шума на реальном изображении
по критерию QBRISQUE
Шум
ROF
20.0455 29.0058
Модиф.
ROF
27.5682
Медиана Винера
фильтр фильтр
30.2719 23.4842
Бельтрами
метод
31.9801
Предложенный
метод
35.7603
Заключение
В работе предложен метод устранения смеси гауссовского и пуассоновского шумов на основе метода полной вариации.
Качество результата устранения шума зависит от значений коэффициентов линейной комбинации λ1 , λ2 и коэффициента сглаживания µ.
Показано, что качество обработки реальных изображений с неизвестным
типом шума оказывается лучше качества обработки методами, специально
предназначенными для устранения шума только одного вида.
Список литературы
1. Chan T.F., Shen J. Image processing and analysis: Variational, PDE, Wavelet and
stochastic methods. SIAM, 2005. 400 p.
2. Level set and PDE based reconstruction methods in imaging / M. Burger [et al.].
Springer, 2008. 319 p.
28
Н. Х. Т. Данг, С. Д. Двоенко
3. An introduction to total variation for image analysis / A. Chambolle [et al.] //
Theoretical foundations and numerical methods for sparse recovery. 2009. V. 9.
P. 263–340.
4. Xu J., Feng X., Hao Y. A coupled variational model for image denoising using
a duality strategy and split Bregman // Multidimensional systems and signal
processing. 2014. V. 25. P. 83–94.
5. Rankovic N., Tuba M . Improved adaptive median filter for denoising ultrasound
images // Advances in computer science, WSEAS ECC’12. 2012. P. 169–174.
6. Lysaker M., Tai X . Iterative image restoration combining total variation
minimization and a second-order functional // International journal of computer
vision. 2006. V. 66. P. 5–18.
7. Li F., Shen C., Pi L . A new diffusion-based variational model for image denoising
and segmentation // Journal mathematical imaging and vision. 2006. V. 26. No. 1-2.
P. 115–125.
8. Noise reduction with low dose CT data based on a modified ROF model / Y. Zhu
[et al.] // Optics express. 2012. V. 20. No. 16. P. 17987–18004.
9. Tran M.P., Peteri R., Bergounioux M. Denoising 3D medical images using a second
order variational model and wavelet shrinkage // Image analysis and recognition.
2012. V. 7325. P. 138–145.
10. Getreuer P. Rudin – Osher – Fatemi total variation denoising using split Bregman //
IPOL 2012. URL: http://www.ipol.im/pub/art/2012/g-tvd/ (дата обращения:
23.07.2015).
11. Caselles V., Chambolle A., Novaga M. Handbook of mathematical methods in
imaging. Springer, 2011. 1607 p.
12. Rudin L.I., Osher S., Fatemi E. Nonlinear total variation based noise removal
algorithms // Physica D. 1992. V. 60. P. 259–268.
13. Chen K. Introduction to variational image processing models and application //
International Journal of computer mathematics. 2013. V. 90. No. 1. P. 1–8.
14. Le T., Chartrand R., Asaki T.J. A variational approach to reconstructing images
corrupted by Poisson noise // Journal of mathematical imaging and vision. 2007.
V. 27. No. 3. P. 257–263.
15. Luisier F., Blu T., Unser M. Image denoising in mixed Poisson-Gaussian noise //
IEEE transaction on Image processing. 2011. V. 20. No. 3. P. 696–708.
16. Poisson-Gaussian noise parameter estimation in fluorescence microscopy imaging /
A. Jezierska [et al.] // IEEE International Symposium on Biomedical Imaging 9th.
2012. P. 1663–1666.
17. An EM approach for Poisson-Gaussian noise modeling / A. Jezierska [et al.] //
EUSIPCO 19th. 2011. V. 62. No. 1. P. 13–30.
18. Wang C., Li T. An improved adaptive median filter for Image denoising // ICCEE.
2012. V. 53. No. 2.64. P. 393–398.
19. Abe C., Shimamura T. Iterative Edge-Preserving adaptive Wiener filter for image
denoising // ICCEE. 2012. V. 4. No. 4. P. 503–506.
20. Zosso D., Bustin A. A Primal-Dual Projected Gradient Algorithm for Efficient
Beltrami Regularization. Computer Vision and Image Understanding, 2014.
http://www.math.ucla.edu/∼zosso/ (дата обращения: 23.07.2015).
О восстановлении функции яркости изображения методом полной вариации
29
21. Mittal A., Moorthy A.K., Bovik A.C. No-Reference Image Quality Assessment in
the Spatial Domain // IEEE Transactions on Image Processing. 2012. V. 21. No. 12.
P. 4695–4708.
22. Zeidler E. Nonlinear functional analysis and its applications: Variational methods
and optimization. Springer, 1985. 662 p.
23. Rubinov A., Yang X. Applied Optimization: Lagrange-type functions in constrained
non-convex optimization. Springer, 2003. 286 p.
24. Gill P.E., Murray W. Numerical methods for constrained optimization. Academic
Press Inc., 1974. 283 p.
25. Immerker J . Fast noise variance estimation // Computer vision and image
understanding. 1996. V. 64. No. 2. P. 300–302.
26. База данных биомедицинских изображений университета Алберта (Alberta
University), Канада: http://www.ualberta.ca/∼mingchen/index.htm (дата обращения: 23.07.2015).
Данг Нгок Хоанг Тхань (myhoangthanh@yahoo.com), аспирант, кафедра
информационной безопасности, Институт прикладной математики и компьютерных наук, Тульский государственный университет.
Двоенко Сергей Данилович (dsd@tsu.tula.ru), д.ф.-м.н., профессор, кафедра информационной безопасности, Институт прикладной математики и
компьютерных наук, Тульский государственный университет.
On restoration of an image intensity function based on the total
variation method
N. H. T. Dang, S. D. Dvoenko
Abstract. we consider the restoration problem of intensity function of an
original noiseless image based on the total variation method. Today, there are
many modern devices to create digital images. These devices have optical defects
to create images. Hence, the quality of images depends on the quality of optical
sensors. Because of the limit of technology, such defects are always included, for
example, as noise. Noise reduces image quality and result of image processing. In
order to reduce noises effectively, many special methods are used. For electronic
microscope images the noise includes Gaussian noise, Poisson noise and also their
linear combination. The solution of the restoration problem of image intensity
function allows us to remove those noises more effectively.
Keywords: restoration problem of intensity function, total variation, ROF
model, Gaussian noise, Poisson noise, image processing, biomedical image, Euler –
Lagrange equation.
30
Н. Х. Т. Данг, С. Д. Двоенко
Dang Ngoc Hoang Thanh (myhoangthanh@yahoo.com), postgraduate
student, department of Information security, Institute of Applied Mathematics
and Computer Sciences, Tula State University.
Dvoenko Sergey (dsd@tsu.tula.ru), doctor of physical and mathematical
sciences, professor, department of Information security, Institute of Applied
Mathematics and Computer Sciences, Tula State University.
Поступила 18.09.2015
Документ
Категория
Без категории
Просмотров
10
Размер файла
878 Кб
Теги
методов, восстановлен, полное, функции, изображение, яркость, вариаций
1/--страниц
Пожаловаться на содержимое документа