بهینه‌سازی فرایند ترکیب تصاویر با استفاده از الگوریتم جستجوی گرانشی

شیخان, منصور; چمن کار, آماج

تعداد نشریات	44
تعداد شماره‌ها	1,877
تعداد مقالات	15,278
تعداد مشاهده مقاله	43,732,220
تعداد دریافت فایل اصل مقاله	17,564,744

بهینه‌سازی فرایند ترکیب تصاویر با استفاده از الگوریتم جستجوی گرانشی

هوش محاسباتی در مهندسی برق

مقاله 5، دوره 6، شماره 3، آبان 1394، صفحه 0-67 اصل مقاله (1.35 M)

نوع مقاله: مقاله پژوهشی فارسی

نویسندگان

منصور شیخان^* ¹؛ آماج چمن کار²

¹- دانشیار، دانشکده مهندسی برق- دانشگاه آزاد اسلامی- واحد تهران جنوب- ایران

²کارشناس ‏ارشد، گروه مهندسی مکاترونیک، دانشگاه آزاد اسلامی- واحد تهران جنوب – ایران

چکیده

از آنجا که به دلیل محدودیتهای عکسبرداری نمیتوان با
گرفتن تنها یک عکس از یک صحنه به تصویری با کیفیت و وضوح مطلوب دست یافت، ترکیب
تصاویر با استفاده از روشهای گوناگون روز به روز در حال گسترش و بهبود است. از طرفی
امروزه از الگوریتمهای جمعیتی مبتنی بر تصادف جهت بهینهیابی استفاده گستردهای میشود
که اغلب با الهام از فرایندهای فیزیکی یا رفتارهای موجودات به وجود آمدهاند. الگوریتم
بهینهیابی جستجوی گرانشی (GSA) از مجموعه
الگوریتمهای بهینهیابی است که با الهام از قانون جاذبه و مفهوم جرم شکل یافته است
و عاملهای جستجوگر مجموعهای از اجرام میباشند. این پژوهش قصد دارد به ترکیب
تصاویر با میزان تمرکزهای مختلف با استفاده از الگوریتم بهینهیابی جستجوی گرانشی
بپردازد. برای این کار ابتدا ترکیب تصاویر با استفاده از اطلاعات فرکانس مکانی
مورد مطالعه قرار گرفته و سپس با بررسی مفهومGSA و تلفیق آن با فرکانس مکانی، روشی
جدید برای ترکیب تصاویر ارائه شده است. عملکرد روش پیشنهادی با روش ترکیب مبتنی بر
اطلاعات پیکسلها بر اساس منطق فازی و روش بهینهسازی فرایند ترکیب تصاویر با بکارگیری
الگوریتم بهینهیابی اجتماع ذرات (PSO)
مورد مقایسه قرار گرفته که نتایج تجربی حاکی از میزان خطای کمتر روش پیشنهادی است.
درضمن، کارآمدی این روش برای ترکیب تصاویر رنگی نیز بررسی شده است.

کلیدواژه‌ها

ترکیب تصاویر؛ فرکانس مکانی؛ بهینه‌یابی؛ الگوریتم جستجوی گرانشی

اصل مقاله

1- مقدمه[1]

ترکیب تصاویر[1] عبارت است از ادغام ویژگی‌های دو یا چند تصویر، به طوری که تصویر حاصل، ویژگی‏های مهم هر یک از تصاویر را شامل شود و تصویری با کیفیت و ویژگی‏های کامل‏تر ارایه نماید [1]. ترکیب تصاویر کاربردهای گسترده‏ایی دارد و با در دسترس بودن اطلاعات حسگرها در زمینه‌های مختلف، مانند عکس‏برداری دیجیتال با میزان تمرکز خودکار [2]، سنجش از دور [3]، تصویربرداری پزشکی و بینایی ماشین [4]، ترکیب حسگری به یک موضوع جدید پژوهشی تبدیل شده است.در عکس‏برداری، هنگامی که عدسی دوربین روی یک نقطه خاص متمرکز است، وضوح تصویر در آن نقطه بسیار خوب است، اما فضاهای دورتر از نقطه تمرکز وضوح کمتری خواهند داشت [1].تمرکز چندگانه[2] روشی قدرتمند برای ترکیب تصاویر با وضوح‏های متفاوت است که در آن‏ها تصاویر با میزان تمرکزهای مختلف گرفته شده و سپس، با یکدیگر ادغام می‏شوند تا تصویری با وضوح کامل حاصل ‌شود. برای ادغام تصاویر با وضوح‏های متفاوت، روش‏های مختلفی وجود دارد که این روش‏ها را می‏توان به سه گروه اصلی دسته‏بندی کرد. نخستین دسته، روش‏هایی هستند که بر پایه اطلاعات پیکسل‏ها می‏باشند و از جمله آن‏ها می‏توان به روش‏های بیشینه‏گیری [5 و 6]، میانگین‏گیری [6 و 7] و به‏کارگیری منطق فازی [8] اشاره کرد. دو روش نخست از ابتدایی‌ترین روش‌های ترکیب تصاویر هستند که از ویژگی‏های آن‌ها می‌توان به سرعت بالا و سادگی الگوریتم اشاره کرد، اما این روش‌ها اصولا دقت پایینی دارند. در روش فازی کارایی بالاست، اما به‏علت سرعت پایین نمی‌توان از آن در سامانه‏های بی‏درنگ استفاده کرد. دسته دوم، روش‏های مبتنی بر توابع تبدیل می‏باشند که بیشتر برگرفته از تبدیل فوریه هستند، مانند ترکیب تصاویر با استفاده از تبدیل لاپلاس [9- 11]، تبدیل موجک [12- 16]، IHS[3] [17] و تبدیل Brovey [18]. این روش‌ها دقت مناسبی دارند، برای انجام کارهای تلفیقی به‏کار گرفته می‌شوند و تبدیل موجک از پرکاربردترین آن‌ها محسوب می‌شود. از آنجا که در IHS از ویژگی‏های رنگ تصاویر استفاده می‏شود، برای تصاویر رنگی مناسب می‏باشند ولی در عمل به‏علت کارایی پایین‏تر روش‏های IHS و Brovey نسبت به سایر روش‏ها، از آن‏ها کمتر استفاده می‏شود. دسته سوم از روش‏ها نیز مبتنی بر استفاده از ترکیب اطلاعات بلوکی و روش‏های پردازش و بهینه‏یابی هوشمند هستند.

شایان ذکر است که در سالیان اخیر، الگوریتم‏های مبتنی بر رفتار موجودات یا پدیده‏های طبیعی در حوزه‏های کاربردی گوناگون با هدف بهینه‏سازی استفاده شده‏اند، مانند ژنتیک [19]، اجتماع زنبورها [20]، بهینه‏سازی مبتنی بر اجتماع ذرات[4] (PSO) [21]، خفاش [22]، کاوش باکتری [23]، و چرخه آب [24]. همچنین، نسخه‏های اصلاحی و آمیختار این الگوریتم‏ها نیز بسیار مورد توجه قرار گرفته‏اند، مانند اجتماع ذرات وفقی [25]، اجتماع ذرات هم‏تکاملی اصلاح‏شده [26]، آمیختار الگوریتم‏های جغرافیای زیستی و تکاملی تفاضلی [27]، و الگوریتم توسعه‏یافته بهینه‏یابی جفت‏گیری زنبورعسل [28].

از الگوریتم‏های بهینه‏یابی فرا ابتکاری در موضوع ترکیب تصاویر نیز با اهداف مختلف استفاده شده است که نمونه‏هایی از این اهداف را می‏توان چنین برشمرد:

الف) بهینه‏یابی ابعاد بلوک‏های تصویر [29]،

ب) بهینه‏سازی ضرایب وزن ترکیب [30 و 31]،

پ) آشکارسازی لبه و بهسازی در سامانه‏های ترکیب تصویر [32]،

ت) بهینه‏سازی ساختار و شاخص‏های شبکه عصبی مصنوعی مورد استفاده برای ترکیب تصاویر [33 و 34].

از نمونه الگوریتم‏ها/ مدل‏های هوش محاسباتی مورد استفاده در موضوع ترکیب تصاویر نیز می‏توان به این موارد اشاره کرد: ژنتیک [35 و 36]، PSO [37- 39]، تکامل تفاضلی [29]، یادگیری- تعلیم[5] [30]، تکاملی چندهدفی [31]، اجتماع مورچگان [32]، اجتماع زنبورها [33]، طوفان مغزی آشوبی[6] [34] و شبکه عصبی [40 و 41].

در این مورد الگوریتم PSO و نسخه‏های اصلاحی آن بیشتر استفاده شده‏اند. در این پژوهش، با استفاده از فرکانس مکانی و بهینه‏سازی اندازه بلوک‏ها مبتنی بر الگوریتم جستجوی گرانشی[7] (GSA)، کوشش شده است روش کارآمدی برای ترکیب تصاویر ارایه شود. لازم به یادآوری است که GSA به‏عنوان الگوریتمی نوین‏تر و نیز با سه مزیت یاد شده در زیر نسبت به الگوریتم PSO مطرح است [42]:

الف) روال به ‏روز رسانی در الگوریتم PSO بدون در نظر گرفتن کیفیت پاسخ‏ها انجام می‏شود (مقادیر برازش در این روال مهم نیستند)، حال آنکه در GSA نیرو متناسب با مقدار برازش است (عامل‏ها فضای جستجوی پیرامون خود را با توجه به نیرو مشاهده می‏کنند).

ب) روال به‏روزرسانی در الگوریتم PSO بدون درنظر گرفتن فاصله بین پاسخ‏ها انجام می‏شود، حال آنکه در GSA نیرو تناسب معکوس با فاصله بین پاسخ‏ها دارد.

پ) در الگوریتم PSO، راستا و جهت یک عامل تنها با توجه به دو موقعیت بهینه محاسبه می‏شود، حال آنکه در GSA این جهت بر اساس نیروی کلی حاصل از تمام عامل‏های دیگر محاسبه می‏شود.

نکته قابل توجه این است که در هر دو الگوریتم PSO و GSA، بهینه‏سازی از طریق حرکت عامل‏ها در فضای جستجو انجام می‏شود، هر چند که روش حرکت در این دو الگوریتم متفاوت است.

البته در این راستا، GSA با هدف افزایش دقت سامانه‏های بازیابی تصویر در مطالعات سالیان اخیر به کار گرفته شده است [43]. همچنین، الگوریتم‏های دیگر مانند بهینه‏سازی مبتنی بر اجتماع گربه‏ها نیز برتری خود را بر PSO نشان داده است [44 و 45]؛ که به تازگی در حوزه مشابه با مقاله حاضر به کار گرفته شده است [46].

ساختار این مقاله چنین است: در بخش دوم الگوریتم جستجوی گرانشی مرور شده و در انتهای این بخش نسخه‏های اصلاحی این الگوریتم نیز معرفی خواهند شد. بخش سوم به موضوع ترکیب تصاویر بر اساس اطلاعات فرکانس مکانی می‏پردازد. در بخش چهارم با تعریف مسأله در قالب یک مسأله بهینه‏یابی و استفاده از فرکانس مکانی به بیان روش پیشنهادی پرداخته شده است. در بخش پنجم نیز نتایج تجربی و مقایسه با سایر روش‏ها ارایه شده‏اند. در بخش ششم نتیجه‏گیری بیان شده است.

2- الگوریتم جستجوی گرانشی

این الگوریتم بر اساس قوانین نیوتن ارایه شده است که در آن عامل‌های جستجوکننده اجرام هستند که به شکل سیاره‏های یک منظومه تصور می‏شوند و منطقه بهینه همچون یک سیاه‏چاله، سیاره‏ها را به سمت خود می‏کشد. تبادل اطلاعات و اثرگذاری اجرام روی یکدیگر تحت نیروی گرانش انجام می‏شود [42].

هر جسمی جسم دیگر را به سمت خود جذب می‏کند و مقدار نیروی جاذبه بین دو جسم با جرم‏های M₁ و M₂ و فاصله R، با حاصل‏ضرب جرم آن دو جسم و عکس مجذور فاصله بین آن‏ها متناسب است رابطه (1):

(1)

در این رابطه G ثابت گرانش نامیده می شود. این رابطه نشان می‏دهد که هر جسم به‏واسطه نیروی جاذبه، محل و جرم سایر اجسام را درک می‏کند و هر جسم به نسبت میزان جرمش و فاصله‏ای که با دیگر اجسام دارد، روی سایر اجسام تأثیر می‏گذارد و به آن‏ها نیرو وارد می‏کند.

بر اساس قانون اول نیوتن، هر جسم حالت سکون یا حرکت یکنواخت خود را بر روی خط راست حفظ می‏کند مگر این‏که تحت تأثیر نیرو یا نیروهایی مجبور به تغییر آن حالت شود. بر اساس قانون دوم نیوتن، وقتی به جسمی نیرویی وارد می‏شود شتابی می‏گیرد که به نیرو و جرم جسم بستگی دارد. رابطه بین شتاب، نیرو و جرم در رابطه (2) ارایه شده است که در آن شتاب با a نشان داده شده است.

(2)

بر اساس قانون گرانش، به هر جسم از یک مجموعه، از جانب سایر اجسام نیروهای گرانشی وارد می‏شود. در نتیجه، جسم به سمت برآیند این نیروها که با F_r نشان داده می‏شود، شتاب می‏گیرد. روابط بالا بیان می‏کنند که هر جسم، جسم دیگر را به سمت خود می‏کشاند اما تأثیر جسم بزرگتر و نزدیکتر بیشتر است. با درنظر گرفتن قانون جاذبه و قوانین حرکت، میزان و جهت حرکت هر جسم، توافقی است بین تأثیر نیروی ثقل وارد بر آن و سرعت فعلی جسم.

در یک سامانه ایزوله با دو جسم i و j، جسم i تحت تأثیر نیروی جاذبه جسم j شتابی برابر a_i می‏گیرد که بر اساس رابطه (3) محاسبه می‏شود. F_ij، مقدار نیروی گرانشی وارد بر جسم i از جانب جسم j است که مطابق رابطه (3) محاسبه می‏شود. در این روابط M_pi و M_ii به‏ترتیب جرم گرانشی غیرفعال و جرم اینرسی جسم i و M_ajجرم گرانشی فعال جسم j هستند.

(3)

در فیزیک ثابت شده است که ضریب گرانشی با آهنگ بسیار کندی در طول زمان، بر اساس رابطه (4) کوچک می‏شود.

(4)

در الگوریتم جستجوی گرانشی، بهینه‏یابی به‏کمک طرح قوانین گرانشی و حرکت در یک سامانه مصنوعی به‏ شکل گسسته در زمان انجام می‏شود. محیط سامانه همان محدوده تعریف مسأله است.

در قدم اول، فضای سامانه مشخص می‏شود. محیط شامل یک دستگاه مختصات چند- بعدی در فضای تعریف مسأله است. هر نقطه از فضا، یک پاسخ مسأله است. عامل‏های جستجوکننده مجموعه‏ای از اجرام هستند. هر جرم چهار مشخصه دارد: موقعیت جرم، جرم گرانشی فعال، جرم گرانشی غیرفعال، و جرم اینرسی [42]. موقعیت جرم، نقطه‏ای در فضا است که پاسخی از مسأله است. مقدار اجرام گرانشی و اینرسی، با توجه به برازندگی هر جرم تعیین می‏شوند.

پس از تشکیل سامانه، قوانین حاکم بر آن مشخص می‏شوند. فرض می‏شود تنها قانون گرانش و قوانین حرکت حاکم هستند. حال سامانه را به شکل مجموعه‏ای از N جرم تصور کنید. موقعیت بعد d از جرم i با x_i^d نشان داده شده است.

(5)

در این سامانه، در زمان t به جرم i از سوی جرم j در جهت بعد d نیرویی به اندازه F_ij^d(t) وارد می‏شود. مقدار این نیرو بر اساس رابطه (6) محاسبه می‏شود. M_aj و M_pi به‏ترتیب جرم گرانشی فعال جرم j و جرم گرانشی غیرفعال جرم i می‏باشند، G(t) ثابت گرانش در زمان t و R_ij فاصله بین دو جرم i و j می‏باشند. ε نیز یک عدد بسیار کوچک است.

(6)

نیروی وارد بر جرم i در جهت بعد d در زمان t (F_i^d(t))، بر اساس رابطه (7) برابر مجموع ضریب‏های تصادفی نیروهایی است که K جرم برتر بر جرم وارد می‏کنند. در این رابطه rand_jیک عدد تصادفی با توزیع یکنواخت در بازه [0,1] است. در رابطه (7) می‏توان مجموع تمام نیروهای وارد بر جسم را درنظر گرفت، اما برای بهبود قدرت کشف الگوریتم، تنها به مجموعه Kbest شامل K عضو برتر، اجازه تأثیرگذاری بر سایر اعضا داده می‏شود. مقدار K به شکل متغیر با زمان تعریف می‏شود. به این شکل که در زمان شروع تمام اجرام روی یکدیگر تأثیر می‏گذارند و با گذشت زمان از تعداد اعضای تأثیرگذار بر جمعیت، به شکل خطی کم می‏شود تا این‏که در انتها، تنها درصد کمی از بهترین‏های جمعیت بر سایر اعضا نیرو وارد می‏کنند.

(7)

شتاب جرم i در جهت بعد d در زمان t با a_i^d(t) و جرم اینرسی جسم i با M_ii نشان داده شده است.

(8)

سرعت بعدی هر جرم برابر مجموع ضریبی از سرعت فعلی جرم و شتاب جرم تعریف می شود (رابطه 9). موقعیت جدید بعد d از جرم i بر اساس رابطه (10) محاسبه می‏شود. v_i^d(t) سرعت بعد d جرم i در زمان t است.

(9)

(10)

rand_i و rand_j اعداد تصادفی با توزیع یکنواخت در بازه [0,1] هستند که برای حفظ ویژگی تصادفی‏ بودن جستجو، استفاده شده‏اند.

برای تنظیم ثابت گرانش، از یک مقدار اولیه شروع کرده، با گذشت زمان مقدار آن کاهش داده می‏شود. ثابت گرانش بر اساس رابطه (11)، تابعی از ثابت گرانش اولیه و زمان است. این موضوع در دنیای واقعی نیز صدق می‏کند و ثابت گرانش با آهنگ بسیار کندی در طول زمان کوچک می‏شود. در نسخه پیوسته این الگوریتم، یک پیشنهاد برای این تابع، استفاده از رابطه نمایی برای کاهش ثابت گرانش است (رابطه 12).

(11)

(12)

در رابطه (12)، G₀ ثابت گرانش اولیه، α یک مقدار ثابت مثبت و T کل تکرارهای الگوریتم و به‏عبارتی طول عمر سامانه است.

در این الگوریتم، اجرام گرانشی و اینرسی بر اساس رابطه (13)، برابر در نظر گرفته شده، برای تنظیم آن‏ها، از مقدار تابع هدف اجرام بر اساس رابطه (14) استفاده می‏شود. مقدار اجرام بر اساس رابطه (15)، هنجار می‏شوند. در این روابط، به اجرام با شایستگی بهتر، جرم بیشتری نسبت داده شود.

(13)

(14)

(15)

در این روابط fit_i(t) بیانگر میزان برازندگی جرم i در زمان t است. در مسائل کمینه‏یابی می‏توان از روابط (16) و (17) برای محاسبه بهترین و بدترین مقدار شایستگی استفاده کرد. در مسائل بیشینه‏یابی نیز بهترین و بدترین مقدار شایستگی بر اساس روابط (18) و (19) تعریف می شوند.

(16)

(17)

(18)

(19)

در ابتدای تشکیل سامانه، هر جسم به شکل تصادفی در یک نقطه از فضا قرار می‏گیرد که پاسخی از مسأله است. در هر لحظه از زمان، اجرام ارزیابی شده، تغییر مکان هر جرم محاسبه شده، در زمان بعدی، جرم در آن موقعیت قرار می‏گیرد. جرم‏های گرانشی، جرم اینرسی و ثابت گرانش نیوتن در هر مرحله به‏روزرسانی می‏شوند. شرط توقف می‏تواند پس از تکرارهای مشخص، تعیین شود.

البته نسخه‏های اصلاحی GSA نیز در سالیان اخیر پیشنهاد شده‏اند که توضیح نمونه‏هایی از آن‏ها در ادامه آورده شده است:

الف) کاهش پیچیدگی محاسباتی GSA استاندارد از طریق کاهش تعداد ارزیابی‏های تابع هدف توسط الگوریتمی به نام GSA- خوشه‏بندی‏شده (مبتنی بر خوشه‏بندی جمعیت و جایگزینی هر خوشه با جرم مرکزی به‏عنوان یک عامل جدید) [47]

ب) به‏کارگیری نظریه مکانیک کوانتومی در GSA استاندارد به منظور جلوگیری از همگرایی زودرس به کمینه محلی [48]

پ) اعمال یک عملگر جستجوی محلی گسسته به بهترین پاسخ GSA استاندارد با یک احتمال متغیر [49]

ت) به‏کارگیری عملگر فروپاشی[8] در GSA استاندارد (با اعمال به بردار موقعیت عامل‏ها) با هدف بهبود توانایی کاوش فضای جستجو [50]

ث) به‏کارگیری ترکیب روش جستجو در الگوریتم PSO و توپ- کشسان[9] برای تسریع همگرایی و نیز استفاده از عملگر جهش آشوبی برای فرار از کمینه‏های محلی [51]

ج) به‏کارگیری یادگیری مبتنی بر تقابل[10] برای آمایش جمعیت و جهش نسل در GSA استاندارد [52]

چ) جلوگیری از همگرایی زودرس و گرفتاری در کمینه‏های محلی با به‏کارگیری عملگر آشوبی [53 و 54] (به‏عنوان نمونه با تغییر در معادله به‏روزرسانی سرعت در GSA استاندارد [55])

ح) جایگزینی معادله نمایی در رابطه تابع گرانش با یک تابع تکه‏ای- خطی به‏منظور بهبود توانایی جستجو در GSA استاندارد [56].

3- ترکیب تصاویر بر اساس اطلاعات فرکانس مکانی

در این قسمت به مرور الگوریتمی پایه برای ترکیب تصاویر بر اساس فرکانس مکانی پرداخته می‏شود. این الگوریتم از نظر محاسباتی ساده بوده و می‏تواند در کاربردهای بی‏درنگ استفاده شود [57]. یادآوری می‏شود که فرکانس مکانی برای اندازه‏گیری سطح فعالیت در یک تصویر استفاده می‏شود. فرکانس مکانی تصویر با استفاده از معادلات ارایه‏شده در رابطه (20) قابل محاسبه است:

(20)

در این رابطه F تصویر نهایی و n×m مبین ابعاد تصویر است. مقدار زیاد فرکانس مکانی سطح فعالیت بالای تصویر را توصیف می‏کند که در نتیجه شفافیت تصویر است.

گام‏های اصلی در این الگوریتم عبارتند از:

1- تجزیه تصاویر اصلی به بلوک‏هایی با ابعاد n×m.

2- محاسبه فرکانس مکانی برای هر بلوک.

3- مقایسه فرکانس مکانی مربوط به دو بلوک متناظر A_iو B_i و تشکیل iامین بلوک (F_i) از تصویر ترکیبی به‏ شکل رابطه (21):

(21)

در این رابطهTh در واقع مقدار آستانه برای ایجاد محدوده‌ای برای میانگین‌گیری بین دو تصویر است. این مقدار خود می‌تواند با استفاده از سایر الگوریتم‌های بهینه‌یابی اصلاح شود، اما با توجه به اندازه تصاویر، مقداری بین 1 تا 2 برای آن مقدار مناسبی است [57]. منظور از نیز میانگین‌گیری پیکسل به پیکسل است.

4- بررسی و اصلاح نتایج ترکیب در مرحله 3. در این مرحله، هدف اصلاح مقادیر در محل‌های اتصال بلوک‌ها به یکدیگر است. انجام این مرحله برای اطمینان از عملکرد مناسب الگوریتم است و در بسیاری موارد این مرحله مد نظر قرار داده نمی‌شود.

4- روش پیشنهادی

با توجه به آنچه گفته شد، در این پژوهش، تلاش شده است با کمک فرکانس مکانی و بهره‌گیری از الگوریتم جستجوی گرانشی، روش کارآمدی برای ترکیب تصاویر با میزان تمرکز‌های مختلف ارایه شود. برای این کار در ابتدا تصاویر اصلی به بلوک‌های مختلف قطعه‌بندی شده و سپس، فرکانس مکانی برای هر قطعه را محاسبه می‌شود. در بلوک‌های متناظر، بلوکی که دارای مقدار فرکانس مکانی بیشتری است، به‏عنوان خروجی محاسبه می‌شود. این قسمت همانند چیزی است که در روش ترکیب تصاویر با استفاده از فرکانس مکانی مطرح شده است. آنچه در این جا دارای اهمیت است، اندازه بلوک‌هاست. در ترکیب تصاویر با وضوح مختلف با استفاده از الگوریتم‌های بهینه‌سازی، به‏دنبال پیدا کردن اندازه بهینه برای قطعه‌بندی تصاویر هستیم. بنابراین، با یک مسأله بهینه‏سازی مواجه‌ایم. در این مورد، GSA می‌تواند به شکل خودکار و تطبیقی، بهترین مقدار را برای اندازه بلوک‌های تصویر پیدا کند. در این روش، فرکانس مکانی برای تعیین تصویر با وضوح بیشتر استفاده می‏شود. شکل (1) بیان‏کننده چگونگی عملکرد این روش است.

شکل (1): چگونگی ترکیب تصاویر با میزان تمرکز‌های مختلف با استفاده از الگوریتم GSA

4-1- تولید جمعیت اولیه

برای شروع باید مسأله را در قالب GSA مدل کرد. نخستین گام تشکیل یک جمعیت اولیه برای شروع کار است. برای این منظور یک جمعیت تصادفی از اجرام ایجاد کرده که هر جرم نشان‏دهنده اندازه بلوک در فضاست. فضای کلی جستجو اندازه هر کدام از تصاویر است. هر جرم دارای دو بعد می‌باشد که بیانگر طول و عرض بلوک است. بنابراین، در رابطه (22) برای موقعیت هر جرم می‏توان نوشت:

(22)

این مقادیر برای جمعیت اولیه به‏ شکل تصادفی ایجاد می‌شوند و به مرور در تکرارهای بعدی مقادیر آن‏ها تغییر کرده تا مقدار بهینه به دست آید. در این روش هر عامل، کل تصویر را بلوک‌بندی می‏کند، به عبارتی در هر تکرار، تعداد i تصویر بلوک‌بندی‏شده بررسی می‏شوند. مقدار تعلق‏گرفته به جمعیت اولیه خود می‌تواند به‌صورتی بهینه تولید شود که در ادامه به آن پرداخته می‌شود. آنچه در این مرحله باید مورد توجه قرار گیرد این است که اندازه آن‏ها نباید از فضای جستجوی تصویر که همان اندازه تصویر اصلی است، بیشتر باشد. در الگوریتم جستجوی گرانشی این مسأله درنظر گرفته می‌شود، اما برای جمعیت اولیه باید این محدودیت اعمال شود تا مقادیر تصادفی انتخاب‏شده از اندازه تصویر بزرگ‌تر نشوند.

4-2- تعیین تابع مناسب برای بهینه‌یابی

انتخاب تابع صحیح برای تعیین میزان وضوح هر تصویر و در حقیقت تعیین تابعی که باید برای رسیدن به تصویر مطلوب، کمینه/بیشینه شود بسیار مهم است و به گونه‌ای طرح مسأله بهینه‌یابی را مشخص می‌کند. با توجه به اینکه میزان وضوح تصویر قرار است بهینه شود، بنابراین، تابع مسأله یا همان تابع برازش[11]، تابع فرکانس مکانی است.

تابع فرکانس مکانی در روش پیشنهادی در دو مرحله بررسی می‏شود. در مرحله اول، برای محاسبه مقدار وضوح هر بلوک در تصویر به‏کار گرفته می‏شود که این مرحله مربوط به تابع بهینه‏یابی نیست و تنها نقش یک انتخاب‏گر بلوک را ایفا می‌کند. مرحله دوم، محاسبه میزان وضوح تصویر نهایی برای هر کدام از جمعیت‌هاست که در اینجا، نقش همان تابع برازش در مسأله بهینه‏یابی را دارد.

بنابراین، تصاویر ورودی با توجه به ابعاد و اندازه جرم‏ها تقسیم‌بندی می‌شوند. سپس، وضوح هر کدام از بلوک‌های هر تصویر محاسبه می‏شود و بلوک با وضوح بالاتر به‏عنوان بلوک متناظر تصویر ترکیبی انتخاب می‌شود و این انتخاب بر اساس تباین[12] تصویر بر روی قسمت‌های مختلف تکرار شده و تصاویر عکس‌های ترکیب‏شده در مرحله اول برای مقادیر اولیه ایجاد می‌شوند. پس از آنکه تصویر ایجاد شد، مقدار برازش هر تصویر با استفاده از فرکانس مکانی محاسبه می‌شود، به این ترتیب که مقدار فرکانس مکانی تصویر نهایی محاسبه شده و مقدار آن به‏عنوان مقدار برازش آن تصویر یا جرم مربوط درنظر گرفته می‌شود. در واقع می‏توان با بیانی هر یک از تصاویر را یک جرم درنظر گرفت و با بیان دقیق‌تر هر بلوک انتخابی را یک جرم درنظر گرفت. بهترین مقدار برازش در بین تصاویر ایجادشده به‏عنوان مقدار بهترین سراسری[13] انتخاب می‏شود. با توجه به اینکه در این مسأله بهینه‌یابی به‏دنبال پیدا کردن فرکانس مکانی با مقدار بیشتر هستیم، از روابط بیشینه‌یابی استفاده می‌کنیم (روابط 18 و 19).

این الگوریتم می‌تواند به دو شکل متوقف شود: با توجه به میزان تکرارهای تعیین‏شده، و یا با توجه به رسیدن به میزان خطای قابل قبول و تعریف‏شده. روند نمای این روش در شکل (2) نشان داده شده است.

شایان ذکر است که برای تصاویر رنگی، در ابتدا تصاویر را خاکستری کرده و سپس، تصاویر خاکستری را بررسی کرده و در پایان، بلوک رنگی متناظر با تصویر با کیفیت بالاتر جایگزین می‏شود. به این ترتیب تصویر ادغامی رنگی حاصل می‌شود.

4-3- گام‏های الگوریتم پیشنهادی

1- دریافت تصاویر و تبدیل آنها به تصاویر خاکستری.

2- آمایش مقادیر جرم‏ها و سرعت‏ها با مقادیر اولیه تصادفی.

3- بلوک‏بندی تصاویر ورودی و محاسبه فرکانس مکانی هر بلوک.

4- انتخاب بلوک با فرکانس مکانی بیشتر به‏عنوان خروجی.

5- ایجاد تصویر نهایی برای تمام جرم‏ها و محاسبه فرکانس مکانی و انتخاب فرکانس مکانی بیشتر به عنوان fbest.

6- انتخاب fbest به‏عنوان fbest global در صورت بیشتر بودن نسبت به تکرار قبل.

شکل (2): روندنمای روش پیشنهادی ترکیب تصاویر با استفاده از الگوریتم جستجوی گرانشی

7- محاسبه سرعت وموقعیت جرم‏ها.

8- محاسبه G، نیرو و شتاب.

9- محاسبه سرعت وموقعیت جدید جرم‏ها.

10- مقداردهی مجدد اجرام خارج‏شده از فضای جستجو.

11- تکرار مراحل قبل تا رسیدن به مقدار خطای مطلوب و یا پایان تعداد تکرارهای تعیین‏شده.

12- انتخاب تصویر دارای مقدار fbest global به‏عنوان تصویر ترکیبی نهایی خروجی.

4-4- مقداردهی اولیه

مقادیر fbest و fbest global اولیه، صفر درنظر گرفته می‏شوند. تعداد جمعیت اولیه پیشنهادی 50 و بیش‏ترین تعداد تکرار 5 بار است که این تعداد تکرار بر اساس بهینه‏بودن زمان و مناسب بودن خطای نهایی تعیین شده است. این مقدار با توجه به اندازه تصویر قابل تغییر است.

با توجه به اینکه در سامانه دو مقدار x و y (طول و عرض هر بلوک) مطرح است، در الگوریتم بُعد برابر 2 می‏باشد (یعنی d=2). باید توجه داشت که این مقادیر باید در بازه فضای جستجو که همان اندازه تصویر می‌باشد، تعریف شوند.

برای محاسبه مقدار G در رابطه 12، مقدار 20=α و 100=G₀ در نظر گرفته شده است.

از آنجا که اندازه بلوک‌ها مقادیر صحیح هستند، پس از ایجاد جمعیت اولیه، ابتدا مقادیر آن‌ها تبدیل به عدد صحیح می‌شود. از آنجا که کل تصویر با اندازه یکسانی بلوک‏بندی می‌شود، باقیمانده حاصل تقسیم طول و عرض تصویر بر مقدار ورودی محاسبه می‌شود و درصورتی‏که برابر صفر نباشد یک واحد از آن کم می‌شود و این کار تا جایی‏که باقیمانده برابر صفر شود، ادامه پیدا می‌کند و به این ترتیب مقادیر جمعیت‌های ایجادشده بر روی تصویر منطبق شده و الگوریتم برای مسأله مورد بحث اجرا می‌شود. برای روشن‌تر شدن موضوع، فرض ‌کنیم که طول و عرض تصویر با ابعاد 325×235 ‌باشد و یکی از مقادیر جمعیت بدین‏ترتیب بدست آمده باشد: (7/13 2/48). بنابراین، ابتدا این مقدار گرد شده و به‏ شکل (14 48) در محاسبات درنظر گرفته می‏شود. سپس، تقسیم ابعاد تصویر به مقادیر فوق انجام شده و چون باقیمانده صفر نیست، یک واحد از آن کم می‌شود و این کار تا جایی‏که باقیمانده صفر شود، ادامه پیدا می‏کند. در این مورد مقادیر انتخابی (13 47) به باقیمانده صفر منجر می‏شوند و بنابراین، تصویر با این اندازه بلوک‏بندی می‌شود.

همان‏طور که پیش‏تر نیز اشاره شد، مقداردهی اولیه اندازه بلوک‌ها نیز می‌تواند بهینه شود. بدیهی است هر چه اندازه بلوک‌های تصویر را کوچک‏تر فرض نماییم نتیجه بهتر خواهد بود، اما آنچه مهم است اندازه بهینه بلوک‌هاست. ممکن است حتی تصویری با دو بلوک نتیجه‌ای کاملا مطلوب را نیز داشته باشد و نیازی به بلوک‌بندی بیش از اندازه نباشد. بلوک‌های بیشتر به معنای عملیات بیشتر و در نتیجه زمان بیشتر عملیات است و هدف از استفاده از الگوریتم‌های بهینه‏یابی در این موارد نیز دقیقا مقابله با همین وضعیت است.

با توجه به بررسی‌های انجام‏شده، مشاهده می‌شود که بلوک‌ها در پایان همگی دارای اندازه‌ا‌ی کمتر از نصف اندازه تصاویر اصلی هستند.

5- نتایج تجربی و مقایسه با سایر روش‏ها

برای بررسی بیشتر موضوع، عملکرد روش پیشنهادی با دو روش دیگر مورد مقایسه قرار گرفته است. این روش‏ها عبارتند از: بهینه‏سازی ‏فرآیند ترکیب تصاویر به روش PSO و روش ترکیب مبتنی بر اطلاعات پیکسل‏ها بر اساس منطق فازی. روش کار در الگوریتم PSO مشابه الگوریتم GSA است، با این تفاوت که مقادیر بهینه اندازه بلوک‌ها به جای الگوریتم GSA، توسط الگوریتم PSO تعیین می‌شوند.

شاخص‏های طراحی در PSO به شکل زیر می‌باشند: ضریب‌های c₁ و c₂ که برای تعیین رفتار جمعی و شناختی است، مقدار 2 در نظر گرفته شده‌اند. مشاهده می‌شود که هنگامی که این مقادیر بزرگتر از 2 باشند، سرعت ذره با مقدار بزرگتری به‏روز می‌شود و در نتیجه ذره پرش زیادی در فضای جستجو خواهد داشت و ممکن است از محدوده جستجو خارج شود. مقدار وزن اولیه w، 1 در نظر گرفته شده است که ضریب w_damp=0.99 می‏باشد [58].

روش فازی پیشنهادی توسط نویسندگان در [8] نیز به‏عنوان روشی مبتنی بر اطلاعات پیکسل‌ها برای مقایسه انتخاب و نتایج آن بررسی شده است. در این روش، پیکسل‌های متناظر دو تصویر با یکدیگر مقایسه شده و بر اساس منطق فازی مقدار خروجی مناسب برای پیکسل انتخاب می‌شود. قوانین فازی پیشنهادی بر این اساس ایجاد شده‌اند که پیکسل‌های با مقادیر دورتر از مقدار 128 وزن بیشتری به خود بگیرند. نحوه عملکرد این روش در شکل (3) نشان داده شده است. سامانه ممدانی فازی دارای دو ورودی (یکی برای تصویر با وضوح سمت راست و دیگری برای تصویر با وضوح سمت چپ) و یک خروجی (برای تصویر خروجی) است. ورودی‏ها و خروجی دارای 7 تابع عضویت از مقدار 0 تا 255 و به‏ شکل مثلثی هستند و تنها در قسمت میانی خروجی از یک تابع ذوزنقه‏ای استفاده شده است تا اطراف مقدار 128 نوعی میانگین‏گیری انجام گیرد. در شکل (4) نیز نتیجه ترکیب تصاویر با سطوح خاکستری برای این سه روش نشان داده شده است. مقادیر بهینه اندازه بلوک‌‌ها برای تصویر شکل (4) (با ابعاد 240×240 پیکسل) با استفاده از الگوریتم PSO، بلوک‌هایی با اندازه‌های 7×7 بدست آمد. همچنین، در روش GSA، بهترین کیفیت تصویر نهایی برای تصویر شکل (4)، بلوک‌هایی با اندازه‌های 60×8 محاسبه شد که مقادیر مربوط به فرکانس مکانی آن‏ها و همچنین، میزان خطا با استفاده از معیار RMSE نیز در جدول (1) ارایه شده است.

شکل (3): روش مبتنی بر منطق فازی برای ترکیب تصاویر با میزان تمرکز‌های مختلف [8]

شکل (4): نتیجه ترکیب تصاویر به سه روش مختلف، a) تصویر با وضوح سمت چپ، b) تصویر با وضوح سمت راست، c) تصویر واضح اصلی، d) نتیجه ترکیب به روش فازی، e) نتیجه ترکیب با به‏کارگیری الگوریتم PSO، f) نتیجه ترکیب با به‏کارگیری GSA

مقادیر RMSE نیز با توجه به رابطه 23 محاسبه شده‏اند:

(23)

که در این رابطه، f و به‏ترتیب مبین تصاویر اصلی و ترکیبی هستند.

در شکل (5) مقادیر بیشینه، کمینه و میانگین تابع هدف در هر تکرار برای تصویر ارایه‏شده در شکل (4) در صورت استفاده از الگوریتم جستجوی گرانشی نشان داده شده است.

جدول (1): مقادیر RMSE تصاویر ترکیبی با سطوح خاکستری برای سه روش مورد بررسی در این پژوهش

معیار	تصویر با وضوح سمت چپ	تصویر با وضوح سمت راست	تصویر ترکیبی به روش فازی	تصویر ترکیبی به روش PSO	تصویر ترکیبی به روش GSA
فرکانس مکانی	3862/9	2540/9	7360/10	7503/10	5964/11
میزان خطا (RMSE)	1685/3	3028/3	961/0	3295/0	2670/0

شکل(5): مقادیر بیشینه، کمینه و میانگین تابع هدف (فرکانس مکانی) در تکرارهای مختلف GSA

آنچه در ترکیب تصاویر تا این مرحله انجام شده است، کار بر روی تصاویر سیاه و سفید بوده است. برای کار بر روی تصاویر رنگی پیشنهادهای گوناگونی داده شده است که در آن‏ها معمولا از ابتدا بر روی تصویر رنگی کار کرده و تصاویر رنگی را با یکدیگر ترکیب می‏کنند [59- 61]. از آنجا که در روش پیشنهادی در این مقاله، بر روی بلوک‌های تصویر کار می‌شود، می‌توان ابتدا تصاویر را به سطح خاکستری انتقال داده و پس از تعیین اینکه بلوک مربوط به کدام تصویر وضوح بیشتری دارد، بلوک مربوط در تصویر رنگی را جایگزین کرد و به این ترتیب ترکیب تصاویر در سطح رنگی انجام می‌شود. به این ترتیب علاوه بر اینکه پیچیدگی کار کمتر می‏شود، تصویر حاصل دارای کیفیت بالاست. این روش، روشی انعطاف‏پذیر است که می‏تواند با سایر روش‏ها نیز تلفیق شود.

در ادامه نمونه‏ای از ترکیب تصاویر رنگی با استفاده از روش پیشنهادی ارایه شده است (شکل 6). برای مشاهده بهتر، قسمتی از نتایج حاصل بزرگ‏نمایی شده و در شکل (7) نشان داده شده است. بخش بزرگ‏نمایی‏شده، مرز بین مات‏شدگی دو تصویر می‏باشد. با دقت در این تصاویر می‌توان مشاهده کرد که ترکیب تصاویر به‏خوبی انجام شده است. البته برای بررسی دقیق‏تر، مجموع مربعات خطا (MSE) نیز در جدول (2) نشان داده شده است. نتایج شهودی خطا نیز در شکل (8) نشان داده شده است. با توجه به جدول (2) و همچنین، با دقت در تصویر قسمت c (تفاضل تصویر اصلی با تصویر ترکیبی)، مشاهده می‏شود که الگوریتم پیشنهادی به‏خوبی دو قسمت راست و چپ تصاویر را با یکدیگر ترکیب کرده است و میزان خطای آن با تصویر اصلی بسیار کم است.

در نهایت، برای بررسی عملکرد روش پیشنهادی، نتایج تجربی دیگری از ترکیب تصاویر با روش‏های مختلف در شکل‏های (9) تا (11) نشان داده شده است. برای مقایسه بهتر، مقادیر خطای تصویر نهایی نسبت به تصویر اصلی نیز در جدول (3) آورده شده است.

شکل(6): ترکیب تصاویر رنگی با استفاده از روش پیشنهادی، a) تصویر با وضوح سمت چپ، b) تصویر با وضوح سمت راست، c) نتیجه ترکیب با به‏کارگیری GSA، d) تصویر اصلی

شکل(7): بزرگ‏نمایی قسمتی از شکل 6 در مرز مات‏شدگی، a) تصویر با وضوح سمت چپ، b) تصویر با وضوح سمت راست، c) نتیجه ترکیب با به‏کارگیری GSA، d) تصویر اصلی

جدول (2): مقادیر MSE برای تصاویر رنگی با روش ترکیب تصاویر پیشنهادی

تصویر ترکیبی به روش GSA	تصویر با وضوح سمت چپ	تصویر با وضوح سمت راست	معیار
5739/0	8945/14	7792/9	MSE

شکل(8): تفاضل تصاویر شکل (6) با تصویر اصلی، a) تفاضل تصویر با وضوح سمت چپ و تصویر اصلی، b) تفاضل تصویر با وضوح سمت راست و تصویر اصلی، c) تفاضل تصویر ترکیبی به روش GSA و تصویر اصلی، d) تفاضل تصویر اصلی و خودش

شکل (9): نتیجه ترکیب تصاویر با روش‏های مختلف، a) تصویر با وضوح سمت چپ، b) تصویر با وضوح سمت راست، c) تصویر واضح اصلی، d) نتیجه ترکیب به روش فازی، e) نتیجه ترکیب با به‏کارگیری الگوریتم PSO، f) نتیجه ترکیب با به‏کارگیری GSA (ابعاد تصویر:325×235، اندازۀ بهینه بلوک 13×47 و fbest=55.5541)

شکل (10): نتیجه ترکیب تصاویر با روش‏های مختلف، a) تصویر با وضوح سمت چپ، b) تصویر با وضوح سمت راست، c) تصویر واضح اصلی، d) نتیجه ترکیب به روش فازی، e) نتیجه ترکیب با به‏کارگیری الگوریتم PSO، f) نتیجه ترکیب با به‏کارگیری GSA (ابعاد تصویر:240×240، اندازۀ بهینه بلوک 120×240 و fbest=13.7351)

شکل (11): نتیجه ترکیب تصاویر با روش‏های مختلف، a) تصویر با وضوح سمت چپ، b) تصویر با وضوح سمت راست، c) تصویر واضح اصلی، d) نتیجه ترکیب به روش فازی، e) نتیجه ترکیب با به‏کارگیری الگوریتم PSO، f) نتیجه ترکیب با به‏کارگیری GSA (ابعاد تصویر:472×631، اندازۀ بهینه بلوک 59×1 و fbest=14.7551)

جدول (3): مقادیر RMSE در روش‏های مختلف ترکیب تصاویر در شکل‏های (9- 11)

شماره شکل	تصویر ترکیبی به روش فازی	تصویر ترکیبی به روش PSO	تصویر ترکیبی به روش GSA
9	9528/0	1118/0	0301/0
10	3707/1	8699/0	2603/0
11	0730/2	8430/0	2287/0

6- نتیجه‏گیری

تاکنون روش‏های مختلفی برای ترکیب تصاویر با هدف رسیدن به تصویری با کیفیت مطلوب ارایه شده است. در این مقاله، روش ترکیب مبتنی بر اطلاعات بلوکی تصاویر و تلفیق آن با الگوریتم جستجوی گرانشی (GSA) استفاده شد. در این راستا برای ترکیب تصاویر با وضوح مختلف به‏دنبال پیدا کردن اندازه بهینه بلوک‏ها برای قطعه‏بندی تصاویر هستیم، و در همین جاست که GSA وظیفه بهینه‏یابی را انجام می‏دهد. بدین‏منظور، تابع فرکانس مکانی به عنوان تابع برازش درنظر گرفته شد.

به منظور ارزیابی روش پیشنهادی، نتایج ترکیب تصاویر به روش GSA با دو روش دیگر مقایسه شد. این روش‏ها عبارت بودند از: بهینه‏سازی ‏فرآیند ترکیب تصاویر به روش PSO و روش ترکیب مبتنی بر اطلاعات پیکسل‏ها بر اساس منطق فازی. نتایج تجربی گویای میزان خطای کمتر روش پیشنهادی در مقایسه با دو روش دیگر است. در ضمن، روش پیشنهادی در ترکیب تصاویر رنگی نیز همچون تصاویر سیاه و سفید به شکل موفقی عمل کرد که نتایج یاد شده در شکل‏های (6) تا (8) نشان داده شد.

برای ادامه پژوهش می‏توان با استفاده از نظریه گراف‌ها به بهبود روش ترکیب تصاویر با به‏کارگیری الگوریتم جستجوی گرانشی پرداخت و از این طریق هر تصویر به بلوک‌هایی با ابعاد متفاوت بلوک‌بندی شده و اندازه بلوک‌ها در یک تصویر نیز بهینه شود. از مشکلاتی که در این راستا پیش‌روی پژوهش‏ها خواهد بود می‌توان به وابستگی شدید تعداد جمعیت جرم‌ها (در GSA) به اندازه تصاویر ورودی اشاره کرد که ممکن است در صورت بزرگ بودن اندازه تصویر، تعداد جمعیت زیادی نیاز شده و باعث پایین آمدن سرعت سامانه شود. همچنین، از آنجا که قرار است یک تصویر با اندازه‌های مختلف بلوک‌بندی شود، اندازه بلوک‌ها به ویژه بلوک‌های نهایی در طول و عرض تصویر بسیار به یکدیگر وابسته بوده و انتخاب محل و اندازه هر بلوک به سادگی امکان‏پذیر نخواهد بود که امید است با استفاده از مفاهیم گراف‌ها این مشکلات از پیش رو برداشته و روشی کاراتر ارایه شود.

[1] تاریخ ارسال مقاله : 3/07/1393

تاریخ پذیرش مقاله : 21/06/1394

نام نویسنده‏ی مسئول : منصور شیخان

نشانی نویسنده‏ مسئول : ایران – تهران – مجتمع فنی و مهندسی دانشگاه آزاد اسلامی (واحد تهران جنوب) – دانشکده‏ برق

1- Image fusion

2- Multi-focus

3- Intensity-Hue-Saturation

4- Particle Swarm Optimization

5- Teaching-learning

6- Chaotic brain storm optimization

7- Gravitational Search Algorithm

8- Disruption

9- Elastic-ball

10- Opposition

11- Fitness

12- Contrast

13- Global

مراجع

[1] Zhang, X., Li, X., Liu, Z., Feng, Y., “Multi-focus image fusion using image-partition-based focus detection”, Signal Processing, Vol. 102, pp. 64-76, 2014.

[2] Saeedi, J., Faez, K., “Fisher classifier and fuzzy logic based multi-focus image fusion”, Proc. IEEE International Conference on Intelligent Computing and Intelligent Systems, Vol. 4, pp. 420-425, Shanghai, Nov. 2009.

[3] Yang, X. H., Jing, Z. L., Liu, G., Hua, L. Z., Ma, D. W., “Fusion of multispectral and panchromatic images using fuzzy rule”, Communications in Nonlinear Science and Numerical Simulation, Vol. 12, No. 7, pp. 1334-1350, 2007.

[4] Garg, S., Kiran, U., Mohan, K., Tiwary, R., “Multilevel medical image fusion using segmented image by level set evolution with region”, Proc. 27th Annual International Conference of Engineering in Medicine and Biology Society, pp. 7680- 7683, Shanghai, Jan. 2006.

[5] Wang, Z., Cui, P., Li, F., Chang, E., Yang, S., “A data-driven study of image feature extraction and fusion”, Information Sciences, Vol. 281, pp. 536-558, Oct. 2014.

[6] Wu, T., Wu, X-J., Luo, X-Q., “A study on fusion of different resolution images”, Procedia Engineering, Vol. 29, pp. 3980-3985, 2012.

[7] Shah, S. K., Shah, D. U., “Comparative study of image fusion techniques based on spatial and transform domain”,International Journal of Innovative Research in Science, Engineering and Technology, Vol. 3, No. 3, pp. 10168-10175, 2014.

[8] Chamankar, A., Sheikhan, M., Razaghian, F., “Multi-focus image fusion using fuzzy logic”, Proc. 13th Iranian Conference on Fuzzy Systems, pp. 1-4, Qazvin, August 2013.

[9] James, A. P., Dasarathy, B. V., “Medical image fusion: A survey of the state of the art”, Information Fusion, Vol. 19, pp. 4-19, Sep. 2014.

[10] Burt, P. J., Adelson, E. H., “The Laplacian pyramid as a compact image code”, IEEE Transactions on Communications, Vol. 3l, No. 4, pp. 532-540, 1983.

[11] Khaleghi, B., Khamis, A., Karray, F. O., Razavi, S. N., “Multisensor data fusion: A review of the state-of-the-art”, Information Fusion, Vol. 14, No. 1, pp. 28-44, 2013.

[12] Tian, J., Chen, L., “Multi-focus image fusion using wavelet-domain statistics”, Proc. 17th International Conference on Image Processing, pp. 1205-1208, Hong Kong, Sep. 2010.

[13] De, I., Chanda, B., “A simple and efficient algorithm for multifocus image fusion using morphological wavelets”, Signal Processing, Vol. 86, No. 5, pp. 924-936, 2006.

[14] Pajares, G., de la Cruz, J. M., “A wavelet-based image fusion tutorial”, Pattern Recognition, Vol. 37, No. 9, pp. 1855-1872, 2004.

[15] Stankovic, R. S., Falkowski, B. J., “The Haar wavelet transform: Its status and achievements”, Computers and Electrical Engineering, Vol. 29, No. 1, pp. 25-44, 2003.

[16] Naidu, V. P. S., Raol, J. R., “Pixel-level image fusion using wavelets and principal component analysis”, Defence Science Journal, Vol. 58, pp. 338-352, 2008.

[17] Hong, G., Zhang, Y., “High resolution image fusion based on wavelet and IHS transformations”, Proc. IEEE/ISPRS Joint Workshop on Remote Sensing and Data Fusion over Urban Areas, pp. 99-104, May 2003.

[18] Nikolakopoulos, K. G., “Comparison of nine fusion techniques for very high resolution data”, Photogrammetric Engineering & Remote Sensing, Vol. 74, No. 5, pp. 647- 659, 2008.

[19] Seyedbarzegar, S.M., Mirzaie, M., “Optimization of spacers placement in metal oxide surge arresters due to the thermal stability increment using genetic algorithm”, Computational Intelligence in Electrical Engineering, Vol. 5, No. 4, pp. 63-70, 2015.

[20] Azarinejadian, F., Mirhosseini Moghaddam, S.M., Marzband, M., Parhizi, N., “Optimal energy management system in grid connected microgrid integrated with distributed generation by using the multi-period artificial bee colony”, Computational Intelligence in Electrical Engineering, Vol. 5, No. 3, pp. 1-14, 2014.

[21] Pourtaheri, Z., Zahiri, S.H., “Design of an integrated phase frequency detector with optimal power consumption and delay by using particle swarm optimization algorithm”, Computational Intelligence in Electrical Engineering, Vol. 5, No. 3, pp. 15-22, 2014.

[22] Razavi Asfali, S.E., Falaghi, H., “Dynamic multi-stage placement of phasor measurement units using bat optimization algorithm”, Computational Intelligence in Electrical Engineering, Vol. 5, No. 2, pp. 91-110, 2014.

[23] Jourabian, M., Afzalan, E., “A hybrid bacterial foraging optimization algorithm and sequential quadratic programming method for dynamic economic dispatch considering the valve-point effects”, Computational Intelligence in Electrical Engineering, Vol. 5, No. 1, pp. 13-26, 2014.

[24] Yazdani-Asrami,M., Khosravi, A.R., A. Ghabeli Juybari, A., “Parameter estimation of induction motors using water cycle optimization”, Computational Intelligence in Electrical Engineering, Vol. 4, No. 3, pp. 71-82, 2013.

[25] Rezaee, Y., Khodabakhshian, A., Moazzami, M., “Optimal UPFC placement in Isfahan-Khozestan transmission network by using APSO algorithm for improving transient stability”,Computational Intelligence in Electrical Engineering, Vol. 3, No. 2, pp. 55-72, 2012.

[26] Abdoli, M., Hatami, A., “Dynamic economic emission dispatch incorporating wind farms using modified co-evolutionary particle swarm optimization meta-heuristic algorithm”, Computational Intelligence in Electrical Engineering, Vol. 5, No. 4, pp. 31-44, 2015.

[27] Abdi Daviran, S., Teshnehlab, M., Aliyari Shouredeli, M., Golahmadi, H., “Multi objective optimization using biogeography based optimization and differentional evolution algorithm”, Computational Intelligence in Electrical Engineering, Vol. 3, No. 3, pp. 11-24, 2012.

[28] Baradaran tavakoli, H., Mozafari, B., “Short term hydrothermal scheduling via improved honey-bee mating optimization algorithm”, Computational Intelligence in Electrical Engineering, Vol. 3, No. 3, pp. 67-80, 2012.

[29] Aslantas, V., Kurban, R., “Fusion of multi-focus images using differential evolution algorithm”, Expert Systems with Applications Vol. 37, No. 12, pp. 8861-8870, 2010.

[30] Jin, H., Wang, Y., “A fusion method for visible and infrared images based on contrast pyramid with teaching learning based optimization”, Infrared Physics & Technology, Vol. 64, pp. 134-142, May 2014.

[31] Jin, H., Xi, Q., Wang, Y., Hei, X., “Fusion of visible and infrared images using multiobjective evolutionary algorithm based on decomposition”, Infrared Physics & Technology, Vol. 71, pp. 151-158, Jul. 2015.

[32] Kavitha, C. T., Chellamuthu, C., “Medical image fusion based on hybrid intelligence”, Applied Soft Computing, Vol. 20, pp. 83-94, Jul. 2014.

[33] Yu, J., Duan, H., “Artificial bee colony approach to information granulation-based fuzzy radial basis function neural networks for image fusion”, Optik, Vol. 124, No. 17, pp. 3103-3111, 2013.

[34] Li, C., Duan, H., “Information granulation-based fuzzy RBFNN for image fusion based on chaotic brain storm optimization”, Optik, Vol. 126, No. 15-16, pp. 1400-1406, 2015.

[35] Zhang, J., Feng, X., Song, B., Li, M., Lu, Y., “Multi-focus image fusion using quality assessment of spatial domain and genetic algorithm”, Proc. IEEE Conference on Human System Interactions, pp. 71-75, Krakow, May 2008.

[36] Taher, G. M., Wahed, M. E., Taweal, G. E., Fouad, A., “Image fusion approach with noise reduction using genetic algorithm”, International Journal of Advanced Computer Science and Applications, Vol. 4, No. 11, pp. 10-16, 2013.

[37] Zhang, L., Zhang, X., Xu, X., Wang, D., Liu, J., Liu, Y., “Quantum-inspired particle swarm optimization algorithm with performance evaluation of fused images”, Optica Applicata, Vol. 43, No. 4, pp. 679-691, 2013.

[38] Siddiqui, A. B., Jaffar, M. A., Hussain, A., Mirz, A. M., “Block-based pixel level multi-focus image fusion using particle swarm optimization”, International Journal of Innovative Computing, Information and Control, Vol. 7, No. 7(A), pp. 3583-3596, 2011.

[39] Zhang, X., Sun, L., Han, J., Chen, G., “An application of swarm intelligence binary particle swarm optimization (BPSO) algorithm to multi-focus image fusion”, Optica Applicata, Vol. 40, No. 4, pp. 949-964, 2010.

[40] Li, S., Kwok, J. T., Wang, Y., “Multi-focus image fusion using artificial neural networks”, Pattern Recognition Letters, Vol. 23, No. 8, pp. 985-997, 2002.

[41] Sheela Rani, C. M., Kumar, V. V., Ramana Reddy, B. V., “Improved block based feature level image fusion technique using multi-wavelet with neural network”, International Journal of Soft Computing and Engineering, Vol. 2, No. 4, pp. 265-271, 2012.

[42] Rashedi, E., Nezamabadi-pour, H., Saryazdi, S., “GSA: A gravitational search algorithm”, Information Sciences, Vol. 179, No. 13, pp. 2232-2248, 2009.

[43] Rashedi, E., Nezamabadi-pour, H., Saryazdi, S., “Improving the precision of CBIR systems by color and texture feature adaptation using GSA”, Computational Intelligence in Electrical Engineering, Vol. 4, No. 3, pp. 71-82, 2013.

[44] Pyari, M. P., Ganapati, P., “Solving multiobjective problems using cat swarm optimization”, Expert Systems with Applications, Vol. 39, No. 3, pp. 2956-2964, 2012.

[45] Wang, Z. H., Chang, C. C., Li, M. C. “Optimizing least-significant-bit substitution using cat swarm optimization strategy”, Information Sciences, Vol. 192, No. 9, pp. 98-108, 2012.

[46] Yang, F., Ding, M., Zhang, X., Hou, W., Zhong, C., “Non-rigid multi-modal medical image registration by combining L-BFGS-B with cat swarm optimization”, Information Sciences, Vol. 316, pp. 440-456, Sep. 2015.

[47] Shams, M., Rashedi, E., Hakimi, A., “Clustered-gravitational search algorithm and its application in parameter optimization of a low noise amplifier”, Applied Mathematics and Computation, Vol. 258, pp. 436-453, May 2015.

[48] Soleimanpour-moghadam, M., Nezamabadi-pour, H., Farsangi, M. M., “A quantum inspired gravitational search algorithm for numerical function optimization”, Information Sciences, Vol. 267, pp. 83-100, May 2014.

[49] Doraghinejad, M., Nezamabadi-pour, H., Mahani, A., “Channel assignment in multi-radio wireless mesh networks using an improved gravitational search algorithm”, Journal of Network and Computer Applications, Vol. 38, pp. 163-171, Feb. 2014.

[50] Sarafrazi, S., Nezamabadi-pour, H., Saryazdi, S., “Disruption: A new operator in gravitational search algorithm”, Scientia Iranica D, Vol. 18, No. 3, pp. 539-548, 2011.

[51] Chen, Z., Yuan, X., Tian, H., Ji, B., “Improved gravitational search algorithm for parameter identification of water turbine regulation system”, Energy Conversion and Management, Vol. 78, pp. 306-315, Feb. 2014.

[52] Shaw, B., Mukherjee, V., Ghoshal, S. P., “Solution of reactive power dispatch of power systems by an opposition-based gravitational search algorithm”, Electrical Power and Energy Systems, Vol. 55, pp. 29-40, Feb. 2014.

[53] Jiang, S., Wang, Y., Ji, Z., “Convergence analysis and performance of an improved gravitational search algorithm”, Applied Soft Computing, Vol. 24, pp. 363-384, Nov. 2014.

[54] Li, C., Zhou, J., Xiao, J., Xiao, H., “Hydraulic turbine governing system identification using T–S fuzzy model optimized by chaotic gravitational search algorithm”, Engineering Applications of Artificial Intelligence, Vol. 26, No. 9, pp. 2073-2082, 2013.

[55] Han, X. H., Chang, X. M., “A chaotic digital secure communication based on a modified gravitational search algorithm filter”, Information Sciences, Vol. 208, pp. 14-27, Nov. 2012.

[56] Li, C., Li, H., Kou, P., “Piecewise function based gravitational search algorithm and its application on parameter identification of AVR system”, Neurocomputing, Vol. 124, pp. 139-148, Jan. 2014.

[57] Kennedy, J., Eberhart, R., “Particle swarm optimization”, Proc. International Conference on Neural Networks, pp. 1942-1948, Perth, Nov. 1995.

[58] Canga, E. F., Image Fusion, MEng. Thesis, Electrical & Electronic Engineering Department, University of Bath, 2002.

[59] Zhao, H., Li, Q., Feng, H., “Multi-focus color image fusion in the HSI space using the sum-modified-laplacian and a coarse edge map”, Image and Vision Computing, Vol. 26, No. 9, pp. 1285-1295, 2008.

[60] Savić, S., Babić, Z., “Color multifocus image fusion using empirical mode decomposition”, Telfor Journal, Vol. 5, No. 2, pp. 128-133, 2013.

[61] Rattanapitak, W., Udomhunsakul, S., “Comparative efficiency of color models for multi-focus color image fusion”, Proc. International MultiConference of Engineers and Computer Scientists, Vol. 2, pp. 1-6, Hong Kong, Mar. 2010.

آمار

تعداد مشاهده مقاله: 1,992

تعداد دریافت فایل اصل مقاله: 916

سامانه مدیریت نشریات علمی. طراحی و پیاده سازی از سیناوب

پیوندهای مفید

اخبار و اعلانات

آمار

بهینه‌سازی فرایند ترکیب تصاویر با استفاده از الگوریتم جستجوی گرانشی