TW496080B - Method and system for normalizing a plurality of signals having a shared component - Google Patents
Method and system for normalizing a plurality of signals having a shared component Download PDFInfo
- Publication number
- TW496080B TW496080B TW089107919A TW89107919A TW496080B TW 496080 B TW496080 B TW 496080B TW 089107919 A TW089107919 A TW 089107919A TW 89107919 A TW89107919 A TW 89107919A TW 496080 B TW496080 B TW 496080B
- Authority
- TW
- Taiwan
- Prior art keywords
- channels
- color channel
- channel
- digital image
- patent application
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 145
- 230000006870 function Effects 0.000 claims description 159
- 238000011049 filling Methods 0.000 claims description 78
- 230000009466 transformation Effects 0.000 claims description 78
- 230000002079 cooperative effect Effects 0.000 claims description 37
- 238000010586 diagram Methods 0.000 claims description 37
- 230000008859 change Effects 0.000 claims description 28
- 238000004364 calculation method Methods 0.000 claims description 23
- 230000004044 response Effects 0.000 claims description 19
- 238000003860 storage Methods 0.000 claims description 17
- 238000004422 calculation algorithm Methods 0.000 claims description 15
- 238000001914 filtration Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 11
- 239000003973 paint Substances 0.000 claims 1
- 239000013589 supplement Substances 0.000 claims 1
- 239000002245 particle Substances 0.000 description 136
- 239000013598 vector Substances 0.000 description 95
- 238000012937 correction Methods 0.000 description 40
- 238000010606 normalization Methods 0.000 description 18
- 230000000875 corresponding effect Effects 0.000 description 17
- 239000011159 matrix material Substances 0.000 description 17
- 238000013459 approach Methods 0.000 description 12
- 238000005259 measurement Methods 0.000 description 11
- 230000009467 reduction Effects 0.000 description 11
- 230000006872 improvement Effects 0.000 description 9
- 230000011218 segmentation Effects 0.000 description 8
- 230000001629 suppression Effects 0.000 description 8
- 238000004891 communication Methods 0.000 description 6
- 238000009434 installation Methods 0.000 description 6
- 239000003086 colorant Substances 0.000 description 5
- 230000001419 dependent effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000001434 glomerular Effects 0.000 description 4
- 230000000717 retained effect Effects 0.000 description 4
- 230000002238 attenuated effect Effects 0.000 description 3
- 230000002596 correlated effect Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000011068 loading method Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000000844 transformation Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 101150108928 CCC1 gene Proteins 0.000 description 2
- 101150050425 CCC2 gene Proteins 0.000 description 2
- 238000007792 addition Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 235000012054 meals Nutrition 0.000 description 2
- BQJCRHHNABKAKU-KBQPJGBKSA-N morphine Chemical compound O([C@H]1[C@H](C=C[C@H]23)O)C4=C5[C@@]12CCN(C)[C@@H]3CC5=CC=C4O BQJCRHHNABKAKU-KBQPJGBKSA-N 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- OVSKIKFHRZPJSS-UHFFFAOYSA-N 2,4-D Chemical compound OC(=O)COC1=CC=C(Cl)C=C1Cl OVSKIKFHRZPJSS-UHFFFAOYSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000013551 empirical research Methods 0.000 description 1
- 239000003292 glue Substances 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 229960005181 morphine Drugs 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009131 signaling function Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration using non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20032—Median filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
Description
496080 A7 B7 經濟部智慧財產局員Η消費合作社印製 五、發明說明(1 ) 發明技術領域 本發明是關於訊號處理,更特別是與複數個具有一共 用的成分之歸一化相關。 相關應用的互相參考 這個應用與由亞伯特•愛德格所提出,在1 9 9 9年8 月2 0日歸檔,代理商的文件摘要號碼第 021971.0127號,且標題爲”自動地改善數位 式影像品質的方法和裝置”的美國應用序號_相 關。 這個應用與由亞伯特•愛德格所提出,在1 9 9 9年8 月2 0日歸檔,代理商的文件摘要號碼第 02 1971 · 0 104號,且標題爲”改善數位式影像 的方法和系統”的美國應用序號__相關。 發明背景 在許多訊號處理應用中,多數個訊號可能存在有一共 用的成分。在一非線性的方式下,一或多個這些相關的訊 號可能變成失真。在某些應用中,每一個訊號可以一不同 的方式使其失真或不失真,且失真的性質可能無法預測。 此外,可能很難決定哪一個相關的訊號變成失真。這樣的 失真正常是不理想的。 例如:在影像處理的領域中,無論影像爲被掃描的照 相影像、人造衛星影像、醫療影像、或其他類型的數位式 (請先閱讀背面之注意事項再填寫本頁) .3丨 裝496080 A7 B7 Printed by a member of the Intellectual Property Bureau of the Ministry of Economic Affairs and a Consumer Cooperative. V. Description of the Invention (1) Technical Field of the Invention The present invention relates to signal processing, and more particularly relates to the normalization of a plurality of components having a common component. Cross Reference to Related Applications This application and the one proposed by Albert Edger were archived on August 20, 1999, with the agent's document abstract number 021971.0127, and the title was "Automatically Improve Digital Method and device for image quality of the United States "US application serial number_related. This application was proposed by Albert Edger, filed on August 20, 1999, with the agent's document summary number 02 1971 · 0 104, and the title is "Improving Digital Imaging Methods and Systems "US Application Serial Number __ related. BACKGROUND OF THE INVENTION In many signal processing applications, most signals may have a common component. In a non-linear manner, one or more of these related signals may become distorted. In some applications, each signal can be distorted or undistorted in a different way, and the nature of the distortion may be unpredictable. In addition, it may be difficult to decide which related signal becomes distorted. Such distortion is not ideal. For example: in the field of image processing, whether the image is a scanned photo image, satellite image, medical image, or other types of digital (please read the precautions on the back before filling this page).
n n n I 一 θ” · flu 1 I n «ϋ an I 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -4- 經濟部智慧財產局員工消費合作社印製 496080 A7 B7__ 五、發明說明(2 ) 影像’都可能遭受非線性的失真。這樣的失真可能存在於 影像以一類比方式表示和/或可能在數位式影像產生時發 生。舉例而言,在照相影像的掃描中,在掃描過程中,一 典型的數位式影像的三種色彩波道(channel )每一個可 能遭受非線性失真,但是,這樣的非線性失真在每一色彩 波道中可能不相同。這樣的非線性失真可能以吾人不想要 的方式,將影像降級(degrade )。對照相影像而言,影 像可能看起來不那麼合意。對人造衛星影像和醫療影像而 言,由影像所提供的資訊,由於這樣的失真,可能不那麼 有用。 發明槪述 本發明的一方面是一種用來使具有一共用的成分的複 數個訊號歸一化的方法,而其中,訊號中的至少一個已經 被以非線性的方式失真。一種失真函數是對複數個訊號中 的至少一個而決定,其中,該訊號函數與相對所剩餘的訊 號中的至少一個的一特別訊號之失真成正比。一種相對反 向的失真函數爲此失真訊號而產生,其中,這函數是爲回 應爲此訊號所決定的失真函數。藉由應用爲此訊號所產生 的相對反向的失真函數而將所討論的訊號歸一化。 本發明具有許多重要的技術優點。本發明的各種不同 的實施例可能不具有、具有一、一些、或所有的這些優點 ,而不悖離本發明的範圍。本發明允許具有一共用的成分 的訊號彼此之間自動地被歸一化,其中’這些訊號中的一 (請先閱讀背面之注意事項再填寫本頁) ·裝--------訂------·— 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -5- 496080 經濟部智慧財產局員工消費合作社印製 A7 B7 五、發明說明(3 ) 或多個已經被以非線性的方式失真。雖然,此失真也許不 會完全地被移除,但是,該訊號可以更正確地代表所要的 訊:號。本發明可以有利地被應用在影像處理的範圍中。例 如:本發明可以被用來修正代表一彩色的照相影像的一數 位式影像之彩色波道中的失真。此彩色波道可以被歸一化 ,以致於整體的影像被改良,且影像的色彩變得看起來更 合意。對於人造衛星影像和醫學影像,本發明可以被用來; 增加影像所提供的資訊的正確度和可用性。因爲本發明可 以使用一電腦和/或掃描器來執行這樣的歸一化,所以, 本發明允許影像的自動歸一化。這樣的自動歸一化可以減 少達到這樣的歸一化所需的時間和力氣,而允許一較大量 的影像之歸一化。例如:在照相影像的情形,可以使用諸 如:Photoshop的軟體,人工地修正非線性失真。本發明 可以在少許的時間、花少許的成本且修正該影像的人士無 須專門的知識之下,讓任何所要的影像的自動彩色修正。 附圖簡述 爲了更完全地了解本發明和其中的優點,吾人參考如 下的說明,連同附圖,其中: 圖1例示根據本發明,可以被使用的一通用電腦之方 frft fgj · 塊圖, 圖2例示包含本發明的一實施例之掃描器之例子; 圖3例示根據本發明的一種方法’描述一數位式影像 的改良的一流程圖; 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) (請先閱讀背面之注意事項再填寫本頁)nnn I a θ ”· flu 1 I n« ϋ an I This paper size is applicable to China National Standard (CNS) A4 (210 X 297 mm) 2. Description of the Invention (2) Images' may suffer from non-linear distortions. Such distortions may exist when the images are represented in an analogous manner and / or may occur when digital images are generated. For example, in scanning of photographic images During the scanning process, each of the three color channels of a typical digital image may suffer from non-linear distortion, but such non-linear distortion may be different in each color channel. Such non-linear distortion May degrade the image in a way we do n’t want. For photographic images, the images may not look so desirable. For satellite images and medical images, the information provided by the images due to such distortion May not be so useful. DESCRIPTION OF THE INVENTION One aspect of the present invention is a method for normalizing a plurality of signals having a common component. Method, wherein at least one of the signals has been distorted in a non-linear manner. A distortion function is determined for at least one of a plurality of signals, wherein the signal function is relative to at least one of the remaining signals. The distortion of a particular signal is proportional. A relatively inverse distortion function is generated for this distortion signal, where this function is a response to the distortion function determined for this signal. By applying the relative inversion generated for this signal Normalizes the signal in question. The invention has many important technical advantages. Various embodiments of the invention may not have, have one, some, or all of these advantages without departing from this Scope of the invention. The present invention allows signals with a common component to be automatically normalized to each other, where 'one of these signals (please read the precautions on the back before filling out this page). ---- Order ------ · — This paper size is applicable to China National Standard (CNS) A4 specification (210 X 297 mm) -5- 496080 Staff Consumption of Intellectual Property Bureau, Ministry of Economic Affairs Printed by Sakusha A7 B7 V. Invention Description (3) or more has been distorted in a non-linear manner. Although this distortion may not be completely removed, the signal can more accurately represent the desired signal : No. The present invention can be advantageously applied in the scope of image processing. For example, the present invention can be used to correct the distortion in the color channel of a digital image representing a color photographic image. This color channel can be Normalize, so that the overall image is improved, and the color of the image becomes more desirable. For artificial satellite images and medical images, the present invention can be used to increase the accuracy and availability of the information provided by the images. Since the present invention can perform such normalization using a computer and / or scanner, the present invention allows automatic normalization of images. Such automatic normalization can reduce the time and effort required to achieve such normalization, while allowing normalization of a larger number of images. For example, in the case of photographic images, software such as: Photoshop can be used to manually correct non-linear distortion. The present invention allows automatic color correction of any desired image with a little time, a little cost, and the person who corrects the image does not need special knowledge. Brief description of the drawings In order to fully understand the present invention and its advantages, I refer to the following description, together with the drawings, in which: Figure 1 illustrates a block diagram of a general-purpose computer frft fgj · according to the present invention, which can be used; FIG. 2 illustrates an example of a scanner including an embodiment of the present invention; FIG. 3 illustrates a flow chart describing a method of improving a digital image according to the present invention; the paper size is applicable to the Chinese National Standard (CNS) A4 specification (210 X 297 mm) (Please read the notes on the back before filling this page)
496080 A7 B7 五、發明說明(4 ) 圖4例示根據本發明,描述一種測量與像素強度( p i X e 1 i n t e n s丨t y )相關的顆粒軌跡(g r a i n t r a c e )的方法之 流程圖; 圖5例不被分割成區段(s e g m e π t )的一影像; 圖6例不一數位式影像的紅色波道的加權(w e i g h t )、 平均強度(intensity)、和顆粒強度之圖; 圖7例示一數位式影像的一種色彩波道的一散佈( scatter )圖; 圖8例示使用曲線擬合演算法(curve fitting a 1 g〇M thm)所晝出經過數位式影像的一種色彩波道的一散 佈圖之線; 圖9例示所畫出經過一數位式影像的紅色資訊的一散 佈圖之線; 圖1 0例示根據本發明,描述一種用來測量與頻率相 關的顆粒軌跡的一種方法之流程圖; 圖1 1例示一數位式影像的一種色彩波道的一散佈圖 9 圖1 2例示所畫出經過一數位式影像的一種色彩波道 的一散佈圖之線; 圖1 3例示從被歸一化的、具有在畫在空間上像素強 度中的改變之數位式影像而來的像素的線上的.穎粒軌跡; 圖1 4例示顆粒軌跡的頻率數値對在跨強度而被歸一 化的一數位式影像中的頻率所畫的圖; 圖1 5例示顆粒軌跡的頻率數値對在跨強度和頻率而 本纸張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) (請先閱讀背面之注意事項再填寫本頁) 裝--------訂------.—· 經濟部智慧財產局員工消費合作社印製 496080 A7 B7 經濟部智慧財產局員工消費合作社印製 五、發明說明(5 ) 被歸一化的一數位式影像中的頻率所晝的圖; 圖1 6例示說明用來從一數位式影像移除諸如:顆粒 軌跡等的高頻雜訊的一種方法的舉例之流程圖; 圖1 7例示說明用來從一數位式影像移除諸如:顆粒 軌跡等的高頻雜訊的一種方法的舉例之流程圖; 圖1 8例示代表與本發明關聯的二維傅立葉變換( Fourier transform)的方法; 圖19例示說明可被用來將具有一共用的成分的複數 個訊號歸一化的方法例子的步驟之流程圖,其中,複數個 訊號中的至少一個被以非線性的方式使之失真; 圖2 0例示說明根據本發明,可被用來改良一數位式 影像的一方法例子之流程圖;和 圖2 1例示可被用來連接例示於圖1 9或圖2 0中的 方法的三個散佈圖例子。 元件對照表 1 0 :通用電腦 1 2 :處理器 1 4 :隨機存取記憶體 1 6 :唯讀記憶體 1 8 :滑鼠 2 0 :鍵盤 2 2 :磁碟機 2 4 :印表機 (請先閱讀背面之注意事項再填寫本頁) 裝 ----訂------.---邊 %· 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -8- 496080 A7 B7 經濟部智慧財產局員工消費合作社印製 五、發明說明(6 ) 2 6 :顯示器 2 8 :通訊連結 ’ 3 0 :掃描器 3 4 :掃描器 3 6 :處理器 3 8 :儲存媒體 4 0 :掃描硬體 4 2 :影像處理軟體 4 4 :控制軟體 較佳實施例之詳細說明 藉由參考附圖的圖1至1 5,如被使用的數字和各種 不同附圖的對應部分,本發明的較佳實施例與其優點最容 易被了解。 圖1例示根據本發明,可以被用作影像改良的一通用 電腦1 0之方塊圖。特別的是,通用電腦1 〇可以包括一 影像改良系統的一部分,也可以被用來執行包含影像改良 軟體的應用程式。通用電腦1 〇可以被應用於執行任何眾 所週知的^13 — 0〇3、?(:一0〇3、〇32、 UN I X、MAC — OS、和視窗(Windows )操作系統 或其他操作系統。通用電腦1 0包括處理器1 2、隨機存 取記憶體(R A Μ ) 1 4、唯讀記憶體(R a Μ ) 1 6、 滑鼠1 8、鍵盤2 0和輸入/輸出裝置,諸如:印表機 2 4 、磁碟機2 2 、顯示器2 6和通訊連結( (請先閱讀背面之注意事項再填寫本頁) 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐〉 496080 Α7 Β7 五、發明說明(7 ) (請先閱讀背面之注意事項再填寫本頁) communication link) 2 8。本發明包括:可以儲存於 RAMI 4、R〇M1 6 、或磁碟機2 2中且可以由處理 器· 1 2所執行的程式。通訊連結2 8被連接至一電腦網路 ,也可被連接至一電3舌線、一天τ泉、一蘭道或任何其他型 式的通訊連結。磁碟機2 2可以包括各種不同的儲存媒體 ,諸如:軟式磁碟機、硬式磁碟機、C D - ROM磁碟機 或磁帶機。雖然這個實施例使用複數個磁碟機2 2,吾人 也可使用單一個磁碟機2 2 ,而不悖離本發明的範圍。圖 1僅只提供可使用於本發明的一電腦的例子。本發明可以 被使用在非通用電腦或不具傳統操作系統的通用電腦。 經濟部智慧財產局員工消費合作社印製 通用電腦1 0更包括掃描器3 0 ,其中該掃描器可被 用來掃描根據本發明的教導而被改良的影像。在這實施例 中,影像改良可以由被儲存在掃描器3 0中且由掃描器 3 0所執行的軟體所實行,且將其結果儲存在包含掃描器 3 0的一部分的一儲存媒體且/或在通用電腦1 〇中的任 何儲存裝置中。替換地,改良影像的軟體可被儲存於與通 用電腦1 0相關聯的任何儲存媒體中,且可以由處理器 1 2所執行,以改良由掃描器3 0所掃描的影像。此外’ 影像改良可以內在地發生在掃描器3 0與通用電腦1 〇中 ,如所述,而不悖離本發明的範圍。掃描器3 0可以包括 一影片掃描器或任何形式的一平台掃描器,而不悖離本發 明的範圍。 圖2例示根據本發明而建構的一掃描器3 4之例子。 掃描器3 4包含··處理器3 6、儲存媒體3 8和掃描硬體 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) ~\〇Ζ 496080 Α7 Β7 經濟部智慧財產局員工消費合作社印製 五、發明說明(8 ) 4〇。處理器3 6藉由執行在儲存媒體3 8中所儲存的控 制軟體4 4而控制掃描硬體4 0的操作。雖然爲了簡單起 見,吾人例示一單一儲存媒體,可是,儲存媒體3 8可以 包括多個儲存媒體,也可以包含不同型式的儲存媒體。因 此,例如:控制軟體4 4可以被儲存於R〇Μ記憶體、 RAM記憶體或一硬碟中。掃描硬體4 0被用來利用某類 型的光學電路系統,將一類比影像轉換成一數位式影像。/ 任何形式的光學電路系統均可被使用於掃描硬體4 0中, 而不悖離本發明的範圍。在掃描硬體4 0掃描影像之後, 那影像可以根據本發明,而使用儲存於儲存媒體3 8中的 影像處理軟體4 2而被改良。同樣地,被掃描的影像可被 儲存於儲存媒體3 8中,改良後的影像也可以。替換地, 掃描器3 4可以不具任何影像處理軟體4 2。這種軟體反 而可以由通用電腦1 0所提供,用來改良從掃描器3 4接 收的一影像。因此,一掃描後的影像和/或一改良後的掃 描影像可以透過一通訊璋(communication port)(未明確 顯示出),由掃描器3 4供應至通用電腦1 〇。雖然吾人 例示舉例的一掃描器3 4之實施例可被用於與本發明相關. 的影像改良,但是,也可使用其他掃描器,而不悖離本發 明的範圍。 圖3例示描述本發明的一實施例所用的一種方法,以 改良一數位式影像的一流程圖。此處所描述的步驟可以使 用電腦軟體而被實行,以下所述的任何方法的步驟也可如 此。此軟體,如上述,可以由掃描器3 4或通用電腦1 0 (請先閱讀背面之注意事項再填寫本頁) 裝 ----訂------.--- ^· 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公餐) -11 - 496080 經濟部智慧財產局員工消費合作社印製 A7 B7 五、發明說明(9 ) 所執行。從掃描器3 0以外所接收的一數位式影像也可以 根據本發明而被改良。在步驟1 0 0 2中,一數位式影像 的:穎粒軌跡相關於像素強度而被測量。在步驟1 〇 〇 4中 ,穎粒軌跡相關於頻率而被測量。在步驟1 〇 〇 6中,校 正陣列(correction array )被計算出。此校正陣列在步 驟1 0 0 8中,被用來將數位式影像與強度相關的顆粒軌 跡歸一化,而後,在步驟1 0 1 0中,將數位式影像相關 於頻率的顆粒軌跡歸一化。此被歸一化的顆粒軌跡資料可 被用來抑制在步驟1 0 1 2中的顆粒軌跡。此外,如步驟 1〇1 4中所示,藉由歸一化顆粒軌跡而得的像素強度値 可被用來平衡所有色彩波道中的這些像素強度。這些步驟 中有某些步驟可能可以被省略,或其他步驟也可以被包括 進來,而不悖離本發明的範圍。 雖然本發明可以被用在由照相影像所產生的數位式影 像,但是,它也可被用於將高頻雜訊歸一化,諸如:在其 他型式的影像中之顆粒軌跡,諸如:人造衛星影像、醫學 影像等等。以下的討論說明黑白和彩色的影像之歸一化。 對彩色影像而言,如下所述的方法可被使用,或者對黑白 影像所述的方法可被使用於彩色影像中的每一個別波道上 。替換地,黑白的方法可被使用於波道的一子集合,而彩 色的方法則被用於其餘的波道上。類似的選項應用於非照 相影像的數位式影像上。 測量相關於強度的顆粒軌跡 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -12- -----^----^----裝--------訂--- (請先閱讀背面之注意事項再填寫本頁) 囈· 496080 經濟部智慧財產局員工消費合作社印製 A7 B7 五、發明說明(i〇) 本發明包括測量相關於像素強度的顆粒軌跡之一種方 法。圖4顯示涉及這方法的步驟。在這應用中’與這個和 其’他方法相關的所有計算可以使用通用電腦1 0、掃描器 3 4或這些裝置之一與其他裝置合倂的組合而執行之。在 步驟1 1 0 2中,此數位式影像可以被分割爲幾個區段’ 接著在步驟1 1 0 4中,進行一諸如:傅立葉變換的變換 ,而後,在步驟1 1 0 6中,移除低和中頻率項。然後’ 在步驟1 1 0 8中,在每一色彩波道中,各種不同的公式 可被用來計算每一區段的顆粒強度(代表因顆粒軌跡引起 的像素強度的大小之偏差的値)。在步驟1 1 1 〇中’爲 了幫忙驗證這些顆粒強度,可以使用一各別的公式’來計 算在每一色彩波道中每一區段的一”加權”値。其次,在 步驟1 1 1 2中,可以由加權値、DC(direct current, 直流電壓,慣例上代表一訊號的平均値)、和每一色彩波道 的顆粒強度決定而產生散佈圖。在步驟1 1 1 4中’吾人 可以使用一曲線擬合演算法,畫出經過散佈圖之線,以顯 示每一像素強度的適當顆粒強度。 區段分割 在步驟1 1 0 2中,數位式影像被分割成資料區段, 每一區段對應數位式影像的一空間區域。區段分割可被應 用至數位式影像的每一波道,或在一黑白影像的情形下的 單一波道。如圖5所表示,這種分割的一代表顯示複數個 區段,如:區段1 2 0 2。雖然吾人可以使用任何區段大 (請先閱讀背面之注意事項再填寫本頁) 裝—------訂----------· 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -13 - 496080 經濟部智慧財產局員工消費合作社印製 B7_____五、發明說明(11 ) 小,但是,每一區段的面積最好大到足夠包含統計上有意 義的數目的像素,諸如:在1 6到1 0 2 4像素之間,而 卻:不能大到包含許多實際上不同的強度之影像的面積(在 2 0 0 0乘3 0 0 0像素的一影像大小)。在大到能夠提 供足夠的資料,又不會太大,與,因此,在一 2 0 0 0像 素乘以3 0 0 0像素級數的一影像的單一區段之中,具有 太多不同的強度成分之間,包含2 5 6 ( 1 6 X 1 6區段 / )個不同的像素的區段可以提供一良好的折衷。使用這種 方法,單一 3 5 m m負片(或幻燈透明片)將典型地形成 大約2 4 K非重疊區段。任何影像大小或任何區段大小均 可被使用,而不悖離本發明的範圍。 圖5將這些區段顯示爲正方形,但是,它們可以是其 他形狀,而不悖離本發明的範圍。不同的區段也無須具相 同的形狀和大小,雖然,使他們一致相同在處理時可能有 用。在一特定的實施例中,區段的數目和他們的形狀與大 小可以由與頻率、強度或頻率內容相關聯的影像之變化所 定義。同時注意:區段可以被重疊和窗格,·如:在 _歸檔,美國專利應用第_號 的“被窗格的影像方塊之混合”中所述’此文已經由參考 而被倂入,宛如在此完全被發表。 傅立葉變換 在將此數位式影像區段分割之後’吾人可以在每一色 彩波道(或在單一黑白波道)上’對每一區段’在像素強 (請先閱讀背面之注意事項再填寫本頁) 裝496080 A7 B7 V. Description of the invention (4) FIG. 4 illustrates a flowchart describing a method for measuring a grain trace related to pixel intensity (pi X e 1 intens 丨 ty) according to the present invention; FIG. 5 does not show An image divided into sections (segme π t); Figure 6 illustrates the weighting, average intensity, and particle intensity of the red channel of different digital images; Figure 7 illustrates a digital image A scatter diagram of a color channel of FIG .; FIG. 8 illustrates a line of a scatter diagram of a color channel that passes through a digital image using a curve fitting algorithm (curve fitting a 1 g0M thm). Figure 9 illustrates a scatter diagram line of red information passing through a digital image; Figure 10 illustrates a flowchart describing a method for measuring frequency-dependent particle trajectories according to the present invention; Figure 1 1 illustrates a scatter of a color channel of a digital image. Figure 9 2 illustrates a line of a scatter diagram of a color channel that passes through a digital image. The morphine particle trace on a line of pixels from a digital image with a change in pixel intensity drawn in space. Figure 14 illustrates the frequency of the particle trajectory normalized across the intensity. A graph of the frequency in a digital image; Figure 15 illustrates the frequency of the particle trajectory. For the intensity and frequency across the paper, this paper scale applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) ( Please read the precautions on the back before filling this page.) -------- Order ------.------ Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs 496080 A7 B7 Printed by a consumer cooperative. 5. Description of the invention (5) A graph of the frequency of day in a digital image normalized; Figure 16 illustrates the removal of high frequencies such as particle trajectories from a digital image. A flowchart of an example of a method of noise; Figure 17 illustrates a flowchart of an example of a method for removing high-frequency noise such as: particle traces from a digital image; Figure 18 illustrates a representative and The two-dimensional Fourier transform associated with the present invention (Fourier t ransform) method; FIG. 19 illustrates a flowchart illustrating the steps of an example method that can be used to normalize a plurality of signals having a common component, in which at least one of the plurality of signals is non-linearly Distort; FIG. 20 illustrates a flowchart illustrating an example of a method that can be used to improve a digital image according to the present invention; and FIG. 21 illustrates an example that can be used to connect an example illustrated in FIG. 19 or FIG. 20 Three scatter diagram examples of the method. Component comparison table 1 0: General-purpose computer 1 2: Processor 1 4: Random access memory 1 6: Read-only memory 1 8: Mouse 2 0: Keyboard 2 2: Disk drive 2 4: Printer ( Please read the precautions on the back before filling in this page) Binding ---- Ordering --------.--- Edge% · This paper size applies to China National Standard (CNS) A4 (210 X 297 mm) -8- 496080 A7 B7 Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs. 5. Description of the invention (6) 2 6: Display 2 8: Communication link '3 0: Scanner 3 4: Scanner 3 6: Processor 3 8 : Storage medium 4 0: Scanning hardware 4 2: Image processing software 4 4: Detailed description of the preferred embodiment of the control software By referring to Figures 1 to 15 of the drawings, such as the numbers used and the various drawings Corresponding parts, the preferred embodiments of the present invention and their advantages are most easily understood. Fig. 1 illustrates a block diagram of a general-purpose computer 10 which can be used for image improvement according to the present invention. In particular, the general-purpose computer 10 may include a part of an image improvement system, or may be used to execute an application program including image improvement software. The general-purpose computer 10 can be used to perform any of the well-known ^ 13 — 0〇3,? (: 003, 032, UN IX, MAC-OS, and Windows (Windows) operating system or other operating systems. General-purpose computer 10 includes processor 1 2. Random access memory (RA M) 1 4 , Read-only memory (R a M) 1 6, mouse 18, keyboard 20, and input / output devices such as printer 2 4, disk drive 2 2, display 26 and communication links ((Please Please read the notes on the back before filling this page) This paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm> 496080 Α7 Β7 V. Description of the invention (7) (Please read the notes on the back before filling (This page) communication link) 2 8. The present invention includes: programs that can be stored in RAMI 4, ROM16, or disk drive 2 2 and can be executed by the processor 1 2. The communication link 2 8 is connected To a computer network, it can also be connected to an electric 3 tongue line, a day τ spring, a Landau or any other type of communication link. The drive 22 can include a variety of different storage media, such as: floppy disk Drive, hard drive, CD-ROM drive, or tape drive. Although this embodiment With multiple disk drives 22, we can also use a single disk drive 22 without departing from the scope of the present invention. Figure 1 provides only an example of a computer that can be used with the present invention. The invention can be used On general-purpose computers or general-purpose computers without traditional operating systems. The general consumer computer 10 printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs further includes a scanner 3 0, which can be used to scan according to the teachings of the present invention. Improved image. In this embodiment, the image improvement may be performed by software stored in the scanner 30 and executed by the scanner 30, and the result may be stored in a file containing a part of the scanner 30. A storage medium and / or in any storage device in the general-purpose computer 10. Alternatively, the software for improving the image may be stored in any storage medium associated with the general-purpose computer 10 and may be executed by the processor 12 To improve the image scanned by the scanner 30. In addition, 'image improvement can occur inherently in the scanner 30 and the general-purpose computer 10, as described, without departing from the scope of the present invention. The scanner 30 may include a film scanner or a platform scanner in any form without departing from the scope of the present invention. Figure 2 illustrates an example of a scanner 34 constructed in accordance with the present invention. The scanner 34 includes ·· Processor 3 6, Storage Media 3 8 and Scanning Hardware The paper size is applicable to Chinese National Standard (CNS) A4 (210 X 297 mm) ~ \ 〇 496080 Α7 Β7 Printed by the Intellectual Property Bureau of the Ministry of Economic Affairs and Consumer Cooperatives The fifth, invention description (8) 40. The processor 36 controls the operation of scanning the hardware 40 by executing the control software 4 4 stored in the storage medium 38. Although for the sake of simplicity, I exemplify a single storage medium, the storage medium 38 may include multiple storage media or different types of storage media. Therefore, for example, the control software 44 can be stored in ROM memory, RAM memory or a hard disk. Scanning hardware 40 is used to convert an analog image into a digital image using a certain type of optical circuit system. / Any form of optical circuit system can be used in the scanning hardware 40 without departing from the scope of the present invention. After scanning the image by the hardware 40, the image can be improved according to the present invention using the image processing software 42 stored in the storage medium 38. Similarly, the scanned image can be stored in the storage medium 38, and the improved image can also be stored. Alternatively, the scanner 34 may be provided without any image processing software 42. Such software can instead be provided by a general-purpose computer 10 for improving an image received from the scanner 34. Therefore, a scanned image and / or an improved scanned image can be supplied to the general-purpose computer 10 by the scanner 34 through a communication port (not explicitly shown). Although the embodiment of a scanner 34 exemplified by us can be used for image improvement related to the present invention, other scanners can be used without departing from the scope of the present invention. FIG. 3 illustrates a flowchart describing a method used to improve a digital image according to an embodiment of the present invention. The steps described here can be performed using computer software, as can any of the methods described below. This software, as mentioned above, can be installed by the scanner 3 4 or the general-purpose computer 10 (please read the precautions on the back before filling in this page). ---- Order ------.--- ^ · This paper Standards are applicable to China National Standard (CNS) A4 specifications (210 X 297 meals) -11-496080 Printed by A7 B7, Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs 5. Implementation of the invention description (9). A digital image received from other than the scanner 30 can also be improved according to the present invention. In step 1 0 02, the glomerular trajectory of a digital image is measured in relation to the pixel intensity. In step 1004, the glomerular locus is measured as a function of frequency. In step 106, a correction array is calculated. The correction array is used to normalize the particle trajectories of the digital image and the intensity in step 10 08, and then, in step 10 10, the particle trajectories of the digital image to frequency are normalized. Into. This normalized particle trajectory data can be used to suppress the particle trajectory in step 1012. In addition, as shown in step 104, the pixel intensity 値 obtained by normalizing the particle trajectory can be used to balance these pixel intensities in all color channels. Some of these steps may be omitted, or other steps may be included without departing from the scope of the invention. Although the present invention can be applied to digital images produced from photographic images, it can also be used to normalize high-frequency noise, such as particle trajectories in other types of images, such as artificial satellites Imaging, medical imaging, and more. The following discussion illustrates the normalization of black and white and color images. For color images, the method described below can be used, or for black and white images, the methods described above can be used for each individual channel in the color image. Alternatively, the black and white method can be used for a subset of the channels, and the color method can be used for the remaining channels. Similar options apply to digital images that are not photographic. Measure the particle trajectory related to the strength. The size of this paper is applicable to China National Standard (CNS) A4 (210 X 297 mm) -12- ----- ^ ---- ^ ---- pack ----- --- Order --- (Please read the precautions on the back before filling out this page) 呓 · 496080 Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs A7 B7 V. Description of the invention (i〇) The invention includes measurement related to pixels A method of particle trajectory of intensity. Figure 4 shows the steps involved in this method. In this application, all calculations related to this and other methods may be performed using a general-purpose computer 10, a scanner 34, or a combination of one of these devices with other devices. In step 1 102, the digital image can be divided into several sections. Then, in step 1 104, a transformation such as a Fourier transform is performed, and then in step 1 106, shifting is performed. Divide low and medium frequency terms. Then, in step 1 108, in each color channel, various formulas can be used to calculate the particle intensity (representing the deviation of the magnitude of the pixel intensity due to the particle trajectory) of each segment. In step 1 1 10, 'to help verify the intensity of these particles, a separate formula' can be used to calculate a "weight" for each segment in each color channel. Secondly, in step 11 12, a scatter diagram can be generated by determining the weighted chirp, DC (direct current, direct voltage, which conventionally represents the average chirp of a signal), and the particle intensity of each color channel. In step 1 1 1 4 we can use a curve fitting algorithm to draw a line through the scatter diagram to show the proper particle intensity for each pixel intensity. Segment segmentation In step 1 102, the digital image is divided into data segments, and each segment corresponds to a spatial region of the digital image. Segmentation can be applied to each channel of a digital image, or a single channel in the case of a black and white image. As shown in Fig. 5, a representative of this segmentation displays a plurality of sections, such as: section 1 2 0 2. Although we can use any section (please read the precautions on the back before filling out this page) Packing ---------- Order ---------- · This paper size applies to Chinese National Standards (CNS ) A4 size (210 X 297 mm) -13-496080 Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs B7_____ V. Invention description (11) Small, but the area of each section is preferably large enough to include statistics A meaningful number of pixels, such as: between 16 and 10 24 pixels, but not large enough to contain an area of many images of virtually different intensities (between 2 0 0 by 3 0 0 0 pixels An image size). In a single segment of an image that is large enough to provide enough data, but not too large, and therefore has a number of 20,000 pixels multiplied by 3,000 pixels, there are too many different Among the intensity components, a segment containing 2 5 6 (16 X 1 6 segments /) different pixels can provide a good compromise. Using this method, a single 35 mm negative (or slide transparencies) will typically form approximately 24 K non-overlapping segments. Any image size or any segment size can be used without departing from the scope of the invention. Figure 5 shows these sections as squares, however, they can be other shapes without departing from the scope of the invention. The different sections need not be the same shape and size, although making them identical may be useful in processing. In a particular embodiment, the number of segments and their shape and size may be defined by changes in the image associated with frequency, intensity, or frequency content. Also note: sections can be overlapped and panes, such as: in "Archives, US Patent Application No. _" Mixed Image Boxes by Pane "" This article has been incorporated by reference, It seems as if it was published here. Fourier Transform After segmenting this digital image segment, 'we can each pixel' on each color channel (or on a single black-and-white channel) in pixel intensity (please read the notes on the back before filling (This page)
i §1 1 flu 一 口f ϋ II —ϋ ϋ— n n ϋ— I %, 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -14 - 496080 A7 B7 經濟部智慧財產局員工消費合作社印製 五、發明說明(12) 度値上執行一傅立葉變換,如圖4中的步驟1 1 0 4中所 示。形成該數位式影像的資料含有在該影像中,每一像素 的·強度値,且對每一色彩波道’在空間上排列。強度是指 一像素的亮度。例如:一白色的像素比~灰色或黑色的像 素具有較大的強度値。不同色彩波道的相同空間面積通常 具有不同的像素強度値。 傅立葉變換將每一區段的像素強度値從空間域轉換成 ; 頻率域。傅立葉變換所計算的其中一値爲每一色彩波道的 每一區段之平均強度,亦即“ D C ” 。傅立葉變換所計算 另一値爲每一色彩波道的每一區段中的每一頻率基元( element )的垂直和水平頻率向量的頻率和大小。雖然這實 施例使用一傅立葉變換,然而,也可使用其他型式的變換 將一訊號從空間域轉換成頻率域,諸如:一哈德瑪變換( Hadamard transfrom )。此外,頻率變換空間可以使用® 通滅波器而被細分。 移除低和中頻率 其次,在圖4中的步驟1 1 0 6中,每一變換的低和· 中頻率項可以從進一步的考量中被移除。替換地’這些項 可以被削弱、或這步驟被省略,而不悖離本發明的範圍。 在這實施例中,每一色彩的波道(或在一黑白影像的單一 波道)之每一區段在本發明的一特殊實施例的有用的範置1 內,接收這移除。具有顯著地高頻之顆粒強度較容易測量 和從影像資訊中區別。因爲在照相軟片上的實質顆粒傾向 (請先閱讀背面之注意事項再填寫本頁) 裝 n ·Βϋ I ^ · ·ϋ «-a— n n ·1 ΜΗ! -v-α- .% 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -15- 496080 Α7 Β7 經濟部智慧財產局員工消費合作社印製 五、發明說明(13) 於小而隨機佈置,顆粒軌跡在所有的頻率中約略相等的存 在,但是,在較高頻率下,具有較少的影像內容。爲了得 到最佳的影像細節和顆粒軌跡的區別’最好集中注意在高 頻率上。 在這實施例中,吾人使用一 1 6 X 1 6區段6肴次強 調:其他區段大小可被使用,而不悖離本發明的範圍。爲 了過濾掉低和中範圍頻率,吾人可以,例如:過濾掉所有 X 一頻率和y -頻率位置的絕對値的總合爲小於等於8的 頻率。替換地,若頻率向量如圖1 8中被配置,則可以選 擇在該D C値周圍的一適合的半徑,例如:7777=8 在這半徑範圍內的頻率可以被移除。其他的半徑或過 濾方法均可被使用’而不悖離本發明的範圍。 計筧顆粒軌跡 其次,在圖4中的步驟1 1 0 8中,吾人可使用互相 關(c r 〇 s s - c 〇 r r e 1 a t i ο η )公式來計算每一區段的顆粒強 度。在一黑白數位式影像情形中,因只存在一波道,故不 使用互相關公式。對於彩色影像和黑白影像而言,此計算 使用在移除低和中頻率項之後,從傅立葉變換得來的數位 式影像之頻率向量。 爲了瞭解這些計算,一數位式影像的解釋將是有用的 。包含一區段的總資訊是由每一色彩的波道的每一像素的 資訊所組成的。例如··在典型的三顏色數位式影像中,一 (請先閱讀背面之注意事項再填寫本頁) 裝i §1 1 flu f f II —ϋ ϋ — nn ϋ — I%, this paper size applies to China National Standard (CNS) A4 (210 X 297 mm) -14-496080 A7 B7 Employees of Intellectual Property Bureau, Ministry of Economic Affairs Printed by the consumer cooperative V. Description of the invention (12) A Fourier transform is performed on the scale, as shown in step 1 104 in FIG. The data forming the digital image contains the intensity of each pixel in the image and is spatially arranged for each color channel. Intensity refers to the brightness of one pixel. For example: a white pixel has greater intensity than a gray or black pixel. The same spatial area of different color channels usually has different pixel intensities. The Fourier transform transforms the pixel intensity of each segment from the spatial domain to the frequency domain. One of the values calculated by the Fourier transform is the average intensity of each segment of each color channel, that is, "D C". Calculated by the Fourier transform The other is the frequency and magnitude of the vertical and horizontal frequency vectors of each frequency element (element) in each section of each color channel. Although this embodiment uses a Fourier transform, other types of transforms can also be used to transform a signal from the spatial domain to the frequency domain, such as a Hadamard transfrom. In addition, the frequency conversion space can be subdivided using the ® pass-through wave canceller. Removing low and medium frequencies Second, in step 1 106 in Figure 4, the low and medium frequency terms of each transformation can be removed from further consideration. Alternatively, these terms may be weakened or this step may be omitted without departing from the scope of the invention. In this embodiment, each segment of each color channel (or a single channel in a black-and-white image) receives this removal within the useful range 1 of a particular embodiment of the present invention. Particle strength with significantly higher frequencies is easier to measure and distinguish from image information. Because of the tendency of substantial particles on photographic film (please read the precautions on the back before filling this page) Loading n · Βϋ I ^ · · ϋ «-a— nn · 1 ΜΗ! -V-α-.% This paper size Applicable to China National Standard (CNS) A4 specification (210 X 297 mm) -15- 496080 Α7 Β7 Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economy There are approximately equal frequencies, but at higher frequencies, there is less image content. In order to get the best difference between image detail and particle trajectory, it is best to focus on high frequencies. In this embodiment, we use a 6 X 1 6 section 6 times emphasis: other section sizes can be used without departing from the scope of the present invention. In order to filter out the low and mid-range frequencies, we can, for example: filter out the sum of all X-frequency and y-frequency absolute absolute chirps to frequencies less than or equal to 8. Alternatively, if the frequency vector is configured as shown in FIG. 18, a suitable radius around the DC can be selected, for example: 7777 = 8 frequencies within this radius can be removed. Other radii or filtering methods can be used 'without departing from the scope of the invention. Calculating the particle trajectory Secondly, in step 1 108 in Fig. 4, we can use the cross-correlation (c r s s-c o r r e 1 a t i ο η) formula to calculate the particle intensity of each section. In the case of a black-and-white digital image, because there is only one channel, the cross-correlation formula is not used. For color and black and white images, this calculation uses the frequency vector of the digital image from the Fourier transform after removing the low and medium frequency terms. To understand these calculations, an explanation of a digital image will be useful. The total information including a section is composed of information of each pixel of each color channel. For example ... In a typical three-color digital image, one (please read the precautions on the back before filling this page)
n an 一OJ0 11 n ϋ ammwm n ϋ I 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -16 - 496080 Α7 Β7 五、發明說明(14) (請先閱讀背面之注意事項再填寫本頁) 區段的總資訊爲來自紅、綠和藍色波道的資訊之組合。“ 色彩波道”這名詞爲方便而使用,且該影像波道不一定必 須·爲色彩波道。本發明可以被應用在一數位式影像的影像 波道,即使這樣的波道並不代表一種色彩。在這實施例中 ,此波道代表一照相影像的三種色彩波道。 紅色資訊,此處表示紅色頻率向量,且可以由公式1 代表之。在這專利中,粗體字代表一複數向量,其中’該 向量慣例上可以被考慮爲具有一實數與虛數分量’或替換 地一大小與相位分量。 R = R s + R II (1) 在公式1中,R代表記錄在整個紅色波道一區段中每 一像素組合的訊號。R s爲該區段來自影像訊號本身的紅色 資訊的部分。r n爲來自”雜訊”的部分,其中該雜訊可實 際地被認爲是顆粒軌跡。 同樣地,綠色資訊和藍色資訊爲: G = G s + G η (2) B = + (3) 經濟部智慧財產局員工消費合作社印製 假設只有r、G和B可以被測量出’ R u、G η和Β η 爲被計算估計的。同樣地,一黑白影像波道的雜訊成分被 計算估計。 因爲顆粒圖樣(pattern)在原有的軟片的每一色彩 層是不同的’吾人可以合理的假設在一數位式影像的區段 內的顆粒軌跡,實際上是雜訊,在整個色彩波道中是不相 關的。換言之,來自一區段中訊號的雜訊部分之強度丨直在 -17 - 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公爱) 經濟部智慧財產局員工消費合作社印製 496080 A7 ____B7__ __ 五、發明說明(15 ) 每一色彩波道中將傾向相異的。另一方面,來自一區段中 訊號的影像部分之強度値在每一色彩波道中將傾向相類似 的,且因此爲相關的。這樣的相關在邏輯上可以起因於自 然影像的亮度成分常常比色彩成分還強的原理。因此,色 彩之間的互相關可被用於R $和包含R的R η的比例混合之 預測値。很弱的互相關暗示訊號非常接近零,因此n an-OJ0 11 n ϋ ammwm n ϋ I This paper size is applicable to China National Standard (CNS) A4 (210 X 297 mm) -16-496080 Α7 Β7 V. Description of the invention (14) (Please read the note on the back first Please fill in this page again) The total information in the section) is a combination of information from the red, green and blue channels. The term “color channel” is used for convenience, and the image channel does not necessarily have to be a color channel. The present invention can be applied to an image channel of a digital image, even if such a channel does not represent a color. In this embodiment, this channel represents three color channels of a photographic image. Red information, which represents the red frequency vector, and can be represented by Equation 1. In this patent, bold text represents a complex vector, where 'the vector can be conventionally considered to have a real and imaginary component' or alternatively a size and phase component. R = R s + R II (1) In Equation 1, R represents the signal of each pixel combination recorded in one segment of the entire red channel. R s is the red information part of the segment from the image signal itself. r n is the part from "noise", where the noise can actually be considered as a particle trajectory. Similarly, the green information and blue information are: G = G s + G η (2) B = + (3) Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs. Assume that only r, G, and B can be measured. R u, G η and B η are calculated and estimated. Similarly, the noise component of a black and white image channel is calculated and estimated. Because the particle pattern is different in each color layer of the original film, 'I can reasonably assume that the particle trajectory in a section of a digital image is actually noise, which is not in the entire color channel. related. In other words, the strength of the noise part of the signal in a section 丨 straight at -17-This paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 public love) Printed by the Intellectual Property Bureau of the Ministry of Economic Affairs Consumer Cooperatives 496080 A7 ____B7__ __ 5. Description of the Invention (15) Each color channel will have a different tendency. On the other hand, the intensity of the image portion of the signal from a segment will tend to be similar in each color channel and is therefore related. Such a correlation can logically result from the principle that the brightness component of a natural image is often stronger than the color component. Therefore, the cross-correlation between colors can be used for the prediction of the ratio mixing of R $ and R η containing R. Weak cross-correlation implies that the signal is very close to zero, so
Rn = R = Rs + Rn (4) 另一方面,強的互相關暗示R S成分比R η佔優勢,因 而Rn = R = Rs + Rn (4) On the other hand, strong cross-correlation implies that the R S component is more dominant than R η, so
Rn«R = Rs+Rn (5) 因此,該雜訊的一合理的估計可變爲: 估計的R = F(互相關)· R ( 6 ) 在公式6中,當互相關趨近於自相關,亦即理想相關 時,F (互相關)趨近於0,且當互相關趨近於〇時’其 趨近於1。 因此,只要給定R,影像波道間的互相關可被用來得 出R “的一改良後的估計値。然而,當一強訊號存在時’該 估計値可能便得較不可靠。因此,互相關的第二個用處爲 得出這可靠度的估計値’且因而得出在與其他區段作平均 時,控制須將R η的估計値作多大的考慮的一“加權”。 來自雨種色彩波道的資訊間的互相關之測量’諸如: 紅色和綠色波道,因此可以代數表示之: i?G = (R· + RJ(G· + GJ = R,G· + R〆· + u 三R,G· ⑺ I--- I I J I I I - I I I I (請先閱讀背面之注意事項再填寫本頁)Rn «R = Rs + Rn (5) Therefore, a reasonable estimate of this noise can be changed to: Estimated R = F (cross-correlation) · R (6) In Equation 6, when the cross-correlation approaches to self Correlation, that is, ideal correlation, F (cross-correlation) approaches 0, and when cross-correlation approaches 0, it approaches 1. Therefore, as long as R is given, the cross-correlation between image channels can be used to arrive at a modified estimate of R ". However, when a strong signal exists, 'the estimate' may be less reliable. Therefore, The second use of cross-correlation is to get an estimate of this reliability 値 'and thus a “weighting” of how much control must be taken into account when averaging with other sectors. Measurement of cross-correlation between information of color channels such as: red and green channels, so it can be expressed algebraically: i? G = (R · + RJ (G · + GJ = R, G · + R〆 ·· + u Three R, G · ⑺ I --- IIJIII-IIII (Please read the precautions on the back before filling this page)
· I I I 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -18- 經濟部智慧財產局員工消費合作社印製 496080 A7 ___B7_ _ 五、發明說明(16 ) 在公式了中’ R G指示出有關來自紅色和綠色波道一 區段之資訊間的互相關。這是在已經移除低和中範圍頻率 之後,對每一非D C頻率項,一個頻率一個頻率地,由紅 色波道的區段頻率値以複數乘積乘上綠色波道的區段頻率 値而完成。例如:在紅色波道中,頻率位置1的頻率向量_ 上在綠色波道中頻率位置1的頻率向量。在一各別的計算 中,在紅色波道中頻率位置2的頻率向量乘上在綠色波道; 中頻率位置2的頻率向量,等等。應注意的是:其中有兩 向量的向量相乘,即”點乘積” ’因此’代表這兩個向量 包含向量方向的相關。這些個別計算的乘積被加總。 公式7也顯示:在數位式影像的大面積上’ R n G s、 G s R n和R n G n的値並無相關。它們的乘積有時候隨機地 爲正,有時候隨機地爲負,故對最後的總和具有較小的影 響。它們的平均値趨近於〇。最後’公式了顯不:一個頻 率一個頻率地,每一區段的總和、相關的紅色和綠色資訊 槪約等於相關的紅色和綠色影像一訊號値。相同地’ B R = BsRs 且 BG = BiiGs。 每一區段的色彩波道間的互相關,一個頻率位置一個 頻率位置地,可以使用公式8 - 1 〇來計算之。 BR = RB= Σ R · B 整個區段 (8 > BG = GB= Σ G · R 整個區段 (9 ) BG = GB= Σ G-B 整個區段 (10) 因此,對一特別的區段而言,在兩色彩波道間的互相關包 含在對應的頻率上’頻率域向量的點乘積之總和(在移除 ·---Μ----裝--------訂------;--- (請先閱讀背面之泛意事項再填寫本頁) 本纸張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -19- 496080 A7 ______B7___ 五、發明說明(17 ) 低和中範圍頻率項之後)。 (請先閱讀背面之注意事項再填寫本頁) 其次,顆粒強度可以對在每一色彩波道中的每一區段 來測量。顆粒強度爲從在原來的軟片中的顆粒所留下的軌 跡數目。具有大量的顆粒軌跡之區段有高顆粒強度。一區 段的紅色顆粒強度(R g s )可以由以下公式計算之: ^ = l|RQ| = SVRn *Ra (⑴ 亦即,一區段中的紅色顆粒強度等於它的紅色雜訊値丨 的絕對値之總和,依次等於紅色雜訊値平方的平方根之總 以下的估算,根據經驗硏究可以被使用: A · Rn = λΙ[Κ (12) 要注意的是:平方根是用來將偶然的強雜訊的衝擊減 至最低。這個估算也可以被用於其他色彩波道,或一黑白 影像的單一波道。 因爲R „ = R - R s,所以紅色波道的顆粒強度之公式 可以變爲如下: 你-Σλ/λλ: = (13) 經濟部智慧財產局員工消費合作社印製 注意:對R的一已知區段,當R s增大且R η變小時, “ R — R s ”的値趨近於〇 ,因此R g s趨近於〇。相反地 ,當當R u增大且R s變小時,“ R — R〆’的値變大,因 而R s s變大。 綠色和藍色波道的顆粒強度之公式可以同樣的方式寫 出: 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -20 - 496080· III This paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) -18- Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economy 496080 A7 ___B7_ _ V. Description of the invention (16) In the formula ' RG indicates cross-correlation between information from a segment of red and green channels. This is after removing the low and mid range frequencies, and for each non-DC frequency term, one frequency at a time, the segment frequency of the red channel is multiplied by the complex frequency product of the segment frequency of the green channel. carry out. For example: in the red channel, the frequency vector at the frequency position 1 on the frequency vector_ in the green channel. In a separate calculation, the frequency vector at frequency position 2 in the red channel is multiplied by the frequency vector at position 2 in the green channel, and so on. It should be noted that there are two vectors in which the vectors are multiplied, that is, the “point product”, so '' represents that these two vectors contain the correlation in the direction of the vector. The products of these individual calculations are summed. Equation 7 also shows that, over a large area of the digital image, there is no correlation between R R n G s, G s R n, and R n G n. Their product is sometimes randomly positive and sometimes randomly negative, so it has a small effect on the final sum. Their average chirp approaches zero. Finally, the formula is shown: one frequency and one frequency, the sum of each segment, the relevant red and green information 槪 is approximately equal to the signal of the relevant red and green images 値. Similarly, B R = BsRs and BG = BiiGs. The cross-correlation between the color channels of each zone, one frequency position and one frequency position, can be calculated using Equation 8-10. BR = RB = Σ R · B Entire Section (8 > BG = GB = Σ G · R Entire Section (9) BG = GB = Σ GB Entire Section (10) Therefore, for a particular section, In other words, the cross-correlation between the two color channels includes the sum of the dot products of the frequency domain vectors at the corresponding frequencies (after removal · · —M ---- 装 -------- Order- -----; --- (Please read the general Italian matter on the back before filling in this page) This paper size applies to China National Standard (CNS) A4 (210 X 297 mm) -19- 496080 A7 ______B7___ 5 (17) After the low and medium frequency terms) (Please read the notes on the back before filling this page) Second, the particle intensity can be measured for each segment in each color channel. Particle intensity Is the number of traces left from the particles in the original film. The segment with a large number of particle traces has high particle intensity. The red particle intensity (R gs) of a segment can be calculated by the following formula: ^ = l | RQ | = SVRn * Ra (⑴, that is, the intensity of the red particle in a section is equal to the sum of its absolute noise 値 丨, which in turn equals The sum of the square roots of the color noise and the square root can be estimated based on experience: A · Rn = λΙ [Κ (12) Note that the square root is used to reduce the impact of occasional strong noise to The lowest. This estimate can also be used for other color channels, or a single channel of a black and white image. Because R „= R-R s, the formula for the particle intensity of the red channel can be changed as follows: You -Σλ / λλ: = (13) Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs Note: For a known section of R, when R s increases and R η becomes small, 値 of “R — R s” approaches 〇 Therefore, R gs approaches 0. Conversely, when R u increases and R s becomes small, the R of “R — R〆 ′” becomes larger, and thus R ss becomes larger. The formula can be written in the same way: The paper size is applicable to the Chinese National Standard (CNS) A4 (210 X 297 mm) -20-496080
五、發明說明(18 )V. Description of the invention (18)
Gg^ = lVG*Ga =lVG*(G~Gs) (14) Bgs =Σλ/Β·(Β-Β5) (15) 經濟部智慧財產局員工消費合作社印製Gg ^ = lVG * Ga = lVG * (G ~ Gs) (14) Bgs = Σλ / Β · (Β-Β5) (15) Printed by the Consumer Cooperative of Intellectual Property Bureau of the Ministry of Economic Affairs
一類似的公式可以被使用在一黑白影像的單一波道15 每一色彩波道的影像訊號存在於另一色彩波道的影f象 訊號中。例如:紅色影像訊號(R s )中的某一些存在於綠 色影像訊號(G s )和藍色影像訊號(B s )中。綠色影像 訊號(G s )中的某一些存在於紅色影像訊號(R s )和^ 色影像訊號(B s )中,且藍色影像訊號(B s )中的某一 些存在於紅色影像訊號(R s )和綠色影像訊號(G s )中 〇 以下的常數(K )値,對辨識顆粒強度是有用的’如 下所述: K R。爲在綠色影像訊號中的紅色影像訊號中的數量。 K。R爲在紅色影像訊號中的綠色影像訊號中的數量° K R B爲在藍色影像訊號中的紅色影像訊號中的數量° K。B爲在藍色影像訊號中的綠色影像訊號中的數量° K B R爲在紅色影像訊號中的藍色影像訊號中的數量° K B。爲在綠色影像訊號中的藍色影像訊號中的數量° 根據經驗硏究,根據這些常數,輸入訊號和雜訊資1只 的頻率大小,可算出每一色彩的影像訊號的估計値。在— 實行中,紅色、綠色和藍色訊號爲: Gs-/2K〇s^^y2KCRR (16) (17) 「靖先閱讀背面之注意事項再填寫本頁) Μ---^--------^------:---. 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -21 - 496080 A7 五、發明說明(19 ) Bs^y2KBRR,/2K3CG (18) :在公式1 6 — 1 8中,R爲一區段的紅色頻率向量的 集合,G爲一區段的綠色頻率向量的集合,B爲一區段的 藍色頻率向量的集合。這些値已經從傅立葉變換而得來° 如前,低和中頻率已經被移除。 對紅色顆粒強度(R g s ) ’求解,則公式爲: (19) (請先閱讀背面之注意事項再填寫本頁) 對綠色和藍色顆粒強度(R g s ) 〜. [G - % 求解 公式爲 (20) Β^ = Σ ylB-[B-{y2KBRR + y2 Ksgq) (21) 另-組的估算可以藉由使用先前所計算的互 相關而被 用對常數(K )値(見公式 可由如下而估算之: 10A similar formula can be used for a single channel of a black and white image. The image signal of each color channel exists in the image signal of another color channel. For example, some of the red image signals (R s) exist in the green image signals (G s) and the blue image signals (B s). Some of the green image signals (G s) exist in the red image signal (R s) and ^ color image signals (B s), and some of the blue image signals (B s) exist in the red image signal ( R s) and a constant (K) 〇 of 0 or less in the green image signal (G s) are useful for distinguishing the intensity of particles' as described below: KR. Is the number of red image signals in the green image signal. K. R is the number of red video signals in the green video signal ° K R B is the number of red video signals in the blue video signal ° K B is the number of green image signals in the blue image signal ° K B R is the number of blue image signals in the red image signal ° K B It is the number of blue image signals in the green image signal. According to empirical research, based on these constants, the frequency of the input signal and noise data can be used to calculate the estimated image signal of each color. In-implementation, the red, green and blue signals are: Gs- / 2K〇s ^^ y2KCRR (16) (17) "Jing first read the precautions on the back before filling out this page) Μ --- ^ --- ----- ^ ------: ---. This paper size applies to China National Standard (CNS) A4 (210 X 297 mm) -21-496080 A7 V. Description of the invention (19) Bs ^ y2KBRR, / 2K3CG (18): In the formulas 16—18, R is a set of red frequency vectors of a section, G is a set of green frequency vectors of a section, and B is a blue frequency of a section The set of vectors. These chirps have been obtained from the Fourier transform. As before, the low and medium frequencies have been removed. Solving for the red particle intensity (R gs) ', the formula is: (19) (Please read the back Please fill in this page again for the intensity of green and blue particles (R gs) ~. The group estimate can be used for the constant (K) 値 by using the previously calculated cross-correlation (see the formula can be estimated as follows: 10
求解。例如:K 裝--------訂------.---邊 %·Solve. For example: K Pack -------- Order ------.--- Side% ·
RG.RBRG.RB
KK
RG 經濟部智慧財產局員工消費合作社印製Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs
GsGs
BGBG
jGB.GR V BRjGB.GR V BR
BG-GB (22)BG-GB (22)
以相同方法求解,其他的K値爲: KGR KBR G5 Jr (23) GR Ίδ (24) BG :^GR (25) 本紙張尺度適用中國國家標準(CNS)A4規格(210 χ 297公釐) 496080 A7 B7 五、發明說明(20 ) KGB: GR z -— BR (26) KBC. BR —__ __ GR (27) 吾人使用已得出之値並使用移除低和中頻率(在這些 公式中的R、G和B )後,從傅立葉變換所得的頻率向量 ,則可計算出每一色彩,在每一區段中的顆粒強度。 例如:本發明透過如下系列的方程式估算紅色顆粒強 度且估算: _ Rgs =z|R„| = Ia/r~r7 (ii) (12) ----···---J----·-裝—— (請先閱讀背面之注意事項再填寫本頁)Solve in the same way, the other K 値 are: KGR KBR G5 Jr (23) GR Ίδ (24) BG: ^ GR (25) This paper size applies the Chinese National Standard (CNS) A4 specification (210 χ 297 mm) 496080 A7 B7 V. Description of the invention (20) KGB: GR z -— BR (26) KBC. BR —__ __ GR (27) We use the obtained 得出 and use to remove the low and medium frequencies (in these formulas R, G, and B), from the frequency vector obtained from the Fourier transform, the intensity of the particles in each color in each segment can be calculated. For example, the present invention estimates the intensity of the red particles through the following series of equations and estimates: _ Rgs = z | R „| = Ia / r ~ r7 (ii) (12) ---- ·· --- J --- -· -Pack—— (Please read the precautions on the back before filling this page)
RgsRn=R-R, (13) ^ = Za/R*(R-rs) (13),RgsRn = R-R, (13) ^ = Za / R * (R-rs) (13),
B (16) 經濟部智慧財產局員工消費合作社印製B (16) Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs
KK
RG K0 BRΊδ GRJr 本紙張尺度適用中國國家標準(CNS)A4規格(21〇 χ 297公S ) -23- (19) (22) (23) 496080 A7 B7 五、發明說明(21 ) R.RG K0 BRΊδ GRJr This paper size is applicable to the Chinese National Standard (CNS) A4 specification (21〇 χ 297 S) -23- (19) (22) (23) 496080 A7 B7 V. Description of the invention (21)
/2BG /2 BR/ 2BG / 2 BR
M = SR*BM = SR * B
5G = ZG -B5G = ZG -B
GR = Σ G · R (28) (8) (10) (9) (請先閱讀背面之注意事項再填寫本頁)GR = Σ G · R (28) (8) (10) (9) (Please read the notes on the back before filling this page)
RnRn
SjR·SjR ·
IRB (29) 估算綠色顆粒強度和藍色顆粒強度的公式爲The formula for IRB (29) to estimate the intensity of green particles and blue particles is
Ggs =Ggs =
GG
ZRBZRB
IR B (30) ^ = Σ· B. 經濟部智慧財產局員工消費合作社印制衣 (31) 在每一情形下,對在空間頻率域中的每一頻率加總。 對黑白影像而言,每一區段的顆粒強度的估算是由將 每一頻率向量的點乘積之平方根和在移除低和中頻率後的 本身相加而得的。 計算區段加權 再根據一組計算估算顆粒強度後,第二組公式可被用 來決定每一色彩波道的每一區段的一”加權”値,如圖4 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) .24 - 496080 A7 _____—__B7 _五、發明說明(22 ) 中的步驟1 1 1 〇所示。再次地,移除低和中頻率後,可 使用從傅立葉變換得來的頻率向量。此加權値提供所計算 出’的穎粒強度(雜訊)實際上爲顆粒軌跡的測量而不是影 像細自P的測星的~信賴水平(c ο n f丨d e n c e 1 e v e 1 )。 爲了計算這加權値,可以使用公式8 - 1 0。替換地 ’也可使用如下求每一區段的頻率向量値的總和所用的公 式。這堅公式(其爲自相關的度量)將在一區段中,針對 所有的頻率位置,將頻率向量的點乘積與它們本身相加。 吾人可找到每一色彩波道的各別的總和。 (請先閱讀背面之注音〉事項再填寫本頁)IR B (30) ^ = Σ · B. Printing of clothing by employees' cooperatives in the Intellectual Property Bureau of the Ministry of Economic Affairs (31) In each case, sum each frequency in the space frequency domain. For black and white images, the estimation of the particle strength of each segment is obtained by adding the square root of the dot product of each frequency vector and itself after removing the low and medium frequencies. After calculating the segment weight and then estimating the particle intensity based on a set of calculations, the second set of formulas can be used to determine a "weight" for each segment of each color channel, as shown in Figure 4. This paper scale applies Chinese national standards (CNS) A4 specification (210 X 297 mm) .24-496080 A7 _____—__ B7 _V. Step 1 1 1 0 in the description of the invention (22). Once again, after removing the low and medium frequencies, a frequency vector derived from the Fourier transform can be used. This weighted 値 provides that the calculated ‘glue intensity (noise) is actually a measure of the particle trajectory rather than the level of trust (c ο n f d e n c e 1 e v e 1) of the star from which the image is fine. To calculate this weighted 値, Equation 8-1 0 can be used. Alternatively, it is also possible to use the formula used to find the total of the frequency vector 値 of each section as follows. This formula (which is a measure of autocorrelation) will add the dot product of the frequency vectors to all frequency positions in themselves in one section. We can find the individual sum of each color channel. (Please read the Zhuyin on the back> Matters before filling out this page)
RR=I R · R GG=E G · G ΒΒ= Σ B · B (32)(33)(34) 在每一情形下,R、G和B代表移除低和中範圍頻率後一 區段的傅立葉變換。 使用從公式8 - 1 0和.3 2 - 3 4所得的結果,可以 根據如下的公式,計算出每一區段的一加權値: 裝--------訂--- Φ 經濟部智慧財產局員工消費合作社印製 加權=1 . 0RR = IR • R GG = EG Transform. Using the results obtained from formulas 8-1 0 and .3 2-3 4, one weighted 値 for each sector can be calculated according to the following formula: 装 -------- Order --- Φ Economy Ministry of Intellectual Property Bureau employee consumer cooperatives print weighting = 1.0
\RG— RB—GB (35) \RRvGG 七 公式3 5是有用的,因爲來自一區段的影像部分的不 同色彩之相同頻率値在各色彩波道之間’比來自顆粒軌跡 者還要相似。如在這應用先前所解釋者’從每一色彩波道 所記錄的顆粒軌跡在一數位式影像中具有各別的圖樣’且 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -25- 496080 經濟部智慧財產局員工消費合作社印製 A7 B7 五、發明說明(23) 因此具有不同組的頻率値。另一方面,影像在所有的色彩 波道中傾向相似’且因而在每一色彩波道中具有相似的頻 率値。 例如:在完全由影像訊號所組成的一區段中’在加權 値方程式的分子中的互相關組合値槪約等於在分母中自相 關的組合値。那區段的加權値將趨近於1 . 〇 一1的平方 根値,即最後爲0。 f 然而,完全由顆粒軌跡訊號所組成的一區段則具有趨 近於1 · 0 - 0的加權値,即最後爲1 · 〇。這是因爲.在 分子中的互相關的組合値會遠小於在分母中自相關的組合 値。 因而,一區段的內容由越多的影像訊號所組成’則加 權値將越趨近於0。一區段的內容由越多的顆粒軌跡(雜 訊)所組成,則加權値將越趨近於1。 對黑白影像而言,此加權通常被設爲一固定値,如: 以加權、D C和顆粒強度作圖 執行這些計算之後,每一區段的加權値可對整個影像 ,被畫在三個圖形上,每一色彩波道畫一個圖形,如圖4 中步驟1 1 1 2中所示。明白地說,一個圖形是爲影像的 每一色彩波道而畫出的,而不是每個區段一個圖形。在此 數位式影像中的每一區段根據它的顆粒強度和D C (其爲 經由傅立葉變換所得該區段的平均強度値)而表示於圖形 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公爱) (請先閱讀背面之注意事項再填寫本頁)\ RG— RB—GB (35) \ RRvGG Seven formulas 3 5 are useful because the same frequencies of different colors from the image part of a section 値 between color channels are more similar than those from particle trajectories . As explained earlier in this application, 'the particle trajectory recorded from each color channel has a separate pattern in a digital image' and this paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 cm) (%) -25- 496080 A7 B7 printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs. 5. Description of the invention (23) Therefore, it has different groups of frequencies. On the other hand, images tend to be similar 'in all color channels and therefore have similar frequencies in each color channel. For example, in a section composed entirely of image signals, the cross-correlation combination 'in the numerator of the weighting' equation is approximately equal to the self-correlation combination in the denominator '. The weighted 値 of that sector will approach the square root 値 of 1.00-1, which is finally 0. f However, a segment composed entirely of particle trajectory signals has a weighted 値 approaching 1 · 0-0, which is finally 1 · 0. This is because the cross-correlation combination 値 in the numerator is much smaller than the auto-correlation combination 分 in the denominator. Therefore, the content of a sector is composed of more video signals', the weighting will be closer to zero. The content of a segment is composed of more particle trajectories (noise), and the weighted chirp will approach 1. For black and white images, this weighting is usually set to a fixed frame, such as: After performing these calculations with weighting, DC, and particle intensity mapping, the weighting frame of each segment can be drawn on the entire image in three graphics For each color channel, draw a pattern, as shown in step 1 1 1 2 in Figure 4. To be clear, a graphic is drawn for each color channel of the image, not a graphic per sector. Each segment in this digital image is represented on the graph according to its particle intensity and DC (which is the average intensity of the segment obtained through Fourier transform 値). This paper's dimensions apply the Chinese National Standard (CNS) A4 specification ( 210 X 297 public love) (Please read the notes on the back before filling this page)
· n «I I— el I n n 一口t · ϋ_ι I— n an an m n I 496080 A7 B7___ 五、發明說明(24 ) (請先閱讀背面之注意事項再填寫本頁) 上。此外,對每一區段所畫出的點被分配到那區段的加權 値。其結果爲一散佈圖,其中X軸代表平均強度,y軸代 表:顆粒強度,且在散佈圖上的一特殊點的大小包含該加權 値。對一黑白影像而言,可產生單一的散佈圖,與上述的 顆粒強度和加權用的修改後公式一致。 圖6例示對於一數位式影像中的一區段,紅色波道的 資訊之圖。y軸代表顆粒強度。X軸代表平均強度。點 1 5 0 2爲此區段的紅色顆粒強度和紅色平均強度的點。 點1 5 0 2的放大試圖1 5 0 4例示該區段的加權値( W V )。 以相同的方式,此區段的一點分別在一綠色散佈圖上 的綠色波道被畫出,且在一藍色散佈圖上的藍色波道被畫 出。畫出一影像的區段之資料後’其結果爲三個散佈圖’ 每一色彩波道一個。圖了例示一數位式影像區段’經一色 彩波道,例如:紅色波道’所記錄一散佈圖。類似的散佈 圖將顯示綠色和藍色’或一黑白影像的單一波道的結果。 經濟部智慧財產局員工消費合作社印製 畫出估計的顆粒強度之線 一旦產生散佈圖之後’曲線擬合技術可被用來’對每 一色彩波道(或一黑白影像的單一波道)’畫出通過霧狀 點的線,如圖4中的步驟1 1 1 4所示。被分配到每一點 位置的加權値可能影響這些線的精確位置’較大的加權値 將線拉向它們的方向。圖8例示經一色彩波道’例如:紅 色波道的一散佈圖所畫出的線1 7 〇 2 °此曲線代表顆粒 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) 〇7 - ""~ 496080 經濟部智慧財產局員工消費合作社印製 A7 B7 五、發明說明(25) 強度對強度所畫的圖。雖然,吾人將在下文中討論曲線擬 合法,但是,任何曲線擬合技術均可被使用,而不悖離本 發明的範圍。此外,不同的技術可被用在不同的色彩波道 上。 這些線因而可被用來決定對應在一特別的色彩波道中 任何像素強度値的該特別的色彩波道中數位式影像中該像 素的估計顆粒強度的有用度量。這度量可被用來以各種不 同的方法,使該數位式影像看起來更合意。 影像資訊存在於所有三個波道的假設大部分時間均爲 真。但是,一數位式影像的明亮色彩區域可能包含具強影 像和色彩之間較小相關的區域,因此導致離平均範圍很遠 的偶然強點。爲了達成估計的顆粒強度之曲線擬合’此曲 線擬合技術可以將偶然強點考慮進來,而不讓此曲線被它 們所扭曲。例如:它可使用點的一中位數(med i an )平均, 或可找出最相合的一原型(prototype )。它也可以做出最 相合的一第一原型,抑制位於該原型線兩倍高的強點,然 後在修改後的資料上作出一第二原型線。這樣的強點,極 可能它們代表影像細節,而不是雜訊,因此’抑制這些點 可能產生一較精確曲線。這樣的最合適操作可以被執行許 多次,直到不拋棄任何點爲止。在這實施例中’吾人使用 一中位數平均技術。 圖9顯示爲一數位式影像的紅色資訊所畫的線 1 8 0 2。線1 8 0 2可以是相同的曲線或和與該散佈圖 關聯而產生的曲線成正比之曲線,諸如:圖8中的線 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -28- (請先閱讀背面之注意事項再填寫本頁)· N «I I— el I n n 口 口 · ϋ_ι I— n an an m n I 496080 A7 B7___ 5. Description of the invention (24) (Please read the precautions on the back before filling this page). In addition, the points drawn for each segment are assigned to the weighted 値 of that segment. The result is a scatter plot, where the X-axis represents the average intensity, and the y-axis represents: particle intensity, and the size of a particular point on the scatter graph contains the weighted 値. For a black and white image, a single scatter diagram can be generated, which is consistent with the modified formula for particle intensity and weighting described above. Fig. 6 illustrates the information of the red channel for a section in a digital image. The y-axis represents particle strength. The X-axis represents the average intensity. The point 1 5 0 2 is the point of the intensity of the red particles and the average intensity of the red in this section. An enlargement of the point 15 0 2 attempts to exemplify the weighted chirp (W V) of the segment. In the same way, one point of this segment is drawn on the green channel on a green scatter diagram, and the blue channel on a blue scatter chart is drawn. After drawing the data of an image segment, the result is three scatter plots, one for each color channel. The figure illustrates a scattergram recorded by a digital image section ′ through a color channel, such as a red channel. A similar scatter plot will show the results for a single channel of green and blue 'or a black and white image. Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs to draw a line of estimated particle intensity. Once a scatter plot is generated, a 'curve fitting technique can be used' for each color channel (or a single channel of a black and white image) Draw a line through the fog point, as shown in step 1 1 1 4 in FIG. 4. The weights assigned to each point position may affect the exact position of these lines' larger weights, pulling the lines towards them. Figure 8 illustrates a line drawn through a color channel. For example, a line drawn by a scatter plot of the red channel 17.0 2 ° This curve represents particles. The paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) ) 〇7-" " ~ 496080 Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs A7 B7 V. Description of the invention (25) Drawing of intensity versus intensity. Although we will discuss curve fitting in the following, any curve fitting technique can be used without departing from the scope of the invention. In addition, different techniques can be used for different color channels. These lines can thus be used to determine a useful measure of the estimated particle intensity of a pixel in a digital image of a particular color channel for any pixel intensity in that particular color channel. This metric can be used in various ways to make the digital image look more desirable. The assumption that image information exists in all three channels is true most of the time. However, the bright color areas of a digital image may contain areas with a small correlation between strong images and colors, thus resulting in occasional strong points far from the average range. In order to achieve a curve fit of the estimated particle strength, this curve fitting technique can take occasional strong points into account without letting the curve be distorted by them. For example, it can use a median (med i an) average of the points, or it can find a prototype that most closely matches (prototype). It can also make the most consistent first prototype, suppress the strong point that is twice as high as the prototype line, and then make a second prototype line on the revised data. Such strong points, most likely they represent image details, not noise, so 'suppressing these points may produce a more accurate curve. This most appropriate operation can be performed many times until no point is discarded. In this embodiment, we use a median averaging technique. Figure 9 shows the line 1 8 0 2 drawn for the red information of a digital image. Line 1 8 0 2 can be the same curve or a curve proportional to the curve associated with the scatter diagram, such as: The paper size of the paper in Figure 8 applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) PCT) -28- (Please read the notes on the back before filling this page)
· taMm i In n 11 ϋ— ϋ 一口、I H ϋ tmmmmm n n ϋ* n I 496080 A7 B7 五、發明說明(26 ) 1 7 0 2。注意:在此X軸代表像素強度値’而非在前面 使用過的整個區段之平均像素強度。對具有點1 8 0 4所 示:強度的一像素,估計的顆粒強度在點1 8 0 6 °換言之 ,雖然曲線1 7 0 2是根據平均強度’從一散佈圖得來的 ,但是,若設定一像素的強度,曲線1 8 0 2 (該曲線可 以與曲線1 7 0 2相同)可被用來提供該像素的顆粒強度 之估計値。 測量相關於頻率的顆粒強度 在完成顆粒軌跡的度量之後,於人可以對每一色彩波 道(或一黑白影像的單一波道)測量相關於頻率的顆粒強 度,如圖3中的步驟1 〇 0 4所示。頻率內容代表影像細 節。具有一較高度影像細節的區段在所有頻率中,將具有 較多內容。 例如:考慮由天空所包圍且豎立於一片草地上的一顆 樹的影像。影像中的草地部分因其許多葉片和色彩的陰影 ,具有在所有頻率中大部分的高能量,亦即’一個大頻率 內容。很多高能量可能具有高頻率內容。另一方面’天空 部分,其缺乏影像細節,故具有很少產生此影像的高頻。 因此,任何被度量的高頻內容因而可以假設爲主要從顆粒 細節產生。 爲了測量相關於頻率的顆粒強度,吾人執行在圖1 〇 中所例示的步驟。在步驟2 0 0 2中,此數位式影像第二 次被區段分割,且在步驟2 0 0 4中,執行第二次傅立葉 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) (請先閱讀背面之注意事項再填寫本頁) 裝 經濟部智慧財產局員工消費合作社印製 -29- 經濟部智慧財產局員工消費合作社印製 496080 A7 B7 ___ 五、發明說明(27 ) 變換,儘管,也可以替換地使用第一次區段分割的結果和 第一次傅立葉變換。若作新的計算,則可使用不同的區段 大小、形狀和重疊選項(這樣的選項在前面已說明),如 同也可使用不同型式的變換(再次地,這樣的選項在前面 已說明)。在步驟2 0 0 6中,爲了辨識相關於頻率的估 計顆粒強度,吾人產生具有被擬合的曲線之散佈圖。其次 ,在步驟2 0 0 8中,吾人根據顆粒強度對強度函數,其 中,該函數爲對每一區段與頻率値相關而得來,畫出由估 計的顆粒強度所分割的大小。在步驟2 0 1 0中,畫出辨 識與頻率相關的估計顆粒軌跡之線。 區段分割 當測量與頻率相關的顆粒強度時’在步驟2 0 1 2中 ,數位式影像可以使用前述的圖4中步驟1 1 0 2的任何 技術和選項而被分割成區段。 傅立葉變換 在步驟2 0 0 4中,吾人使用前述與圖4中步驟 1 1 0 4相關聯的任何技術和選項’對在每一色彩波道( 或在一黑白影像的單一波道)中的每一區段的像素強度値 上執行一第二傅立葉變換。然而,低和中頻率項的範圍並 未被移除,如先前圖4的步驟1 1 0 6中所示。 辨識相關於頻率的估計題粒軌1 本紙張尺度適用令國國家標準(CNS)A4規格(210 X 297公爱) -30 - (請先閱讀背面之注意事項再填寫本頁) 裝 . 496080 A7 B7 五、發明說明(28) (請先閱讀背面之注音?事項再填寫本頁) 使用從第二次傅立葉變換得來的D C (平均強度)値 ,吾人可以識別在每一色彩波道(或在一黑白影像的單一 波道)中,每一區段與頻率相關的顆粒強度,如步驟 2 0 0 6中所示者。先前經由每一色彩波道(或黑白影像 的單一波道)的散佈圖所畫的線(諸如:線1 7 〇 2和 1802),可以被應用在這個目的上。對每一區段,從 第二次傅立葉變換得來的D C (平均強度)値可以被用來 / 對曲線的強度値,在X軸上定位一點,且因此,找出在y 軸上估計的顆粒強度。 例如:若一區段具有一與圖9中在點1 8 0 8上的値 相當的D C値,它的估計顆粒強度將是在線1 8 0 2的點 1 8 1 0上的値。因此’在一區段的每一頻率被認爲具有 從該區段的平均強度得來的一估計的顆粒強度。 書M/G S對頻率之圖 經濟部智慧財產局員工消費合作社印製 在圖1 0的步驟2 0 0 8中,在識別相關於頻率的顆 粒強度之後,在每一色彩波道(或在一黑白影像的單一波 道)中,每一區段在每一頻率的大小可以被在每一色彩波. 道(或在一黑白影像的單一波道)中’每一區段的估計顆 粒強度所分割。每一頻率的大小是由在步驟2 0 1 〇中的 傅立葉變換所決定的。 於是可以產生新的散佈圖,每一色彩波道一個圖形( 或一黑白影像一個圖形)。對所討論的波道而言,散佈圖 的y軸代表每一區段在每一頻道上大小的値’其中該値是 本纸張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -31 - 496080 A7 ____ B7 _ 五、發明說明(29 ) (請先閱讀背面之注意事項再填寫本頁) 由那區度的估計顆粒強度所分割。對所討論的波道而言, 散佈圖的X軸代表呈現於數位式影像中的頻率値。因此, 對在數位式影像中的每一區段之每一頻率晝出在散佈圖上 的一點。在上文的舉例中,2 8 9個頻率被呈現,因此, 該例中新散佈圖將具有2 8 9 S個被畫出的點,其中S代 表區段的數目。在這散佈圖中,並無被附加至各點上的加 權値,除了被畫在同一位置上的點之數目。明確地說,每 一點具有加權1 ,但是,若第二點被畫在同一位置上,則 所存在的點可被分配一加權2。替換地,一個點的表列可 以被簡單的維護且包含兩個點。 圖1 1例示一色彩波道的這樣一散佈圖。較低的“霧 狀點”代表那些已經被以上的加權所隔離的區段’而具有 較高的機率代表主要的顆粒強度而非影像訊號。那些從霧 狀帶被移除的點具有較高的機率代表影像訊號。 書線來辨識相關於頻率的顆粒軌跡 經濟部智慧財產局員工消費合作社印製 其次,使用任何曲線擬合演算法,畫出經由在每一散 佈圖中的較低霧狀帶中的所有點之線,如圖1 〇中的步驟 2 0 1 0所示。在這實施例中,吾人使用一中位數平均型 式的演算法。當畫出這線時’在這霧狀帶外面’錯誤的較 高圖點可以被忽略,因而不會使結果失真。爲了做到這樣 ,可以使用如下的對每一點之”加權”公式,其中y是在 y軸上一點的値: 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐〉 -32- 經濟部智慧財產局員工消費合作社印製 496080 A7 B7_ 五、發明說明(3〇 ) 任何其他的適合的加權可以被使用,而不悻離本發明 的範圍。圖1 2顯示經由一色彩波道的霧狀帶所畫出的一 線Έ 2 0 2 ,作爲一例子。注意:當畫此線時,錯誤的高 點2 2 0 4和2 2 0 6被忽略。同時,僅僅使用中和高頻 率圖點來決定此線可能是有利的。以此方式,強的低頻內 容將不會不利地影響被估計的D C。 計算校正陣列 然後可以計算每一色彩波道的一校正陣列,如圖3中 的步驟1 0 0 6所示。 找出由圖形所估計的D C項 爲了計算校正陣列,吾人使用經由在步驟2 0 1 0 ( 圖1 0 )中所畫出的散佈圖來計算在每一色彩波道(或在 一黑白影像的單一波道)中的一估計D C ( E D C )。例 如:圖1 2顯示由在線2 2 0 2上的點2 2 0 8所指示的 E D C値。然後,在這例子中,這三個散佈圖提供紅色、 綠色和藍色的一估計的D C値。 計算陣列 決定E D C値之後,使用公式3 7,計算每一色彩波 道(或在一黑白影像的單一波道)的一陣列,其中〔x〕 涵蓋像素強度値的範圍: 陣列U ]=(复·顆粒強度[,]) (3 了) 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) - 33: (請先閱讀背面之注意事項再填寫本頁) · in a—· n I «ϋ ai« i·— J·· n «ϋ I n n _ 片Ύ 言 4 496080 A7 B7 五、發明說明(31) (請先閱讀背面之沒意事項再填寫本頁) 使用從步驟2 0 1 0中計算的曲線所得的估計的D C 値和從步驟1 1 1 4中計算的曲線所得的每一強度的估§十 顆粒強度,對每一波道計算一陣列。 計筧校正陣列 陣列〔X〕値可以被用來計算每一波道(或單一黑白 波道)的一校正陣列。吾人可以使用公式3 8來計算這陣 列。 相關[〇]=〇 相關[小έΣ(陣列陣列⑽ (38) /7=»! 2 將相關於強度的顆粒軌跡歸一化 經濟部智慧財產局員工消費合作社印製 計算校正陣列之後,數位式影像的顆粒軌跡就可以相 關於強度地被歸一化,如圖3中的步驟1 0 0 8所示。經 校正陣列而得的値可以被應用在此數位式影像的資料,而 達成歸一化。在數位式影像中的每一像素強度値” X ”可 以被每一色彩波道(或一黑白影像的單一波道),在校正 〔X〕中的値所取代。這程序使得在數位式影像中的顆粒 軌跡在每一波道的灰階中,具有槪約相同的可見衝擊( V 1 s u a 1 1 mp a c t ),這使得該數位式影像看起來更爲合意。例 如:當其被以一對比係數(gamma )約2,以一正片顯示在一 典型的顯示器上時,更爲合意。 圖1 3爲來自一被歸一化的、具有在畫在空間上像素 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公爱):34 - ' 經濟部智慧財產局員工消費合作社印製 496080 A7 B7____ 五、發明說明(32 ) 強度中的改變之數位式影像的像素的線之顆粒軌跡的圖形 表示。應注意的是:在整個圖形上的像素強度的範圍是平 均的,且並不會像在一未歸一化的數位式影像中像素強度 一樣地劇烈地波動。在圖1 3中,點2 3 0 2上的波動大 約與點2 3 0 4上的相同。 將相關於頻率的顆粒軌跡歸一化 數位式影像的顆粒軌跡現在可以與頻率相關的被歸一 化,如圖3中的步驟1 0 1 0所示。應用校正陣列之後, 吾人可以執行另一個頻率變換,而將已被強度上歸一化的 影像之每一波道(或該單一波道)從時間域變換乘頻率域 。替換地,校正陣列可被轉換成一頻率域的函數,用來應 用在代表此影像的原有空間頻率域。在這實施例中,校正 陣列被應用在時間域且在每一波道上執行轉換,以實行將 要說明的頻率歸一化程序。 圖1 4爲以圖形表示大小對在軟體已經跨強度地歸一 化顆粒軌跡之後的一數位式影像中的頻率的一 2 D圖。在 此平均(DC)値2402爲1 · 0。 爲了跨頻率而將顆粒軌跡歸一化,此軟體將每一區段 中的每一頻率的大小(Μ )乘上一因數(F )。此因數指 示出所畫出經由散佈圖的線2 4 0 2偏離1 . 0的程度。 偏差因數(F ) = (39) 本紙張尺度適用中國國家標準(CMS)A4規格(210 X 297公釐) .35 - (請先閱讀背面之注意事項再填寫本頁) 裝--------訂------·—· 496080 A7 B7 五、發明說明(33) (請先閱讀背面之注意事項再填寫本頁) 一偏差因數的實際圖將顯示一頻率升高曲線,其中, 該曲線爲圖1 4中的線2 4 0 2的倒置。在此相位並未受 影響且被保留的。 圖1 5以圖形表示頻率大小對在軟體已乘上偏差因數 之後的一數位式影像中的頻率的一 2 D圖。以慣例方法’ 經由這圖所畫的一線2 5 0 2顯示:顆粒軌跡的大小大槪 與高和低頻的面積相等。它們不會像在一未歸一化的數位 式影像一樣地劇烈地波動。 抑制顆粒軌跡 計算校正陣列之後,由校正陣列所得的被歸一化的値 可以被用來抑制數位式影像中的顆粒軌跡,如圖3中的步 驟1 0 1 2所示。也可使用其他顆粒抑制方法’而不悖離 本發明的教導,諸如以下所述的第二種方法。 在根據如上的教導之數位式影像歸一化中,在每一色 彩波道中,每一區段和在每一頻率中的穎粒強度已經被調 整以致大小分量爲1 · 0 ,如圖1 5中所示。相位不受影 經濟部智慧財產局員工消費合作社印製 響。因此,一分量的大小是一有力指示器,指示:從在 1 . 0範圍內的影像細節中隔離顆粒軌跡(雜訊)是一分 水線數量。大於1 · 0的大小較有可能是從一影像細節所 產生的,小於1 · 0的大小較有可能是從一顆粒軌跡所產 生的。因此,爲了抑制在一數位式影像中的顆粒軌跡,最 好抑制大小接近和比1 · 〇小的成分。另一方面’比 1 . 0大的成分可以被認爲代表影像細節,而不應被抑制 -36- 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) 496080 經濟部智慧財產局員工消費合作社印製 Α7 Β7 五、發明說明(34) 的。 這與未被歸一化的數位式影像相反’其中’在穎粒軌 跡(雜訊)與影像之間的分水線數量隨著強度和頻率而變 ,使得幾乎不可能精確地將顆粒軌跡與影像隔離。如先前 所述,影像可以被歸一化’或者用強度和頻率描述雜訊的 分佈的資料可以被包含於一未被歸一化的影像中’以引導 顆粒的減少,而不悖離本發明的範圍。 爲了抑制顆粒軌跡,如下的方法可以被應用在被歸一 化的數位式影像區段中的每一頻率之大小(M): F 〇 r在歸一化的數位式影像中的每一區段 維持原有的相位資訊 F 〇 r在該區段(包括D C項)中的每一Μ I f ( Μ < = 0 · 5 ) Μ = Ο else 抑制因數=1 .〇一 e 一( Μ - 0 . 5 ) Μ = Μ *抑制因數 這公式藉由將大小爲或比Ο · 5小的被歸一化之頻率 由0替代而抑制其値。因爲有些穎粒軌跡具有稍微比 〇.5大的被歸一化的大小,所以這公式也可使用一抑制 因數來將它們的値用較小的値替代。此抑制因數對具有一 較高被歸一化的大小之頻率値減小情形較緩和’因爲’這 些値趨向於對應影像細節。 以上例子具有一數位式影像,其中該數位式影像由爲 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -37 - -----ί---:----_裝 i — (請先閱讀背面之注意事項再填寫本頁) . 經濟部智慧財產局員工消費合作社印製 496080 A7 ____B7__ 五、發明說明(35 ) 頻率和強度的函數之已知的顆粒強度對每一區段分割而將 之歸一化,然後應用以一固定爲1 . 〇的分水線的減少量 。一相當的方法也可將具有等於爲頻率和強度的函數之期 望的雜訊的一分水線點的顆粒減少量,直接應用至未被歸 一化的成分上。爲了進行這方法,吾人可以使用上列的穎 粒抑制方法,而在抑制因數中,以Μ ’代替Μ ’其中 ,__Μ 所期望的顆粒軌跡 有利地使用來自校正陣列的被歸一化値,則其他抑制 穎粒軌跡的方法都是可能的。例如:諸如PHOTOSHOP的影 像改良軟體可以依人工方式,與該歸一化後的値合用,雖然 這樣的軟體需要廣泛的使用者的回饋,以找出顆粒軌跡的 適當抑制水平。 現在說明另一種,自動地或在使用者的選擇下,可以 抑制一影像的某些或實質上所有顆粒之方法。 圖1 6例示說明從一數位式影像中移除諸如:顆粒軌 跡的高頻雜訊的一方法之流程圖。在圖1 6中所晝出的方 法可被應用在數位式影像的一或多個波道或,若該數位式 影像只有一波道,應用在數位式影像的單一波道上。對一 具有多個波道的影像而言,圖1 6所述的方法可被應用在 總波道中的子集合上,而不悖離本發明的範圍。同樣地, 圖1 6所述的方法可被應用在一多波道數位式影像中的多 個波道的每一個上,而不悖離本發明的範圍。 在步驟1 6 1 〇中’一特別的影像波道的雜訊被歸一 化。此雜訊可使用上述的方法或任何其他方法,而被歸一 ---,---»(*----^裝—— (請先閱讀背面之注意事項再填寫本頁) · 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公爱) -38- 496080 A7 B7___ 五、發明說明(36 ) 化。替換地,雖然在進行圖1 6所述的方法的其餘部分之 前,歸一化此雜訊較好,但是,此方法也可被應用在一未 經歸一化的數位式影像。 圖1 6所述的方法可被應用在曾具有諸如:顆粒軌跡 的高頻雜訊,在所有頻率上被歸一化至大小爲1的一數位 式影像。若此顆粒已被歸一化至某個不是1的値,則在應 用圖1 6所述的方法的其餘步驟之前,影像波道的強度値 可被縮放至1。例如:若此顆粒在所有頻率上被歸一化至 大小爲1 0 0 ,則在執行步驟1 6 1 2後且在執行其他步 驟之前,每一頻率大小可被1 0 0除。在一更複雜的例子 中,若每一頻率已對顆粒或其他高頻雜訊,而被歸一化成 一唯一値,則在此程序的其餘部分開始之前,可先除以適 當値。這除法也可發生在步驟1 6 1 2的變換之後。 如上述,圖1 6所例示的方法可被用在具一單一波道 或多個波道的影像。只要存在有多個波道,圖1 6所述的 方法可被使用。替換地,利用波道間的關係的一改良方法 可被使用。這樣的方法關聯的圖1 7而被說明如下。 圖1 6的說明假設此方法被應用至一單一波道,Μ ’ ,其中Μ ’代表一數位式影像的每一像素以空間排列的強 度値。對多波道的影像而言,存在這樣的多個波道。 在步驟1 6 1 2中,影像波道的多數個區段從空間域 被變換至頻率域。這樣的區段分割和變換可以根據上述與 雜訊歸一化相關聯的方法發生。吾人可以使用任何將影像 區段分割和從空間域變換至頻率域的變換。在這實施例中 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) (請先閱讀背面之注意事項再填寫本頁) 裝--------訂-------—. 經濟部智慧財產局員工消費合作社印製 -39- 經濟部智慧財產局員工消費合作社印製 496080 A7 B7 _ 五、發明說明(37 ) ’吾人使用一 1 6像素乘以1 6像素的區段大小。其他區 段大小可被使用,而不悖離本發明的範圍。此外,如上述 ,區段可以重疊,而不悖離本發明的範圍。同樣地,其他 區段形狀和大小可被使用,而不悖離本發明的範圍。對一 特別的區段而言,所產生的變換可被標示作Μ ’ Τ。如以 上所注意者,這樣的變換一般包括一系列的二維向量,其 中’此向量的一分量代表一特定頻率成分的大小,而向量 的另一分量代表在該特定頻率的相位。在這實施例中,每 一變換的向量可方便地被配置於兩個方向,而簡化在圖 1 6中所執行的方法之說明。圖1 8例示組織的這兩個方 向。 在圖1 8中,向量被配置在1 7x 1 7矩陣中。DC 頻率項被置於矩陣中心且被指定一座標爲(0 ,0)。此 外’如圖1 8中所不,這向量被指定一指標(index) 1 4 4。當此頻率向量被配置於一向量的一維陣列時可使 用此指標。在一維陣列中,DC向量在指標1 44上。在 二維陣列中,頻率向量被置於一指標,其中該指標對應在 頻率上和在x和y方向上於該位置上的適當改變。在一維 陣列中,此一維陣列在指標〇上的第一頻率向量爲圖1 8 所例示的二維陣列左上角(座標一 8 ,- 8 )所出現的向 量。在一維陣列中的向量之放置依數字順序繼續,其中, 該二維陣列的讀取是從第一列由左至右然後第二列由左至 右’依此類推,直到最後一列爲止。這配置是由圖1 8中 括弧起來的一維陣列所標示的,其中,該標示爲一特別的 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -40 - (請先閱讀背面之注意事項再填寫本頁) 裝 ----訂------;---' 496080 Α7 Β7 經濟部智慧財產局員工消費合作社印製 五、發明說明(38) 頻率向量,沿著二維陣列座標所列出的。 在步驟]L 6 χ 4中,於步驟1 6 1 2所得的每一區段 的:變換可被過濾。替換地’這步驟可被省略’而不淳離本 發明的範圍。爲了方便和與圖1 7關聯的說明的一貫性’ 在步驟1 6 1 2所得的變換可被標示爲Μ ’ T 2。這變換 可以能強調訊號內容和衰減影像中的雜訊的方式’而被過 濾。在此方法的這步驟中’可以使用一低通濾波器(丨0^ pas s f 1 1 t e r )較有利。適合的濾波器趨於在隨機雜訊中’ 得到平順的衰減與降低。因此,雖然可以使用任何低通濾 波器,但是,一適合之濾波器可能是用來提高存在於變換 中的有效訊號和抑制存在於變換中的雜訊的理想濾波器。 一RMS (方均根)計算可被用來進一步加重有效的訊號 。此濾波運算趨於保留在頻率中強的變化且抑制在頻率中 小的變化。換言之,可加重極可能代表有效的訊號的強頻 率成分。 在這實施例中,濾波器包含一 3 X 3矩陣’其中該矩 陣在中心的値爲4、在中心的水平左右或垂直上下點的値 爲2、且在矩陣的4個角的値爲1。其他大小的濾波器或 其他濾波器値可被使用,而不悖離本發明的範圍。若w ( a ,b )代表上述的濾波器,則可根據公式4 0計算出被 過濾的變換Μ ’ T 3。 、,y]:· TaMm i In n 11 ϋ— ϋ One sip, I H ϋ tmmmmm n n ϋ * n I 496080 A7 B7 5. Description of the invention (26) 1 7 0 2. Note: Here the X-axis represents the pixel intensity 値 'and not the average pixel intensity of the entire segment used previously. For a pixel with point 1 0 0: the intensity of the particle is estimated at the point 18 0 6 ° In other words, although the curve 1 70 2 is obtained from a scatter plot based on the average intensity, but if To set the intensity of a pixel, curve 1 0 2 (the curve can be the same as curve 1 7 2) can be used to provide an estimate of the particle intensity of the pixel. Measure the intensity of particles related to frequency After measuring the particle trajectory, you can measure the intensity of particles related to frequency for each color channel (or a single channel of a black and white image), as shown in step 1 in Figure 3. 0 4 shown. Frequency content represents image detail. A section with a higher degree of image detail will have more content at all frequencies. Example: Consider an image of a tree surrounded by the sky and standing on a meadow. The grassy part of the image, due to its many leaves and shades of color, has most of the high energy at all frequencies, that is, 'a large frequency content. Many high energies may have high frequency content. On the other hand, the sky part lacks image details, so it has high frequencies that rarely produce this image. Therefore, any measured high-frequency content can therefore be assumed to be generated primarily from grainy details. To measure the intensity of the particles in relation to frequency, we performed the steps illustrated in Figure 10. In step 2 0 2, the digital image is segmented for the second time by the segment, and in step 2 0 4, the second Fourier is performed. The paper size is in accordance with Chinese National Standard (CNS) A4 (210 X 297). (Mm) (Please read the notes on the back before filling out this page) Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs -29- Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs 496080 A7 B7 ___ V. Description of Invention (27 ) Transform, although, the results of the first segmentation and the first Fourier transform may alternatively be used. For new calculations, different section size, shape, and overlap options can be used (such options have been described previously), and different types of transformations can be used (again, such options have been described previously). In step 2 0 0, in order to identify the estimated particle intensity as a function of frequency, we generate a scatter plot with a fitted curve. Secondly, in step 2008, we calculate the intensity function according to the particle intensity. Among them, the function is obtained by correlating each segment with the frequency 値, and draws the size divided by the estimated particle intensity. In step 2 0 0, draw a line identifying the frequency-dependent estimated particle trajectory. Segment Segmentation When measuring frequency-dependent particle intensities ' In step 2012, the digital image can be segmented into segments using any of the techniques and options of step 1102 in FIG. 4 previously described. Fourier Transform In step 2 0 4 we use any of the techniques and options previously described in connection with step 1 1 4 in Fig. 4 'for each color channel (or a single channel in a black and white image) A second Fourier transform is performed on the pixel intensity of each segment. However, the range of low and medium frequency terms has not been removed, as shown previously in step 1 106 of FIG. Identification of frequency-related estimation questions. Grain track 1 This paper applies the national standard (CNS) A4 specification (210 X 297 public love) -30-(Please read the precautions on the back before filling this page) Pack. 496080 A7 B7 V. Description of the invention (28) (Please read the note on the back? Matters before filling out this page) Using the DC (average intensity) obtained from the second Fourier transform 値, we can identify each color channel (or In a single channel of a black-and-white image, the frequency-dependent particle intensity of each segment is as shown in step 2006. Lines previously drawn through the scatter plot of each color channel (or a single channel of a black and white image) (such as lines 1702 and 1802) can be used for this purpose. For each segment, the DC (average intensity) 値 obtained from the second Fourier transform can be used / for the intensity of the curve 値 to locate a point on the X axis, and therefore, find the estimated on the y axis Particle strength. For example, if a segment has a D C 値 which is equivalent to 値 at point 1 8 0 in FIG. 9, its estimated particle intensity will be 値 at point 1 8 1 0 on line 1 8 2. Therefore, each frequency of a segment is considered to have an estimated particle intensity derived from the average intensity of the segment. Book M / GS vs. Frequency Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs in step 10 of Figure 10, after identifying the intensity of particles related to frequency, in each color channel (or In a single channel of a black and white image), the size of each segment at each frequency can be determined by the estimated particle intensity of each segment in each color wave channel (or a single channel of a black and white image). segmentation. The magnitude of each frequency is determined by the Fourier transform in step 2010. A new scattergram can then be generated, one for each color channel (or one for a black and white image). For the channel in question, the y-axis of the scatter diagram represents the size of each section on each channel, where 'the 値 is the paper size applicable to the Chinese National Standard (CNS) A4 (210 X 297 mm -31-496080 A7 ____ B7 _ V. Description of the invention (29) (Please read the notes on the back before filling this page) It is divided by the estimated particle intensity of that area. For the channel in question, the X-axis of the scatter plot represents the frequency chirp present in the digital image. Therefore, for each frequency in each segment in the digital image, a point on the scatter plot appears during the day. In the example above, 289 frequencies are presented, so the new scatter plot in this example will have 2 889 S drawn points, where S represents the number of segments. In this scatter diagram, there is no weighting added to each point except for the number of points drawn at the same position. Specifically, each point has a weight of 1, but if the second point is drawn at the same position, the points that exist can be assigned a weight of 2. Alternatively, a list of points can be easily maintained and contain two points. FIG. 11 illustrates such a scatter diagram of a color channel. A lower "fog point" represents those segments that have been isolated by the above weights, and a higher probability represents the main particle intensity rather than the image signal. Those points removed from the fog zone have a higher probability of representing the image signal. Book lines to identify the particle trajectories related to frequency. Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs. Second, use any curve-fitting algorithm to draw all points in the lower misty bands in each scatter plot. Line, as shown in step 2 10 in Figure 10. In this embodiment, we use a median average type algorithm. When this line is drawn, the erroneously higher picture points 'outside the foggy band' can be ignored, and thus the results will not be distorted. In order to do this, you can use the following "weighting" formula for each point, where y is a point on the y-axis: This paper size applies the Chinese National Standard (CNS) A4 (210 X 297 mm> -32 -Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs 496080 A7 B7_ V. Description of the invention (30) Any other suitable weighting can be used without departing from the scope of the present invention. Figure 12 shows a color channel The line Έ 2 0 2 drawn by the foggy band is used as an example. Note: When drawing this line, the wrong high points 2 2 0 4 and 2 2 0 6 are ignored. At the same time, only medium and high frequencies are used It may be advantageous to determine this line by plot points. In this way, strong low-frequency content will not adversely affect the estimated DC. Calculate the correction array and then calculate a correction array for each color channel, as shown in Figure 3 Step 1 0 0 6. Find the DC term estimated by the graph. To calculate the correction array, we use the scatter plot drawn in step 2 1 0 (Figure 1 0) to calculate each color. Channel (or single in a black and white image Channel) in an estimated DC (EDC). For example: Figure 12 shows the EDC 値 indicated by the point 2 2 0 8 on the line 2 2 0 2. Then, in this example, these three scatter plots provide An estimated DC 値 for red, green, and blue. After calculating the EDC 値 for the array, use Equation 37 to calculate an array for each color channel (or a single channel in a black and white image), where [x] Cover the range of pixel intensity 値: Array U] = (complex · grain intensity [,]) (3 times) This paper size applies to China National Standard (CNS) A4 (210 X 297 mm)-33: (Please read first Note on the back then fill in this page) · in a— · n I «ϋ ai« i · — J ·· n «ϋ I nn _ Ύ Ύ 4 496080 A7 B7 V. Description of the invention (31) (Please read first Please fill in this page again on the unintended side) Use the estimated DC 所得 obtained from the curve calculated in step 2 1 10 and the intensity of each particle obtained from the curve calculated in step 1 1 4 § ten particle strength, An array is calculated for each channel. The calculation and correction array [X] can be used to calculate each channel (or A black and white channel) is a correction array. We can use Equation 3 8 to calculate this array. Correlation [〇] = 〇 Correlation [小 έΣ (Array Array ⑽ (38) / 7 = »! 2 Particle trajectory that will be related to intensity After the consumer cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs prints and calculates the correction array, the particle trajectory of the digital image can be normalized with respect to intensity, as shown in step 10 of FIG. 3. The chirp obtained from the corrected array can be applied to the data of this digital image to achieve normalization. The intensity 値 "X" of each pixel in the digital image can be replaced by 色彩 in the correction [X] for each color channel (or a single channel of a black and white image). This procedure makes the particle trajectory in the digital image have approximately the same visible impact (V 1 s u a 1 1 mp a c t) in the gray level of each channel, which makes the digital image look more desirable. For example, it is more desirable when it is displayed on a typical display with a contrast factor (gamma) of about two and a positive film. Figure 13 is from a normalized paper with pixel dimensions in space. The paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 public love): 34-'Employee Cooperative of Intellectual Property Bureau of the Ministry of Economic Affairs Print 496080 A7 B7____ V. Description of the invention (32) Graphical representation of the particle trajectories of the lines of the pixels of the digital image with a change in intensity. It should be noted that the range of pixel intensities over the entire graphic is uniform and does not fluctuate as sharply as the pixel intensities in an unnormalized digital image. In Figure 13 the fluctuations at point 2 3 0 2 are approximately the same as at point 2 3 0 4. Normalizing the particle trajectories related to frequency The particle trajectories of digital images can now be normalized in relation to frequency, as shown in step 1 0 1 0 in FIG. 3. After applying the correction array, we can perform another frequency transformation and transform each channel (or the single channel) of the image that has been normalized in intensity from the time domain to the frequency domain. Alternatively, the correction array can be converted into a function in the frequency domain for use in the original spatial frequency domain representing the image. In this embodiment, a correction array is applied in the time domain and conversion is performed on each channel to perform a frequency normalization procedure to be described. Figure 14 is a graphical 2D plot of the size versus frequency in a digital image after the software has normalized the particle trajectory across the intensity. Here, the average (DC) 値 2402 is 1 · 0. To normalize the particle trajectory across frequencies, this software multiplies the magnitude (M) of each frequency in each segment by a factor (F). This factor refers to the extent to which the line 2 4 0 2 drawn via the scatter diagram deviates from 1.0. Deviation factor (F) = (39) This paper size is applicable to the Chinese National Standard (CMS) A4 specification (210 X 297 mm). 35-(Please read the precautions on the back before filling this page) Installation ----- --- Order ------... 496080 A7 B7 V. Description of Invention (33) (Please read the notes on the back before filling this page) The actual graph of the deviation factor will show a frequency increase curve, The curve is an inversion of the line 2 4 0 2 in FIG. 14. This phase is not affected and is retained. Figure 15 graphically represents a 2D plot of frequency versus frequency in a digital image after the software has multiplied the deviation factor. In the conventional way, a line 2 5 0 2 drawn through this figure shows that the size of the particle trajectory is large and the area of the high and low frequencies is equal. They do not fluctuate as sharply as in an unnormalized digital image. Suppressing particle trajectories After calculating the correction array, the normalized 値 obtained from the correction array can be used to suppress the particle trajectories in the digital image, as shown in step 1 0 1 2 in Figure 3. Other methods of particle suppression may be used without departing from the teachings of the present invention, such as the second method described below. In the digital image normalization according to the teachings above, in each color channel, the intensity of the particles in each segment and in each frequency has been adjusted so that the size component is 1 · 0, as shown in Figure 1 5 As shown. The aspect is not affected by the printing of employee cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs. Therefore, the size of a component is a powerful indicator that isolates the particle trajectory (noise) from the image details in the 1.0 range is a number of watermarks. Sizes greater than 1 · 0 are more likely to be generated from an image detail, and sizes less than 1 · 0 are more likely to be generated from a particle trajectory. Therefore, in order to suppress particle trajectories in a digital image, it is best to suppress components that are close in size and smaller than 1.0. On the other hand, components larger than 1.0 can be considered to represent image details and should not be suppressed. -36- This paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) 496080 Ministry of Economic Intelligence Printed by the Consumer Cooperative of the Property Bureau A7 B7 V. Invention Description (34). This is in contrast to unnormalized digital images, where 'the' number of watershed lines between glume trajectory (noise) and the image varies with intensity and frequency, making it almost impossible to accurately match the particle trajectory with Image isolation. As mentioned earlier, images can be normalized 'or data describing the distribution of noise with intensity and frequency can be included in an unnormalized image' to guide particle reduction without departing from the invention Range. In order to suppress particle trajectory, the following method can be applied to the size (M) of each frequency in the normalized digital image segment: F οr each segment in the normalized digital image segment Maintain the original phase information F οr each M I f (Μ < = 0 · 5) in this section (including the DC term) Μ = Ο else suppression factor = 1.0. E-(Μ- 0.5) Μ = Μ * The suppression factor formula suppresses 値 by replacing the normalized frequency of 0 or less with 0. Because some glume trajectories have a normalized size slightly larger than 0.5, this formula can also use a suppression factor to replace their 値 with smaller 値. This suppression factor eases the reduction of frequencies with a higher normalized size, 'because' these tend to correspond to the image details. The above example has a digital image, where the digital image is based on the Chinese paper standard (CNS) A4 (210 X 297 mm) -37------ ί ---: --- -_ 装 i — (Please read the notes on the back before filling out this page). Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs 496080 A7 ____B7__ 5. Description of the invention (35) Known particle strength as a function of frequency and intensity Each segment is divided to normalize it, and then the reduction is applied with a watershed fixed at 1.0. A comparable method can also apply the particle reduction of a watershed point with expected noise equal to a function of frequency and intensity directly to the unnormalized components. In order to carry out this method, we can use the above-mentioned glume suppression method, and in the suppression factor, M ′ is used instead of M ′, where the desired particle trajectory of __M advantageously uses the normalized 値 from the calibration array, Other methods of suppressing glomerular trajectory are possible. For example, image-improving software such as PHOTOSHOP can be combined with the normalized combination manually, although such software requires extensive user feedback to find the appropriate level of suppression of particle trajectories. Now, another method for suppressing some or substantially all particles of an image, either automatically or at the user's choice, will now be described. FIG. 16 illustrates a flowchart illustrating a method for removing high-frequency noise such as grain tracks from a digital image. The method shown in Figure 16 can be applied to one or more channels of the digital image or, if the digital image has only one channel, it is applied to a single channel of the digital image. For an image with multiple channels, the method described in Figure 16 can be applied to a subset of the total channels without departing from the scope of the present invention. Likewise, the method described in FIG. 16 can be applied to each of a plurality of channels in a multi-channel digital image without departing from the scope of the present invention. In step 1610, the noise of a particular image channel is normalized. This noise can be normalized using the above method or any other method ---, --- »(* ---- ^ 装 —— (Please read the precautions on the back before filling this page) · This The paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 public love) -38- 496080 A7 B7___ V. Description of the invention (36). Alternatively, although before the rest of the method described in Figure 16 is performed It is better to normalize this noise, but this method can also be applied to an unnormalized digital image. The method described in Figure 16 can be applied to high frequencies that used to have particle traces such as: Noise is normalized to a digital image of size 1 at all frequencies. If the particle has been normalized to a chirp other than 1, then the remaining steps of the method described in Figure 16 are applied Previously, the intensity 値 of the image channel could be scaled to 1. For example: if this particle is normalized to a size of 0 0 at all frequencies, then after performing step 16 2 and before performing other steps, Each frequency can be divided by 100. In a more complex example, if each frequency High-frequency noise, which is normalized to a unique chirp, can be divided by the appropriate chirp before the rest of this procedure begins. This division can also occur after the transformation of step 16 2 as shown above. The method exemplified in 16 can be used for images with a single channel or multiple channels. As long as there are multiple channels, the method described in Figure 16 can be used. Alternatively, the inter-channel An improved method of relationship can be used. Figure 17 related to this method is illustrated as follows. The description of Figure 16 assumes that this method is applied to a single channel, M ', where M' represents a digital image. The intensity of each pixel is arranged in space. For a multi-channel image, there are such multiple channels. In step 16 2, most sections of the image channel are transformed from the spatial domain to the frequency. Domain. Such segment segmentation and transformation can occur according to the above-mentioned methods associated with noise normalization. We can use any transformation that segments and transforms image segments from the spatial domain to the frequency domain. In this embodiment, this Paper size applies to Chinese national standard Standard (CNS) A4 (210 X 297 mm) (Please read the precautions on the back before filling out this page) Installation -------- Order -----------. Intellectual Property Bureau, Ministry of Economic Affairs Printed by employee consumer cooperatives-39- Printed by employee consumer cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs 496080 A7 B7 _ V. Description of the invention (37) 'I use a size of 16 pixels by 16 pixels. The size of other segments It can be used without departing from the scope of the invention. In addition, as described above, the sections can overlap without departing from the scope of the invention. Similarly, other section shapes and sizes can be used without departing from the scope of the invention. The scope of the invention. For a particular segment, the resulting transformation can be labeled as M'T. As noted above, such transformations generally include a series of two-dimensional vectors, where one component of this vector represents the magnitude of a particular frequency component, and the other component of the vector represents the phase at that particular frequency. In this embodiment, each transformed vector can be conveniently configured in two directions, simplifying the description of the method performed in FIG. 16. Figure 18 illustrates these two directions of the organization. In Figure 18, the vectors are arranged in a 17x17 matrix. The DC frequency term is placed in the center of the matrix and is assigned a label (0, 0). In addition, as shown in Figure 18, this vector is assigned an index 1 4 4. This index can be used when the frequency vector is configured in a one-dimensional array of vectors. In a one-dimensional array, the DC vector is on index 1 44. In a two-dimensional array, the frequency vector is placed on an index, where the index corresponds to a suitable change in frequency and in that position in the x and y directions. In a one-dimensional array, the first frequency vector of the one-dimensional array on the index 0 is the vector appearing in the upper left corner of the two-dimensional array (coordinates 8,-8) illustrated in FIG. 18. The placement of the vectors in the one-dimensional array continues in numerical order, where the reading of the two-dimensional array is from the first column from left to right and the second column from left to right, and so on, until the last column. This configuration is marked by the one-dimensional array enclosed by the brackets in Figure 18, where the marking is a special paper size applicable to the Chinese National Standard (CNS) A4 (210 X 297 mm) -40-(Please (Please read the notes on the back before filling this page) Binding ---- Order ------; --- '496080 Α7 Β7 Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economy , Listed along the two-dimensional array coordinates. In step] L 6 χ 4, the: transformation of each section obtained in step 16 12 can be filtered. Alternatively, 'this step may be omitted' without departing from the scope of the present invention. For the sake of convenience and consistency of the description associated with FIG. 17 ', the transformation obtained at step 1612 may be designated as M'T2. This transformation may be filtered to emphasize signal content and attenuate noise in the image. In this step of the method, it may be advantageous to use a low-pass filter (0 0 pas s f 1 1 t e r). Suitable filters tend to get smooth attenuation and reduction in random noise. Therefore, although any low-pass filter can be used, a suitable filter may be an ideal filter for improving the effective signal present in the transform and suppressing the noise present in the transform. An RMS (root mean square) calculation can be used to further emphasize the effective signal. This filtering operation tends to retain strong changes in frequency and suppress small changes in frequency. In other words, the strong frequency components that are likely to represent a valid signal can be exacerbated. In this embodiment, the filter contains a 3 X 3 matrix, where 値 in the center is 4; 在 in the horizontal left and right or vertical points in the center is 2; and 値 in the 4 corners of the matrix is 1. . Other size filters or other filters may be used without departing from the scope of the invention. If w (a, b) represents the above-mentioned filter, the filtered transformation M'T3 can be calculated according to Formula 40. ,, y]:
Σ Σ M.a,b)(Mr2(x^ayy^b))1 ·—_丨 ‘__I Σ Σ Κα,ό) (40) f請先閱讀背面之注音〗事項再填寫本頁} M,T3[DC] 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -41 - 496080 經濟部智慧財產局員工消費合作社印製 A7 B7 五、發明說明(39) 如所注意到的:在被過濾的變換中,D C項被設爲〇 。此外,在此變換的邊緣之外的元件被視爲其値爲0。在 這公式中,這個處置並不明顯的表示出。此濾波運算只有 影響頻率向量的大小。頻率向量的相位維持不變。在公式 4 0中,吾人使用變換的二維表示來計算。爲了記號法的 目的,使用變換的大小之事實由在公式4 0中使用平常字 母(與粗體字母相反)作爲變換Μ ’ Τ2和Μ” Τ3的事 $ 實所指示出。這記號法將在說明圖1 6和1 7中一致地被 使用。當執行一向量運算時’將使用粗體字母在所討論的 特別變換上。 在步驟1 6 1 6中,每一區段的被過濾的變換可被頻 率展開(frequency spread )。替換地’若步驟1 6 1 4 被省略,則在步驟1 6 1 2中所得的原始的變換可在步驟 1 6 1 6中被頻率展開。步驟1 6 1 6也可以被省略’而 不悖離本發明的範圍。 因爲真實世界是由類似的重複圖樣所組成的’所以, 一照相影像的變換空間趨於在D C頻率周圍的任何放大’ 具有相同的外表,除非存在重製的人工製品。對重複的圖 樣而言,一影像的頻率特性趨於以2次方重複。在這步驟 中,在較低頻的頻率特性被用來以雙倍頻率,預先使門濫 細節易於感光。 爲了執行這樣的頻率展開,若吾人使用圖1 8所闡述 的變換之二維版本’則此變換資料的中心一半被擴大成變 換資料的全尺寸區塊。換言之,頻率的中心9 x 9區塊被 (請先閱讀背面之注意事項再填寫本頁) 裝--------訂----- 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -42 - 496080 A7 B7 五、發明說明(4〇 ) 展開爲在1 7 X 1 7變換中的每一其他位置。當一頻半在 變換的剩餘部份中’並未明確地被塡入時’則可使用鄰近 頻率的平均値。再次地,這運算只有影響頻率大小。 在得出被過濾的變換之頻率展開之後’可取決於頻率 展開變換,和在步驟1 6 1 4所得的被過濾的變換’在步 驟1 6 1 8中計算一新的變換。替換地,若步驟1 6 1 4 被省略,則新的變換取決於在步驟1 6 1 2所得的原始的 變換和頻率展開變換。公式4 1可被使用來計算此新的變 換 Μ ’ T 4。 M,T4[x] = 73[x]V(0.2^ BT[x]) 2阱 4 (請先閱讀背面之注意事項再填寫本頁) 裝 經濟部智慧財產局員工消費合作社印製 在公式4 1中,變換Μ’ ΕΤ代表由頻率展開所得的 估計的變換。再次地,這運算只有影響頻率向量的大小° 在步驟1620中’吾人可產生一遮罩(n^sk) ’使 得在步驟1 6 1 8所得的變換之高頻項減幅。因爲顆粒軌 跡趨於具有主要的高頻成分,所以,該遮罩可被用來使顆 粒軌跡和/或其他高頻雜訊減幅。 可使用由公式4 2所代表的選通函數(2at inS f u n c t i ο η )得到此遮罩。Σ Σ Ma, b) (Mr2 (x ^ ayy ^ b)) 1 · -_ 丨 '__I Σ Σ Κα, ό) (40) f Please read the note on the back first and then fill out this page} M , T3 [ DC] This paper size is in accordance with China National Standard (CNS) A4 (210 X 297 mm) -41-496080 Printed by A7 B7, Consumer Cooperative of Intellectual Property Bureau of the Ministry of Economic Affairs 5. Description of Invention (39) As noted: In the filtered transform, the DC term is set to zero. In addition, elements outside the edges of this transformation are considered to have zero. In this formula, this treatment is not clearly shown. This filtering operation only affects the size of the frequency vector. The phase of the frequency vector remains unchanged. In Equation 40, we use the transformed two-dimensional representation for calculations. For the purposes of notation, the fact that the size of the transformation is used is indicated by the use of ordinary letters (as opposed to bold letters) in Equation 40 as a matter of transforming M'T2 and M "T3. This notation will be in Illustrate that Figures 16 and 17 are used consistently. When performing a vector operation 'will use bold letters on the particular transformation in question. In steps 1 6 1 6 the filtered transformation of each segment Can be frequency spread. Alternatively, if step 1 6 1 4 is omitted, the original transform obtained in step 1 6 1 2 can be frequency spread in step 1 6 1 6. Step 1 6 1 6 can also be omitted without departing from the scope of the present invention. Because the real world is composed of similar repeating patterns, so the transformation space of a photographic image tends to be any magnification around the DC frequency, and has the same appearance Unless there are duplicate artefacts. For repeated patterns, the frequency characteristics of an image tend to repeat in the second power. In this step, the frequency characteristics at the lower frequency are used to double the frequency in advance. Make door detail Easily photosensitive. In order to perform such frequency expansion, if we use a two-dimensional version of the transformation illustrated in Figure 18, then the center half of this transformation data is enlarged to a full-size block of the transformation data. In other words, the center of the frequency is 9 x Block 9 (Please read the precautions on the back before filling out this page) Loading -------- Order ----- This paper size applies to China National Standard (CNS) A4 (210 X 297 mm) ) -42-496080 A7 B7 V. Description of the invention (4〇) expands to every other position in the 1 7 X 1 7 transformation. When a frequency and a half are in the remainder of the transformation, it is not explicitly entered You can use the mean 値 of adjacent frequencies. Again, this operation only affects the frequency. After the frequency expansion of the filtered transform is obtained, it can depend on the frequency expansion transform and the result obtained in step 1 6 1 4 The filtered transform 'calculates a new transform in step 16 1 8. Alternatively, if step 1 6 1 4 is omitted, the new transform depends on the original transform and frequency expansion transform obtained in step 1 6 1 2 Equation 4 1 can be used to calculate this new transformation M T 4. M, T4 [x] = 73 [x] V (0.2 ^ BT [x]) 2 well 4 (Please read the precautions on the back before filling this page) Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs In Equation 41, the transformation M ′ ET represents the estimated transformation obtained by the frequency expansion. Again, this operation only affects the size of the frequency vector ° In step 1620, 'we can generate a mask (n ^ sk)' such that The high-frequency term of the transform obtained in step 16 18 is attenuated. Since the particle trajectory tends to have the main high-frequency component, this mask can be used to attenuate the particle trajectory and / or other high-frequency noise . This mask can be obtained using a gating function (2at inS f u n c t i ο η) represented by Equation 42.
ϋ ϋ «ϋ n —Bi 一< n i n an n «n n I %· 64 (42) 公式4 2是根據圖1 8所例示的傅立葉變換的二維表 示。在公式4 2中的X和y項當變換被配置再如圖1 8所 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -43 - 496080 A7 B7 五、發明說明(41 ) 例示的一二維矩陣時,對應在變換中特別的頻率向量之座 標。因此,此選通函數對高頻,趨近於2 ,且對低頻,趨 近於0。此遮罩於是可使用公式4 3而被計算,其中,B (X,y )包括在X,y的遮罩,也是指圖1 8所例示的 二維變換中的座標。 B(x,y)=l-Fgatc (x,y) Fgatc <0.4 (43) B(x,y)=0.6ϋ ϋ «ϋ n —Bi i < n i n an n« n n I% · 64 (42) Equation 4 2 is a two-dimensional representation of the Fourier transform illustrated in FIG. 18. The X and y terms in Equation 4 2 are configured as shown in Figure 18. The paper size applies the Chinese National Standard (CNS) A4 (210 X 297 mm) -43-496080 A7 B7 V. Description of the invention ( 41) When the two-dimensional matrix is exemplified, it corresponds to the coordinates of the special frequency vector in the transformation. Therefore, this gating function approaches 2 for high frequencies and approaches 0 for low frequencies. This mask can then be calculated using Equation 43, where B (X, y) is included in the mask of X, y, which also refers to the coordinates in the two-dimensional transformation illustrated in Figure 18. B (x, y) = l-Fgatc (x, y) Fgatc < 0.4 (43) B (x, y) = 0.6
Fgacc > 0.4 在這實作中,此遮罩函數被設計爲具有最小値爲 〇.6。也可使用其他遮罩函數,而不悖離本發明的範圍 。此外,這步驟可被省略,而不悖離本發明的範圍。 在步驟1 6 2 2中,可以對每一區段計算一雜訊校正 函數。在這實施例中,爲回應遮罩和先前所計算出的影像 的頻率域表示之一,計算此雜訊校正函數。例如:在步驟 16 18所得的變換Μ’ Τ 4可被用來得出此雜訊校正函 數。可由公式4 4和4 5得出雜訊校正函數。 (請先閱讀背面之注意事項再填寫本頁) 一裝 II 一 r · n ·Βϋ I 11 n ·Βϋ I iso· 經濟部智慧財產局員工消費合作社印製 (44) C(x,y)=l + FKatc (x,y) c > 1 + F_(x,y) M'TM[x,y] = GR + (1-GR) C(x,y) (45) 在公式4 4中,雜訊校正函數被限制在1加上在變換 矩陣中一特別的位置之選通函數的値。再次地’公式4 4 只有影響頻率向量的大小分量。本發明的使用者可能想要 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -44- 496080 Α7 Β7 經濟部智慧財產局員工消費合作社印製 五、發明說明(42 ) 在根據本發明的改良後的影像上存留有殘餘的顆粒性。因 此,本發明允許使用者指定在影像被改良後,其所欲留下 的顆粒殘餘(g r a 1 n r e s i d u e ) G r。公式4 5在計算雜訊 校正函數Μ ’ T Μ時,將所要的殘留顆粒考慮進去。 在步驟1 6 2 4中,每個區段的雜訊校正函數被應用 至該區段的原始變換。其結果爲如:顆粒軌跡之高頻雜訊 被抑制後的影像之頻率域表示。使用公式4 6可得出被抑 制顆粒的頻率表示。 M’TSU]=M’TU]M,TM[x] (46) 公式46產生包括向量的Μ’ Τ系列和在Μ’ ΤΜ中 對應的無向量的値之乘積的一系列向量。 在步驟1 6 2 6中,可以對每一區段算出雜訊校正變 換Μ ’ T S的反變換,以得出已經將雜訊抑制後的改良影 像的空間域表示。在圖1 7中所描繪的方法可以被使用在 具多數個波道的數位式影像。若一數位式影像僅具有一單 一波道,則與圖1 6關聯所述的方法較爲可用。吾人也可 在數位式影像的多數個波道中的某些波道上使用圖1 7所 述的方法,而在其他的波道上應用與圖1 6關聯所述的方 法。 與圖1 7關聯所述的方法可被應用至任何多波道數位 式影像,而不悖離本發明的範圍。實例包括:彩色照片( 印出的照片、負片或幻燈片)的數位式影像表示、從人造 衛星所接收的數位式影像和醫學影像。本發明在從數位式 影像移除諸如:由顆粒軌跡所產生的非相關高頻雜訊上很 (請先閱讀背面之注意事項再填寫本頁) 裝Fgacc > 0.4 In this implementation, this mask function is designed to have a minimum 値 of 0.6. Other mask functions may be used without departing from the scope of the invention. Moreover, this step can be omitted without departing from the scope of the invention. In step 162, a noise correction function can be calculated for each segment. In this embodiment, the noise correction function is calculated in response to one of the mask and one of the frequency domain representations of the previously calculated image. For example, the transform M ′ T 4 obtained in step 16 18 can be used to derive the noise correction function. The noise correction function can be obtained from Equations 4 4 and 4 5. (Please read the precautions on the back before filling out this page) One Pack II One r · n · Βϋ I 11 n · Βϋ I iso · Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economy l + FKatc (x, y) c > 1 + F_ (x, y) M'TM [x, y] = GR + (1-GR) C (x, y) (45) In Equation 4 4 The noise correction function is limited to 1 plus the gating function at a particular position in the transformation matrix. Again 'Formula 4 4 affects only the magnitude components of the frequency vector. Users of the present invention may want this paper size to comply with the Chinese National Standard (CNS) A4 specification (210 X 297 mm) -44-496080 Α7 Β7 Residual graininess remains on the improved image according to the invention. Therefore, the present invention allows the user to specify the grain residue (g r a 1 n r e s i d u e) G r that he wants to leave after the image is improved. Equation 4 5 takes into account the desired residual particles when calculating the noise correction function M'TM. In step 16 2 4 the noise correction function of each zone is applied to the original transformation of that zone. The result is, for example, the frequency domain representation of the image after the high-frequency noise of the particle trajectory is suppressed. Using Equation 4 6 gives a frequency representation of the suppressed particles. M'TSU] = M'TU] M, TM [x] (46) Equation 46 produces a series of vectors including the product of the vector's M 'T series and the corresponding vectorless unitary unitary in M' TM. In step 16 26, the inverse transform of the noise correction transform M'TS can be calculated for each sector to obtain the spatial domain representation of the improved image after the noise has been suppressed. The method depicted in Figure 17 can be used in digital images with multiple channels. If a digital image has only a single channel, the method described in association with FIG. 16 is more applicable. We can also use the method described in Figure 17 on some of the many channels of the digital image, and apply the method described in association with Figure 16 on other channels. The method described in association with FIG. 17 can be applied to any multi-channel digital image without departing from the scope of the present invention. Examples include digital image representations of color photos (printed photos, negatives, or slides), digital images received from satellites, and medical images. The present invention is very effective in removing uncorrelated high-frequency noises such as those caused by particle trajectories from digital images (please read the precautions on the back before filling this page).
l #11 I 一 口,· ϋ in I i n· Ί·— 11 I %- 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公爱) 45 496080 A7 B7 ____ 五、發明說明(43 ) 有用。 在步驟1 7 1 〇中,每一影像波道的雜訊被歸一化。 可‘以使用上述的或任何其他的方法而將雜訊歸一化。替換 地,雖然在實行圖1 7所述的方法的剩下部分之前,將雜 訊歸一化較好,但是,此方法也可被應用在一未經歸一化 的數位式影像。 圖1 7所述的方法可被應用至具諸如:顆粒軌跡的高 頻雜訊,在所有的頻率上被歸一化至大小爲1的一數位式 影像。若此顆粒已經被歸一化至某一非1的値,則在應用 圖1 7所述的方法的剩下步驟之前,此影像波道的強度値 可被縮放爲1。例如··若此顆粒在所有頻率上,已被歸一 化至1 0 0的値,則在執行步驟1 7 1 2之後且在執行其 他步驟之前,每一頻率大小可被除以1 0 0。在一更複雜 的例子中,若每一頻率已對顆粒或其他高頻雜訊而被歸一 化至一唯一値,則在進行剩下的步驟前,可以除之以一適 當的値。這除法也在步驟1 7 1 2中發生的變換中進行。 由圖1 7所例示的方法利用一數位式影像的多個波道 間的關係。爲方便起見,由圖1 7所例示的方法將在下文 中與代表一彩色照相影像的一數位式影像關聯地被說明。 這樣的數位式影像通常具有三個波道一一紅色波道、綠色 波道和藍色波道。爲了延伸圖1 7所述的方法至具有更多 波道的影像,以下所計算出的互相關函數可必要地被延伸 至將所有有用的互相關考慮進來。同樣地’若圖1 7被應 用至一雙波道數位式影像,則以下所述的互相關可被調整 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -46 - (請先閱讀背面之注意事項再填寫本頁) 裝 % 經濟部智慧財產局員工消費合作社印製 496080 A7 B7___ 五、發明說明(44 ) 至只考慮這兩個波道。 圖1 7應用此方法至三個波道,R ’ ’ G ’和B ’ 。 (請先閱讀背面之注意事項再填寫本頁) R ’代表一數位式影像的紅色波道的每一像素,在空間上 配置的強度値。G ’代表一數位式影像的綠色波道的每一 像素,在空間上配置的強度値。B ’代表一數位式影像的 藍色波道的每一像素,在空間上配置的強度値。在此所用 的標示法可避免與先前的公式混淆。 經濟部智慧財產局員工消費合作社印製 在步驟1 7 1 2中,每個影像波道的多數個區段從空 間域被變換至頻率域。若高頻雜訊只有從複數個影像波道 的一子集合移除,則只有那些影像波道可被變換或所有的 影像波道可被變換。換言之,若雜訊只有從一單一波道被 移除,則其他波道的資訊可被使用或可不必使用,來決定 將被移除雜訊的波道的雜訊內容。每一波道的區段分割和 變換可以根據與雜訊歸一化關聯的上述方法而發生。任何 分割影像波道的方法或任何從空間域被變換至頻率域的方 法均可被使用。在這實施例中,吾人選用1 6像素X 1 6 像素的區段大小。其他區段大小可被使用’而不悖離本發 明的範圍。此外,如上述,區段可以重疊,而不悖離本發 明的範圍。同樣地,其他區段形狀和大小可被使用’而+ 悖離本發明的範圍。可變的形狀和大小的區段也可被使用 〇 對一特別的區段而言,用做本發明的一例的照相影像 ,所產生的變換可以被標示爲R’ T、G’ τ和B’ 如以上所注意者,這樣的變換包括一系列的二維向量’其 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) 717^ 496080 A7 B7 經濟部智慧財產局員工消費合作社印製 五、發明說明(45) 中,此向量的一分量代表一特定頻率成分的大小,而向量 的第二分量代表在該特定頻率的相位。在這實施例中,每 一變換的向量可方便地被配置於兩個方向,而簡化在圖 1 了中所執行的方法之說明。圖1 8例示組織的這兩個方 向。圖1 8在前文中與圖1 6的說明相關聯地被描述。 在步驟1 7 1 4中,複數影像波道的對應區段間的互 相關,可以針對每一區段已在步驟1 7 1 2中計算出的變 換,做計算。替換地,可以使用區段的一子集,而不悖離 本發明的範圍。因爲一數位式影像的複數個波道趨於顯示 相同的影像圖樣,不同波道的互相關之對應區段將趨於爲 影像-細節,而非如顆粒的隨機雜訊。相關不足的頻率將 趨於代表雜訊,諸如:顆粒。因此,本發明利用對影像細 節,在波道間的相關關係。本發明的各種不同的實施例可 以各種不同的方式利用這關係。以下的數學敘述提供互相 關如何被用來移除雜訊的一例。然而,藉由執行與一數位 式影像的波道間互相關成比例的計算,吾人可使用任何數 學方法,於數位式影像中的雜訊之減少,包括顆粒的減少 ,而不悖離本發明的範圍。 公式4 7 - 4 9提供得到該數位式影像的波道間的互 相關的量測的方法。 RT[x]-GT[x] RO,X[x] = 1 + |RT [x]||GT[x| 其中DC = 0 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) (請先閱讀背面之沒意事項再填寫本頁) 以 « tt tm§ ammmmw I n« ϋ— n 一e,t t n I n < (47) -48- 496080 A7 B7 五、發明說明(46) RT[x]-BT[x] R,B,X[x]: (48) 其中DC = 0 (49) (請先閱讀背面之注意事項再填寫本頁) ⑽呵个—P,TW,B,TW 1 十丨 G,T_B,TW| 其中DC = 0 經濟部智慧財產局員工消費合作社印製 這些公式可被應用在數位式影像的每一區段的每一頻 率向量。公式4 7 - 4 9的計算結果爲互相關的數量値之 矩陣,其中,每一値與在不同波道中兩個頻率向量間的互 相關成正比,其中,該頻率向量代表此在二波道的空間對 應區段中的相同頻率。爲方便起見,公式47 - 49使用 與圖1 8相關聯而被討論的影像區段的空間頻率域變換的 一維矩陣之向量表示法。因此,每個公式的分子爲兩向量 的點乘積,而每個公式的分母將這些向量的每一大小考慮 進來。公式4 7產生的矩陣與一特別的區段在特別的頻率 ,紅色和綠色波道之間的互相關成正比。公式4 8產生的 結果與數位式影像的一特別的區段在特別的頻率,紅色和 藍色波道之間的互相關成正比。公式4 9產生的矩陣與一 特別的區段在不同的頻率,綠色和藍色波道之間的互相關 成正比。吾人可以依此三個公式,對數位式影像的每一區 段,計算出一矩陣。此互相關的D C項被設爲〇。吾人使 用R, G’ 、R, B’和G’ B’記號,避免與上述的顆 粒歸一化之公式相混淆。爲避免混淆,吾人使用R ’ T、 B ’ T和G ’ T記號代表步驟1 7 1 2所得的變換。 -49- 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) 496080 B7 經濟部智慧財產局員工消費合作社印製 五、發明說明(47 ) 對大的同相頻率向量而言,公式4 7 一 4 9的結果趨 於趨近値1 ,且對大的異相頻率向量而言’趨近値一 1 ° 很小的向量趨於得到趨近於0的結果,而同相且具約爲1 的> 大小之向量產生約爲0 · 5的結果。 在步驟1 7 1 6中,涉及紅色波道的互相關高頻成分 可被衰減。替換地,這步驟可被省略,而不悖離本發明的 範圍。通常,用來產生照相影像的彩色軟片具有一綠色層 ,其中,光線穿透該層,接著藍色層,而後是紅色層。因 爲紅色層在光線穿過綠色和藍色層後接收光線,所以,在 紅色層的影像比起在綠色和藍色層,稍微模糊不淸。紅色 層趨於缺乏構成影像資訊的高頻資訊。反而’來自紅色波 道的大部分高頻資訊趨於爲雜訊,且可被拋棄。右使用 不同類型的彩色軟片’則可能要衰減此數位式影像的不同 波道。此外,若其他影像,諸如:人造衛星影像或醫學影 像具有趨於缺乏高頻影像內容的的波道’則可以使用以下 所述的方法或類似的方法’衰減那些影像波道的高頻成分 〇 公式5 0說明可被使用於本發明的一紅色頻率衰丨咸S 數l # 11 I bite, · ϋ in I in · Ί · — 11 I%-This paper size applies to China National Standard (CNS) A4 (210 X 297 public love) 45 496080 A7 B7 ____ V. Description of the invention (43) it works. In step 1710, the noise of each image channel is normalized. Noise may be normalized using the above or any other method. Alternatively, although the noise is better normalized before implementing the rest of the method described in FIG. 17, this method can also be applied to an unnormalized digital image. The method described in Figure 17 can be applied to high-frequency noise with particle trajectories, which are normalized to a digital image of size 1 at all frequencies. If this particle has been normalized to a non-1 chirp, the intensity 値 of this image channel can be scaled to 1 before the remaining steps of the method described in Figure 17 are applied. For example ... If this particle has been normalized to 値 of 1 0 at all frequencies, then after performing step 1 7 1 2 and before performing other steps, each frequency can be divided by 1 0 0 . In a more complex example, if each frequency has been normalized to a unique chirp for particles or other high-frequency noise, then it can be divided by a suitable chirp before proceeding to the remaining steps. This division is also performed in the transformation that occurs in step 1712. The method exemplified by Fig. 17 makes use of the relationship between multiple channels of a digital image. For convenience, the method illustrated by FIG. 17 will be described below in association with a digital image representing a color photographic image. Such digital images usually have three channels-a red channel, a green channel, and a blue channel. To extend the method described in Figure 17 to images with more channels, the cross-correlation functions calculated below may be extended to take into account all useful cross-correlation. Similarly 'If Figure 17 is applied to a dual-channel digital image, the cross-correlation described below can be adjusted. This paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) -46- (Please read the precautions on the back before filling out this page) Installation% Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs 496080 A7 B7___ V. Description of Invention (44) To consider only these two channels. Figure 17 applies this method to three channels, R''G 'and B'. (Please read the notes on the back before filling in this page) R ′ represents the intensity of each pixel in the red channel of a digital image. G 'represents the intensity of spatial arrangement of each pixel of the green channel of a digital image. B ′ represents the spatially arranged intensity of each pixel of the blue channel of a digital image. The notation used here can avoid confusion with the previous formula. Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs. In step 17 12, most sections of each image channel are transformed from the spatial domain to the frequency domain. If high-frequency noise is removed from only a subset of the plurality of image channels, only those image channels can be transformed or all image channels can be transformed. In other words, if the noise is only removed from a single channel, the information of other channels may or may not be used to determine the noise content of the channel from which the noise is to be removed. The segmentation and transformation of each channel can occur according to the above-mentioned methods associated with noise normalization. Any method of segmenting the image channels or any method of transforming from the spatial domain to the frequency domain can be used. In this embodiment, we choose a segment size of 16 pixels X 16 pixels. Other segment sizes can be used 'without departing from the scope of the invention. Further, as described above, the sections may overlap without departing from the scope of the present invention. Likewise, other segment shapes and sizes can be used 'and + deviates from the scope of the invention. Variable shape and size segments can also be used. For a particular segment, used as a photographic image of the present invention, the resulting transformation can be labeled R 'T, G' τ, and B. 'As noted above, such a transformation includes a series of two-dimensional vectors.' The paper size is applicable to the Chinese National Standard (CNS) A4 specification (210 X 297 mm) 717 ^ 496080 A7 B7 Staff Consumption of Intellectual Property Bureau, Ministry of Economic Affairs Cooperative printed 5. Invention Description (45), one component of this vector represents the magnitude of a specific frequency component, and the second component of the vector represents the phase at that specific frequency. In this embodiment, each transformed vector can be conveniently configured in two directions, and the description of the method performed in FIG. 1 is simplified. Figure 18 illustrates these two directions of the organization. Figure 18 was previously described in association with the description of Figure 16. In step 1714, the cross-correlation between the corresponding sections of the complex image channel can be calculated for the transformation that has been calculated in step 1712 for each section. Alternatively, a subset of the segments may be used without departing from the scope of the invention. Because a plurality of channels of a digital image tend to display the same image pattern, the corresponding sections of cross-correlation of different channels will tend to be image-details, rather than random noise such as particles. Inadequately correlated frequencies will tend to represent noise, such as particles. Therefore, the present invention makes use of the correlation between the channels of the image details. Various embodiments of the invention may utilize this relationship in a variety of different ways. The following mathematical description provides an example of how correlation can be used to remove noise. However, by performing calculations that are proportional to the cross-correlation between channels of a digital image, we can use any mathematical method to reduce the noise in the digital image, including the reduction of particles, without departing from the present invention. Range. Equations 4 7-4 9 provide a method of measuring the cross-correlation between the channels of the digital image. RT [x] -GT [x] RO, X [x] = 1 + | RT [x] || GT [x | where DC = 0 This paper size applies to China National Standard (CNS) A4 (210 X 297 mm) (%) (Please read the unintentional matter on the back before filling this page) «tt tm§ ammmmw I n« ϋ — n a e, ttn I n < (47) -48- 496080 A7 B7 V. Description of the invention ( 46) RT [x] -BT [x] R, B, X [x]: (48) where DC = 0 (49) (Please read the notes on the back before filling this page) 个 呵 一 —P, TW , B, TW 1 10 丨 G, T_B, TW | Where DC = 0 These formulas printed by the Employees' Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs can be applied to each frequency vector in each section of the digital image. The calculation result of Equation 4 7-4 9 is the matrix of the number of cross-correlation, in which each unit is proportional to the cross-correlation between two frequency vectors in different channels, where the frequency vector represents the two channels in this channel. Space corresponds to the same frequency in the segment. For convenience, Equations 47-49 use a vector representation of a one-dimensional matrix of the spatial frequency domain transformation of the image section discussed in connection with FIG. Therefore, the numerator of each formula is the dot product of two vectors, and the denominator of each formula takes into account each size of these vectors. The matrix generated by Equation 4 7 is proportional to the cross-correlation between a particular segment at a particular frequency, red and green channels. Equation 4 8 produces a result that is proportional to the cross-correlation between a particular segment of the digital image at a particular frequency, the red and blue channels. Equation 4 9 produces a matrix that is proportional to the cross-correlation of a particular segment at different frequencies, between green and blue channels. We can use these three formulas to calculate a matrix for each segment of the digital image. The cross-correlation DC term is set to zero. I use the R, G ', R, B' and G 'B' notation to avoid confusion with the particle normalization formula described above. To avoid confusion, I use R ′ T, B ′ T and G ′ T notation to represent the transformation obtained in step 17 12. -49- This paper size applies to China National Standard (CNS) A4 (210 X 297 mm) 496080 B7 Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs 5. Description of the invention (47) For large in-phase frequency vectors, The results of Equations 4 7 to 4 9 tend to approach 値 1, and for large out-of-phase frequency vectors, 'approaching to 値 1 °, small vectors tend to get results close to 0, while in-phase and approximately A vector of size > of 1 produces a result of approximately 0 · 5. In step 1 7 1 6, the cross-correlation high-frequency components related to the red channel may be attenuated. Alternatively, this step can be omitted without departing from the scope of the invention. Generally, the color film used to produce a photographic image has a green layer, in which light passes through the layer, then the blue layer, and then the red layer. Because the red layer receives light after it passes through the green and blue layers, the image in the red layer is slightly blurred compared to the green and blue layers. The red layer tends to lack the high-frequency information that makes up the image information. Instead, most of the high-frequency information from the red channel tends to be noisy and can be discarded. Using different types of color film to the right 'may attenuate different channels of this digital image. In addition, if other images, such as satellite images or medical images have channels that tend to lack high-frequency image content, then the methods described below or similar methods can be used to attenuate the high-frequency components of those image channels. Equation 5 0 illustrates a red frequency attenuation that can be used in the present invention.
RFA(0,0)=0 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -50- (請先閱讀背面之注意事項再填寫本頁)RFA (0,0) = 0 This paper size applies to China National Standard (CNS) A4 (210 X 297 mm) -50- (Please read the precautions on the back before filling this page)
· ·ϋ ϋ m I n I ϋ 一°J0 i — n m —a 1 n I · 496080 A7 B7 五、發明說明(48) (請先閱讀背面之注意事項再填寫本頁) 在公式5 0中,X對應例示於圖1 8中頻率向量的二 維陣列表示中一特別的點的X座標,而y對應其y座標。 在座標(0 ’ 0 )的D C値被設爲〇。雖然,吾人可使用 對應在例示於圖1 8中頻率變換的二維配置中,每一頻率 向量位置之座標,方便地計算出紅色頻率衰減値,但是, 公式5 0的計算結果可被儲存於多個値的矩陣中,其中, 矩陣中的每一値包括:公式5 0的結果和在矩陣中的那値 的位置,其中,該位置對應在頻率變換的向量表示的一維 陣列中,對應的頻率的適當位置。例如:使用圖1 8所例 示的範例變換,R F A〔 2 8 8〕的値應是應用公式5〇 至座標(8,8)所得,而RFA〔144〕是RFA ( 0,0 )的値。 使用公式5 0計算出紅色頻率衰減値之後,這些値可 使用公式5 1而被應用至紅色和綠色波道之間的互相關與 使用公式5 2而被應用至紅色和藍色波道之間的互相關。 R,G,X2[X]=RFA[X]R,G,X[X] (51) 經濟部智慧財產局員工消費合作社印製 R,B,X2[X]=RFA[X]R,B,X[X] ( 5 2 ) 這些公式包括:每一頻率的、多數個雙數量的値。雖 然吾人已舉例一種衰減高頻紅色資訊的方法’但是’也可 使用另一種方法,而不悖離本發明的範圍1 ° 在步驟1 7 1 8中,可以計算出波道間的總互相關之 量測。這計算是對此數位式影像的每一區段而執行的’且 本纸張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -51 - 496080 A7 B7 經濟部智慧財產局員工消費合作社印製 五、發明說明(49) 是在那區段中各種不同影像波道之間的總互相關的表示。 雖然,吾人可使用任何與波道間的總互相關成正比的量測 ,但是,公式5 3 — 5 6提供一特別的區段的總互相關的 一種量測。 288 CCC1 = X(RPA[X](RT[X]GTPC])) + (BT[X]GT[X]) (53) χ-0 但限制條件爲:cccm 288 CCCN = J]RFAPC] χ-0 CCC2 = CCC1/ - /1.2 VCCCN CCC3 = —0.2-_5CCC2 1 + CCC2(0.25) 在這實施例中,公式5 3應用只考慮紅色和綠色波道 之間與綠色和藍色波道之間的互相關的一種互相關的量測 。因爲,影像通過照相軟片的其他層而到達紅色和藍色層 ,所以,紅色和藍色層之間的互相關通常並不加進任何與 類似的影像細節有關的有價値資訊。因此,該項已從公式 5 3中被省略。然而,這樣的一項可以與或不與頻率衰減 一起使用,而不悖離本發明的範圍。此外’爲了與紅色頻 率衰減器關聯的上述之理由,紅色和綠色波道之間的互相 關在公式5 3中被衰減。在這情形下’公式5 3使用在圖 1 8的例子中所用的2 8 9向量例,計算整個頻率域內的 總互相關。公式5 4反映由紅色頻率衰減所加進的累積加 櫂之量測。最後,公式5 5和5 6可被應用因而所產生的 互相關量測的範圍爲0 - 1且當影像細節不存在時,大約 (54) (55) (56) 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -52 - (請先閱讀背面之注意事項再填寫本頁) 裝· · Ϋ ϋ m I n I ϋ 1 ° J0 i — nm —a 1 n I · 496080 A7 B7 V. Description of the invention (48) (Please read the precautions on the back before filling this page) In formula 5 0, X corresponds to the X coordinate exemplified in the two-dimensional array of frequency vectors shown in FIG. 18 for a particular point, and y corresponds to its y coordinate. D C 値 at coordinates (0 '0) is set to zero. Although we can easily calculate the red frequency attenuation 値 using the coordinates corresponding to the position of each frequency vector in the two-dimensional configuration of the frequency transformation illustrated in Figure 18, the calculation result of Equation 50 can be stored in In a matrix of multiple chirps, each chirp in the matrix includes: the result of formula 50 and the position of the chirp in the matrix, where the position corresponds to a one-dimensional array represented by a vector of frequency transformation, corresponding to The proper position of the frequency. For example: using the example transformation illustrated in Figure 18, the 値 of R F A [2 8 8] should be obtained by applying formula 50 to coordinates (8, 8), and RFA [144] is the 値 of RFA (0, 0). After calculating the red frequency attenuation 値 using Equation 5 0, these 値 can be applied to the cross-correlation between the red and green channels using Equation 5 1 and between the red and blue channels using Equation 5 2 Cross-correlation. R, G, X2 [X] = RFA [X] R, G, X [X] (51) Printed by R, B, X2 [X] = RFA [X] R, B , X [X] (5 2) These formulas include: a majority of a double number of chirps per frequency. Although I have exemplified a method of attenuating high-frequency red information 'but' another method can be used without departing from the scope of the present invention 1 ° In step 1 7 1 8 the total cross-correlation between channels can be calculated Measurement. This calculation is performed for each segment of this digital image 'and this paper size applies Chinese National Standard (CNS) A4 (210 X 297 mm) -51-496080 A7 B7 Employees of the Intellectual Property Bureau, Ministry of Economic Affairs Printed by the Consumer Cooperative V. Invention Description (49) is a representation of the total cross-correlation between the various image channels in that sector. Although we can use any measurement that is proportional to the total cross-correlation between channels, Equations 5 3-5 6 provide a measure of the total cross-correlation of a particular segment. 288 CCC1 = X (RPA [X] (RT [X] GTPC])) + (BT [X] GT [X]) (53) χ-0 but the limitation is: cccm 288 CCCN = J] RFAPC] χ- 0 CCC2 = CCC1 /-/1.2 VCCCN CCC3 = --0.2-_5CCC2 1 + CCC2 (0.25) In this embodiment, Equation 5 3 is applied to consider only between the red and green channels and between the green and blue channels. A cross-correlation measurement. Because the image passes through the other layers of the photographic film to reach the red and blue layers, the cross-correlation between the red and blue layers usually does not add any valuable information about similar image details. Therefore, this term has been omitted from Equation 53. However, such an item may be used with or without frequency attenuation without departing from the scope of the invention. In addition, for the above reasons associated with the red frequency attenuator, the correlation between the red and green channels is attenuated in Equation 53. In this case, Equation 5 3 uses the 2 8 9 vector example used in the example in Fig. 18 to calculate the total cross-correlation in the entire frequency domain. Equation 5 4 reflects the measurement of the cumulative plus 棹 added by the red frequency attenuation. Finally, equations 5 5 and 5 6 can be applied. The resulting cross-correlation measurement range is 0-1 and when the image details are not present, about (54) (55) (56) This paper scale applies Chinese national standards (CNS) A4 size (210 X 297 mm) -52-(Please read the precautions on the back before filling this page)
I n n n^eJI βϋ n I n n ψί i I 496080 A7 B7 五、發明說明(50) 爲〇.2 。 在步驟1 7 2 0中,根據在步驟1 7 1 2中所算出的 原始變換和每個波道與其他波道間的互相關’可計算出該 波道的每個區段的一加權後的變換。吾人可使用任何加權 後的變換,以回應此數位式影像的一或多個波道間之互相 關的任何量測,而不悖離本發明的範圍。在這實施例中’ 可使用公式5 7 — 5 9來計算加權後的變換。 RT2[X] = |RT:[Xj + 〇·9 R’GX2[X]+ OJRBXZfX] 但限制條件爲R’T2[XP0 (57) DC = 0 GT2[X] = |GT[X]卜 〇.9G'B'X[X] + OJR'GpfX] 但限制條件爲G’T2[X]>〇 (58) DC = 0 BT2PC] = |BT[X]| + O^GIB^PC] + OJRO^pC] 但限制條件爲B,T2[XP〇 (59) DC = 0 在每個情形中,加權後的變換的D C項被設爲0。此 外,加權後變換的每一成分也被限制在某個大於等於〇的 値。注意··公式5 7 - 5 9只對頻率變換的大小部份運算 ,且對每一·變換產生一修改後的大小。 在步驟1 7 2 2中,於步驟1 7 2 0所得的每個波道 的每一區段的加權後變換可被過濾。替換地’這步驟可被 省略,而不悖離本發明的範圍。此外’若步驟1 7 2 0被 省略,則步驟1 7 1 2所得的變換可對每一波道而被過濾 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) (請先閱讀背面之注意事項再填寫本頁) 裝--------訂------,---遍 經濟部智慧財產局員工消費合作社印製 -53- 496080 A7 ____Β7____五、發明說明(51 ) 。每一變換可以能強調訊號內容和衰減影像中的雜訊的方 式,而被過濾。在此方法的這步驟中,可以使用一低通濾 波器較有利。適合的濾波器趨於在隨機雜訊中,得到平順 的衰減與降低。因此,雖然可以使用任何低通濾波器’但 是,一適合之濾波器可能是用來提高存在於變換中的有效 訊號和抑制存在於變換中的雜訊的理想濾波器。一 R M S 計算可被用來進一步加重有效的訊號。此濾波運算趨於保 留在頻率中強的變化且抑制在頻率中小的變化。換言之, 可加重極可能代表有效的訊號的強頻率成分。 在這實施例中,濾波器包含一 3 X 3矩陣’其中該矩 陣在中心的値爲4、在中心的水平左右或垂直上下點的値 爲2、且在矩陣的4個角的値爲1。其他大小的濾波器或 其他濾波器値可被使用,而不悖離本發明的範圍。若w ( a ,b )代表上述的濾波器,則可根據公式5 8 — 6 2 ’ 計算出被過濾的變換R ’ T 3、G ’ T 3和B ’ T 3。 --------γ----裝--------訂--------- (請先閱讀背面之注意事項再填寫本頁) RT3[x,y]=亡亡 W(a,b)(RT2[X + a,y + bj)2 •-丨 (60) 經濟部智慧財產局員工消費合作社印製I n n n ^ eJI βϋ n I n n ψί i 496080 A7 B7 5. The description of the invention (50) is 0.2. In step 17 2 0, a weighted value for each section of the channel can be calculated based on the original transformation calculated in step 1 7 1 2 and the cross-correlation between each channel and other channels. Transform. We can use any weighted transformation in response to any measurement related to one or more channels of this digital image without departing from the scope of the invention. In this embodiment, equations 5 7-5 9 can be used to calculate the weighted transformation. RT2 [X] = | RT: [Xj + 〇 · 9 R'GX2 [X] + OJRBXZfX], but the restriction is R'T2 [XP0 (57) DC = 0 GT2 [X] = | GT [X]. .9G'B'X [X] + OJR'GpfX] but the restriction is G'T2 [X] > 〇 (58) DC = 0 BT2PC] = | BT [X] | + O ^ GIB ^ PC] + OJRO ^ pC] But the restriction is B, T2 [XP〇 (59) DC = 0 In each case, the weighted transformed DC term is set to 0. In addition, each component of the weighted transformation is also limited to a value 値 greater than or equal to zero. Note ·· Equation 5 7-5 9 calculates only the size part of the frequency transform, and generates a modified size for each transform. In step 17 2 2, a weighted post-transformation of each section of each channel obtained in step 17 20 may be filtered. Alternatively, this step may be omitted without departing from the scope of the present invention. In addition, 'If step 1 7 2 0 is omitted, the transformation obtained in step 1 7 1 2 can be filtered for each channel. This paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) (please Please read the notes on the back before filling this page.) -------- Order ------, --- Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs -53- 496080 A7 ____ Β7 ____ Description of the invention (51). Each transform may be filtered to emphasize signal content and attenuate noise in the image. In this step of the method, it may be advantageous to use a low-pass filter. Suitable filters tend to get smooth attenuation and reduction in random noise. Therefore, although any low-pass filter can be used, a suitable filter may be an ideal filter for improving the effective signal existing in the transformation and suppressing the noise existing in the transformation. An R M S calculation can be used to further emphasize the effective signal. This filtering operation tends to keep strong changes in frequency and suppress small changes in frequency. In other words, the aggravable strong frequency components that are likely to represent an effective signal. In this embodiment, the filter contains a 3 X 3 matrix, where 値 in the center is 4; 在 in the horizontal left and right or vertical points in the center is 2; and 値 in the 4 corners of the matrix is 1. . Other size filters or other filters may be used without departing from the scope of the invention. If w (a, b) represents the above-mentioned filter, the filtered transforms R'T3, G'T3, and B'T3 can be calculated according to the formula 5 8-6 2 '. -------- γ ---- Installation -------- Order --------- (Please read the precautions on the back before filling this page) RT3 [x, y ] = 亡 死 W (a, b) (RT2 [X + a, y + bj) 2 •-丨 (60) Printed by the Consumer Cooperative of Intellectual Property Bureau, Ministry of Economic Affairs
1 I RT3[DC] = 0 〇T3[x,y]: I;XW(a,b)GT2[xH.a,y^b]2 a·*·! b*—1 a*·—I b»— GT3PC]: (61) 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -54 - 496080 Α7 Β7 經濟部智慧財產局員工消費合作社印製 五、發明說明(52 ) ZZW(a,b)BT2[x + a,y + b]2 B'T3[x, y]=卜1 ( v- (62)| ZZW(^b) I a —lb—l BT3[DC] = 0 如所注意到的:在每一被過濾的變換中’ D c項被設爲〇 。此外,在此變換的邊緣之外的元件被視爲其値爲0 °在 這公式中,這個處置並不明顯的表示出。此濾波運算只有 影響頻率向量的大小。每一變換的頻率向量之相位維持不 變。在公式6 0 — 6 2中,吾人使用變換的二維表示來計 算。爲了記號法的目的,使用變換的大小之事實由在公式 6 0 - 6 2中使用平常字母(與粗體字母相反)作爲變換 r, T2、R, T3、G, T2、G’ 丁 3、B’ T2 和 Β ’ Τ 3的事實所指示出。 在步驟1 7 2 4中,每個波道的每一區段的被過濾的 變換可被頻率展開。替換地,若步驟1 7 2 2被省略’則 在步驟1 7 2 0中所得的變換或在步驟1 7 1 2中所得的 變換,可在步驟1 7 2 4中被頻率展開。步驟1 了 2 4也 可以被省略,而不悖離本發明的範圍。 因爲真實世界是由類似的重複圖樣所組成的’所以’ 一照相影像的變換空間趨於在D C頻率周圍的任何放大’ 具有相同的外表,除非存在重製的人工製品。對重複的圖 樣而言,一影像的頻率特性趨於以2次方重複。在這步驟 中,在較低頻的頻率特性被用來以雙倍頻率,預先使門檻 細節易於感光。爲了執行這樣的頻率展開,若吾人使用圖 (請先閱讀背面之注意事項再填寫本頁) 裝 訂------^--- %- 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公爱) -55- 496080 A7 B7 五、發明說明(53) 1 8所闡述的變換之二維版本,則此變換資料的中心〜半 被擴大成變換資料的全尺寸區塊。換言之,頻率的中心9 X Ό區塊被展開爲在1 7 X 1 7變換中的每一其他位置。 當一頻率在變換的剩餘部份中,並未明確地被塡入時,M[J 可使用鄰近頻率的平均値。再次地,這運算只有影響頻$ 大小。甚至沒有在變換的中央一半的,大小已被改變的頻 率之既有的相位資訊也被保留。 在得出被過濾的變換之頻率展開之後,可在步驟 1 7 2 6中,對每一波道,計算一新的變換。此新的變換 可以取決於該波道的頻率展開變換和在步驟1 7 2 2所得 的該波道的過濾後變換。替換地,若步驟1 7 2 2被省略 ,則新的變換取決於頻率展開變換和原來在步驟1 7 1 2 所得的變換。公式6 3 - 6 5可被使用來計算此新的變換 R’ T4、G, T4 和 B’ 丁 4。 RT4[X] = V〇.8RT3[X]2 +0.2R^[X]2 (63) GTT4[X] = V〇.8GT3[X]2 +0.2G'ET[X]2 (64) B,T4[X] = VoJBTSpq2 +0·2ΒΈΤ[Χ]2 (65) 在公式63-65中,變換R’ ET、G, ET和 B ’ E T代表在步驟1 7 2 4中,由頻率展開所得的估計 的變換。再次地,公式6 3 - 6 5所執行的運算只有影響 頻率向量的大小。相位資訊被保留。 在步驟1 7 2 8中,吾人可產生一遮罩’使得在步驟 1 7 2 6所得的變換之高頻項減幅。因爲顆粒軌跡趨於具 有主要的高頻成分,所以,該遮罩可被用來使穎粒軌跡和 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) ----、----1----裝— (請先閱讀背面之注音?事項再填寫本頁) 訂: 經濟部智慧財產局員工消費合作社印製 -56- 496080 A7 B7 五、發明說明(54) /或其他高頻雜訊減幅。 可使用由公式4 2所表示的選通函數得到此遮罩。然 後便用公式4 3 ,以得出該遮罩。公式4 2和4 3在上文 中,與圖1 6相關聯地被說明。其他遮罩函數和/或選通 函數可被使用,而不悖離本發明的範圍。此外,這步驟可 被省略,而不1孛離本發明的範圍。 在步驟1 7 3 0中’可以㈢每個波道的每一區段計算 一雜訊校正函數。在這實施例中,爲回應遮罩和先前所計 算出的影像的頻率域表示之一,計算此雜訊校正函數。例 如:在步驟1 7 2 6所得的變換可被用來得出此雜訊校正 函數。可由公式6 6 - 7 1得出雜訊校正函數。在這些公 式中的B ( X,y )得自於公式4 3。 CR(x,y):—M〇._f C3)B(x,y) c“1 +‘㈣ (請先閱讀背面之注意事項再填寫本頁) (66) CR(x, y) = l + F.atc(x,y) CR >1 + F邮(x,y) CG(x,y) = 2Ifei3㈣c〇,1+, (x,y) 2 (67) C〇(x,y) = l + Fgale(x,y) cg >1 + Fgatc(x,y) 經濟部智慧財產局員工消費合作社印製 CB (X, y) = (1 : 0-4CCC3) Β(χ〇〇 d(x,y) CB(x,y) = l + F.atc(x,y) 2 cb >l+Wx,y) RTMIXyl^G^ +(1-GR)CR 〇c,y) GTM[x,y] = GR +(l-GR)C〇(x,y) BTM[x,y] = GR +(KR)C“x,y) (68) (69) (70) (71) -57- 本纸張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) 496080 A7 B7 經濟部智慧財產局員工消費合作社印製 五、發明說明(55 ) 在公式6 6 — 6 8中’雜訊校正函數被限制在1加上 在變換矩陣中一特別的位置之選通函數的値。再次地,公 式;6 6 — 6 8和6 9 — 了 1只有影響頻率向量的大小分量 。原始的相位資訊可被保留。本發明的使用者可能想要在 根據本發明的改良後的影像上存留有殘餘的顆粒性。因此 ,本發明允許使用者指定在影像被改良後’其所欲留下的 顆粒殘餘GR。公式6 9 - 7 1在計算雜訊校正函數 R’ TM、G’ TNI和B’ TM時,將所要的殘留顆粒考 慮進去。 在步驟1 7 3 2中,每個波道的每一區段之雜訊校正 函數被應用至適當的波道的適當的區段之原始變換。其結 果爲如··顆粒軌跡之高頻雜訊被抑制後的影像之頻率域表 示。使用公式7 2 - 7 4可得出被抑制顆粒的頻率表示。 R,TS[X]=R’T[X]R,TM[X] (72) G,TS[X]:G’T[X]G’TM[X] (73) B’TS[X]=B,T[X]B,TM[X] (74) 再次地,公式7 2 - 7 4所執行的運算只有影響頻率 向量的大小,而相位資訊維持不變。 在步驟1 7 3 4中,影像細節可被加入至影像的紅色 波道的每一區段的頻率域表示中。如先前所注意者,光線 通常在光線穿過彩色軟片的綠色和藍色層後,才到達軟片 的紅色層。因此,影像的紅色部分中的細節趨於變成模糊 (請先閱讀背面之注意事項再填寫本頁) 裝 —訂------;— 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) 58 496080 A7 __B7 ____ 五、發明說明(56 ) (請先閱讀背面之注意事項再填寫本頁) 不淸。在步驟1 7 3 4中,吾人可使用綠色和藍色波道的 資訊,將這樣的細節回復至影像的紅色波道。注意:步驟 1 7 3 4可在任何數位式影像的的紅色波道上被執行,不 僅止於應用圖1 7的步驟後的紅色波道上。此外,步驟 1 7 3 4可被省略,而不悖離本發明的範圍。此外,單獨 的步驟1 7 3 4可包括本發明的一實施例。因爲影像的綠 色波道趨於具有最多的影像細節,故其可被用來將細節加 至其他波道上。步驟1 7 3 4所述的技術也可被應用在: 將細節加至不同型式的一不同的波道上或結構上不同於典 型的彩色軟片的彩色軟片所得的數位式影像上。也可使用 其他改良紅色影像細節的方法,而不悖離本發明的範圍。 在這實施例中,公式7 5可被用來對紅色波道的每一區段 計算該區段的修改後的紅色波道變換。 (RTSPC]RFAPC> RTSPC]= 1 1-RFA[X]+RFA[X] (75) ,··· (1 - RFA[X]) (0.7GTS[X]+〇.4B'TS[X])) 經濟部智慧財產局員工消費合作社印製 公式7 5使用與步驟1 7 1 6關聯的上述之紅色頻率 衰減器。 在步驟1 7 3 6中,雜訊校正變換的反變換爲 R’ TS、G’ TS和B’ TS,其中,這些變換可對每 一區段取得,以得出已被抑制雑訊後的改良影像的空間域 表示。若在步驟1 7 3 4中’影像細節已被加進來’則使 用公式7 5所得的紅色頻率變換可以被應用。 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) _ 59 - 496080 A7 B7 五、發明說明(57 ) 平衡彩色波道 (請先閱讀背面之注意事項再填寫本頁) 歸一化顆粒軌跡後,經由校正陣列所得被歸一化的値 可被使用在平衡數位式影像的色彩內容,如圖3的步驟 1 0 1 4所示。替換地,也可以不使用被歸一化的値而平 衡色彩。 一數位式影像的原始値常常包含來自掃描器嘗試改良 該數位式影像的失真和由於軟片因放置時間太久在原有軟 $ 片上影像的失真。這些失真可以由不再反映原有軟片的彩 色層中的顆粒的隨機和平均分佈的色彩圖樣所組成。藉由 使用圖3中步驟1 0 0 8和1 0 1 0最後所得的穎粒軌跡 之被歸一化的値,(或藉其他方法)本發明也可產生原始 色彩値的一估計的重顯(r e c ο n s t r u c t i ο η )。這些被重顯 的値,在它們經由歸一化的數位式影像比原有數位式影像 更平均地分佈的意義上,可被稱爲平衡的。換言之,平衡 色彩內容可以更加改良數位式影像。 經濟部智慧財產局員工消費合作社印製 爲了平衡色彩內容,吾人找出每一波道,在數位式影 像的歸一化的値中最暗的像素強度値,即最低値。也找出 每一波道,在數位式影像的歸一化的値中最明亮的像素強. 度値,即最高値。在這技術上的改良爲:使用來自數位式 影像的一低通版本之最低和最高値。所找出對這方法有利 的一低通過濾器爲具有通常超過3個像素的半徑的一中位 數過濾器。替換地,以像素爲單位,該半徑可以與數位式 影像的最大線性尺寸的平方根成正比。此軟體於是重分佈 像素強度値,以配合特殊的比特深度(Ιμ t depth)(強度 -60- 本纸張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) 496080 A71 I RT3 [DC] = 0 〇T3 [x, y]: I; XW (a, b) GT2 [xH.a, y ^ b] 2 a · * ·! B * —1 a * · —I b »— GT3PC]: (61) This paper size is in accordance with Chinese National Standard (CNS) A4 (210 X 297 mm) -54-496080 Α7 Β7 Printed by the Consumers’ Cooperative of Intellectual Property Bureau of the Ministry of Economic Affairs 5. Description of the invention (52) ZZW (a, b) BT2 [x + a, y + b] 2 B'T3 [x, y] = Bu 1 (v- (62) | ZZW (^ b) I a —lb—l BT3 [DC] = 0 As noted: In each filtered transformation, the 'D c term is set to 0. In addition, components outside the edges of this transformation are considered as having 値 0 °. In this formula, this The treatment is not clearly shown. This filtering operation only affects the size of the frequency vector. The phase of each transformed frequency vector remains unchanged. In Equations 60-62, we use the transformed two-dimensional representation to calculate. To For the purpose of the notation, the fact that the size of the transformation is used is derived from the use of ordinary letters (as opposed to bold letters) as the transformations r, T2, R, T3, G, T2, G 'D3, B in Equations 60-62. 'T2 and B' T 3. This is indicated by the fact that in step 1 7 2 4 The filtered transformation of each section of the channel may be expanded by frequency. Alternatively, if step 17 2 2 is omitted, then the transformation obtained in step 17 20 or the transformation obtained in step 17 12 The transformation can be frequency expanded in steps 1 7 2 4. Steps 1 2 and 4 can also be omitted without departing from the scope of the present invention. Because the real world is composed of similar repeated patterns 'so' a photo The transformation space of an image tends to have any appearance around the DC frequency ', with the same appearance, unless there are duplicate artifacts. For repeated patterns, the frequency characteristics of an image tend to be repeated to the second power. In this step In the lower frequency, the frequency characteristics are used to double the frequency, to make the threshold details easy to be sensitive. In order to perform such frequency expansion, if I use the figure (please read the precautions on the back before filling this page) Binding ----- ^ ---%-This paper size is in accordance with Chinese National Standard (CNS) A4 (210 X 297 public love) -55- 496080 A7 B7 V. The transformation described in (53) 1 8 Two-dimensional version, then the center of this transformation data ~ It is half-expanded into a full-size block of transform data. In other words, the center of the frequency 9 X Ό block is expanded to every other position in the 1 7 X 1 7 transform. When a frequency is in the remainder of the transform, When not explicitly entered, M [J can use the average chirp of adjacent frequencies. Again, this operation only affects the frequency $ size. Existing phase information that is not even at the center half of the transformation, whose frequency has been changed, is retained. After the frequency expansion of the filtered transform is obtained, a new transform can be calculated for each channel in step 1726. This new transformation may depend on the frequency expansion transformation of the channel and the filtered transformation of the channel obtained in step 172. Alternatively, if step 17 2 is omitted, the new transformation depends on the frequency expansion transformation and the transformation obtained in step 1 7 1 2. Equations 6 3-6 5 can be used to calculate this new transformation R 'T4, G, T4 and B' Ding 4. RT4 [X] = V〇.8 RT3 [X] 2 + 0.2R ^ [X] 2 (63) GTT4 [X] = V〇.8GT3 [X] 2 + 0.2G'ET [X] 2 (64) B , T4 [X] = VoJBTSpq2 + 0 · 2ΒΈΤ [Χ] 2 (65) In Equation 63-65, transform R 'ET, G, ET, and B' ET to represent the frequency expansion obtained in step 1 7 2 4 The estimated transformation. Again, the operations performed by Equations 6 3-6 5 only affect the size of the frequency vector. Phase information is retained. In step 1728, we can generate a mask 'to reduce the high-frequency terms of the transformation obtained in step 1726. Because the particle trajectory tends to have the main high-frequency component, this mask can be used to make the particle trajectory and the paper size apply the Chinese National Standard (CNS) A4 specification (210 X 297 mm) ----, ---- 1 ---- Equipment— (Please read the phonetic on the back? Matters before filling out this page) Order: Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs-56- 496080 A7 B7 V. Description of Invention / Or other high frequency noise reduction. This mask can be obtained using a gating function represented by Equation 42. Then use Equation 4 3 to get the mask. Equations 4 2 and 4 3 are explained above in connection with FIG. 16. Other mask functions and / or gating functions may be used without departing from the scope of the invention. In addition, this step can be omitted without departing from the scope of the present invention. In step 1730, a noise correction function can be calculated for each section of each channel. In this embodiment, the noise correction function is calculated in response to one of the mask and one of the frequency domain representations of the previously calculated image. For example, the transform obtained in step 1 7 2 6 can be used to derive this noise correction function. The noise correction function can be obtained from Equation 6 6-7 1. B (X, y) in these formulas is derived from Formula 43. CR (x, y): —M〇._f C3) B (x, y) c “1 + '㈣ (Please read the notes on the back before filling this page) (66) CR (x, y) = l + F.atc (x, y) CR > 1 + Fpost (x, y) CG (x, y) = 2Ifei3㈣c〇, 1 +, (x, y) 2 (67) C〇 (x, y) = l + Fgale (x, y) cg > 1 + Fgatc (x, y) Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economy CB (X, y) = (1: 0-4CCC3) Β (χ〇〇d (x, y) CB (x, y) = l + F. atc (x, y) 2 cb > l + Wx, y) RTMIXyl ^ G ^ + (1-GR) CR 〇c, y) GTM [ x, y] = GR + (l-GR) C〇 (x, y) BTM [x, y] = GR + (KR) C "x, y) (68) (69) (70) (71)- 57- This paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) 496080 A7 B7 Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs 5. Description of the invention (55) In formulas 6 6 — 6 8 'The noise correction function is limited to 1 plus the gating function at a particular position in the transformation matrix. Again, the formulas; 6 6 — 6 8 and 6 9 — 1 only affect the magnitude components of the frequency vector. The original phase information can be retained. The user of the present invention may want to leave residual graininess on the improved image according to the present invention. Therefore, the present invention allows the user to specify the grain residual GR that he wants to leave after the image is improved. Equation 6 9-7 1 In calculating the noise correction functions R 'TM, G' TNI, and B 'TM, the desired residual particles are taken into account. In step 172, the noise correction function of each section of each channel is applied to the original transformation of the appropriate section of the appropriate channel. The result is the frequency domain representation of the image after the high-frequency noise of the particle trajectory is suppressed. Using equations 7 2-7 4 gives a frequency representation of the suppressed particles. R, TS [X] = R'T [X] R, TM [X] (72) G, TS [X]: G'T [X] G'TM [X] (73) B'TS [X] = B, T [X] B, TM [X] (74) Again, the operations performed by Equations 7 2-7 4 only affect the size of the frequency vector, while the phase information remains unchanged. In step 1 7 3 4 the image details can be added to the frequency domain representation of each segment of the red channel of the image. As previously noted, light usually reaches the red layer of the film after it passes through the green and blue layers of the color film. Therefore, the details in the red part of the image tend to become blurred (please read the precautions on the back before filling out this page). Binding — Binding ------; — This paper size applies the Chinese National Standard (CNS) A4 specification 210 X 297 mm) 58 496080 A7 __B7 ____ 5. Description of the invention (56) (Please read the notes on the back before filling this page) Yes. In step 1 7 3 4 we can use the information of the green and blue channels to restore such details to the red channel of the image. Note: Steps 1 7 3 4 can be performed on the red channel of any digital image, not just the red channel after applying the steps of Figure 17. Furthermore, steps 1 7 3 4 can be omitted without departing from the scope of the invention. Furthermore, the individual steps 1 7 3 4 may include an embodiment of the invention. Because the green channel of an image tends to have the most image detail, it can be used to add detail to other channels. The technique described in step 1 7 3 4 can also be applied to: digital images obtained by adding detail to a different channel of a different type or a color film that is structurally different from a typical color film. Other methods of improving the details of the red image can be used without departing from the scope of the invention. In this embodiment, Equation 75 can be used to calculate, for each section of the red channel, the modified red channel transform for that section. (RTSPC] RFAPC > RTSPC] = 1 1-RFA [X] + RFA [X] (75), (1-RFA [X]) (0.7GTS [X] + 〇.4B'TS [X] )) The Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs printed the formula 7 5 using the above-mentioned red frequency attenuator associated with step 1 7 1 6. In step 1 7 3 6, the inverse transform of the noise correction transform is R 'TS, G' TS, and B 'TS, where these transforms can be obtained for each segment to obtain the Improved spatial domain representation of images. If 'Image detail has been added' in step 1 7 3 4 then the red frequency transform obtained using Equation 7 5 can be applied. This paper size applies to China National Standard (CNS) A4 (210 X 297 mm) _ 59-496080 A7 B7 V. Description of the invention (57) Balanced color channel (please read the precautions on the back before filling this page) After the particle trajectories are normalized, the normalized chirp obtained through the correction array can be used to balance the color content of the digital image, as shown in step 1 0 1 4 of FIG. 3. Alternatively, the colors can be balanced without using the normalized chirp. The original frame of a digital image often contains distortion from the scanner's attempt to improve the digital image and distortion of the original image on the film due to the film being left too long. These distortions can consist of random and evenly distributed color patterns that no longer reflect particles in the color layer of the original film. By using the normalized 颖 of the glomerular trajectory finally obtained in steps 1 0 8 and 1 0 1 0 in FIG. 3, (or by other methods) the present invention can also generate an estimated re-display of the original color 値(Rec ο nstructi ο η). These reappeared puppets can be said to be balanced in the sense that they are more evenly distributed through the normalized digital image than the original digital image. In other words, balancing color content can further improve digital images. Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs In order to balance the color content, we find each channel, the darkest pixel intensity in the normalized frame of digital images, that is, the lowest. Also find the brightest pixel intensity of each channel in the normalized frame of a digital image. The degree is the highest frame. The technical improvement is to use the lowest and highest thresholds from a low-pass version of the digital image. A low-pass filter found to be advantageous for this method is a median filter with a radius of usually more than 3 pixels. Alternatively, the radius may be proportional to the square root of the maximum linear dimension of the digital image in pixels. This software then redistributes the pixel intensity 値 to match the special bit depth (Ιμ t depth) (intensity -60- this paper size applies to the Chinese National Standard (CNS) A4 specification (210 X 297 mm) 496080 A7
五、發明說明(58 ) 範圍),其中該深度適合本發明的特殊實施例之硬體。 例如:下表代表來自一數位式影像的一估計重顯之歸 一化的像素強度値之範圍:R 3 · 5 ____2 0 · 3 42.6 _ 2 7.910.1 H5. The scope of the invention description (58)), wherein the depth is suitable for the hardware of the special embodiment of the present invention. For example: The following table represents the range of normalized pixel intensities of an estimated reappearance from a digital image: R 3 · 5 ____2 0 · 3 42.6 _ 2 7.910.1 H
G R 4 爲了改良這些値的可見顯示,可使用利用非低通資料 每一波道的如下方程式: 平衡後的値養的4直最=像素素), 整個舉例的平衡後的値於是爲: R 〇__0 __ 〇GR 4 In order to improve the visible display of these chirps, the following equations can be used for each channel using non-low-pass data: balanced 4 値 maximum = pixel pixels). The balanced example for the entire example is then: R 〇__0 __ 〇
G R .2 5 5 .2 5 5 2 5 5 (請先閱讀背面之注意事項再填寫本頁) 一裝 經濟部智慧財產局員工消費合作钍印製 也可使用色彩改良或校正在相關訊號之間的非線性失 真的其他方法,其可具或不具有顆粒歸一化和/或衰減。 任何這樣的方法可被使用,而不悖離本發明的範圍。 圖19例示說明可被用來將具有一共用的成分的複數 個訊號歸一化的方法的步驟之流程圖,其中,複數個訊號 中的至少一個被以非線性的方式使之失真。雖然,本發明 的這方面可被用於改良包括一被掃描的照相影像之數位式 影像的色彩,圖1 9所述的方法可被使用在許多額外的應 用上。例如:它可被用於校正具有共用的成分的複數個訊 號中的非線性失真。因此,它也可被用於諸如:聲頻應用 、醫學影像應用、衛星影像應用等等。圖1 9所述的方法 —訂------Λ---. 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公爱) -61 - 496080 A7 B7 _ 五、發明說明(59 ) 可被用於這些應用的任一種’而不悖離本發明的範圍。 (請先閱讀背面之注意事項再填寫本頁) 在步驟1 9 1 0中,複數個訊號之一的失真相對於其 他訊號而被測量。一失真函數可對一訊號而產生’其中’ 該訊號與那訊號相對剩餘訊號的至少一個和此複數個訊號 的失真成正比。雖然,可相對於剩餘相關的訊號的僅只一 個測量失真,但是,相對於多個或所有剩餘的訊號量測失 真可能較有利。如上述,包含三個色彩波道的一數位式影 像,一波道的失真可相對其他兩色彩波道而測量。 在步驟1 9 1 2中,可根據在步驟1 9 1 0中所量測 的失真,產生一反失真函數。這樣的反失真函數可以對具 有一共用的成分的複數個訊號中的一個、某些、或全部訊 號而產生。 在步驟1 9 1 4中,可使用在步驟1 9 1 2中所計算 出的反失真函數,歸一化一訊號。再次地,可根據在步驟 1 9 1 2中所計算出的反失真函數歸一化複數個訊號中的 一個、某些、或全部訊號。 經濟部智慧財產局員工消費合作社印製 在本發明的某些實施例中,可能需要在決定步驟 1 9 1 0中的失真函數之前過濾此複數個訊號中的一或多 個。這樣的過濾在雜訊和/或其他效應可能和一正確的失 真函數的計算干涉時,是需要的。例如:移除失真時,一 中位數型的濾波器可能有用。 可能需要被使用在校正複數個訊號中的非線性失真, 其中,每一訊號或多個訊號已經以非線性方式使之失真, 但是,其中,決定一個、某些、或所有訊號相對於彼此的 -62- 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) 經濟部智慧財產局員工消費合作社印製 496080 A7 B7___ 五、發明說明(6〇 ) 失真是困難的。對這樣相關的訊號而言,可藉由使此複數 個訊號在訊號中的小步相對地有等量的變化,而校正失真 。對其爲一單一變數的函數之訊號而言,步驟1 9 1 0中 所計算出的失真函數可取決於複數個訊號中的一個 '某些 、或所有訊號在影像中多個點上之斜率。在這情形下’在 步驟1 9 1 2中所產生的反失真函數將試圖使複數個訊號 的斜率在訊號中小增量上相等。對其爲多個變數的函數之 訊號而言,失真函數可取決於複數個訊號中的一個、某些 、或所有訊號在訊號中多個點上多個方向之變化率。例如 :在步驟1 9 1 0中所計算出的失真函數可取決於複數個 訊號中的一個、某些、或所有訊號多個點上之梯度( gradient ) 0 其中,一梯度量測被用來決定一訊號相對其他訊號的 失真,一散佈圖可被用來產生此失真函數。在散佈圖上的 每一點可藉由計算一散佈圖指標而得出,以回應在此訊號 中的共同點上,複數個訊號中的某些訊號的梯度大小。一 加權可被計算,以回應在這些梯度之間的角度(例如:在 梯度之間的角度之餘弦)。此外,這樣的加權可以與梯度 的大小成正比。可使用曲線擬合演算法將一曲線擬合至散 佈圖上的點,且使用此曲線計算出一反失真函數。例如: 此反失真函數可能與曲線下的面積成正比。 圖2 0例示說明根據本發明,可被用來改良一數位式 影像的一方法之流程圖。本發明的這方面可被用於改良一 數位式影像的色彩。數位式影像,無論其爲被掃描的照相 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -63- —1 ----•裝--------訂------,— (請先閱讀背面之注意事項再填寫本頁) 496080 A7 B7 五、發明說明(61 ) (請先閱讀背面之注意事項再填寫本頁) 影像、人造衛星影像、醫學影像或其他型式的數位式影像 ,由於在類比影像中的失真和/或由於在產生數位式影像 時發生的失真,可能被以非線性方式使之失真。此失真可 能關於強度和/或頻率而變化,且可能對每一波道而不同 地變化。這樣的非線性失真可能以吾人不想要的方式,將 影像降級。對照相影像而言,影像可能看起來不那麼合意 。對人造衛星影像和醫療影像而言,由影像所提供的資訊 ,由於這樣的失真,可能不那麼有用。 經濟部智慧財產局員Η消費合作社印製 如上所述,許多數位式影像包括多個波道。例如:以 一彩色相片爲根據的數位式影像一般包含一紅色波道、一 綠色波道和一藍色波道。其中,一影像是由多個波道所組 成的,包含在每一波道的影像資訊在此三個影像波道之間 產生一共用的成分。因此,在此數位式影像中的小區域上 ,強度趨於在此影像中跨小區域,以等量而變化。雖然, 強度的大小可能相當不同,對一未失真的影像而言,這樣 大小的變化率趨於相似。因此,本發明嘗試根據對一影像 的每一波道(或波道的子集合),在多個方向上強度的變 化,來決定影像的失真。一反失真函數可對一影像中的一 個、某些、或所有波道而產生,以致,當應用這樣的反失 真函數時,每一波道的強度趨於等量的變化,尤當與其他 波道在影像的小區域上比較時。 在某些影像的區域中,當影像未失真時,與其他波道 比較,一波道的強度値趨於線性地變化,然而,在其他影 像中或同一影像的區域中,當影像未失真時,與其他波道 -64- 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公爱) 經濟部智慧財產局員工消費合作社印製 496080 A7 ____B7 _ 五、發明說明(62 ) 比較,一波道的強度値趨於在一小區域上維持大小之間的 一常數比例。爲了自動化彩色數位式影像的改良,本發明 利用這兩種變化型式的平均値,其中,該變化對應至一彩 色影像趨於被數位式地儲存之方式。 在步驟1 9 1 6中,一影i象被掃描以產生一數位式影 像。這樣的掃描可使用上述的掃描器3 0或掃描器3 4而 發生。替換地,步驟1 9 1 6可被省略且本發明被應用至 先前已被掃描出的一影像。此外,本發明可被應用至具有 多個波道的數位式影像,其中,每一波道具有以一非線性 方式,使一或多個波道失真的一共同訊號成分。”共同成 分”一詞廣泛地指以某種方式,互相關的訊號。 可選擇地,在執行圖2 0所述的方法中之剩餘步驟前 ,以上所討論的本發明之顆粒歸一化方面可被應用至一數 位式影像。可選擇地,在執行圖2 0所述的方法中之剩餘 步驟前,上述的一或多種顆粒減少技術也可被應用至影像 。替換地,圖2 0所述的方法可被應用至影像,無須任何 顆粒歸一化和/或顆粒減少,或可使用不同的技術,在歸 一化或減少之後行之。 在步驟1 9 1 8中,影像的每一波道(或影像的波道 之子集合)可被過濾,因而,各種不同的效應,諸如:雜 訊和影像解析度效應,不會與非線性失真的移除相干涉。 在一小數目像素中,諸如··由穎粒所引起的雜訊可能與非 線性失真的移除相干涉。此外,因爲用來產生一彩色照相 影像的光線,通過軟片的多數層,所以’數位式影像的一 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -65 · " (請先閱讀背面之注意事項再填寫本頁) 裝GR .2 5 5 .2 5 5 2 5 5 (Please read the precautions on the back before filling out this page) The consumer cooperation of the Intellectual Property Bureau of the Ministry of Economic Affairs can be used to print or use color improvement or correction between related signals. Other methods of non-linear distortion, which may or may not have particle normalization and / or attenuation. Any such method can be used without departing from the scope of the invention. Figure 19 illustrates a flowchart illustrating the steps of a method that can be used to normalize a plurality of signals having a common component, wherein at least one of the plurality of signals is distorted in a non-linear manner. Although this aspect of the invention can be used to improve the color of a digital image including a scanned photographic image, the method described in Figure 19 can be used in many additional applications. For example, it can be used to correct non-linear distortion in a plurality of signals having a common component. Therefore, it can also be used in applications such as: audio applications, medical imaging applications, satellite imaging applications, and so on. The method described in Figure 19—Order ------ Λ ---. This paper size applies to China National Standard (CNS) A4 (210 X 297 public love) -61-496080 A7 B7 _ V. Description of the invention (59) Can be used for any of these applications' without departing from the scope of the invention. (Please read the precautions on the back before filling this page) In step 1910, the distortion of one of the multiple signals is measured relative to the other signals. A distortion function can generate 'where' the signal is proportional to the distortion of at least one of the signal relative to the remaining signal and the plurality of signals. Although only one measurement can be distorted relative to the remaining correlated signals, it may be advantageous to measure the distortion relative to multiple or all remaining signals. As described above, for a digital image containing three color channels, the distortion of one channel can be measured relative to the other two color channels. In step 1912, an inverse distortion function can be generated based on the distortion measured in step 1910. Such an inverse distortion function can be generated for one, some, or all of a plurality of signals having a common component. In step 1914, a signal can be normalized using the inverse distortion function calculated in step 1912. Again, one, some, or all of the plurality of signals can be normalized according to the inverse distortion function calculated in step 1912. Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs In some embodiments of the present invention, it may be necessary to filter one or more of the plurality of signals before determining the distortion function in step 1910. Such filtering is needed when noise and / or other effects may interfere with the calculation of a correct distortion function. For example, a median filter may be useful when removing distortion. May need to be used to correct non-linear distortion in multiple signals, where each signal or multiple signals have been distorted in a non-linear manner, but where the determination of one, some, or all of the signals relative to each other -62- This paper size is in accordance with Chinese National Standard (CNS) A4 (210 X 297 mm) Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs 496080 A7 B7___ 5. Description of the invention (60) Distortion is difficult. For such a correlated signal, the distortion can be corrected by making the small steps of the plurality of signals relatively change by the same amount. For signals that are a function of a single variable, the distortion function calculated in step 19 10 may depend on the slope of one or more of the signals at some point in the image . In this case, the anti-distortion function generated in step 1912 will attempt to equalize the slopes of the plurality of signals in small increments in the signal. For signals whose functions are multiple variables, the distortion function may depend on the rate of change of one, some, or all of the signals in multiple directions at multiple points in the signal. For example, the distortion function calculated in step 19 10 may depend on the gradient at one, some, or all points of the plurality of signals. Among them, a gradient measurement is used To determine the distortion of a signal relative to other signals, a scattergram can be used to generate this distortion function. Each point on the scatter diagram can be obtained by calculating a scatter diagram indicator in response to the magnitude of the gradient of some of the signals on the common point in the signal. A weight can be calculated in response to the angle between these gradients (for example: the cosine of the angle between the gradients). In addition, such weighting can be proportional to the magnitude of the gradient. A curve fitting algorithm can be used to fit a curve to the points on the scatter plot and use this curve to calculate an inverse distortion function. For example: This inverse distortion function may be proportional to the area under the curve. FIG. 20 illustrates a flowchart illustrating a method that can be used to improve a digital image according to the present invention. This aspect of the invention can be used to improve the color of a digital image. Digital image, regardless of whether it is a scanned photograph. The paper size applies the Chinese National Standard (CNS) A4 (210 X 297 mm) -63- —1 ---- • Installation -------- Order ------, — (Please read the notes on the back before filling this page) 496080 A7 B7 V. Description of the invention (61) (Please read the notes on the back before filling this page) Images, satellite images, Medical or other types of digital images may be distorted in a non-linear manner due to distortions in analog images and / or due to distortions that occur when digital images are generated. This distortion may vary with respect to intensity and / or frequency, and may vary differently for each channel. Such non-linear distortion may degrade the image in ways we don't want. For photographic images, the images may look less desirable. For satellite images and medical images, the information provided by the images may be less useful due to such distortions. Printed by a member of the Intellectual Property Bureau of the Ministry of Economic Affairs and a Consumer Cooperative. As mentioned above, many digital images include multiple channels. For example, digital images based on a color photo generally include a red channel, a green channel, and a blue channel. Among them, an image is composed of multiple channels, and the image information contained in each channel generates a common component among the three image channels. Therefore, over a small area in this digital image, the intensity tends to change across the small area in this image by an equal amount. Although the magnitude of the intensity may be quite different, the rate of change of such magnitude tends to be similar for an undistorted image. Therefore, the present invention attempts to determine the distortion of an image according to the intensity change in multiple directions for each channel (or a subset of channels) of an image. An inverse distortion function can be generated for one, some, or all channels in an image, so that when such an inverse distortion function is applied, the intensity of each channel tends to change by the same amount, especially when compared with other When the channels are compared over a small area of the image. In some image areas, when the image is not distorted, compared with other channels, the intensity of one channel tends to change linearly. However, in other images or in the same image area, when the image is not distorted , Compared with other channels -64- This paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 public love) Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs 496080 A7 ____B7 _ V. Description of the invention (62) The intensity of a channel tends to maintain a constant ratio between magnitudes over a small area. In order to automate the improvement of color digital images, the present invention uses the average of these two variations, where the change corresponds to the way a color image tends to be stored digitally. In step 1916, an image i is scanned to produce a digital image. Such scanning can occur using the scanner 30 or the scanner 34 described above. Alternatively, steps 1 9 1 6 can be omitted and the invention is applied to an image that has been previously scanned. In addition, the present invention can be applied to digital images having a plurality of channels, wherein each channel has a common signal component that distorts one or more channels in a non-linear manner. The term “common components” broadly refers to signals that are in some way interrelated. Alternatively, before performing the remaining steps in the method described in FIG. 20, the particle normalization aspect of the present invention discussed above may be applied to a digital image. Alternatively, before performing the remaining steps in the method described in FIG. 20, the one or more particle reduction techniques described above may also be applied to the image. Alternatively, the method described in FIG. 20 may be applied to the image without any particle normalization and / or particle reduction, or different techniques may be used after normalization or reduction. In step 119, each channel of the image (or a subset of the channels of the image) can be filtered. Therefore, various effects, such as noise and image resolution effects, will not be distorted with nonlinearity. The removal of the interference. In a small number of pixels, noise caused by particles such as ... may interfere with the removal of non-linear distortion. In addition, because the light used to produce a color photographic image passes through most layers of the film, 'the paper size of a digital image applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) -65 · " (Please read the notes on the back before filling this page)
n n n n 一:°4· ·ϋ n an ft n n I 496080 A7 B7 經濟部智慧財產局員工消費合作社印製 五、發明說明(63 ) 或多個波道可能有失真。如上述,彩色軟片的紅色層趨於 在光線穿過軟片的綠色和藍色層後,才接收光線。因此’ 在軟片的紅色層的影像比起在綠色和藍色層,趨於具稍微 多的失真。例如:影像可能在紅色波道上較不淸晰。這樣 的光學失真可能與在影像的色彩中的非線性失真的移除相 千涉。 因此,在步驟1 9 1 8中,可藉由過濾影像,而將影 像的每一波道預先模糊(p r e - b 1 u r )。在這實施例中,半 徑爲1 0像素的一中位數濾波器可被使用在每一色彩波道 在2000dp i下被掃描的軟片。可使用不同的半徑’ 而不悖離本發明的範圍。此外,可使用不同的掃描解析度 ,而不悖離本發明的範圍。除了應用一中位數濾波器之外 ,其他方法可被用來將影像預先模糊。基本上’影像的每 一波道可被改變,因而可將雜訊和如上所述的解析度效應 減幅或消除。 在步驟1 9 2 0中’對數位式影像的至少兩個波道’ 吾人可根據在影像的多數的點上,多個方向上的強度變化 率(斜率)而產生向量。在這實施例中,雖然這樣的向量 是對多數個波道的每一個產生的,但是,也可對多數個波 道的一子集合產生,而不悖離本發明的範圍。在這實施例 中,可以在影像的所有點上或影像的點之子集合上計算斜 率而產生向量,而不悖離本發明的範圍。 在這實施例中,所產生的向量代表一特別的影像波道 ,在所討論的點上的強度之梯度。換言之’此向量包括: (請先閱讀背面之注意事項再填寫本頁) 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -66 - 496080 A7 B7 經濟部智慧財產局員工消費合作社印製 五、發明說明(64 ) 在那影像波道中一特殊點上,代表在X 一方向的強度斜率 的第一分量和代表在y -方向的強度斜率的第二分量。與 在影像的點,於其他方向上的強度斜率成比例的向量也可 被選擇,而不悖離本發明的範圍。選擇在X和y方向的斜 率只是爲了方便起見。 吾人可使用公式7 6 - 7 8,計算在一特殊點(X, y )上紅色、綠色和藍色波道的向量,其分別爲:R ” 2 (X,y ) 、G ” 2 ( X,y )和 B ” 2 ( x,y ),其 中i和j分別代表在X和y方向上的單位向量。 R,2〇c,y) = 〇r(x + l,y)-IT〇c-l,y))i + (B/(x,y + l)-F/(#x,y-l))j (76) G,2〇c,y) = (G,(x +1, y) - G\x -1,y)) i + (G,(x, y +1) - G,〇c,y -1)) j (77) B^2(x, y) =: (B\x +1, y) - Bn(x -1, y)) i + (B\xy y +1) - B^(x, y -1)) j (78) 在這些公式中,R ” ·、G ”和B ”代表此三影像波道 被過濾後的強度値。因此’可以對一數位式影像的每一波 道中的每一點計算一向量。對於具有較多或較少的波道之 影像而言,類似的公式可用在每一波道上。吾人也可只對 每一波道的點之子集合計算向量,而不悖離本發明的範圍 。公式7 6 - 7 8並不I兌明在影像波道的邊緣上的斜率計 算。在這樣的邊緣上’可使用不同的斜率量測或者這些點 可以被忽略,而不彳孛離本發明的範圍。同樣地’可以使用 不同的計算斜率方法’而不悖離本發明的範圍。這實施例 利用討論的像素的鄰左和鄰右點’而不是此像素本身’來 計算在該點上$ 一方向的槪約斜率°可使用類似的方法在 y-方向上。替換地,可以使用在χ一或y一方向上更多 n n H n ϋ n«l— n ϋ · n n n n n n —1 一 01 I n· ·ϋ ·ϋ n n n n I (請先閱讀背面之注意事項再填冩本頁) 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -67- 經濟部智慧財產局員工消費合作社印製 496080 A7 B7____ 五、發明說明(65 ) 的像素,來計算X -方向或y —方向上斜率,而不悖離本 發明的範圍。此外,也可對一特別像素使用槪約斜率的某 平均値,而不悖離本發明的範圍。 在步驟1 9 2 2中,吾人可從多數個波道比較向量, 以產生散佈圖資料,其中,該資料對產生代表影像的一波 道與影像的其他波道比較的失真之失真函·數有用。再次地 ,這樣的散佈圖資料可以對一個、某些或所有的波道而產 生,而不悖離本發明的範圍。此外,完成向量的比較以產 生散佈圖資料可牽涉某些或所有的影像的其他波道。 在這實施例中,散佈圖資料取決於在步驟1 9 2 0所 計算出的向量。因爲未失真的影像趨於具有具類似方向的 梯度向量,所以,可對散佈圖上的每一點計算一加權,其 中,該加櫂與在步驟1 9 2 0中所算出的向量之間的點乘 積成正比。在這實施例中,對一特別的波道,對其他波道 的每一個計算這樣的點乘積,且每一這樣的點乘積可晝一 點。例如:對紅色散佈圖而言,吾人畫上具一加權的一點 ,其中,該加權與在步驟1 9 2 0中對紅色波道和綠色波 道所算出的向量的點乘積成正比。此外,對紅色波道在散 佈圖上將被放置的每一點計算一加權,其中,該加權與在 步驟1 9 2 0中對紅色波道和藍色波道所算出在該點的向 量的點乘積成正比。因此,在散佈圖上的一特殊點的加權 可以與在特殊點上兩波道的強度斜率之間的相似性成正比 。因爲波道相對彼此而失真,所以,本發明的這實施例決 定一特別波道的散佈圖點,其中,該點將所討論的波道的 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公爱) 68 (請先閱讀背面之注意事項再填寫本頁) -------I ^ ---------- %· 496080 A7 B7 五、發明說明(66 ) 特殊點上之斜率與在每一額外的波道對應點上強度的斜率 做比較。替換地,可對波道的子集合做比較,而不悖離本 發明的範圍。 本發明利用亦將在步驟1 9 2 0中對一特別波道產生 ,在一特殊點上斜率向量的相對大小,與在步驟1 9 2 0 中對影像的另一波道在一對應點上所計算出的斜率向量的 大小做比較,而決定散佈圖的一指標。 散佈圖的另一指標包括在影像波道中一點的強度。 散佈圖的點可如下而被產生。吾人對與將失真 因此 校正的影 對紅色波 子中,吾人 圖。散佈圖 軸指標爲在 像的每一波道產生一散佈圖。在這例 道、綠色波道和藍色波道產生一散佈 的X軸指標爲所討論的像素強度。散佈圖的y 1與- 1間的値,其中該値爲對應該特別波道的一點上斜 率向量的相對大小與影像的另一波道在一對應 向量的大小做比較而得。被分配至一特殊點的 述。因此,可使用公式7 9 — 9 3將在紅色、 散佈圖上的特殊點定位。在公式8 5 — 9 0中 的情形下,該加權値可被設爲一如0之預定値。 點上的斜率 加權如上所 綠色或藍色 的負數結果 (請先閱讀背面之注意事項再填寫本頁)n n n n One: ° 4 ·· ϋ n an ft n n I 496080 A7 B7 Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs. 5. Description of the invention (63) or multiple channels may be distorted. As mentioned above, the red layer of a color film tends to receive light after it passes through the green and blue layers of the film. So the image in the red layer of the film tends to be slightly more distorted than in the green and blue layers. For example: the image may be less clear on the red channel. Such optical distortions may be related to the removal of non-linear distortions in the colors of the image. Therefore, in step 1918, each channel of the image can be blurred in advance by filtering the image (p r e-b 1 u r). In this embodiment, a median filter with a radius of 10 pixels can be used for a film scanned at 2000 dpi per color channel. Different radii 'can be used without departing from the scope of the invention. In addition, different scanning resolutions can be used without departing from the scope of the invention. In addition to applying a median filter, other methods can be used to pre-blur the image. Basically every channel of the image can be changed, thereby reducing or eliminating the noise and resolution effects as described above. In step 1920, 'at least two channels of the logarithmic image', we can generate a vector based on the intensity change rates (slopes) in multiple directions at most points of the image. In this embodiment, although such a vector is generated for each of the plurality of channels, it may be generated for a subset of the plurality of channels without departing from the scope of the present invention. In this embodiment, the slope can be calculated on all points of the image or on a subset of the points of the image to generate a vector without departing from the scope of the present invention. In this embodiment, the generated vector represents a particular image channel, a gradient in intensity at the point in question. In other words, this vector includes: (Please read the precautions on the back before filling out this page) This paper size applies to China National Standard (CNS) A4 (210 X 297 mm) -66-496080 A7 B7 Employees of Intellectual Property Bureau, Ministry of Economic Affairs Printed by the Consumer Cooperative. 5. Description of the Invention (64) At a particular point in the image channel, the first component representing the intensity slope in the X-direction and the second component representing the intensity slope in the y-direction. Vectors that are proportional to the intensity slopes at other points in the image can also be selected without departing from the scope of the invention. The slopes in the X and y directions are chosen for convenience only. We can use the formulas 7 6-7 8 to calculate the vectors of the red, green and blue channels at a particular point (X, y), which are: R ”2 (X, y), G” 2 (X , Y) and B ″ 2 (x, y), where i and j represent unit vectors in the X and y directions, respectively. R, 2〇c, y) = 〇r (x + 1, y) -IT〇 cl, y)) i + (B / (x, y + l) -F / (# x, yl)) j (76) G, 2〇c, y) = (G, (x +1, y) -G \ x -1, y)) i + (G, (x, y +1)-G, 〇c, y -1)) j (77) B ^ 2 (x, y) =: (B \ x +1, y)-Bn (x -1, y)) i + (B \ xy y +1)-B ^ (x, y -1)) j (78) In these formulas, R ”·, G ”and B” represent the filtered intensity 値 of the three image channels. Therefore, 'a vector can be calculated for each point in each channel of a digital image. For images with more or fewer channels, a similar formula can be used for each channel. We can also calculate the vector only for the subset of points in each channel without departing from the scope of the present invention. Equations 7 6-7 8 do not calculate the slope calculated at the edges of the image channel. On such edges, different slope measurements can be used or these points can be ignored without departing from the scope of the invention. Likewise, different slope calculation methods can be used without departing from the scope of the present invention. This embodiment uses the adjacent left and right points of the pixel in question 'instead of the pixel itself' to calculate the approximate slope of the $ -direction at that point. A similar method can be used in the y-direction. Alternatively, you can use more in the χ one or y direction. Nn H n ϋ n «l— n ϋ · nnnnnn — 1 one 01 I n ·· ϋ · 注意 nnnn I (Please read the notes on the back before filling in 冩(This page) This paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) -67- Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs 496080 A7 B7____ 5. The pixel of invention description (65) X-direction or y-direction slope without departing from the scope of the present invention. In addition, a certain average of the approximate slope can be used for a particular pixel without departing from the scope of the present invention. In step 192, we can compare vectors from a plurality of channels to generate scatter plot data, where the data produces a distortion function that compares the distortion of one channel representing the image with the other channels of the image. it works. Again, such scatter plot data can be generated for one, some, or all channels without departing from the scope of the present invention. In addition, the comparison of vectors is done to generate other channels of scatter plot data that may involve some or all of the imagery. In this embodiment, the scatter plot data depends on the vector calculated in step 1920. Because undistorted images tend to have gradient vectors with similar directions, a weight can be calculated for each point on the scatter plot, where the point between this addition and the vector calculated in step 1920 The product is proportional. In this embodiment, such a point product is calculated for a particular channel, and for each of the other channels, and each such point product may be one point in the day. For example, for the red scatter diagram, we draw a point with a weight, where the weight is proportional to the point product of the vector calculated for the red channel and the green channel in step 1920. In addition, a weight is calculated for each point where the red channel is to be placed on the scatter diagram, where the weight is equal to the point of the vector at that point calculated for the red channel and the blue channel in step 1920. The product is proportional. Therefore, the weighting of a particular point on the scatter plot can be proportional to the similarity between the intensity slopes of the two channels at the particular point. Because the channels are distorted relative to each other, this embodiment of the present invention determines a scatter plot point for a particular channel, where this point applies the Chinese paper standard of the channel in question to the Chinese National Standard (CNS) A4 specification (210 X 297 public love) 68 (Please read the notes on the back before filling out this page) ------- I ^ ----------% · 496080 A7 B7 V. Description of the invention (66) The slope at a particular point is compared to the slope of the intensity at the corresponding point of each additional channel. Alternatively, a subset of the channels may be compared without departing from the scope of the invention. The present invention utilizes that the relative magnitude of the slope vector at a particular point will also be generated for a particular channel in step 1920, and at a corresponding point with another channel of the image in step 1920. The magnitudes of the calculated slope vectors are compared to determine an index of the scatter plot. Another indicator of the scatter plot includes the intensity of a point in the image channel. The points of the scatter plot can be generated as follows. Our image will be distorted and therefore the correction will be in the red wave, our image. Scatter chart The axis index produces a scatter chart for each channel of the image. In this example, a scattered X-axis index is generated for the channel, the green channel, and the blue channel for the pixel intensity in question. The 値 between y 1 and -1 of the scatter plot, where the 値 is obtained by comparing the relative magnitude of the slope vector at a point corresponding to a particular channel with the size of a corresponding vector at another channel of the image. A description assigned to a particular point. Therefore, you can use formulas 7 9-9 3 to locate special points on the red, scatter plot. In the case of Equations 8 5-9 0, the weighted chirp can be set to a predetermined chirp like zero. The slope of the points is weighted as above. Negative results of green or blue (please read the precautions on the back before filling this page)
· «II I I in ·ϋ n I 一-eJ· --« ϋ n n I - - I 經濟部智慧財產局員工消費合作社印製 Υι k 2R-2(x.y) G*2(x.y) |R"2) + |G"2| j· 2|R^2| 2|G^2| 、 (79) (80) (81) 本紙張尺度適用中國國家標準(CNS)A4規格(210 x 297公釐) .69 - 496080 Α7 Β7 五、發明說明(67) ( 1 2丨 G,2| 1 Y2Gi(x.y> 一 L \ |G#2|+ / 1 2|B,2| ) Y丨 B.2(x.y) 一 1 \ |R,2 卜 |B2|〉 f t 2|B"2| ) V — 1 2B^2(x.y) 1 |B,2 +|G,2| \ 111 \J IT2G.2 ^lR*2(x.y) 秦 〃2||G〃2| W_______ R,2B'2 (82) (83) (84) (85) (85) (請先閱讀背面之注意事項再填寫本頁)· «II II in · ϋ n I 一 -eJ ·-« ϋ nn I--I Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs 2 k-2R-2 (xy) G * 2 (xy) R " 2 ) + | G " 2 | j · 2 | R ^ 2 | 2 | G ^ 2 |, (79) (80) (81) This paper size applies to China National Standard (CNS) A4 (210 x 297 mm) .69-496080 Α7 Β7 V. Description of the invention (67) (1 2 丨 G, 2 | 1 Y2Gi (x.y >-L \ | G # 2 | + / 1 2 | B, 2 |) Y 丨 B. 2 (xy)-1 \ | R, 2 Bu | B2 |〉 ft 2 | B " 2 |) V — 1 2B ^ 2 (xy) 1 | B, 2 + | G, 2 | \ 111 \ J IT2G. 2 ^ lR * 2 (xy) 秦 〃2 || G〃2 | W_______ R, 2B'2 (82) (83) (84) (85) (85) (Please read the notes on the back before filling this page )
WlG-2(Xty) W- 2G*2(x*y) GH2VF2p2| B,2.G,2V|G#2||B^2| BH2V「B,2"R,2 丨 B"2 (87) (88) (89) 經濟部智慧財產局員工消費合作社印製 2B-2(x,y) >嘗2| ^Rix,y) ’ = R'x,y) (91) Χ(3(χ·χ) = G"(x,y) (92) ^B(x.y) = B'(x,y) (93) 公式9 1指示出:藉由決定在一特殊點上的紅:色波道 ^ g而決定該特殊點的散佈圖的X軸指標。對紅色波道的 該特殊點畫出兩點。每一點將具有在紅色散佈圖上相同的 χ 一指標。對晝出的這兩點而言,一點對應在該特殊點上 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -70 - 496080 A7 B7 五、發明說明(68 ) 紅色和綠色波道的斜率向量的比較,而另一點對應在該特 殊點上紅色和藍色向量間的比較。對代表紅色和綠色向量 的比較的點,公式7 9可被用來決定散佈圖的y —指標。 對紅色和藍色向量在特殊點上的比較,可使用公式8 0來 決定散佈圖上該點的y -指標。比較紅色波道對綠色波道 在散佈圖上的點之加權可使用公式8 5計算之,而代表紅 色與藍色波道間的比較的點之加權可使用公式8 6決定之 0 對在綠色波道上的任何點而言,也可產生綠色散佈圖 的兩點。在每一情形下,綠色散佈圖的X -軸指標是使用 公式9 2而決定的。對對應綠色波道和紅色波道的向量之 間的比較的第一點,公式8 1被用來決定在綠色散佈圖上 該點的y -軸指標,而公式8 7代表在該點上被畫出的加 權。對對應在該特殊點上綠色波道和藍色波道的向量之間 的比較之點,公式8 2被用來決定y -軸指標’而公式 8 8被用來決定該特殊點的加權。 對在藍色波道上的任何點而言,也可產生藍色散佈圖 的兩點。在每一情形下,藍色散佈圖的X -軸指標是使用· 公式9 3而決定的。對對應藍色波道和紅色波道的向量之 間的比較的第一點,公式8 3被用來決定在藍色散佈圖上 該點的y -軸指標,而公式8 9代表在該點上被晝出的加 權。對對應在該特殊點上藍色波道和綠色波道的向量之間 的比較之點,公式8 4被用來決定y -軸指標’而公式 9 0被用來決定該特殊點的加權。 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) (請先閱讀背面之注意事項再填寫本頁) -裝--------訂------J---, 經濟部智慧財產局員工消費合作社印製 -71 - 496080 A7 B7 一丨丨· * 丨丨 " — 五、發明說明(69 ) f請先閱讀背面之注音》事項再填寫本頁) 在步驟1 9 2 4中,吾人根據在步驟1 9 2 2中所產 生的散佈圖資料對多數個波道的每一個產生一散佈圖。這 些步驟可被合倂,而不悖離本發明的範圍。此外,若只有 多數個波道的中的一個或某些波道將被校正其失真,則只 有對那些將被校正的波道產生一散佈圖。再次地,如上述 ,在本發明的這實施例中所產生的散佈圖包括與波道的失 真成正比的點,其中,散佈圖是對其他波道的每一個而產 生的。然而,可使用相對於額外波道的一子集合的失真, 而不悖離本發明的範圍。 經濟部智慧財產局員工消費合作社印製 圖2 1例示對應在影像中一單一空間位置而可被畫出 的6個點。對紅色散佈圖而言,所討論的像素強度決定兩 點的每一點之X —軸指標1946。點1934的y —軸 指標由公式7 9所決定。被分配至點1 9 3 4的加權由公 式85所決定。此處,公式79產生一値爲〇.5的y-軸指標。因爲公式8 1其竇是公式7 9的負數,所以,對 應在綠色散佈圖上點1 9 3 6的y -指標爲一 〇 . 5的y 一指標。點1 9 3 6的X —指標是使用公式9 2而決定, 且代表在綠色波道中所討論的像素的強度。由公式8 7產 生點1 9 3 6的加權將爲由公式8 5產生點1 9 3 4的加 權之負數。因此,本發明的各種不同實行可僅僅使用一套 公式,如公式8 5、8 6和8 8來計算加權,而不悖離本 發明的範圍。 再次地,點1 9 3 8具有由公式9 1所決定,在所討 論的像素上的X —指標。其y -指標是由公式8 〇所決定 72 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) 496080WlG-2 (Xty) W- 2G * 2 (x * y) GH2VF2p2 | B, 2.G, 2V | G # 2 || B ^ 2 | BH2V 「B , 2 " R, 2 丨 B " 2 (87 ) (88) (89) 2B-2 (x, y) > taste 2 | ^ Rix, y) '= R'x, y) (91) χ (3 ( χ · χ) = G " (x, y) (92) ^ B (xy) = B '(x, y) (93) Formula 9 1 indicates: by determining red at a particular point: color wave Channel ^ g to determine the X-axis index of the special point scatter plot. Draw two points for the special point of the red channel. Each point will have the same χ-one index on the red scatter plot. For the two In terms of points, one point corresponds to this special point. The paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) -70-496080 A7 B7. 5. Description of the invention (68) The slope of the red and green channels Comparison of vectors, and another point corresponds to the comparison between the red and blue vectors at this particular point. For points representing the comparison of red and green vectors, Equation 7 9 can be used to determine the y-index of the scatter plot. For red Compared with the blue vector at special points, you can use the formula 8 0 to determine the scatter plot. The y-index of this point. Comparing the weight of the red channel to the points of the green channel on the scatter plot can be calculated using Equation 8 5 and the weight of the points representing the comparison between the red and blue channels can be calculated using Equation 8. 6 of decision 0 For any point on the green channel, two points of the green scatter diagram can also be generated. In each case, the X-axis index of the green scatter diagram is determined using Equation 9 2. For the corresponding The first point of the comparison between the vectors of the green channel and the red channel, Equation 8 1 is used to determine the y-axis index of the point on the green scatter plot, while Equation 8 7 represents the point at which the point is drawn. For the point of comparison between the vectors corresponding to the green and blue channels at this particular point, Equation 8 2 is used to determine the y-axis index 'and Equation 8 8 is used to determine the special point. For any point on the blue channel, two points of the blue scatter diagram can also be generated. In each case, the X-axis index of the blue scatter diagram is determined using Equation 9 3 The first point of the comparison between the vectors corresponding to the blue and red channels is 8 3 is used to determine the y-axis index of the point on the blue scatter chart, and formula 8 9 represents the weighting by day out at this point. The corresponding blue and green channels at this particular point The point of comparison between the vectors, formula 8 4 is used to determine the y-axis index 'and formula 9 0 is used to determine the weighting of this particular point. This paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) (Please read the precautions on the back before filling out this page) -Installation -------- Order ------ J ---, Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs- 71-496080 A7 B7 I 丨 丨 · * 丨 丨 " — V. Description of the invention (69) f Please read the note on the back first and then fill out this page) In step 1 9 2 4 we according to step 1 9 The scatter plot data produced in 2 produces a scatter plot for each of the plurality of channels. These steps can be combined without departing from the scope of the invention. In addition, if only one or some of the plurality of channels will be corrected for distortion, only a scattergram will be generated for those channels to be corrected. Again, as described above, the scatter plot generated in this embodiment of the present invention includes points that are directly proportional to the distortion of the channel, where the scatter plot is generated for each of the other channels. However, distortions relative to a subset of the extra channels can be used without departing from the scope of the invention. Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs Figure 21 illustrates six points that can be drawn corresponding to a single spatial location in the image. For the red scatter plot, the pixel intensity in question determines the X-axis index of each of the two points, 1946. The y-axis index of the point 1934 is determined by Equation 79. The weight assigned to point 1 9 3 4 is determined by Formula 85. Here, Formula 79 yields a y-axis index of 0.5. Since the sinus of formula 81 is the negative of formula 79, the y-index corresponding to the point 1 9 3 6 on the green scatter chart is a y-index of 0.5. The X-point at 1 9 3 6 is determined using Equation 9 2 and represents the intensity of the pixel in question in the green channel. The weighting of points 1 9 3 6 generated by Equation 8 7 will be the negative weighting of points 1 9 3 4 generated by Equation 8 5. Therefore, various implementations of the present invention may use only one set of formulas, such as formulas 8, 5, 8 and 88, to calculate the weighting without departing from the scope of the invention. Again, point 1 9 3 8 has an X-index on the pixel in question, determined by Equation 91. The y-index is determined by the formula 80. 72 The paper size is applicable to the Chinese National Standard (CNS) A4 specification (210 X 297 mm) 496080
經濟部智慧財產局員工消費合作社印製 五、發明說明(π ) 且其加權是由公式8 6所決定。這個點提供當與藍色訊號 比較時,紅色訊號的失真的量測。在藍色散佈圖上’點 1 ·9 4 0的y —指標爲點1 9 3 8的y —指標之負數’正 如公式8 3僅是公式8 0的負數。點1 9 4 0的X -指標 1 9 5 0對應所討論的像素在藍色波道中的強度値,且是 使用公式9 3所決定的。 同樣地,在綠色散佈圖上的點1 9 4 2和在監色散佈 圖上的點1 9 4 4具有彼此互爲負數的y -指標(加法上 的相反),正如公式84是公式82的負數。 在步驟1 9 2 4中產生散佈圖之後,可使用一曲線擬 合演算法,經由散佈圖,產生一曲線,如步驟1 9 2 6中 所例示者。在這實施例中,吾人使用曲線擬合演算法的一 平均型。可使用任何曲線擬合演算法,而不悖離本發明的 範圍。可對每一散佈圖產生一曲線。若失真校正被應用至 散佈圖所產生的波道上,則假設,將只在步驟1 9 2 4中 產生一散佈圖。若不對特別波道產生任何失真校正,則在 步驟1 9 2 4中散佈圖的產生和在步驟1 9 2 6中曲線的 產生均可對該波道而被省略。 在步驟1 9 2 8中,可根據在步驟1 9 2 6中所產生 的曲線,產生一反失真函數。這反失真函數可爲與步驟 1 9 2 6中所產生的曲線成正比的任何函數。步驟 1 9 2 6中所產生的曲線與一特別波道對具有一共同成分 的一或多個額外訊號之相對失真成正比。在這實施例中, 每一波道的反失真函數可使用公式9 4 一 9 6而產生。 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -73 - (請先閱讀背面之注意事項再填寫本頁) 裝--------訂------*---· 496080 A7 ________B7___ 五、發明說明(71 ) ΓΚ[€] = Σήί^ψ (94) «·-〇 y a — -Kj (請先閱讀背面之注意事項再填寫本頁) fg,!餐·· 剛(⑹ 公式9 4可被用在紅色波道’公式9 5被用在綠色波 道,而公式9 6被用在藍色波道i: °因爲所涉及數位 式影像,所以’可以使用強度値的一分離數目。例如:這 實施例使用一 8 —位元強度値’其可產生2 5 6個不同可 能的強度値。因此’此反失真函數包括2 5 6個値的一查 表,其範圍爲0 — 2 5 5。因此,在公式94 一 9 6中, “C”具一最大値爲255 ,其中,數位式影像具256 個可能強度値。其他數目的位元可被使用’而不悖離本發 明的範圍。在公式9 4 一 9 6中的” y ”値代表組成此數 位式影像的波道數目。公式9 4中的R 1代表在紅色波道 上,在s上,X -指標爲” 1 ” ’於步驟1 9 2 6中所產 經濟部智慧財產局員工消費合作社印製 生曲線的y値。換言之’在步驟1 9 2 6中,經由散佈圖 所產生的曲線爲代表回應像素強度所畫的失真。公式9 4 對每一所討論的像素強度計算一反失真函數,其中,該函 數與在所討論的強度之前’在步驟1 9 2 6中所產生的曲 線下面積成正比。公式9 5和9 6與公式9 4相類似。 在公式9 4 一 9 6中使用y/y — 1方根,因爲不用 此則可能會導致過度校正。每一波道的校正是基於相對其 他波道的它的失真。若進行一校正而不考慮此因素,假設 ^紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -74 - ' 496080 A7 B7 經濟部智慧財產局員工消費合作社印製 五、發明說明(72 ) 只要校正一波道,則可對每一波道做相對校正。在這例子 中,因爲每一波道都要被校正,所以,若不取適當的方根 ,則會導致淨過度校正。 在步驟1 9 3 0中,可使用在步驟1 9 2 8中所決定 的反失真函數,將每一影像波道歸一化。爲了歸一化,將 處理影像中每一像素強度値。對一特別的空間位置而言, 藉由使用在那像素的紅色波道的強度作爲一指標,而將紅 色強度改變爲使用公式9 4所產生的反失真函數。此紅色 強度値被由反失真函數所供應的強度所取代。同樣地,在 一特殊像素的綠色強度値被在步驟1 9 2 8中,對此綠色 波道所產生的反失真函數所產生的綠色強度値所取代,其 中,此綠色波道的反失真函數之指標包括此綠色波道的目 前像素強度値。最後,使用在步驟1 9 2 8中所產生的藍 色反失真函數,類似地校正像素的藍色波道強度値。 雖然此處所描述的本發明渉及在空間域和頻率域的計 算,但是,也可能相當地使用相反域中的計算,而不悖離 本發明的範圍。 雖然吾人已詳細地描述本發明,可理解的是:可至此 作各種不同的改變、代替和替換,而不悖離如所附加的申 請專利範圍定義的本發明的精神和範圍。 爲了幫助專利局和在這應用上發表專利的任一讀者解 釋至此所附加的申請專利範圍,專利申請人希望注意的是 :他們並不意欲所附的申請專利範圍之任何項懇求美國專 利法(U · S . C · ) § 1 1 2中3 5之6節當作它在歸檔 (請先閱讀背面之注意事項再填冩本頁) 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐) -75- 496080 A7 _B7_ 五、發明說明(73 ) 日期存在,除非“…的裝置” (’’means for")或“…的步驟 ” "s t e p f 〇 r ")的字眼被使用在特殊的申請專利範圍中。 (請先閱讀背面之注意事項再填寫本頁) 經濟部智慧財產局員工消費合作社印製Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs. 5. Description of invention (π) and its weight is determined by formula 86. This point provides a measure of the distortion of the red signal when compared to the blue signal. On the blue scatter chart, y of point 1 · 9 4 0-the index is y of point 1 9 3 8-the negative number of the index is positive, as Equation 8 3 is only a negative number of Equation 80. The X-index 1 950 of the point 1 9 0 0 corresponds to the intensity 値 of the pixel in question in the blue channel, and is determined using Equation 9 3. Similarly, points 1 9 4 2 on the green scatter diagram and points 1 9 4 4 on the monitor scatter diagram have y-indexes (negative additions) that are negative to each other, as Equation 84 is Equation 82 negative number. After generating the scatter diagram in step 19 2 4, a curve fitting algorithm can be used to generate a curve through the scatter diagram, as illustrated in step 1 9 2 6. In this embodiment, we use an average type of curve fitting algorithm. Any curve-fitting algorithm can be used without departing from the scope of the invention. A curve can be generated for each scatter plot. If distortion correction is applied to the channels produced by the scattergram, it is assumed that only a scattergram will be generated in step 192 2. If no distortion correction is generated for a particular channel, both the generation of the scatter diagram in step 19 24 and the generation of the curve in step 19 26 can be omitted for this channel. In step 1928, an inverse distortion function can be generated based on the curve generated in step 1926. This inverse distortion function can be any function that is proportional to the curve produced in step 1 9 2 6. The curve produced in step 1 9 2 6 is proportional to the relative distortion of a particular channel to one or more additional signals having a common component. In this embodiment, the inverse distortion function of each channel can be generated using equations 9 4-9 6. This paper size applies to China National Standard (CNS) A4 specifications (210 X 297 mm) -73-(Please read the precautions on the back before filling this page) Loading -------- Order ----- -* --- · 496080 A7 ________B7___ 5. Description of the invention (71) ΓΚ [€] = Σήί ^ ψ (94) «· -〇ya — -Kj (Please read the notes on the back before filling this page) fg, Meal ... Gang (⑹ Formula 9 4 can be used in the red channel 'Formula 9 5 is used in the green channel, while Formula 9 6 is used in the blue channel i: ° Because the digital image is involved, so 'A separate number of intensity 可以 can be used. For example: this embodiment uses an 8-bit intensity 値' which can produce 2 5 6 different possible strengths 因此. Therefore 'this anti-distortion function includes 2 5 6 値A look-up table with a range of 0-2 5 5. Therefore, in formulas 94 to 96, "C" has a maximum value of 255, of which digital images have 256 possible intensities. Other numbers of bits Can be used 'without departing from the scope of the present invention. "Y" in Equations 9 4-9 6 represents the number of channels that make up this digital image. Equation 9 4 R 1 represents on the red channel, on s, the X-indicator is “1” 'y' printed in step 1 9 2 6 by the consumer property cooperative of the Intellectual Property Bureau of the Ministry of Economics to print the production curve. In other words, in the step In 1 9 2 6, the curve generated by the scatter plot is the distortion drawn in response to the intensity of the pixel. Equation 9 4 calculates an inverse distortion function for each pixel intensity in question, where the function is the same as the intensity in question. Previously, the area under the curve generated in step 1 9 2 6 is directly proportional. Equations 9 5 and 9 6 are similar to equation 9 4. Use y / y — 1 square root in equation 9 4-9 6 because it is not used This may lead to over-correction. The correction of each channel is based on its distortion relative to other channels. If a correction is made without considering this factor, it is assumed that the paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) -74-'496080 A7 B7 Printed by the Consumer Cooperative of the Intellectual Property Bureau of the Ministry of Economic Affairs V. Invention Description (72) As long as one channel is corrected, relative correction can be made for each channel. In this example Because every channel has to be corrected Therefore, if the proper square root is not taken, it will cause a net overcorrection. In step 1930, the inverse distortion function determined in step 1928 can be used to normalize each image channel. For normalization, the intensity of each pixel in the image is processed. For a particular spatial location, the intensity of the red channel at that pixel is used as an indicator, and the intensity of red is changed to the formula 9 The resulting inverse distortion function. This red intensity 値 is replaced by the intensity supplied by the inverse distortion function. Similarly, the green intensity 値 of a particular pixel is replaced by the green intensity 产生 generated by the inverse distortion function generated for this green channel in step 192, where the inverse distortion function of this green channel is The indicators include the current pixel intensity of this green channel. Finally, using the blue inverse distortion function generated in step 192, the pixel's blue channel intensity 値 is similarly corrected. Although the invention described herein deals with calculations in the spatial and frequency domains, it is also possible to use calculations in the opposite domain fairly without departing from the scope of the invention. Although I have described the present invention in detail, it is understood that various changes, substitutions and substitutions can be made hereto without departing from the spirit and scope of the present invention as defined by the appended claims. In order to help the patent office and any readers who have published patents on this application explain the scope of the patent application attached to this point, patent applicants wish to note that they do not intend that any of the attached patent application scopes solicit US patent law U · S. C ·) § 1 1 2 3 3 6 6 is archived (Please read the precautions on the back before filling this page) This paper size applies the Chinese National Standard (CNS) A4 specification (210 X 297 mm) -75- 496080 A7 _B7_ V. Description of the invention (73) Date exists, unless the word "means for" ("means for ") or" step of "" stepf 〇r ") Used in special patent applications. (Please read the notes on the back before filling out this page) Printed by the Consumer Cooperatives of the Intellectual Property Bureau of the Ministry of Economic Affairs
-76- 本紙張尺度適用中國國家標準(CNS)A4規格(210 X 297公釐)-76- This paper size applies to China National Standard (CNS) A4 (210 X 297 mm)
Claims (1)
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15036899P | 1999-08-20 | 1999-08-20 |
Publications (1)
Publication Number | Publication Date |
---|---|
TW496080B true TW496080B (en) | 2002-07-21 |
Family
ID=22534209
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW089107919A TW496080B (en) | 1999-08-20 | 2000-04-26 | Method and system for normalizing a plurality of signals having a shared component |
Country Status (5)
Country | Link |
---|---|
EP (1) | EP1208533A1 (en) |
CN (1) | CN1379891A (en) |
AU (1) | AU4822300A (en) |
TW (1) | TW496080B (en) |
WO (1) | WO2001015085A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7852800B2 (en) | 2002-07-23 | 2010-12-14 | Qualcomm Incorporated | Reducing interference between users in a communications system through time scattering |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8416986B2 (en) | 2009-10-29 | 2013-04-09 | Raytheon Company | Methods and systems for processing data using non-linear slope compensation |
US8738678B2 (en) | 2009-10-29 | 2014-05-27 | Raytheon Company | Methods and systems for determining an enhanced rank order value of a data set |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4677465A (en) * | 1985-11-01 | 1987-06-30 | Eastman Kodak Company | Digital color image processing method with shape correction of histograms used to produce color reproduction functions |
US6442301B1 (en) * | 1997-01-06 | 2002-08-27 | Applied Science Fiction, Inc. | Apparatus and method for defect channel nulling |
-
2000
- 2000-04-26 TW TW089107919A patent/TW496080B/en not_active IP Right Cessation
- 2000-05-03 CN CN00814540.7A patent/CN1379891A/en active Pending
- 2000-05-03 EP EP00930395A patent/EP1208533A1/en not_active Withdrawn
- 2000-05-03 WO PCT/US2000/012312 patent/WO2001015085A1/en active Application Filing
- 2000-05-03 AU AU48223/00A patent/AU4822300A/en not_active Abandoned
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7852800B2 (en) | 2002-07-23 | 2010-12-14 | Qualcomm Incorporated | Reducing interference between users in a communications system through time scattering |
Also Published As
Publication number | Publication date |
---|---|
EP1208533A1 (en) | 2002-05-29 |
WO2001015085A1 (en) | 2001-03-01 |
AU4822300A (en) | 2001-03-19 |
CN1379891A (en) | 2002-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1222631B1 (en) | Method and system for enhancing digital images | |
US7356194B2 (en) | Image processing method, apparatus, and program | |
US6792162B1 (en) | Method and apparatus to automatically enhance the quality of digital images by measuring grain trace magnitudes | |
US7146059B1 (en) | Method of performing fast bilateral filtering and using the same for the display of high-dynamic-range images | |
Velde | Multi-scale color image enhancement | |
EP0971314B1 (en) | A method for preserving image detail when adjusting the contrast of a digital image | |
US7260266B2 (en) | Image processing method, image processing apparatus, image processing program and image recording apparatus | |
CN106157273B (en) | Method and device for generating composite picture | |
US8249341B2 (en) | Image processing method for removing noise contained in an image | |
JP2000050081A (en) | Automatic tone control method by contrast gain control of edge | |
JPH10187966A (en) | Method and device for filtering image | |
US6687414B1 (en) | Method and system for normalizing a plurality of signals having a shared component | |
Trentacoste et al. | Unsharp masking, countershading and halos: enhancements or artifacts? | |
JP2000187728A (en) | Method for adjusting picture detail for correcting applied tone scale | |
US20180307024A1 (en) | Method for generating a reflection-reduced contrast image and corresponding device | |
Li et al. | Content adaptive guided image filtering | |
TW496080B (en) | Method and system for normalizing a plurality of signals having a shared component | |
Krawczyk et al. | Contrast restoration by adaptive countershading | |
US6549678B1 (en) | Method for preserving spatial detail when applying a multidimensional tonal transform to a digital color image | |
KR101462440B1 (en) | Apparatus and Method for Correcting Color | |
US6453075B1 (en) | Method for preserving image detail when adjusting the tone scale of a digital color image | |
TW478279B (en) | Method and apparatus to automatically enhance the quality of digital images | |
Kok et al. | Digital Image Denoising in MATLAB | |
Lee et al. | Multi-image high dynamic range algorithm using a hybrid camera | |
Aon et al. | Scene adaptive color compensation and multi-weight fusion of underwater image |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
GD4A | Issue of patent certificate for granted invention patent | ||
MM4A | Annulment or lapse of patent due to non-payment of fees |