RU2292060C1 - Способ наблюдения за воздушными объектами и поверхностью на базе бортовой рлс - Google Patents
Способ наблюдения за воздушными объектами и поверхностью на базе бортовой рлс Download PDFInfo
- Publication number
- RU2292060C1 RU2292060C1 RU2005120107/09A RU2005120107A RU2292060C1 RU 2292060 C1 RU2292060 C1 RU 2292060C1 RU 2005120107/09 A RU2005120107/09 A RU 2005120107/09A RU 2005120107 A RU2005120107 A RU 2005120107A RU 2292060 C1 RU2292060 C1 RU 2292060C1
- Authority
- RU
- Russia
- Prior art keywords
- matrix
- radar
- amplitudes
- azimuth
- elements
- Prior art date
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
Изобретение относится к радиолокации, а именно к радиолокационным системам наблюдения за воздушной обстановкой и поверхностью на базе бортовой РЛС, работающей в режиме "реального луча", с электронным сканированием. Технический результат заключается в повышении разрешающей способности по азимуту и углу места в заданных сечениях дальности с сохранением области обзора РЛС по азимуту и углу места и увеличение точности оценивания амплитуд сигналов в синтезированных элементах разрешения. Технический результат достигается тем, что в известном способе, основанном на работе в режиме реального луча с электронным сканированием, заключающийся в формировании матрицы радиолокационного изображения воздушной обстановки или поверхности в срезах дальности, при этом за счет быстрого электронного переключения луча РЛС смещают луч по азимуту и углу места построчно соответственно на величину n-й и m-й части ширины ДНА в зоне обзора, измеряют амплитуды сигналов отражения при каждом i,j-м положении луча и формируют из этих амплитуд матрицу измерений y(i,j), , , суммарного канала, которую далее обрабатывают, дополнительно формируют матрицу измерений y'(i,j), разностного канала, затем обрабатывают полученные матрицы для каждого i, j-го положения луча, при этом элементы матриц y(i+k, j+1) и y'(i+k, j+1), суммируют с весами h(k,l) и h'(k,l) и оценивают амплитуду x(i,j), указанные операции повторяют для всех i, j в зоне обзора и тем самым получают матрицу оценок амплитуд
Description
Изобретение относится к радиолокации, а именно к радиолокационным системам наблюдения за воздушной обстановкой и поверхностью на базе бортовой РЛС, работающей в режиме "реального луча" с электронным сканированием.
При наблюдении бортовой РЛС (БРЛС) за группой воздушных объектов или поверхностью в режиме реального луча осуществляется построчное сканирование лучом РЛС заданного сектора пространства путем последовательного смещения луча по азимуту и углу места на ширину диаграммы направленности антенны (ДНА). Четкость радиолокационного изображения (РЛИ) и точность определения угловых координат объектов в срезах дальности при таком способе наблюдения ограничена шириной ДНА. При наблюдении за одиночными воздушными объектами известны методы пеленгации для определения угловых координат (Финкельштейн М.И. Основы радиолокации: Учебник для вузов. М.: Радио и связь, 1983. 536 с.). Однако при наличии группы объектов в одном сечении дальности в пределах одной ДНА (тем более при наблюдении за поверхностью) такие методы не работают. Возникает проблема повышения разрешения РЛС совместно по азимуту и углу места в режиме "реального луча" путем синтезирования новых элементов разрешения с угловыми размерами, меньшими ширины ДНА, в заданных сечениях дальности.
Наиболее близким по технической сущности является способ синтезирования новых элементов разрешения по азимуту при переднем обзоре в режиме "реального луча" (Пат.2249832 РФ. Способ наблюдения за поверхностью на базе бортовой РЛС / В.К.Клочко, Г.Н.Колодько, В.И.Мойбенко, А.А.Ермаков (РФ). Заявка №2003126516. Приоритет 02.09.03), который заключается в следующем. Повышение разрешающей способности с расширением зоны обзора РЛС по азимуту и формирование матрицы радиолокационного изображения поверхности в режиме реального луча с электронным сканированием достигается за счет быстрого электронного переключения (смещения) луча РЛС по азимуту на величину n-й части ширины ДНА и обработки полученных при каждом положении луча амплитуд отраженных сигналов РЛС, которая заключается в следующем.
1. Амплитуды отраженных сигналов на выходе суммарного канала РЛС y1, y2, ..., yn, полученные при n первых положениях луча РЛС в данном i-м элементе разрешения по дальности суммируются с весами h1, h2, ..., hn, которые вычисляются по определенной методике. Результатом такой обработки является оценка амплитуды хn отраженного сигнала, соответствующего n-й части ДНА при первом положении луча:
2. При последующих смещениях луча на n-ю часть ДНА амплитуды сигналов, полученные при n последних положениях луча, суммируются с теми же весами, в результате чего последовательно находятся оценки xn+1, xn+2, ..., xN:
3. Оценки амплитуд xj (j=n, n+1, ..., N), найденные независимо в каждом i-м (i=1, 2, ..., M) элементе разрешения по дальности, располагают в М строк и N-n+1 столбцов и тем самым формируют матрицу радиолокационного изображения поверхности в виде совокупности амплитуд A(i,j), , сигналов, отраженных от соответствующих i-, j-x элементов поверхности.
Однако такой способ обладает следующими недостатками.
1. В указанном способе повышение разрешающей способности достигается только по одной угловой координате - азимуту.
2. Указанный способ не применим при наблюдении за группой воздушных объектов, расположенных в одном сечении (элементе разрешения) дальности и в пределах одной ДНА, т.е. не различимых по угловым координатам.
3. Использование данных только одного суммарного канала РЛС дает оценки амплитуд с ограниченной точностью.
Технический результат направлен на совместное повышение разрешающей способности по азимуту и углу места в заданных сечениях дальности с сохранением области обзора РЛС по азимуту и углу места и увеличение точности оценивания амплитуд сигналов в синтезированных элементах разрешения.
Технический результат предлагаемого технического решения достигается тем, что при наблюдении за воздушной обстановкой или поверхностью с помощью бортовой РЛС в режиме реального луча с электронным сканированием формируют матрицу радиолокационного изображения воздушной обстановки или поверхности в срезах дальности, при этом за счет быстрого электронного переключения луча РЛС смещают луч по азимуту и углу места построчно соответственно на величину n-й и m-й части ширины ДНА (на уровне 0,5 мощности) в зоне обзора, измеряют амплитуды сигналов отражения при каждом i, j-м положении луча и формируют из этих амплитуд матрицу измерений y(i,j), , суммарного канала, которую далее обрабатывают, отличающийся тем, что дополнительно формируют матрицу измерений y'(i,j), разностного канала, затем обрабатывают полученные матрицы для каждого i, j-го положения луча, при этом элементы матриц y(i+k, j+l) и y'(i+k, j+l), , , взятые относительно i, j в окне размера M×N, суммируют с весами h(k,l) и h'(k,l), найденными заранее, и оценивают амплитуду x(i,j), соответствующую центральной m, n-й части ДНА при i, j-м положении луча
указанные операции повторяют для всех i, j в зоне обзора и тем самым получают матрицу оценок амплитуд , , представляющую восстановленное радиолокационное изображение воздушной обстановки или поверхности в заданных элементах дальности с повышенным в несколько раз разрешением по угловым координатам.
Способ осуществляется следующим образом.
Повышение разрешающей способности РЛС по азимуту и углу места в заданной угломерной области воздушного пространства или поверхности в данном элементе разрешения дальности и формирование матрицы радиолокационного изображения поверхности или воздушной обстановки в режиме реального луча с электронным сканированием достигается за счет быстрого электронного переключения (смещения) луча РЛС по азимуту и углу места построчно соответственно на величину n-й и m-й части ширины ДНА и обработки полученных при сканировании лучом зоны обзора амплитуд отраженных сигналов РЛС на выходе не только разностного, но и суммарного канала РЛС, которая заключается в следующем.
1. Амплитуды отраженных сигналов на выходе суммарного и разностного каналов РЛС y(i,j) и y'(i,j), полученные при каждом i, j-м положении луча, формируют в матрицы y(i,j) и у'(i,j), , .
2. Относительно каждого i, j-го элемента искомой матрицы изображения строится окно размера M×N и амплитуды суммарного у(i+k, j+l) и разностного y'(i+k, j+l) каналов РЛС, взятые в этом окне при , , суммируют с весовыми коэффициентами h(k,l) и h'(k,l), расчет которых дан ниже, в результате чего получают оценку амплитуды i, j-го элемента дискретизации поля отражения:
3. Полученная матрица , , представляет восстановленное в данном сечении дальности радиолокационное изображение в i, j-x синтезированных элементах разрешения по азимуту и углу места, размеры которых в несколько раз меньше ширины ДНА. По совокупности всех элементов разрешения дальности формируется трехмерное изображение области воздушного пространства или поверхности.
Расчет весовых коэффициентов сводится к следующему.
Сканирование лучом РЛС осуществляется построчно: фиксируется i-я строка и луч последовательно перемещается по строке (меняется j) на n-ю часть ширины ДНА по азимуту. Затем осуществляется переход к следующей (i+1)-й строке путем изменения (смещения) положения луча на m-ю часть ширины ДНА по углу места (меняется i), после чего сканирование по азимуту (по j) повторяется. Амплитуда y(i,j) сигнала на выходе суммарного канала в k-м фиксированном элементе разрешения дальности при i, j-м положении луча по азимуту и углу места после некоторых предварительных преобразований представляет сумму амплитуд x(i+k, j+l) сигналов поля отражения, взятых с коэффициентами ДНА α(k,l) по всем k, l-м элементам дискретизации в пределах ДНА, с учетом помехи p(i,j):
Модель (1) аналогична модели пространственного смазывания РЛИ по строкам и столбцам. Коэффициенты α(k,l) ДНА суммарного канала представляют функцию с разделенными переменными: α(k,l)=α1(k)·α2(l) и модель измерения (1) допускает разделение
и последующую двухэтапную обработку: вначале по строкам, затем по столбцам, направленную на восстановление искомого поля отражения x(i,j). Это позволяет использовать известные одномерные алгоритмы оценивания параметров поля по интегральным (суммарным) наблюдениям (например, Клочко В.К., Чураков Е.П., Фатьянов С.О. Калмановский алгоритм восстановления смазанного радиолокационного изображения // Известия вузов. Радиоэлектроника. 2004. Том 47. №9-10. С.54-59). В общем случае оценка параметрах x(i,j) имеет вид
где h1(k), h2(l) - весовые коэффициенты, вычисляемые заранее по определенной методике; М≥m, N≥n. Так как число MN измерений, участвующих в оценивании x(i,j), меньше числа (M+m-1)(N+n-1) оцениваемых параметров, входящих в состав этих измерений, то точность оценки (3) недостаточна для различения нескольких близко расположенных точечных объектов (на расстоянии размера элемента дискретизации).
Для увеличения точности оценивания поля x(i,j) в элементах дискретизации i, j, приводящему к повышению разрешающей способности по угловым координатам, предлагается дополнительно с суммарным каналом использовать данные разностного канала. Модель амплитуды сигнала на выходе разностного канала аналогична (1):
но отличается тем, что коэффициенты β(k,l) ДНА разностного канала не описываются функцией с разделенными переменными: β(k,l)≠β1(k)/β2(l), и двухэтапная процедура обработки для разностного канала оказывается неприменимой.
Для получения оценки по данным суммарного и разностного каналов РЛС рассматриваются измерения y(i,j), y'(i,j) в M×N-окрестности i, j-го элемента, которые удобно пронумеровать в сквозном порядке (построчно) и представить в составе 2МN-вектора
где М≥m, N≥n. Параметры x(i,j), входящие в состав 2MN измерений, представляются вектором
где S=M+m-1, L=N+n-1. Помехи p(i,j) и p'(i,j) собираются в 2MN-векторе P(i,j)=(p1, р2, ..., р2MN)T. Тогда модель измерений (1), (4) представляется в следующей матричной форме:
где А - матрица размерности 2MN×SL, сформированная в соответствии с (1), (4) из двумерных коэффициентов ДНА первого канала αij, , и второго канала βij, , :
Уравнение (5) в развернутом виде представляет систему 2MN уравнений с (M+m-1)(N+n-l) неизвестными Х. При определенных значениях М и N получается избыточность измерений по отношению к оцениваемым параметрам: 2MN>(M+m-1)(N+n-1), что позволяет находить оценки искомых параметров Х методом наименьших квадратов (МНК). В матричной форме минимизация функции по приводит к регуляризованным МНК-оценкам
где δ - параметр регуляризации, необходимый для обращения плохо обусловленной матрицы АTА; Е - единичная матрица; Н=(АTА)-1АT - матрица весовых коэффициентов; Т - символ транспонирования.
При обращении АTА удобно использовать рекуррентную процедуру:
где ak - k-я строка матрицы А коэффициентов ДНА.
В качестве оценки элемента поля x(i.j), расположенного в центре M×N-окна, принимается центральный элемент вектора с наименьшей дисперсией ошибки оценивания, который вычисляется по формуле
Оценка представляет амплитуду сигнала отражения в m, n-й части ДНА (элементе дискретизации), соответствующей центру ij-го положения луча РЛС (синтезированного элемента разрешения по азимуту и углу места). При вычислении используется центральная строка матрицы H, соответствующая минимуму дисперсии ошибки оценивания ковариационной матрицы . Алгоритм оценивания представляет линейную комбинацию измерений, взятых с заранее просчитанными весовыми коэффициентами и его удобно представить в виде
где h1, h2, ..., h2MN - соответствующая строка матрицы Н; y1, y2, ..., y2MN - измерения суммарного и разностного каналов, причем М≥m, N≥n и принимается , если . Возможно получение сразу нескольких оценок , , для одного вектора измерений Y путем использования соответствующих строк матрицы Н, однако для обеспечения одинаковой точности этих оценок следует увеличить число измерений 2MN: М>m, N>n.
Возвращаясь к первоначальной нумерации (по строкам и столбцам), алгоритм (9) принимает вид
Точность оценок (10) описывается корреляционной матрицей ошибок оценивания . Свойство матрицы КX при заданных коэффициентах αij и βij полностью определяется параметрами m, n, М, N. В табл.1 представлены следующие расчетные величины: среднеквадратическое отклонение (СКО) ошибки оценивания σ=σ[Δх], взятое из корреляционной матрицы КX, и оценка СКО найденная путем моделирования, при σP=1 в зависимости от размеров M×N окна измерений для двух и одного каналов. Моделировались 4 близко расположенных точечных объекта в пределах ДНА размера в m×n=3×3 элементов дискретизации. В табл.2 показано слева направо в окне m×n: искомое поле x(i,j), измерения первого канала y(i,j) и оценки , полученные с помощью (9) при m×n=3×3 в окне измерений M×N=5×5 для двухканальной обработки с округлением до целого. Для сравнения в табл.3 даны результаты, полученные для одноканальной обработки. Использовались следующие аппроксимации коэффициентов ДНА суммарного и разностного каналов: αij=exp[-(i2+j2)], βij=0,7ехр[-1.5(ρ-1)2], . Видно, что точность двухканальной обработки существенно выше одноканальной.
Алгоритм (9) или (10) при наблюдении за воздушными объектами реализуется в тех элементах разрешения по дальности, в которых зафиксирован отраженный сигнал. При наблюдении за поверхностью восстанавливается трехмерное изображение по совокупности всех k-x срезов дальности в области обзора или обрабатываются по предложенной методике амплитуды сигналов, принятых в заданном диапазоне дальности (imin, imax), накрывающем наблюдаемый участок поверхности. Пространственные элементы дискретизации (i, j, k) представляют синтезированные элементы разрешения, угловые размеры которых в несколько раз меньше ширины ДНА.
Предлагаемый способ позволяет в несколько раз повысить разрешающую способность РЛС по азимуту и углу места в сечениях дальности в режиме "реального луча" с сохранением области обзора РЛС по азимуту и углу места и сформировать матрицу радиолокационного изображения воздушной обстановки или поверхности в виде совокупности амплитуд , , сигналов, отраженных от соответствующих i, j-x элементов угломерного пространства в k-x элементах разрешения дальности, которая позволяет наблюдать на экране индикатора группу воздушных объектов или поверхность (объекты на поверхности) в условиях отсутствия оптической видимости, что повышает безопасность полетов и эффективность решения поставленных перед летчиком задач.
Claims (1)
- Способ наблюдения за воздушными объектами и поверхностью на базе бортовой РЛС, основанный на работе в режиме реального луча с электронным сканированием, заключающийся в формировании матрицы радиолокационного изображения воздушной обстановки или поверхности в срезах дальности, при этом за счет быстрого электронного переключения луча РЛС смещают луч по азимуту и углу места построчно соответственно на величину n-й и m-й части ширины ДНА в зоне обзора, измеряют амплитуды сигналов отражения при каждом i-м, j-м положении луча и формируют из этих амплитуд матрицу измерений y'(i,j), , суммарного канала, которую далее обрабатывают, отличающийся тем, что дополнительно формируют матрицу измерений y'(i,j), , разностного канала, затем обрабатывают полученные матрицы для каждого i,j-го положения луча, при этом элементы матриц y(i+k,j+1) и y'(i+k,j+1), , взятые относительно i,j в окне размера M×N, суммируют с весами h(k,l) и h'(k,l), найденными заранее, и оценивают амплитуду x(i,j), соответствующую центральной m-й, n-й части ДНА при i-й, j-м положении луча
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2005120107/09A RU2292060C1 (ru) | 2005-06-28 | 2005-06-28 | Способ наблюдения за воздушными объектами и поверхностью на базе бортовой рлс |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2005120107/09A RU2292060C1 (ru) | 2005-06-28 | 2005-06-28 | Способ наблюдения за воздушными объектами и поверхностью на базе бортовой рлс |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2292060C1 true RU2292060C1 (ru) | 2007-01-20 |
Family
ID=37774762
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2005120107/09A RU2292060C1 (ru) | 2005-06-28 | 2005-06-28 | Способ наблюдения за воздушными объектами и поверхностью на базе бортовой рлс |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2292060C1 (ru) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2558674C1 (ru) * | 2014-03-04 | 2015-08-10 | Акционерное общество "Концерн радиостроения "Вега" | Способ автоматизированного выявления компактных групп взаимодействующих воздушных объектов |
RU2600573C1 (ru) * | 2015-11-23 | 2016-10-27 | Общество С Ограниченной Ответственностью "Научно-Технический Центр "Версия" (Ооо "Нтц "Версия") | Способ восстановления изображений объектов по разреженной матрице радиометрических наблюдений |
RU2604720C1 (ru) * | 2015-12-28 | 2016-12-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Рязанский государственный радиотехнический университет" | Способ восстановления изображений при неизвестной аппаратной функции |
RU2618088C1 (ru) * | 2016-01-27 | 2017-05-02 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Военно-космическая академия имени А.Ф. Можайского" Министерства обороны Российской Федерации | Способ оптимального восстановления изображений в радиолокационных системах дистанционного зондирования Земли в телескопическом режиме |
RU2622899C1 (ru) * | 2016-05-25 | 2017-06-21 | Закрытое акционерное общество "Рязанская радиоэлектронная компания" (ЗАО "РРК") | Способ определения аппаратной функции радиометра |
RU2624460C1 (ru) * | 2016-01-27 | 2017-07-04 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Военно-космическая академия имени А.Ф. Можайского" Министерства обороны Российской Федерации | Способ оптимального восстановления изображений в радиолокационных системах дистанционного зондирования Земли |
-
2005
- 2005-06-28 RU RU2005120107/09A patent/RU2292060C1/ru not_active IP Right Cessation
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2558674C1 (ru) * | 2014-03-04 | 2015-08-10 | Акционерное общество "Концерн радиостроения "Вега" | Способ автоматизированного выявления компактных групп взаимодействующих воздушных объектов |
RU2600573C1 (ru) * | 2015-11-23 | 2016-10-27 | Общество С Ограниченной Ответственностью "Научно-Технический Центр "Версия" (Ооо "Нтц "Версия") | Способ восстановления изображений объектов по разреженной матрице радиометрических наблюдений |
RU2604720C1 (ru) * | 2015-12-28 | 2016-12-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Рязанский государственный радиотехнический университет" | Способ восстановления изображений при неизвестной аппаратной функции |
RU2618088C1 (ru) * | 2016-01-27 | 2017-05-02 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Военно-космическая академия имени А.Ф. Можайского" Министерства обороны Российской Федерации | Способ оптимального восстановления изображений в радиолокационных системах дистанционного зондирования Земли в телескопическом режиме |
RU2624460C1 (ru) * | 2016-01-27 | 2017-07-04 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Военно-космическая академия имени А.Ф. Можайского" Министерства обороны Российской Федерации | Способ оптимального восстановления изображений в радиолокационных системах дистанционного зондирования Земли |
RU2622899C1 (ru) * | 2016-05-25 | 2017-06-21 | Закрытое акционерное общество "Рязанская радиоэлектронная компания" (ЗАО "РРК") | Способ определения аппаратной функции радиометра |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2292060C1 (ru) | Способ наблюдения за воздушными объектами и поверхностью на базе бортовой рлс | |
EP1286180B1 (en) | Periodic repetition interval staggered post-doppler adaptive monopulse processing for detection and location of a moving target in ground clutter | |
JP6660857B2 (ja) | レーダーシステムを用いて関心領域の画像を生成する方法 | |
US6226321B1 (en) | Multichannel parametric adaptive matched filter receiver | |
CN105137425B (zh) | 基于卷积反演原理的扫描雷达前视角超分辨方法 | |
JP5134878B2 (ja) | レーダ装置 | |
CN107748362A (zh) | 一种基于最大锐度的线阵sar快速自聚焦成像方法 | |
CN113238225B (zh) | 一种用于动目标检测的雷达 | |
RU2368917C1 (ru) | Способ формирования изображений в многоканальных ртлс и рлс | |
CN113655478A (zh) | 一种成像方法及装置 | |
RU2379706C2 (ru) | Способ повышения разрешающей способности радиотепловых изображений | |
CN110378896A (zh) | 基于极化相干性的TomoSAR植被病虫害监测方法及装置 | |
KR20190036325A (ko) | 압축 센싱을 이용한 isar 영상의 자동 초점 및 수직 거리 스케일링 장치와 그 방법 | |
CN104297735B (zh) | 基于先验道路信息的杂波抑制方法 | |
RU2284548C1 (ru) | Способ наблюдения за поверхностью и воздушной обстановкой на базе бортовой рлс | |
RU2379705C2 (ru) | Способ двухэтапного восстановления изображений в многоканальных радиолокационных и радиотеплолокационных станциях | |
CN113608192B (zh) | 一种探地雷达远场定位方法、装置及计算机可读存储介质 | |
RU2373552C1 (ru) | Способ формирования радиолокационного изображения в бортовых системах радиовидения | |
RU2539558C1 (ru) | Способ формирования трехмерного изображения земной поверхности и воздушной обстановки с помощью антенной решетки | |
CN114265058A (zh) | Mimo雷达目标测角方法、装置、电子设备及存储介质 | |
RU2300780C1 (ru) | Способ измерения высоты объектов на поверхности при получении трехмерного радиолокационного изображения поверхности с объектами на базе бортовой рлс маловысотного полета | |
RU2612323C1 (ru) | Способ восстановления изображений в двухканальной сканирующей системе | |
CN110399832A (zh) | 基于相干性的TomoSAR植被病虫害监测方法及装置 | |
Bahçeci et al. | Transmitter and Receiver Design in Compressive Sensing based Direction of Arrival Estimation | |
CN110378894B (zh) | 基于相关性的TomoSAR植被病虫害监测方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20070629 |