Способ оптической обработки изображений и устройство для его осуществления
Использование: в оптических процессорах при обработке изображений методом пространственно-частотной фильтрации. Сущность изобретения: способ включает освещение полихроматической волной транспаранта с обрабатываемым изображением, первое оптическое преобразование Фурье над волной, промодулированной транспарантом , фильтрацию диапазонов пространственных частот путем их цветной селекции в зависимости от фильтруемых диапазонов частот по представленному закону. Устройство состоит из последовательно расположенных источника полихроматического излучения, осветительную систему, транспарант , фурье-объектов, пространственночастотный фильтр, выполненный в виде одного или нескольких кольцевых фильтров из различных цветных стекол с функцией пропускания, определяемый по формуле, фурье-преобразующий объектив и плоскость регистрации изображения транспаранта . 9 ил. ел
СОКОВ СОВЕТСКИХ
СОЦИАЛИСТИЧЕСКИХ
РЕСПУБЛИК (я)з G 02 В 27/46
ГОСУДАРСТВЕННОЕ ПАТЕНТНОЕ
ВЕДОМСТВО СССР (госпАтент ссса) ОПИСАНИЕ ИЗОБРЕТЕНИф
К ПАТЕНТУ (21) 4861480/10 (22) 31.09.90 (46) 23.04.93. Бюл. hh 15 (71) Научно-исследовательский институт радиоэлектроники и лазерной техники МГТУ им. Баумана (72) Д.B.Çóáêoâ, А.M.Kèëèí, О,В.Рожков и
Л.Н.Тимашова (73) Д. В;Зубков (56) Зль-Сам П, Считывание информации, записанной в фазомодулирующей среде (В кн.: Оптическая обработка информации)
Под ред. С.П.Ерковича. — М.; Мир, 1966.
Патент CLUA N 3729252, кл. G 02 В 27/38, 1971. (54) СПОСОБ ОПТИЧЕСКОЙ ОБРАБОТКИ
ИЗОБРАЖЕНИЙ И УСТРОЙСТВО ДЛЯ ЕГО
ОСУЩЕСТВЛЕНИЯ (57) Использование; в оптических процессорах при обработке изображений методом
Изобретение относится к области оптического приборостроения и может быть использовано в оптических процессорах при обработке изображений методом пространственноо-частотной фильтрации.
Целью изобретения является повышение информационной емкости обрабатываемого изображения эа счет повышения освещенности крупных объектов и сокращение времени обнаружения объектов эа счет получения большого цветового различия объектов разных форм и размеров.
„„5Q „„1811628 А3 пространственно-частотной фильтрации.
Сущность изобретения. способ включает освещение полихроматической волной транспаранта с обрабатываемым изображением, первое оптическое преобразование Фурье над волной, промодулированной транспарантом, фильтрацию диапазонов пространственных частот путем их цветной селекции в зависимости от фильтруемых диапазонов частот по представленному закону. Устройство состоит иэ последовательно расположенных источника полихроматического
° излучения, осветительную систему, транспарант, фурье-объектов, пространственночастотный фильтр, выполненный в виде одного или нескольких кольцевых фильтров из различных цветных стекол с функцией пропускания, определяемый по формуле, фурье-преобразующий объект <в и плоскость регистрации иэображения транспаранта. 9 ил. ееавЪ
Поставленная цель достигается тем, что С в способе, включающем формирование по- Я лихроматической освещающей волны, осве- (ф щение этой волной транспаранта с обрабатываемым изображением, первое оптическое преобразование Фурье над вол- Ы ной, промодулированной транспарантом, фильтрацию диапазонов пространственных частот, второе оптическое преобразование
Фурье, в результате которого формируется обработанное изображение, регистрацию обработанного изображения, осуществляют фильтрацию диапазонов пространствен1811628
PN — макс макс 40
40 ных частот путем их цветовой селекции по закону
"«9) = iip9) g Jclrc ф) — circ i )i Л
k =1 У".к уа — 1 У 5 где Uq(p) - амплитуда поля на длине волны
А до фильтрации, I
Uq p)- амплитуда поля на длине волны
Х после фильтрации, p — радиальная полярная координата в 10 плоскости пространственно-частотного спектра;
p< — внешний радиус k-й кольцевой зоны пространственно- частотного спектра, г „- коэффициент пропускания по интенсивйости k-й кольцевой зоны частотного спектра для длины волы А, N — количество фильтруемых диапазонов частот, причем ро = О.
В устройстве, реализующем предложенный способ и содержащем расположенные последовательно на одной оптической оси источник полихромэтического излучения, осветительную систему, транспарант, фурье-объектив, пространственно-частотный фильтр, фурье-преобразующий объек тив, плоскость регистрации изображения .транспаранта, пространственно-частотный фильтр выполнен в виде одного или более цветных кольцевых фильтров из различных цветных стекол и имеет функцию пропускания а максимальный внешний радиус M-й (последней) кольцевой зоны фильтра:
ГДЕ 1м,кс — МаКСИМаЛЬНаЯ ПРОСтРаНСтВЕНная частота обрабатываемоro иэображения;
Лмакс максимальная длина волны;
fe.o — фокусное расстояние фурье-объ- 45 ективэ.
На фиг. 1 изображена функциональная оптическая схема устройства для осуществления предложенного способа; на фиг. 2— цветной кольцевой фильтр с рабочей зоной 50 р«р <р2, выполненный из цветного стекла марки СЗС22; на фиг. 3 — зависимость времени обнаружения объектов от их цветового различия; на фиг. 4 — входной сигнал— функция пропускания объекта типа 55 г
cIrc(); на фиг. 5 — цветовой график с гоб результатами проверки способа методом численного моделирования; на фиг. 6-9— графики зависимостей яркости Яот координаты r поля обработанного изображения, 1 построенные по данным таблиц 1-4 — соответственно. В таблицах 1-4 приведены расчеты координат цветности Х, Y, Z и яркости
Ь для объектов с радиусами r
0 < r < 1,5 rpQ. для случаев; табл. 1—
1 объект r
Устройство (фиг. 1) представляет собой проекционную систему, которая содержит последовательно установленные вдоль оптической оси источник полихроматического излучения 1, конденсор 2, транспарант 3, фурье-объектив (ФО) 4, цветной кольцевой фильтр 5, фурье-преобразующий объектив (ФПО) 6, плоскость регистрации (экран) 7, Транспарант 3 установлен в передней фокальной плоскости ФО 4; фильтр 5 — в задней фокальной плоскости ФО 4, совмещенной с передней фокальной плоскостью ФПО 6; экран 7 — в задней фокальной плоскости ФПО 6.
Устройство работает следующим образом.
Обрабатываемый фотоснимок — транспарант 3 с помощью осветительной системы — источника 1 и конденсора 2 — освещают полихроматическим светом. Рассмотрим осесимметричный объект фотоснимка, функция пропускания которого имеет вид (см. фиг. 4): тоб (ь,rr) = тоб (r) = CirC () (1) де (Q — координаты входной плоскости (плоскости транспаранта);
roe — радиус объекта.
Распределение поля в задней фокальной плоскости ФО 4 на длине волны k
2а
Оф(Эр) = ф я Я - т — 1 а (1) 3,(è 1 3 ) 1 (= г,а, S,(eа,1гй.г, Х,i 2 p "as) (2) где P — символ прямого преобразования
Фурье-Бесселя; м = — Š— -радиальная пространствен Q.o ная частота;
1811628
Фк
I Iq(I ))IqdS
Х(т ) ))
5 )12
1 ()" )Х ))4Ъ+ I (I Jt)«Xdr()+ ) I () «)Я )It
1 1
10 (3) U P) =гоб
«
Из выражения (3) видно, что масштаб спектра существенно зависит от размера объекта г0б: для малоразмерных объектов масштаб больше (спектр "растянут" сильнее), для крупноразмерных — меньше, Таким образом, кольцевой фильтр 5 будет оказывать различное влияние на спектры объектов разных размеров. В соответствии 20 с изобретением, функция пропускания цветного кольцевого фильтра, состоящего из N кольцевых зон:
«ф(««) =2, )circ (+) — с1гс(— — )1 rr«2«, 25
Pk . P< — 1 (4)
Тогда распределение поля непосредственно за фильтром будет иметь вид:
1 Т(, "оЕР1 () (р) (),„(гк) t () Г, . ф.о
9 0b Р У н
; «,0„((с««с — )-c)rc(),) (j 5)
После обратного преобразования Фурье с помощью ФПО 6 получим распределение поля в плоскости регистрации 7 (задней фо- 40 кальной плоскости ФПО 6):
III(,) ° J рк,(р)1. tI," " Р)гр.
=г;2:",J 1(, С "«р).) 4 ).(кт.—;. р)4р 45
" J)x-< (Ф )
С учетом спектральной мощности источника света P(il.) распределение интенсивности в плоскости регистрации 7:
)„((г ) = l 0 (г ) j Р(Л) В международной колориметрической системе XYZ координаты цвета изображе- 55 ния вычисляются по формулам: а) при сплошном спектре источника света в диапазоне длин волн (Л1; Л2):
Jo(X), J1(X) — функции Бесселя 1-го рода
0-ro и 1-го порядков соответственно (Х-аргумент). ф.0 — фокусное расстояние фурье-объ1 ектива 4.
Переходя к координатам в частотной плоскости спектра, получим:
«)«и
Iq(r ltd 4) )(J I>()"))c5Ide+J Г (и ))
r)I I «Х1
%
Iq(1" ) Е,И () )-.,„
I ())()")х4дэ+ J I (r ) 4)d ))4 (yq() ) de
I «))1
1 (8)а а относительная визуальная яркость; hp
j Z>(r ) r)t,d4r, Lv(61)1 р(.) >, ««1 где Х Уд, Ля- удельные цветовые координаты, зависящие от длины волны излучения k б) при дискретном спектре источника света с набором длин волн Л,..., Л():
Z 4„()х,„ Х(п 1 . xi«
I. и I
)2„(c)II, ()„("1 к,4Z t)Д(г 1Я„ к к к
Q (<„(t ) g4I
+ I4, l)(yew+ lac(clg4Is tX 1)к(г1 Кк
° К.1
11 « кк("1 кк
c«(at) s
;>„ 15„(п)ха„+ 1 к(г1 лккй ккк(")x x к «Xs ° и
Z I (")3% „
Ък("1- "*.
Я а)к 32 к (5á) .
Таким образом, цветность изображения осесимметричного объекта типа c)rc(— ) опГ (об ределяется его размером, спектральной мощностью излучения источника и спектральным коэффициентом пропускания цветного кольцевого фильтра.
Рассмотрим работу устройства, фильтр
5 которого выполнен (согласно изобретению) в виде кольца и изображен на фиг. 2.
Его функция пропускания получается из формулы (4) при N 2:
«««(о) circф)+ (circ (g) — circ ô)) r r«
f 9) 1811628 где p, pz — внутренний и внешний радиусы фильтра соответственно; тг- зависимость коэффициента пропускания фильтра по интенсивности от длины лопни излучения. Формула (5J примет еид. (У1 1 1 oS P
Ц(р) О (р) (р) = ".5
1 с(гс — t с (с — -с1 C— (7О)
Соответственно формула (6) преобразуется к виду: («((=3 (р(.(Ä," "р}ь,"я pIзя(л(з)Jp+
Ра (11)
Таким образом, как следует из полученных выражений, при большом гм спектр в плоскости фильтра 5 "сжат" настолько, что почти весь находится в области 0 (p
Работа предлагаемого устройства с источником полихроматического излучения
ДРШ-100W и с фильтром, изображенным на фиг. 2, моделировалась на ПЭВМ. Результаты расчета координат цветности для объектов г<>е7 = 0,3 мм и гуе2 = 0,03 мм при
Ъ толщине фильтра L1 - 0 мм (т.е. при отсутствии фильтра) и Lz - 3 мм, выполненного из стекла СЗС22, в интервале О < r < 1,5гоб представлены в таблицах фиг. 6-9 соответственно.
На цветовом графике (фиг. 5) проиллюстрированоизменениецветностималораэмерНОГО ОбЪЕКта(ГуВ2еи 0,03 ММ) ПО СраВНЕНИЮ с крупноразмерным объектом (гое1- 0,3 мм). Откуда видно, что цвет первого смещается в фиолетовую область спектра (точка
Ц), в то время как цвет последнего не выходит за область, выделенную пунктиром, т.е. совпадает с цветом источника — точка И, Цветовое различие Л Е объектов вычисляется по известным в цветоведении формулам (6):
Л E = ((U>* Uz ) + (Vi -Vz*) +
+(« «)2)1" где U* = 13L*(U-Uo); V* = 13L*(V — Vo); (13) 10
L* =25!7 з — 16; U—
Для исследуемых объектов Л Е - 9 порогов. На фиг. 2 приведен график зависимости
20 времени обнаружения от цветового различия ЛЕ, откуда видно, что для данных объектов цветовое различие ЛЕ = 9 определяет время обнаружения около 3,5 с.
На графиках — фиг. 10, 11, 72, 13 — пока25 заны распределения по радиусу яркости Ly для объектов rpe7= 0,3 мм и гоб2 = 0,03 мм при толщине фильтра L< = 0 мм (т.е. при отсутствии фильтра) и Lz = 3 мм, Незначительное уменьшение яркости для объекта гоб2= 0,03 мм при наличии фильтра обусловлено поглощением части энергии фильтром
СЗС22 толщиной 2 = 3 мм.
В прототипе относительная освещенность крупного объекта составляет 0,5, мелкого — 1, Их цветовое различие за счет контраста освещенностей равно Л Е =- 4 и соответствующее этому различию время обнаружения составляет около 8 с.
Таким образом, значительно сокращается время обнаружения объектов выделенных размеров при повышенной информационной емкости изображения, так как не происходит падения освещенности
45 крупных объектов.
Формула изобретения
1, Способ оптической обработки изображений, включающий формирование полихроматической освещающей волны, освещение этой волной транспаранта с обрабатываемым изображением, первое оптическое преобразование Фурье над волной, промодулированной транспарантом, фильтрацию диапазонов пространственных частот, второе оптическое преобразование
Фурье, в результате которого формируется обработанное изображение. регистрацию обработанного изображения, отличающийся тем, что, с целью повышения информационной емкости обработанного изобра181162Я
10 жения и сокращения времени обнэруже::ия обьектов иэображения, осуществляют фильтрацию диапазонов пространственных частот путем их цветовой селекции по закону
N 5
Щр) = 0 (р) g fclcc (" — ) — circ(— )1 k
Р Р ) транспаранта, о т л и ч а ю щ е е с я тем, ч;и пространственно-частотный фильтр Выполнен в виде одного или нескольких кольцевых фильтров иэ различных цветных стекол и имеет функцию пропускания.гф ф = Jcirc (>) — circ(-" -}1 +<4 — р где Оа(р) — амплитуда поля на длине волны до фильтрации;
U (p) — амплитуда поля на длине волны 1О после фильтрации; р- радиальная полярная координата в плоскости пространственно-частотного спектра; рк — внешний радиус к-й кольцевой зо- 15 ны пространственно-частотного спектра; ф — коэффициент пропускания по интенсивности к-й кольцевой зоны частотного спектра для длины волны 1, К вЂ” количество фильтруемых диапазо- 2О нов частот, причем ро = О.
2. Устройство для оптической обработки изображений, содержащее расположенные последовательно на одной оптической оси источник полихроматического излучения, осветительную систему, транспарант, фурье-объектив, и ро стра нствен но-частотный фильтр, фурье-преобразующий обьектив, плоскость регистрации изображения
Раайчс оА ъейта ),О«.ОЗ ми
Кааоаннаты ыве7тн<эсти Источника где p — радиальная полярная координата в плоскости пространственно-частотного фильтра; р,— внешний радиус к-й кольцевой зоны фильтра, — коэффициент пропускания по интенсивности k-й кольцевой эоны фильтра для длины волны;
N — количество кольцевых зон фильтра, причем po = О, а максимальный внешний радиус К-й (последней) кольцевой зоны фильтра
PN = макс 1 макс тф.о, ГДЕ 1 акс — МЭКСИМЭЛЬНЭЯ ПРОСтРЭНСтВЕННаЯ частота обрабатываемого изображения;
Яуакс МЭКСИМЭЛЬНЭЯ ДЛИНЭ ВОЛНЫ;
f y.о — фокусное расстояние фурье-объектива.
Таблица!
Х-0,2664, Y=0,2983, Z=0,4353.
Х Y о. 2649«2 0.-29985
О. 26492 О. 29980
О."26490 О ° 29968
О 264SS О ° 29947
О. .26484 О. 2991 9
0;26481 о.29834
0;26477 0.29843 .О 26473 0.29796
- 0.26469. 0.2974* 0 26467. 0.29694
:0»26467. О ° 29640 0«26469 0.29587
0.26474 0.29538 о 264аз о 29494
О .26498 0.29453
0«к26519 0 ° 29434
О, 26548 . 0.29425
О. 26586 0.29435
О, i 26637.. О ° 29469
О. 26702 .: . О. 29535
О"26786 i 0.29639
О 26893 " 0«29793
0;27029 0.30011 0;27204 0.30315
О ;27432, 0.30735
О.ДT739 0.31326
0-28168 0.32182
О;2881 1 0.33496
О . 29889 О. 35730
0 321 17 0.40318
Точка
2
3 ,4.
6
7 о, Ф
1,1
12 .13
14
16
17, 13 ,19
2021
22
23
24
26 .7
23
29
3Ь
0 ° 4 523
О. 43523
О. 4 542
О. 43565
0.43597
0.436 5
0.43681
О.4. 731
О. 4 734.
0.43339
О. 4 393, О. 43944 о. 43933
О. 44023
О. 44О44
0.44047
0.44023
0.4 979
О. 4 394
0.43763
0.43575
О.43314 .0.42960
0.42432
0.41332
0.40934
0.39649
0.37693 о.34381
О ° 27565
l . 34797 ооое
1 «««4««u J
1.8.1174
1.76735
i 10631
1.63167 1 ° 543322
1 ° 44540
1 ° 33375
1.22632 1.11056
О.99337
37852 О„7665
0.65982
0.55977
О;46759,о. 33413
0.30939
0.24505 о.13952
О.14294
0.10473
0 ° 07417
0.05042, 0.03255 о. о1965
0.01081
О;00515
0.00190
1811628
Табпица2
Радиус объекта: 0,30 мн
Координаты цветности источника: К 0.2884: Y 0,2983; 70,4353.
2, L>
О. 43515 1. 03414
0 ° 43976 1.00091
0.43496 0.98640
0.42954., 1.01490
0 ° 43768 1.01351
0.43957 0.97586
0. 43314 0;99392
0.43263 1.03266
0.43459 0 ° 99813.
0.44204 0.96197
0.43605 1.01507
О ° 42211 1.04277
0.44304 0.9Ь982
0.45060 0.95671
0.41431 1° . 05412
О ° 43058 1 ° 04280:
0 ° 46488 0.91082
0.42056 0 ° 97078.
0.41445 1.18975.
0.46550 0.92772:
0.44043 0.24441.
0.32590." 0 00170
0.167111 0.00641
0.76224 0.00026
0 ° 24608 О ° 00136
О ° 17350 О ° 00032
О. 46783. О. 00040
О. 46561 О; 00028 О. 03450, О. 00013 .0.4482» 0.00019
ТаблицаЗ
Раличс,объекта i .,0»;b3 ин:.
Коооаинаты иветмости. источника ф Х=0.2664: Y0,2983, Z0,4353.
)(Y Z
О" 2б253 О;.17863 О. 61885
О, 20253 О. 17861 О ° .61886
О"20254 О 17856" 0.61890
О. .20256 0.17849 :0.61895
ОГ20259 О. 17839 0.61902
О 20263 0.17828 0.61909
О 202Ь9 0.17817 О.Ü1914
0. 0277 О 17807 О. 61916, О 20287 О 17800 О. 61913. ,.0.20301 0.17797 0.61902 0320319 . 0.17802 О.Ь1879
Оь 20341 О. 17817 i О. 61841 ,О "20370 ; О. 17846 О. 61784
О 20406 О. 17893 0. 61700
". О."20452 . О. 17963 О. 61584 .О. 20510 О. 18063 . О. 61427
"О,20583 . О 18200 О. 61217
О.".20Ь75 0.18384 0.60941 ,О;.20792 О. 18629 О. 60580
O» 20940 0.18952 0»60107
О 21132 0.193ТВ 0.59490
;О;"21382, О. 19941 . 0.58677
О. 217 i 4 О. 2ОЬ89 О. 57597 О, 22164 О. 21697 . О. 5Ь139 О 22791 О. 23077 О» 54132 ,0.„23694 О. 25013 О. 51293
О» 25046 g О. 27804 О. 4Т150
;О; 27161 О ° 31938 О. 40901
:0.30581 0.38127 0.31293
: .О 35971 0.4*743 О. 17286
Точка
2., 4
6,,10
11,:
12 13 . .14
i5 16.;
17.
18.
19
20 21 .
22;
Я4 .
27
29
Точка
1. °
4, Ь
?.
8 ,9
1Î. .1,.1
:12
1,3 14;
: 1 5 .}6.
;17; 18
19,:
20, 21
;22
23 .24
25"
26
27, .28
Х
0.25503
0.26055
О ° 27058
0.26973
026249
0,26450
О ° 26977
О ° 26711
0.26409 О ° 26509
0.26822
0.26890
О. 26270
0. 26354
0 ° 27308
0.2Ь628
0.25790
О ° 27160
0.27223
:0.25729
0.264Ь4
С..33313
0 33702 О. 17557
0. 32865
О., 25482
0.28184
0.21743
0.46963
0.22255
0» 30982
0.29969
0.29445
0.30073
0.29983
0.29593
0 ° 29708
О. 30027
О. 30132
0.29286
О"29573
0.30900
Î;29426
0.31261
0 30313
О ° 27722
0.30784
0.31332 . 0.27721
0..29494
0.34098
0.49587
0.06219
0.42528
0, 57 168
0.25033 .О .31696
0..49586
О 32924
Ь
0.73750 О.?3397
0.72349
О ° 70631
0.68288
0.65378
0 ° Ь1975
О. 58158
О ° 54019
О. 49651
О. 45146
0.40599
0.36094
О ° 31712
0.27520
0,23578
0.19930
0.16610
0.13638
0.11021
0.08757
0.06834
0.05231 ,0.03921
О 02875
0.02059
0»01438
0»00980
0»00653
0»00427
Таблица4
Радиус объекта 2 О. 30 >1м
Кс>оодинаты иветнпсти источника: X=0 2 54;
Y=P,2983; Z=0,4353.
Фиг.1
Точка
2
4
5, 6
7::, 8
1 0
11:
12
13
14 ,15 .16, 17
18:
19. 20
21
2 Р
2«
24, 255 с6 ,л7
28 29
ЗО
О. 24397
О. 24535
О. 24981
0.25294
0.25406
0.25938
0 ° 265>4
О. 26787 .
О. 27176
О. 27453
О. 27599
0; 27842
0 27482
О. 27158 . 0.27
О. 26449
0.25451
0.25620
0.24441 .0 ° 23109 .0.26046
0.44954 .0.46762
0.37692
0.39608
0.41208
0 ..2. 5 6 1
0.32218
О 46652
0.40292
0.27423
0.26062
0.27192
28008
0.28540
0.29242
О.«0175
0 ° "1189
О. «»101.>
0.31201
0.52454
О. 1»" >В
О. -0154
С>. 1505
О. 29871
Г>. 270:",0
0 28280 г> 265"6
О. 22745
С>. 28608
О. 53620
0.45898
0.42052
0.50997
О ° 402 8
0.19329
0.24401 0.51675
0.51621
0 . 48 i 80, 0. 8 2666
0.48967 0,.79454.
О ° 48958 0.77556
0.47515. 0.82875. 0 ° 4658* 0. 88753
О ° 45521 . 0 ° 90693, О. 44224 О. 95620
0. 4=.038 1. 04073
0.41635 1 ° 077ВД,, 0.41535 1.0 498
0 41201 1.13508 О.79703, 1. l6934
Î. 41160 1 .,11774
О. 42689 . 1. 06968,О. 41 14 " . 1 ° 0754>.
О. 43*80; 1. 01 442,.C>. 47519 О. 86644
0:46099, 0.86988. О;4902 0.83 30 ;
0.54146 0;*4065
0.45346.. 0.23019, 0.01426 . 0.0 496
0.07БГ9 0.00894,.
0.20256 0.007>>
0.09="95 О.ОО:586
;. О. 1855 0.0002
:0.57310. 0 .00020.. 0..,4 =-81 0„0O028
0.0167Г "0.00053;
О. 08087 . О. 00 1 28
1811b28
Фиг,2, \
1811628
1811628
1811628
Составитель Д.Зубков
Техред M.Ìoðãåíòàë
Корректор Л;Ливринц
Редактор А.Купряков
Заказ 1467 Тираж Подписное
ВНИИПИ Государственного комитета по изобретениям и открытиям при ГКНТ СССР
113035, Москва, Ж 35, Раушская наб„4/5
Производственно-издательский комбинат "Патент", г. Ужгород, ул.Гагарина. 101
