CC BY
12Zhilyakov E.G., Likholob P.G., Lichogodina E.S. Steganographic method encoding of geodata in the images of the earth's surface// Сетевой журнал «Научный результат».
Серия «Информационные технологии». - Т.1, №2,2016.
UDC 621.397
DOI 10.18413/2518-1092-2016-1-2-64-70
Zhilyakov E.G. Likholob P.G. Lichogodina E.S.
STEGANOGRAPHIC METHOD ENCODING OF GEODATA IN THE IMAGES OF THE EARTH'S SURFACE
1) Doctor of Technical Sciences, Professor, Head of Department of Information and Telecommunication Systems and Technologies, Belgorod State National Research University, 85 Pobedy St., Belgorod, 308015, Russia
e-mail: Zhilyakov@bsu.edu.ru 2) Senior Lecturer. Department of Information and Telecommunication Systems and Technologies Belgorod State National Research University, 85 Pobeda St., Belgorod, 308015, Russia. e-mail: likholob@bsu.edu.ru
3) Student. Department of Information and Telecommunication Systems and Technologies Belgorod State National Research University, 85 Pobeda St., Belgorod, 308015, Russia. e-mail: elza9313@gmail.ru
Abstract
This article presents the possibility of implementing geodata in images of the earth surface using the methods of steganography. The authors considered the steganographic modified method of Koch-Zhao and proposed a subband method that allows to put geodata in the image, introducing fewer distortions.
Keywords: geodata; steganography; Koch-Zhao method; subband method.
УДК 621.397
Жиляков Е.Г. Лихолоб П.Г. Лихогодина Е.С.
СТЕГАНОГРАФИЧЕСКИЙ МЕТОД КОДИРОВАНИЯ ГЕОДАННЫХ В СНИМКИ ЗЕМНОЙ ПОВЕРХНОСТИ
1) заведующий кафедрой информационно-телекоммуникационных систем и технологий, доктор технических наук, профессор. Белгородский государственный национальный исследовательский университет, ул. Победы д.85,
г. Белгород, 308015, Россия. e-mail: Zhilyakov@bsu.edu.ru 2) старший преподаватель кафедры информационно-телекоммуникационных систем и технологий Белгородский государственный национальный исследовательский университет, ул. Победы д.85, г. Белгород, 308015, Россия. e-mail: likholob@bsu.edu.ru 3) студент кафедры информационно-телекоммуникационных систем и технологий Белгородский государственный национальный исследовательский университет, ул. Победы д.85, г. Белгород, 308015, Россия. e-mail: elza9313@gmail.ru
Аннотация
В данной статье рассматривается возможность внедрения геоданных в снимки земной поверхности при помощи методов стеганографии. Рассмотрен стеганографический модифицированный метод Коха-Жао, а также предложен субполосный метод кодирования информации, позволяющий помещать геоданные в изображение с внесением меньшего числа искажений.
Ключевые слова: геоданные; стеганография; метод Коха-Жао; субполосный метод.
Images, obtained with the help of aerial photographs and artificial satellites, play a significant role in modern cartography. Increasingly, they are used to find the optimal route, analyzing the state of the objects and their identification. Currently, the use of imagery of the earth surface is typical for different spheres of human activities: nature management, agriculture, environmental safety, the land registry, navigation, modeling, cartography, monitoring, energy, etc. Images, obtained with the help of aerial
photographs and artificial satellites for analysis added to the database, comparing them with geo-information. Geo-information is the information about geographic location and condition of features, as well as about changes in their condition.
Existing technologies of storage geo-information using database and the corresponding image noticeably increase the traffic to mobile devices when using geo-information. Thus, it seems reasonable to store the image information (geodata) in the image,
Zhilyakov E.G., Likholob P.G., Lichogodina E.S. Steganographic method encoding of geodata in the images of the earth's surface// Сетевой журнал «Научный результат».
Серия «Информационные технологии». - Т.1, №2,2016.
using secretive coding, without introducing any user-visible distortion and without affecting the system of
automatic identification. Methods to carry out covert encoding are called steganography [1, 2].
Fig. 1. The scheme of encoding information
Fig. 2. The scheme of decoding information
Figures 1 and 2 represent a scheme of encoding and decoding geodata in the image.
For concealing information for coding in the image will be used matrix notation, the model is presented below.
а)
X - image size [/ x K ]
Ф - fragment of image in the frequency do:
Ф1,1 Ф1, 2 ф1,(а-1) Ф1,а
Фц ф2 .2 Ф2,а
Фи. я
Ф(П-1), 1
фи, 1 фи, 2 Фп,а
Fig. 3. Matrix representation of image: a) image; b) the fragment of image in frequency representation;
c) fragment of image size
Zhilyakov E.G., Likholob P.G., Lichogodina E.S. Steganographic method encoding of geodata in the images of the earth's surface// Сетевой журнал «Научный результат».
Серия «Информационные технологии». - Т.1, №2,2016.
Figure 3a presents in numerical form, the result of the registration of energy [3] from the visual image of an object or phenomenon, stored, transmitted or used for further reproduction in the form of pixels. The pixel means the smallest indivisible element of the registered image in the spatial domain, which is a quantitative estimate of the energy in the recorded frequency band.
In practice of digital image processing [4-7] there are methods based not on the image itself, but on its discrete cosine transform (DCT). There is a direct discrete cosine transform and an inverse discrete cosine transform [8-11].
Direct DCT:
$ = D-F • DT (1)
where F - a block image size [n x M], with pixel values f ;
J n,m '
Π- a block of DCT coefficients size [u x q], with values ^ ;
D - a DCT matrix, d = ^uai}, u = \,2,...,U, c = \,2,..., Q .
Inverse DCT:
f = dr • o^d, (2)
where F - a block of image retrieved from the frequency coefficients of DCT.
The matrix D stores a set of cosine functions in the discrete form, which are necessary to carry out an orthogonal transformation image F (1) and (2) to DCT.
d (и, a>) =
N M
■Jn x M
C(и)c(a) 2 2 cos
n = 1m = 1
n(2n + l)u cos n(2m +
2N 2M
(3)
where u , C - step selectable frequencies
u = \,2,...,U C = \,2,..., Q ; C(l) - coefficients: '=\
C(l) = w 2 , where l corresponds u or C . \, l > \
(4)
The method of steganography of Koch-Zhao is based on this mathematical technique. The method consists in encoding bits of information by replacing a valid value of the DCT coefficients selected after preliminary analysis of the frequency domain of the image. The preliminary analysis allowing selecting frequency coordinates is known before a decisive rule. Encoding and decoding are performed in the DCT coefficients.
Figure 4 represents a known frequency image model, which carried out the separation of the matrix of coefficients of DCT on the frequency domain. Thus, the lowest frequency is located in the upper left corner of the image, and the highest - in the lower right. It is known that the most effective is the encoding of information in the field of mid frequencies, because low frequencies contain the most part of the energy of the image, and high frequencies - are most susceptible to the distortion of a communications channel and more easily influenced during processing.
Evaluation of the distortions introduced in the image after secretive encoding are determined by the correlation equation (5) that evaluates the similarity measure between the initial image fragment and the fragment with the covertly encoded information.
P = -
N M f 1 N M \ N M f 1 N M
22\f ——22F 122IF ——22F
n=1 m=1 V N ■ M n=1 m=1 ) n=1 m=1 V N • M n=1 m=1
i
N M í 1 N M V N M í 1 N M
221 f —— ^2F ^I F —— 22^
22\ n,m ЛТ , , 22 n,m 22\ n,m Лт , , 22 n
(5)
where Fn m - a pixel of the original image fragment F ;
Fn,m a pixel with the covertly encoded information in the image fragment F ;
p - the degree of correlation between the image F and the image F .
1
2
n=1 m=1
n=1 m=1
n=1 m=1
n=1 m=1
Zhilyakov E.G., Likholob P.G., Lichogodina E.S. Steganographic method encoding of geodata in the images of the earth's surface// Сетевой журнал «Научный результат».
Серия «Информационные технологии». - Т.1, №2,2016.
Low frequency domain The constant component
ii--» -----
Фи ф1,2 Ф1,3 ... Ф1, 63 Ф,64 'I
Ф2 , 1 Ф2,2 ... •<7 Ф2,64
Фз,1 ... / , / / /" /
у' ... Ф5,6 ... У ' /
У У ... Ф6,5 У У ' . . . , /...
/ ... У к' Ф7,6 у ... у
... У У У ... Ф10,30 л
... / ... ( ^30,10
Ф 63,1 ... Z-
Ф64, 1 Ф64,2 ... ** Ф64, 64
t Mid frequency domain
High frequency domain
Fig. 4. Frequency model of image
The results of calculations are presented in table 1. The main drawback of the method and, in general, the methods used to encode the frequency region, is the use of DCT for all pixels of the image that does not allow spatial selectivity. As can be seen in figure 5, when changing the coefficients located in the low frequency region, and with increasing frequency
encoding (Fig. 5 b-d) in the image there appears a significant distortion, such as blemishes and ripples. Since the ratio (1; 1) contains information about the constant component of the signal, as it changes the image into a completely shaded one (Fig. 5b). This feature is connected with the energy distribution in the image.
Fig. 5. The fragment of image: a) input; b) by inverting the DCT coefficients with coordinates (1,1); c) by inverting the DCT coefficients with coordinates (3,1); d) by inverting the DCT coefficients with coordinates (6,7); e) by inverting the DCT
coefficients with coordinates (45,46)
In the process of the preliminary research, it was noted that the information encoded in the high frequency and mid-frequency regions is less resistant to external influences compared with information encoded in the low frequency region. However, as can be seen in table 1 and the previously shown Fig. 5 b-d, distortion in the high frequency and mid-frequency regions is significant. Therefore, it is necessary to develop a more selective method of encoding covert, in which all transformations and
analysis will occur in a spatial region with minimal losses.
The mathematical basis of the proposed method is a subband transformation [3, 4].
The direct subband transformation is performed by using a receiving subband projection (6), wherein the core of the subband matrix has a form (7). The reverse transformation takes place by a line-by-line scalar product of the transpose of the matrix of eigenvectors for the matrix of subband projection (8).
Zhilyakov E.G., Likholob P.G., Lichogodina E.S. Steganographic method encoding of geodata in the images of the earth's surface // Сетевой журнал «Научный результат».
Серия «Информационные технологии». - Т.1, №2,2016.
Table 1
Evaluation of the distortion caused by embedding
sin
Frequency domain № of experiment Coefficients Correlation, P
u a
1 2 3 4 5
LF 1 1 2 0.9864
2 1
LF 2 1 1 œ
2 2
LF 3 1 3 0.8072
3 1
LF 4 5 6 0.9699
6 5
LF 5 6 5 0.9865
7 6
LF 7 9 10 0.9978
10 9
MF 6 25 12 0.9981
26 13
MF 8 10 30 0.9998
30 10
MF 9 37 38 0.9999
38 37
HF 10 45 46 Close to 1
46 45
HF 11 62 63 Close to 1
63 62
- 0, / \ 2,r 1 ,r /,, ^ I
—,-— (u — a)
B = 2—^-7-г-^ cos
n(u — a)
v, + v,
2,r 1,r
2
(u — a)
(7)
where u, c - the step of the discrete frequencies u = 1,2,...,U, c = 1 ,2...., Q.
l q = b q ,
r ^ r r '
where Lr - the matrix of eigenvalues. Reverse subband transformation:
R
F = Z( (Q r , Ar):
(8)
The direct subband transformation (generate subband projection):
Ar =<Qr,F , r = 1,2,3,^,R (6)
where F - a block of image size [n x M] with the
pixel values f ;
A J n,m 7
Ar - a the vector-matrix containing the values of subband projections;
Qr - the matrix, the rows of which are eigenvectors of the subband matrix [4, 5].
The core of the subband matrix has a form:
where F - a block of image restored from the values of subband projections.
( )r - the transpose operation; - the row-wise dot product.
Figure 6 proposes a model image in subband domain, given for one line of the image. To reduce the seepage in the adjacent frequency domain eigenvalues of the subband matrix it is advisable to pick them up in such a way that they would be close to one. Therefore, for the image shown earlier, and the implementation of the encoding in a single string, the frequency band was divided into R=12 frequency intervals, with the first two numbers were close to one, respectively, the coding was carried out in the corresponding projections.
The subband method of secretive coding of an image is performed in the spatial domain of the image without transition frequency in (9), but the number of decision rules allows to choose the frequency band for encoding, the degree of infiltration in an adjacent lane, the amount of introduced distortion, as well as to determine the stability of geodata to distortions (e.g., noise).
Fig. 6. Model of image in frequency domain
Zhilyakov E.G., Likholob P.G., Lichogodina E.S. Steganographic method encoding of geodata in the images of the earth's surface // Сетевой журнал «Научный результат».
Серия «Информационные технологии». - Т.1, №2,2016.
Decoding is the same in the spatial domain.
It is worth noting that the proposed method of subband coding allows recovering the original pixel values, knowing only decision rules that have been coded. The difference in energy between the restored and original image, according to the experiment is 10-16.
The encoding algorithm:
1. The image represented in a digital form that is split into blocks F [nxM] pixels.
2. To represent the geodata in the form of symbols ew e {l,-l} .
3. To choose the number of frequency regions
R.
4. To calculate the subband matrix l q = b q .
r ^ r r
5. To find decomposition of the subband matrix B on the eigenvectors q = {5^} and eigenvalues L , r = 2,3,...,R, c = 1.2.. ...Q.
the projection of the line:
6. To find
a = (qr, f ,
( = 1,2...,О..
7. According
A,
r = 2, 3,...,R,
to several decision rules, to
choose subband projection: агш .
8. According to the encoded symbol
encode:
F = F + К• -lal-qr
to
(9)
The decoding algorithm:
1. The image represented in a digital form that is split into blocks F [nxM] pixels.
2. To represent the geodata in the form of
symbols ew e
{1,-1}-
3. To choose the number of frequency regions R (must be matched with the encoder).
4. To calculate the subband matrix L q = b q .
5. To find decomposition of the subband matrix on the eigenvectors Q = {q^} and eigenvalues r = 2, 3,.,R, c = 1,2..,Q.
6. To find the projection of the line:
r = 2,3,., R,
B r
L
A,
Ar ={ar,J'
=( qr, f, c = 1.2.... Q.
7. According to several decision rules, to choose subband projection: arm .
8. To carry out decoding, ew .
The proposed algorithm allows carrying out the adaptive subband encoding in the specified frequency-spatial domain. It should be noted that for coding information there must be used one-dimensional decomposition. Secretive coding was carried out for all lines, therefore, the throughput of the method above is not less than 64 times more effective in comparison with the method of Koch-Zhao.
Figure 7 and table 2 show the results of a computational experiment for steganographic encoding of information in the image using subband projections. One can see that the distortion is manifested in the form of "ripples" (spatial harmonic signal). The distortion in the form of the ripples is less visible on the fragments having the same intensity, in the case of change of the intensity the ripples increase (figure 7 c, d), and correlation decreases. To minimize these distortions it is necessary to develop additional decision rules, which can realize the selectivity in the spatial domain.
а)
Fig. 7. The fragment of image: a) input; b) by inverting the subband coefficients with coordinates (1,1) - low frequency domain; c) by inverting the subband coefficients with coordinates (10,1) - high frequency domain; d) by inverting the subband coefficients with coordinates (2,1) - low frequency domain; e) by inverting the subband coefficients with coordinates (4,1) -
mid-frequency domain
e
w
Zhilyakov E.G., Likholob P.G., Lichogodina E.S. Steganographic method encoding of geodata in the images of the earth's surface // Сетевой журнал «Научный результат».
Серия «Информационные технологии». - Т.1, №2,2016.
Table 2 presents the results of the experiment for the method of subband projections. By encoding the same portion of an image, distortion in all frequency regions becomes much fewer and, consequently, the correlation increases.
Table 2
Evaluation of the distortion caused by embedding the information using the method of subband
Thus, the proposed steganographic method in digital image processing allows for encoding geodata directly in the image without noticeable deterioration of visual quality.
The advantage of the proposed method compared to the existing method is in the increased bandwidth when using fragments greater than 64 pixels, the smaller degree of distortion, the possibility of selectivity in the spatial and frequency domain, inherent in the mathematical model. This is achieved through a mathematical apparatus and the proposed method of encoding.
Работа выполнена при поддержке гранта РФФИ 15-07-01570 (Субполосная скрытная интеграция/извлечение дополнительной
информации в аудио или видео контентах) / The work is executed with the support of RFBR grant 15-07-01570 (Subband secretive
integration/extraction of additional information in audio or video content).
References
1. Zhilyakov E.G. O About the Subband Properties of Images / E.G. Zhilyakov, S.P. Belov, A.A. Chernomorets, E.V. Bolgova // Nauchnye vedomosti BelGU. Ser. Istorija. Politologija. Jekonomika. Informatika. 2013. №8 (151). V.26/1. Pp. 175-182.
2. Zhilyakov E.G. The Secretive Technology of Encoding Data in Images of the Earth's Surface / E.G. Zhilyakov, T.N. Balabanova, E.S. Likhogodina, P.G. Lykholob // Nauchnye vedomosti BelGU. Ser. Istorija. Politologija. Jekonomika. Informatika. 2016. №2 (223). Pp. 182-190.
3. Zhilyakov, E.G. Optimal Subband Methods for Analysis and Synthesis of Finite-duration Signals // Automation and remote control. 2015. №4 (76). Pp. 589-602.
4. Zhilyakov E.G., Konstantinov I.S., Chernomorets A.A., Bolgova E.V., The Image Compression Subband Method // International Journal of Soft Computing. 2015. № 10. Pp. 442-447.
5. Fridrich, J. Steganography in digital media: Principles, algorithms, and applications, Steganography in Digital Media. 2012. Pp. 1-441.
6. Kryzhevich L.S., Hahaeva E.V. A method of hiding text information in images using biorthogonal splash transformation/ Auditorium: Electronic scientific journal of the Kursk state University. 2014. № 2. URL: http://auditorium.kursksu.ru/index.php?page=6&new=2 (дата обращения 22.04.2016).
7. Barni M. Capacity of the watermark-channel: how many bits can be hidden within a digital image? / Barni M, Bartolini F, De Rosa A, Piva A // Security and Watermarking of Multimedia Contents, San Jose, Calif, USA, Proceedings of SPIE. January 1999. Pp. 437-448.
8. Christoph Loeffler, Adriaan Lieenberg, and George S. Moschytz. Practical fast 1-d DCT algorithms with 11 multiplications, Acoustics, Speech, and Signal Processing // ICASSP-89, International Conference on IEEE. 1989. pp. 988-991.
9. Edmund Y. Lam, Joseph W. Goodman. A Mathematical Analysis of the DCT Coefficient Distributions for Images // IEEE TRANSACTIONS ON IMAGE PROCESSING, №10 (9). OCTOBER 2000. pp. 1661-1666.
10. Mauro Barni. A DCT-domain system for robust image watermarking / Mauro Barni, Franco Bartolini, Vito Cappellini, Alessandro Piva //Signal Processing 1998, №66 (3). pp. 357-372.
11. Piva A., Barni M., Bartolini F., Cappellini V. DCT-based watermark recovering without resorting to the uncorrupted original image // ICIP '97 Proceedings of the 1997 International Conference on Image Processing (ICIP '97). 1997. №3. p. 520.
12. Stankovic, S., Orovic, I., Sejdic, E. Multimedia signals and systems // Multimedia Signals and Systems. 2012. pp. 1-349.
projections
Frequency domain № of experiment Coefficients Correlation, P
r j
1 2 3 4 5
LF 1 1 1 0.9164
LF 2 1 2 0.9999
MF 3 4 1 0.9989
MF 4 4 5 0.9987
MF 5 4 10 Close to 1
MF 7 4 30 0.9999
MF 6 6 15 Close to 1
MF 8 6 30 Close to 1
MF 9 6 50 0.9999
HF 10 11 50 Close to 1
HF 11 12 60 Close to 1
