International Journal of Biomedical Imaging

International Journal of Biomedical Imaging / 2016 / Article

Research Article | Open Access

Volume 2016 |Article ID 9583727 | https://doi.org/10.1155/2016/9583727

Hui Liang Khor, Siau-Chuin Liew, Jasni Mohd. Zain, "Parallel Digital Watermarking Process on Ultrasound Medical Images in Multicores Environment", International Journal of Biomedical Imaging, vol. 2016, Article ID 9583727, 14 pages, 2016. https://doi.org/10.1155/2016/9583727

Parallel Digital Watermarking Process on Ultrasound Medical Images in Multicores Environment

Academic Editor: Shuihua Wang
Received20 Nov 2015
Revised02 Jan 2016
Accepted04 Jan 2016
Published11 Feb 2016

Abstract

With the advancement of technology in communication network, it facilitated digital medical images transmitted to healthcare professionals via internal network or public network (e.g., Internet), but it also exposes the transmitted digital medical images to the security threats, such as images tampering or inserting false data in the images, which may cause an inaccurate diagnosis and treatment. Medical image distortion is not to be tolerated for diagnosis purposes; thus a digital watermarking on medical image is introduced. So far most of the watermarking research has been done on single frame medical image which is impractical in the real environment. In this paper, a digital watermarking on multiframes medical images is proposed. In order to speed up multiframes watermarking processing time, a parallel watermarking processing on medical images processing by utilizing multicores technology is introduced. An experiment result has shown that elapsed time on parallel watermarking processing is much shorter than sequential watermarking processing.

1. Introduction

The technological advancement in communication network has facilitated healthcare professionals across the world in accessing electronic medical records, such as medical images, and obtaining second opinions for high-quality diagnosis. As a consequence, medical images are exposed to security threat such as tampering of images, which may lead to wrong diagnosis and treatment [1]; thus a digital watermarking on medical image is introduced. So far most of the watermarking research has been done on single frame medical image which is impractical in the real environment, where most of the ultrasound medical images consist of multiframes; thus a digital watermarking on multiframes medical images is proposed. Watermarking could be applied on ultrasound medical images frame by frame sequentially but it would be time consuming for a large dataset; thus a parallel environment is necessary for speeding up multiframes watermarking process by utilizing multicore technology. The performance improvements gained by the use of a multicore processor largely depended on the software algorithms used and their implementation. According to Amdahl’s law, the possible gains are only limited on software portion that can be run in parallel simultaneously on multiple cores. Our work explores efficient parallel implementations of the digital watermarking scheme in multicore environment.

2. Theory: Concept of Digital Watermarking

Digital watermarking is the technology that imperceptibly modifies original data and embeds them directly into image. Digital watermarking in general is comprised of three major components [2]:(1)Watermark generator: a desired watermark(s) is generated for particular applications, which are optionally dependent on some keys.(2)Watermark embedder: watermark(s) are embedded into the object, sometimes based on an embedding key.(3)Watermark detector: detecting the existence of some predefined watermark in the object. It is sometimes desirable to extract a message as well.

The purpose of medical image security is to maintain privacy of the patient information in the image and to assure data integrity that prevents the image from tampering [3]. Watermarking can be used in medical images to prevent unauthorized modification by authenticating the content of the image. Tamper localization capable watermarking scheme can detect and locate modification of pixel values on the image [4]. The tampered area can be recovered by retrieving the original pixel values that were stored on the image itself as a watermark. Tamper localization is useful for deducing the motive of the tampering and whether any modification is legitimate [5]. Liew and Zain proposed a reversible watermarking scheme (TALLOR watermarking scheme) by dividing image into ROI (Region of Interest) and RONI (Region of Noninterest) [6]. ROI is the significant part of the medical images that is used by doctors to diagnose the patient, and RONI is the area outside the ROI. Watermarking for tamper detection and recovery is done in the ROI area based on Jasni’s scheme [7]. The original Least Significant Digits (LSBs) that are removed in watermark embedding process are stored in RONI after compression. The stored LSBs later can be used to restore the image to its original bits value so the watermarking scheme can be reversible [8]. The research was conducted on a single frame of ultrasound which is impractical in the real world; thus a digital watermarking on multiframes medical images is proposed. In order to speed up multiframes watermarking processing time, a parallel computation in multicores watermarking processing on medical images is introduced.

3. Theory: Parallel Computing in Multicores

In recent years there has been a surge of interest in running application in parallel to take advantage of multiprocessor and multicore systems. Developments in microprocessor technologies have resulted in most processors having multiple computing cores in a single chip [9]. Parallel computing is a concept of performing tasks simultaneously by partitioning a large and complex problem into smaller tasks and solving each of them concurrently. There are two forms of parallelism: task parallelism and data parallelism (as shown in Table 1).


Task parallelismData parallelism

It is a form of parallelization of computer code across multiple processors in parallel computing environments [22]. It is a form of parallelization of computing across multiple processors in parallel computing environments [23].

Task parallelism focuses on distributing execution processes (threads) across different parallel computing nodes [23].Data parallelism focuses on distributing the data across different parallel computing nodes.

Data parallelism emphasizes the distributed (parallelized) nature of the data, as opposed to the processing (task parallelism) [10]. Data parallelism is adopted in this experiment since each processor performs the same code (watermarking code) on different pieces of distributed ultrasound frames.

Optimal speedup from parallelization should be linear if the number of processing elements is inversely proportional to its run time. However, not many parallel algorithms achieve optimal speedup. Most of them have a near-linear speedup for small numbers of processing elements, which flattens out into a constant value for large numbers of processing elements [11].

According to Amdahl’s law, the overall speedup from parallelization would be restricted by a small portion of the program which cannot be parallelized. A large and complex program usually consists of several parallelizable parts and several nonparallelizable (sequential) parts. If is the fraction of running time a program spends on nonparallelizable parts [12], then

It is the maximum speedup with parallelization of the program, with being the number of processors used. If the sequential portion of a program accounts for 10% of the runtime ( = 0.1), then a 10x speedup will be the maximum, regardless of how many processors are added. This generates an upper limit on the usefulness of adding more parallel execution units [13].

Gustafson’s law is another law in computing, closely related to Amdahl’s law [9]. It states that the speedup with processors isAmdahl’s law assumes that the total amount of work to be done in parallel is also independent of the number of processors, whereas Gustafson’s law assumes that the total amount of work to be done in parallel varies linearly with the number of processors [12].

Applications are often classified based on the frequency of synchronization and communication needs between their subtasks. Fine-grained parallelism is where an application has a high rate of communication among subtasks; coarse-grained parallelism is where an application does not communicate many times per second, and it is embarrassingly parallel if an application seldom or never has to communicate. Embarrassingly parallel applications are considered the easiest to parallelize [14]. In the embarrassingly parallel problems, speedup factors could be achieved near the number of cores, or even more if the problem is partitioned enough to fit within each core’s cache(s), avoiding use of much slower main system memory.

4. Parallel Computing with MATLAB

4.1. How Parallel Computing Runs a Job

The MATLAB job scheduler (MJS) is the process that coordinates the jobs execution and their tasks evaluation. The MJS can be run on any machine on the network. The MJS runs the submitted jobs in queue order, unless any jobs in its queue are promoted, demoted, cancelled, or deleted. MJS assigns task from the running job to each worker for execution and fetches result from workers upon the task completion. The cycle is repeated with another task. When all tasks for a running job have been assigned to workers, the MJS starts running the next job on the next available worker. Tasks were executed simultaneously by all workers in order to speed up execution of large MATLAB jobs. The MJS then returns the results of all the tasks in the job to the client session (as shown in Figure 1) [15]. In this research, client and MJS are located in a single computer where a client (end user) has sent a job (a stack of ultrasound medical images) into MJS where it segregates the ultrasound medical images into multiple tasks (subset of ultrasound medical images frames) into workers (cores) and executes watermarking process in parallel. For example, in quad cores computer, MJS will segregate 15 frames of ultrasound medical images into 4,4,4,3 frames to each core/worker, respectively. The details process will be discussed in Section 5.

4.2. Life Cycle of a Job

Job progresses through a number of stages upon its creation. In the MJS (or other schedulers), each stage of a job is categorized by their state, such as pending, queued, running, or finished (refer to Table 2). Figure 2 illustrates the stages in the life cycle of a job. Functions used in job management are createJob, submit, and fetchOutputs [14].


Job stageDescription

PendingA job is created on the scheduler with the createJob function in client session of Parallel Computing Toolbox software. The job’s first state is pending. This is when the job was defined by adding tasks to it.

QueuedWhen the submit function is executed on a job, the MJS or scheduler places the job in the queue, and the job’s state is queued. The scheduler executes jobs in the queue in the sequence in which they are submitted, all jobs moving up the queue as the jobs before them are finished. The sequence of the jobs in the queue can be changed with the promote and demote functions.

RunningWhen a job reaches the top of the queue, the scheduler distributes the job’s tasks to worker sessions for evaluation. The job’s state is now running. If more workers are available than are required for a job’s tasks, the scheduler begins executing the next job. In this way, there can be more than one job running at a time.

FinishedWhen all of a job’s tasks have been evaluated, the job is moved to the finished state. At this time, the results can be retrieved from all the tasks in the job with the function fetchOutputs.

FailedWhen using a third-party scheduler, a job might fail if the scheduler encounters an error when attempting to execute its commands or access necessary files.

DeletedWhen a job’s data has been removed from its data location or from the MJS with the delete function, the state of the job in the client is deleted. This state is available only as long as the job object remains in the client.

Table 2 describes each stage in the life cycle of a job [15].

5. Research Methodology: Sequential and Parallel Watermarking Embedding and Authentication Process

In watermarking process on a single frame of ultrasound medical image, watermark is embedded into ultrasound medical images and becomes an input file for watermarking authentication process (as shown in Figure 3). The purpose of authentication process is to localize and recover the tamper region in medical images. In other words, the prerequisite of watermarking authentication process is watermarked ultrasound medical images, in which it is the output file generated by watermarking embedding process. Process watermarking frame by frame sequentially was time consuming; thus parallel computing on multicores was introduced to solve the problem concern. This is accomplished by dividing ultrasounds frames into tasks to each core, respectively, and performing watermarking process simultaneously with others.

Two modes of watermarking process in multiframes will be developed and compared; there is sequential (by using for loop) versus parallel watermarking process (watermarking on multicores); the purpose is to prove that the parallel watermarking will have a significant improvement on elapsed time. In both sequential watermarking embedding and authentication process, watermarking is processed frame by frame sequentially by using a control loop. Therefore the elapsed time is proportional to the number of frames processed. The elapsed time could be reduced by using parallel computing on multicore processing technique and this technique enables ultrasound frames to be divided and distributed to multicore for parallel watermarking processing.

5.1. Sequential Watermarking Process

Sequential watermarking embedding and authentication process are sharing a common framework, where ultrasound medical images in DICOM format are read and perform watermarking process frame by frame sequentially by using a for loop. Processed frames will then concatenate into a variable named as “A” which will convert into DICOM format at the end of the watermarking process. The relationship between both processes is that the output file of watermarking embedding process is the input file for watermarking authentication process. The difference between them is that authentication process has an additional step in identifying the tampered frames; if it is 0, it means nontampered; else it is tampered and then image recovery is performed and tampered frame number is recorded and will be displayed upon the completion of watermarking process (as demonstrated in Figure 4). The main algorithms of parallel watermarking process are dividing volumetric ultrasound medical images and distributing them into a number of cores and executing sequential watermarking processes on each core in parallel; therefore a successful sequential watermarking process is a prerequisite in parallel watermarking process. The details of parallel watermarking process will be discussed in Section 5.2.

5.2. Parallel Watermarking Process

In parallel watermarking process (as illustrated in Figures 5 and 6), ultrasound multiframes medical images were loaded into a quad core microprocessor/cluster and create a job on the scheduler; the job is then divided into tasks according to the number of cores in the microprocessor.

The code implemented enables cluster to autodetect the number of cores available in the processor; if the processor used is a quad core, then the job is divided into 4 tasks, where ultrasound frames are equally divided by 4; for example, if the total number of ultrasound frames is 30, then it will be divided into 8,8,7,7 frames, and if the total number of frames is 15, then it will divided into 4,4,4,3 frames. Those divided frames will then distribute to 4 cores, respectively. In each core, watermarking process is carried out sequentially on the divided frames and at the same time it runs concurrently with other cores (as illustrated in Figure 7).

It is important to ensure that the generated frames output is in order after parallel watermarking process; therefore the frame number has to be assigned and keep track before the frames are sent to the cores, respectively. Upon the tasks completion, tasks were reassembled into a job, whereby all the frames will concatenate into an array and submit back to the cluster. The result is retrieved from all the tasks in the job with the function fetch output. All the frames will concatenate according to the frame number order and write into a DICOM file. A job is deleted on two circumstances:(1)When the scheduler encounters an error.(2)When the job is finished.

Both parallel watermarking embedding and authentication process have a similar process flow as described above, except the process applied on each task, input and output files as listed in Table 3. The input file of authentication process is the tampered output file of embedding process.


Parallel watermarkingInput fileProcess applied on each taskOutput file

Embedding processUltrasound (US) multiframe medical images/raw fileWatermarking embedding processWatermarked US multiframe medical images

Authentication processTampered watermarked US multiframe medical imagesWatermarking authentication processRecovered US medical images and a message of tampered frame numbers

6. Experimental Design and Set-Up

TALLOR watermarking scheme, developed by Liew and Zain [4], will be executed in sequential and parallel modes. The elapsed time obtained from both modes will be compared and speedup factor of parallel relative to sequential watermarking process will be measured. It is to verify the efficiency of parallel framework.

Three important performance metrics were studied. These are(1)imperceptibility: testing the quality of medical images in terms of invisibility of watermarking in multiframes environment;(2)elapsed time: the time taken to perform watermarking embedding and authentication process on medical images in multiframes environment;(3)robustness to tampering: testing the effectiveness and efficiency of the tamper detection, localization, and recovery function in multiframes environment.

The evaluation was performed by running MATLAB (The Mathworks, Inc., Natick, MA, USA) program on a laptop with quad core CPU of Intel® Core i7-3630QM CPU @ 2.4 GHz, 2401 MHz, 4 Core(s), 8 Logical Processor(s), and RAM of 8 GB. Three samples of ultrasound medical images in DICOM format were used to test the system (as shown in Table 4).


Ultrasound medical imagesImage dimension in pixelsBits per pixelNumber of frames

Ultraound_Sample_1.dcm830
Ultraound_Sample_2.dcm815
Ultraound_Sample_3.dcm815

7. Experimental Result

7.1. Imperceptibility

The perceptibility of a watermarked image can be judged according to its fidelity and quality. Fidelity measures the similarity between images before and after watermarking [16]. A high fidelity means that watermarked image is very similar to the original image. The mean-squared-error (MSE) and peak signal-to-noise ratio (PSNR) were calculated by comparing the watermarked image and original image. Watermarked images may bear visible or invisible distortion due to the embedding process. One way to quantify distortion is the mean-square error. If is a vector of predictions, and is the vector of observed values corresponding to the inputs to the function which generated the predictions, then the MSE of the predictor can be estimated by

MSE is the average term by term difference between the original image, , and the watermarked image, . If and are identical, then . A related distortion measure is the peak signal-to-noise ratio (PSNR), measured in decibels (dB). The problem with mean-square error is that it depends strongly on the image intensity scaling while PSNR rectifies this problem by scaling the mean-square error according to the image range [17]. PSNR is defined as follows:where is the peak value of the original image. If the signals are identical, then PSNR is equal to infinity. A high PSNR represents a high fidelity of a watermarked image. In this thesis, PSNR is used as a measurement for image fidelity. A high-quality watermarked image does not have any obvious noticeable distortion caused by the watermark embedding process. The assessment of quality is usually evaluated by human observers and is influenced by personal preferences which are subjective in nature.

Three different sets of ultrasound medical images which contain thirty and fifteen frames have been watermarked in two different ways: (1) sequentially and (2) parallel. It is important to ensure that the quality and fidelity of images were not affected by the way watermarking embedding process performed. Both sequential and parallel watermarking embedding processes have produced the same MSE and PSNR result for each frame as indicated in Table 5 except some negligible differences in the highlighted areas. This means that the operation either in sequential or in parallel mode does not affect the image quality and its fidelity. The PSNR values reflect the medical image integrity and high PSNR values indicate lesser distortion on medical image after watermarking process. PSNR value reflected medical image fidelity; ideally the watermarked medical image should be visually indistinguishable as original image. The PSNR values are calculated for all images ranging within 48.29~48.74 dB, which are within the acceptable range for diagnosis purposes and it has achieved imperceptibility as shown in Figure 9 where the images before and after watermarking embedding process are visually indistinguishable as the original images.

(a) Watermarking embedding process on Ultrasound_Sample_1.dcm

Ultrasound frame numberSequentialParallel
MSEPSNRMSEPSNR

11.044447.9421.044447.942
21.011948.0791.011948.079
31.024148.0271.024148.027
41.020048.0451.020048.045
50.982748.2070.982748.207
61.005548.1071.005548.107
71.024048.0281.024048.028
80.9905818.1720.9905818.172
90.9970448.1440.9970448.144
101.073347.8241.073347.824
111.022448.0351.017848.054
121.021848.0371.021848.037
131.073247.8241.073247.824
141.022448.0351.022448.035
151.021648.0381.021648.038
161.079647.7981.079647.798
171.023748.0291.023748.029
181.005148.1091.005148.109
191.017548.0551.039347.963
200.9734148.2480.9734148.248
210.9913448.1690.9913448.169
221.006748.1021.006748.102
230.9756848.2380.9756848.238
240.9903248.1730.9903248.173
251.000848.1271.000848.127
260.9728348.250.9728348.25
270.9817848.2110.9817848.211
281.002948.1181.002948.118
290.9732348.2490.9732348.249
300.9929948.1610.9929948.161

(b) Watermarking embedding process on Ultrasound_Sample_2.dcm

Ultrasound frame numberSequentialParallel
MSEPSNRMSEPSNR

11.033847.9861.033847.986
21.032947.991.032947.99
31.049147.9231.049147.923
41.050647.9161.050647.916
51.046547.9341.059747.879
61.061647.811.061647.81
71.050247.9181.050247.918
81.045847.9361.045847.936
91.062147.8691.059547.88
101.063247.8651.063247.865
111.061147.8731.061147.873
120.9394648.4020.9394648.402
131.065847.8541.059147.882
141.05547.8981.05547.898
150.9395548.4020.9395548.402

(c) Watermarking embedding process on Ultrasound_Sample_3.dcm

Ultrasound frame numberSequentialParallel
MSEPSNRMSEPSNR

10.8692148.740.8692148.74
20.8771348.70.8771348.7
30.8835848.6680.8835848.668
40.900348.5870.900348.587
50.9027348.5750.9027348.575
60.9196748.4940.9196748.494
70.9119848.5310.9119848.531
80.910848.5370.910848.537
90.9135548.5230.9135548.523
100.9325848.4340.9325848.434
110.954848.3320.954848.332
120.9530748.340.9530748.34
130.9623948.2970.9623948.297
140.9385648.4060.9385648.406
150.9271648.4590.9271648.459

7.2. Elapsed Time

Elapsed time is the time taken to perform watermarking embedding (Figure 8) and authentication process (Figure 12) on medical images in multiframes environment. This section is to test the speedup factor in parallel mode as compared to sequential mode in watermarking process. The formula of speedup factor is defined as follows:

Speedup factor is to measure the speed of parallel mode by factors relative to sequential mode. For example, if the speedup is 3, this means parallel process is three times faster than sequential mode.

7.2.1. Watermarking Embedding Process

Table 6 shows that watermarking embedding process in parallel has achieved a significant speedup (14.13~19.29) relative to sequential process. In Figure 10, the elapsed time in sequential watermarking embedding process on sample_2 and sample_3 is similar but increases by double in sample_1; this means that the elapsed time in sequential process is proportional to the number of frames, whereas in parallel watermarking embedding process, the elapsed time results are consistent despite the number of frames processed. In conclusion, the number of frames does not have much impact on the parallel process as compared to sequential process. The watermarking embedding scheme is pixel oriented, which means different ultrasound sample with the same frame size will produce a similar result.


Input file
(ultrasound medical images)
Number of framesElapsed time in watermarking embedding process (seconds)Speedup Output file
(watermarked ultrasound medical images)
SequentialParallel

Ultraound_Sample_130738.7038.3019.29Watermarked_US1
Ultraound_Sample_215336.4023.8014.13Watermarked_US2
Ultraound_Sample_315332.6022.8014.59Watermarked_US3

The elapsed time is actually depending on many factors such as processor type, clock rate, memory speeds and use of memory catches, location of code in memory, compiler efficiency, and compiler optimization technique.

7.2.2. Watermarking Authentication Process

Different to watermarking embedding process, watermarking authentication process is to verify whether there is any tampering that occurred in the watermarked ultrasound medical images and then recovered the tampered frame to its original state. In Table 7, three watermarked ultrasound medical images have been fully tampered and the speedup factor ranges within 2.92~3.64. The speedup factor is lesser as compared to watermarking embedding process. This is because, firstly, watermarking authentication algorithm has more if-else branches as compared to watermarking embedding algorithm and, secondly, watermarking authentication process needs to return two results (a string of tampered frames and concatenated ultrasound frames), whereas watermarking embedding process just returns one result, that is, concatenated frames of watermarked ultrasound medical images.


Input file (watermarked ultrasound medical images that have been fully tampered) Tampered frame/total frameElapsed time in watermarking authentication process (seconds)Speedup Output file (tampered
ultrasound file recovered as
original ultrasound file)
SequentialParallel

Tampered_Watermarked_US130/30597.1164.23.64Recovered_US1.dcm
Tampered_Watermarked_US215/15270.092.62.92Recovered_US2.dcm
Tampered_Watermarked_US3_a15/15241.273.33.29Recovered_US3_a.dcm
Tampered_Watermarked_US3_b15/15259.171.03.64Recovered_US3_b.dcm

Tampered_Watermarked_US3_a and Tampered_watermarked_US3_b are from the same source but have tampered differently.

In watermarking authentication process also is pixel oriented; therefore the different sources have shown a little impact on elapsed time. It could be observed in Figure 11 that watermarking authentication process in sequential mode is proportional to frame size but it is not the case in parallel mode; the elapsed time does not have many changes in parallel mode at different frame size.

7.2.3. Overall Performance of Watermarking Embedding and Authentication Process

The whole package of watermarking process is involving two steps: (1) watermarking embedding process and (2) watermarking authentication process. Therefore it is necessary to test the overall elapsed time involved in both processes. The high speedup in watermarking embedding process is compromised by the low speedup in watermarking authentication process. Table 8 has shown that the overall speedup factor for three ultrasound medical images samples ranges within 5.21~6.60.

(a) Overall elapsed time taken in watermarking process on Ultrasound_Sample_1.dcm

Watermarking processSequentialParallelSpeedup

Embedding738.7038.3019.29
Authentication597.10164.203.64
Overall time taken1335.80202.506.60

(b) Overall elapsed time taken in watermarking process on Ultrasound_Sample_2.dcm

Watermarking processSequentialParallelSpeedup

Embedding336.4023.8014.13
Authentication270.0092.602.92
Overall time taken606.40116.405.21

(c) Overall elapsed time taken in watermarking process on Ultrasound_Sample_3.dcm

Watermarking processSequentialParallelSpeedup

Embedding332.6022.8014.59
Authentication241.273.33.29
Overall time taken573.896.15.97

7.3. Robustness to Tampering

In order to demonstrate the tamper localization function in detecting forgery, counterfeited images were created by manually modifying the pixel values in the watermarked images using image processing software—ImageJ 1.46r. Figure 13 shows an example of tampering on three frames (frames numbers 2, 4, and 6) in ultrasound watermarked medical images.

Different set of tampered watermarked ultrasound medical images has been used to test the effectiveness and efficiency of the tamper detection, localization, and recovery function in multiframes environment (as shown in Table 9). The function’s effectiveness was measured by checking whether it could detect the tampered frame number in multiframes environment (as shown in Figures 12 and 13) and is able to recover to its original form (as shown in Figure 14). The function’s efficiency was measured by comparing the elapsed time taken for both sequential and parallel watermarking authentication process (as shown in Figure 15). Both effectiveness and efficiency testing were performed while testing the function’s robustness to tampering.


Tampered watermarked US imagesElapsed time for watermarking authentication process
SeqDRParDR

Tampered_0_frames18.2YesYes23.1YesYes
Tampered_3_frames56.9YesYes35.3YesYes
Tampered_6_frames97.0YesYes39.1YesYes
Tampered_9_frames139.9YesYes53.8YesYes
Tampered_12_frames241.6YesYes58.2YesYes
Tampered_15_frames259.1YesYes71.0YesYes

Seq: sequential version (seconds).
Par: parallel version (seconds).
D: able to detect and display the tampered frame number?
R: able to recover to its original form after authentication process?

For tampered free watermarked ultrasound medical images, the initial setup of parallel authentication process such as determining number of cores and dividing and distributing frames to each core was time consuming; therefore it takes longer processing time than sequential version.

For both sequential and parallel watermarking authentication process, the elapsed time is proportional to the total number of tampered frames. Relative to sequential watermarking authentication, the efficiency of parallel watermarking authentication showed remarkably when there are more tampered frames. In summary, the elapsed time for parallel watermarking authentication process is varied based on the total number of tampered frames in ultrasound medical images.

With 100% tampered watermarked ultrasound medical images, the efficiency of parallel version can be calculated based on elapsed time, such as 259.1/71 = 3.65, which means parallel version can perform 3.65 times faster than sequential version. In other words, sequential watermarking authentication has performed a job all by itself, whereas parallel watermarking authentication has delegated a job to four workers/cores, in which parallel version has performed approximately 4 times faster than sequential version. In conclusion, the speedup factor was approximate to the number of cores since there is little communication between subtasks.

The frame order of tampered frame also is a factor that affects the performance in parallel process (as shown in Figure 16). The black rectangle boxes with bold font represent a tampered frame. Both (a) and (b) have distributed 15 frames (containing 9 tampered frames) of watermarked ultrasound medical images to 4 cores, respectively. The main difference between (a) and (b) was that the tampered frames were organized in different order. The elapsed time for (a) is 5 seconds lesser than (b). It is because (a) has an even distribution of tampered frames and hence fairer workloads among 4 cores as compared to (b).

8. Conclusion and Future Work

The approach used to overcome the performance constraint in sequential watermarking process is by distributing ultrasound frames over multiple cores. The performance constraint for sequential watermarking process can be categorized into two problems: capacity and capability. The capacity problem occurred when the existing hardware and software are unable to perform the anticipated computations in an estimated time [18]. For example, it may not be feasible to conduct watermarking process on a large data size of DICOM files in any reasonable manner. Therefore, to run the entire computation may become impractical even though the existing hardware and software are capable of performing the watermarking process as required. The actual physical constraint, such as the processor speeds or total memory on a system, may cause the problem of capability, in which it will restrict the amount of watermarking processes performed. System upgrades may be a solution for the problem concern, but it is bounded by technology and cost constraint. In this case, the sequential watermarking problem may be partitioned into smaller and manageable parts that can be performed in parallel. Parallel computing with MATLAB is the simplest approach to leveraging multicores processor. However, the maximum number of parallel threads cannot exceed the number of cores available on the system. The performance gain obtained by using multiple cores on a single system is also limited and varies depending on the specific computation and the data size [9]. The proposed parallel watermarking scheme required little effort to separate the ultrasound medical images into a number of parallel tasks due to the little dependency (or communication) between those parallel tasks and achieved a speedup factor that is almost equivalent to the number of cores; for example, quad cores microprocessor will result in a speedup factor of approximately four, which means parallel watermarking process is approximately four times faster than sequential mode in quad cores microprocessor. Hence it could be classified as “embarrassingly parallel problems,” a phrase that comments on the ease of parallelizing such applications and the fact that it would be embarrassing for the programmer or compiler to not take advantage of such an obvious opportunity to improve performance [19]. The performance of proposed parallel digital watermarking scheme could be further improved by using graphics processing unit (GPU). Both CPU and GPU could run thousands of threads concurrently, but GPU will have a better performance than CPU due to its larger number of cores possessed relative to CPU; the difference in theoretical performance can differ by a factor ten in favour of the GPU; therefore a parallel digital watermarking operation on GPU is recommended. In many cases a hybrid CPU-GPU implementation yields the best performance. A good example is image registration algorithms, where the GPU can be used to calculate a chosen similarity measure in parallel, while the CPU can run a serial optimization algorithm. So et al. made a comparison between CPUs and GPUs for ultrasound systems, in terms of power efficiency and cost effectiveness. The conclusion was that a hybrid CPU-GPU system performed best [20]. The adoption of either GPU or hybrid CPU-GPU is largely dependent on the parallel adaption in an algorithm; an algorithm that exhibited “embarrassingly parallel problem” will be suitably used in GPU whereas hybrid CPU-GPU is suitably applied in an algorithm that exhibited “fine-grained parallelism.” For future work, the proposed method could be applied on Magnetic Resonance Images (MRI) where the ROI could be classified by using weighted-type fractional Fourier transform approach [21] prior to watermarking process. Since the watermarking is pixel oriented, thus it could be also applied on nature images with lesser restriction on image fidelity requirement as compared to medical images.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

  1. C. K. Tan, J. C. Ng, X. Xu, C. L. Poh, Y. L. Guan, and K. Sheah, “Security protection of DICOM medical images using dual-layer reversible watermarking with tamper detection capability,” Journal of Digital Imaging, vol. 24, no. 3, pp. 528–540, 2011. View at: Publisher Site | Google Scholar
  2. Q. Li and N. Memon, “Security models of digital watermarking,” in Multimedia Content Analysis and Mining: International Workshop, MCAM 2007, Weihai, China, June 30-July 1, 2007. Proceedings, vol. 4577 of Lecture Notes in Computer Science, pp. 60–64, Springer, Berlin, Germany, 2007. View at: Publisher Site | Google Scholar
  3. G. Wang and N.-N. Rao, “A fragile watermarking scheme for medical image,” in Proceedings of the 27th Annual International Conference of the Engineering in Medicine and Biology Society (IEEE-EMBS '05), pp. 3406–3409, IEEE, Shanghai, China, January 2005. View at: Publisher Site | Google Scholar
  4. S.-C. Liew and J. M. Zain, “The usage of block average intensity in tamper localization for image watermarking,” in Proceedings of the 4th International Congress on Image and Signal Processing (CISP '11), pp. 1044–1048, Shanghai, China, October 2011. View at: Publisher Site | Google Scholar
  5. S.-C. Liew, S.-W. Liew, and J. M. Zain, “Tamper localization and lossless recovery watermarking scheme with ROI segmentation and multilevel authentication,” Journal of Digital Imaging, vol. 26, no. 2, pp. 316–325, 2013. View at: Publisher Site | Google Scholar
  6. S.-C. Liew and J. M. Zain, “Reversible medical image watermarking for tamper detection and recovery,” in Proceedings of the 3rd IEEE International Conference on Computer Science and Information Technology (ICCSIT '10), pp. 417–420, Chengdu, China, July 2010. View at: Publisher Site | Google Scholar
  7. J. M. Zain and A. R. M. Fauzi, “Medical image watermarking with tamper detection and recovery,” in Proceedings of the 28th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '06), pp. 3270–3273, New York, NY, USA, September 2006. View at: Publisher Site | Google Scholar
  8. B. W. R. Agung and P. F. P. Adiwijaya, “Medical image watermarking with tamper detection and recovery using reversible watermarking with LSB modification and run length encoding (RLE) compression,” in Proceedings of the IEEE International Conference on Communication, Networks and Satellite (ComNetSat '12), pp. 167–171, Bali, Indonesia, July 2012. View at: Publisher Site | Google Scholar
  9. S. Samsi, V. Gadepally, and A. Krishnamurthy, “MATLAB for signal processing on multiprocessors and multicores,” IEEE Signal Processing Magazine, vol. 27, no. 2, pp. 40–49, 2010. View at: Publisher Site | Google Scholar
  10. S. P. Wang and R. S. Ledley, “Advanced computer architecture,” in Computer Architecture and Security Fundamentals of Designing Secure Computer Systems, John Wiley & Sons, New York, NY, USA, 2013. View at: Google Scholar
  11. J. Degroote and J. Vierendeels, “Multi-solver algorithms for the partitioned simulation of fluid-structure interaction,” Computer Methods in Applied Mechanics and Engineering, vol. 200, no. 25–28, pp. 2195–2210, 2011. View at: Publisher Site | Google Scholar
  12. R. Cieszewski, M. Linczuk, K. Pozniak, and R. Romaniuk, “Review of parallel computing methods and tools for FPGA technology,” in Photonics Applications in Astronomy, Communications, Industry, and High-Energy Physics Experiments 2013, vol. 8903 of Proceedings of SPIE, Wilga, Poland, October 2013. View at: Publisher Site | Google Scholar
  13. F. P. Brooks, The Mythical Man-Month: Essays on Software Engineering, Addison-Wesley, Reading, Mass, USA, 1996.
  14. P. Vicat-Blanc, S. Soudan, R. Guillier, and B. Goglin, “Utilization of network computing technologies,” in Computing Networks: From Cluster to Cloud Computing, chapter 2, John Wiley & Sons, Hoboken, NJ, USA, 2013. View at: Publisher Site | Google Scholar
  15. The MathWorks, How Parallel Computing Products Run a Job, 2015, http://www.mathworks.com/help/distcomp/how-parallel-computing-products-run-a-job.html.
  16. I. J. Cox, M. L. Miller, and J. A. Bloom, Digital Watermarking, Morgan Kaufmann, San Francisco, Calif, USA, 2002.
  17. Y.-D. Zhang, S.-H. Wang, X.-J. Yang et al., “Pathological brain detection in MRI scanning by wavelet packet Tsallis entropy and fuzzy support vector machine,” SpringerPlus, vol. 4, article 716, 2015. View at: Publisher Site | Google Scholar
  18. B. V. Krishna and K. Baskaran, “Parallel computing for efficient time-frequency feature extraction of power quality disturbances,” IET Signal Processing, vol. 7, no. 4, pp. 312–326, 2013. View at: Publisher Site | Google Scholar
  19. A. Leykin, J. Verschelde, and Y. Zhuang, “Parallel homotopy algorithms to solve polynomial systems,” in Mathematical Software—ICMS 2006: Second International Congress on Mathematical Software, Castro Urdiales, Spain, September 1–3, 2006. Proceedings, vol. 4151 of Lecture Notes in Computer Science, pp. 225–234, Springer, Berlin, Germany, 2006. View at: Publisher Site | Google Scholar
  20. H. K.-H. So, J. Chen, B. Y. S. Yiu, and A. C. H. Yu, “Medical ultrasound imaging: to GPU or not to GPU?” IEEE Micro, vol. 31, no. 5, pp. 54–65, 2011. View at: Publisher Site | Google Scholar
  21. Y.-D. Zhang, S. Chen, S.-H. Wang, J.-F. Yang, and P. Phillips, “Magnetic resonance brain image classification based on weighted-type fractional Fourier transform and nonparallel support vector machine,” International Journal of Imaging Systems and Technology, vol. 25, no. 4, pp. 317–327, 2015. View at: Publisher Site | Google Scholar
  22. R. A. Jain and D. V. Padole, “Scalable and flexible heterogeneous multi-core system,” International Journal of Advanced Computer Science and Applications, vol. 3, no. 12, pp. 174–179, 2012. View at: Publisher Site | Google Scholar
  23. G. Luo and D. Zhang, “Efficiency improvement for data-processing of partial discharge signals using parallel computing,” in Proceedings of the 10th IEEE International Conference on Solid Dielectrics (ICSD '10), pp. 1–4, IEEE, Potsdam, Germany, July 2010. View at: Publisher Site | Google Scholar

Copyright © 2016 Hui Liang Khor et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views1242
Downloads570
Citations

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.