BRPI0808283A2 - "método implementado em computador para registrar imagens de radar de abertura sintética (sar) e sistema computadorizado para registrar imagens de radar de abertura sintética (sar)" - Google Patents

Description

MÉTODO IMPLEMENTADO EM COMPUTADOR PARA REGISTRAR IMAGENS DE RADAR DE ABERTURA SINTÉTICA (SAR) E SISTEMA COMPUTADORIZADO PARA REGISTRAR IMAGENS DE RADAR DE ABERTURA SINTÉTICA (SAR)

A presente invenção refere-se ao campo de

processamento de imagens, e mais particularmente ao processamento de imagens de radar de abertura sintética (SAR).

A resolução de dados SAR não é comparável à 10 resolução dos dados eletro-ópticos (EO) . Os sensores EO incluem dispositivos fotográficos e outros dispositivos ópticos de formação de imagem, tais como coletores de detecção e variação da luz (LIDAR) . Os sensores EO são passivos, uma vez que eles capturam a refletividade da luz 15 das cenas para formar imagens fotográficas das mesmas. No entanto, os sensores EO são limitados pelo horário do dia e pelas condições atmosféricas.

Um radar de abertura sintético (SAR) é vantajoso, uma vez que as imagens podem ser adquiridas de dia ou à

2 0 noite, bem como sob condições meteorológicas inclementes. Um

SAR é ativo, uma vez que ele registra a radiação para retrodifusão dos sinais de radiofreqüência (RF) para gerar imagens SAR. Cada célula de resolução SAR tem geralmente muitos difusores. As fases dos sinais de retorno destes difusores 25 são distribuídas aleatoriamente, e a interferência resultante causa salpicos.

Os salpicos conferem uma aparência granulosa na imagem detectada que é vista no final, e desse modo uma resolução mais baixa quando comparada a uma imagem EO. Os

3 0 salpicos impõem uma limitação significativa na exatidão das

medições que podem ser feitas. Por exemplo, a mensuração é frequentemente inclusiva nos dados SAR. A interferência de lobo lateral também cria uma aparência ruidosa nos dados SAR. Além disso, o mau funcionamento do hardware ou a interferência de rádio podem diminuir a fidelidade dos dados SAR.

Os dados SAR estão sendo tratados atualmente com alguma forma de apodização em que os lobos principal e laterais são removidos. No entanto, a apodização faz com que os dados SAR pareçam binários. Isto também resulta em uma imagem detectada com uma aparência granulosa. Os dados SAR também estão sendo tratados com filtros de baixa passagem, tal como a ponderação de Taylor. No entanto, os difusores podem ficar borrados conjuntamente, resultando uma resolução reduzida. Em conseqüência das abordagens atuais utilizadas para tratar os dados SAR, determinadas aplicações de análise podem ser inclusivas, incluindo o registro, a detecção de estrada, a detecção de mudança, a extração de elevação e a mensuração.

Para as imagens SAR que contêm salpicos, um objetivo de realce consiste em remover os salpicos sem destruir características importantes da imagem. 0 brilho de

2 0 um pixel é determinado não somente pelas propriedades dos

difusores na célula de resolução, mas também pelas relações de fases entre os retornos destes difusores. Em determinadas aplicações, no entanto, a remoção dos salpicos pode ser contra produtiva. Um exemplo em que a preservação dos 25 salpicos é importante é onde a detecção das características é da mesma escala que os padrões de salpico. Uma técnica conhecida para remover salpicos dos dados SAR, bem com realçar a resolução, é a aplicação de algoritmos de difusão anisotrópica.

3 0 Uma abordagem para remover salpicos dos dados SAR é

apresentada no artigo intitulado "Difusão Anisotrópica de Redução de Salpico", da autoria de Yu et al. Uma abordagem de equação diferencial parcial (PDE) é utilizada para a remoção do salpico. Particularmente, um espaço em escala de imagem é gerado, o qual é um conjunto de imagens filtradas que variam de finas a grossas. Uma outra abordagem é apresentada no artigo intitulado "Remoção de Salpico por Difusão 5 Anisotrópica Para Imagens SAR de Alta Resolução", da autoria de Xi et al. Um algoritmo de filtragem de difusão não-linear baseado em um esquema de distinção, isto é, um esquema de separação de operador aditivo (AOS), é aplicado nos dados de imagens distintas. Embora ambas estas abordagens resultem em 10 uma melhora na resolução dos dados SAR pela redução do ruído e preservação das bordas, ainda existe uma demanda para fazer com que os dados SAR se pareçam mais como dados EO de alta resolução.

Em vista dos antecedentes expostos acima, portanto, um objetivo da presente invenção consiste em melhorar a resolução de dados SAR para que se pareçam mais como dados EO.

Este e outros objetivos, características e vantagens de acordo com a presente invenção são providos por

2 0 um método implementado em computador para processamento das

imagens de radar de abertura sintética (SAR), o qual compreende a determinação do ruído em uma imagem SAR a ser processada, a seleção de um limite de ruído para a imagem SAR com base no ruído determinado, e o ajuste matemático de um 25 algoritmo de difusão anisotrópica com base no limite de ruído selecionado. 0 algoritmo de difusão anisotrópica ajustado é aplicado à imagem SAR.

0 ruído pode ser determinado com base na análise estatística dos valores do gradiente da imagem SAR. A análise

3 0 estatística pode ser baseada em um desvio padrão dos valores

do gradiente mais uma constante, por exemplo. Alternativamente, o ruído pode ser determinado com base em um esquema de janelas de Fourier ou em uma decomposição de ondulações.

0 algoritmo de difusão anisotrópica pode ser baseado em uma equação de calor que compreende um termo nãoconstante. 0 ajuste matemático do algoritmo de difusão 5 anisotrópica pode compreender o ajuste do termo não-constante com base no limite de ruído selecionado. Ao ajustar o termo não-constante, é vantajosamente possibilitado para que a equação de calor seja adaptada aos dados SAR que estão sendo processados. Em conseqüência disto, uma classe de funções

pode ser criada para obter os resultados desejados, em que cada função corresponde aos dados específicos SAR que estão sendo processados. As séries de dados SAR discrepantes podem ser melhor processadas ao melhorar a resolução da imagem SAR visualizada.

Um outro aspecto da invenção refere-se a um

algoritmo de difusão anisotrópica complexo. As equações para o algoritmo de difusão anisotrópica acima descrita foram rederivadas para dados complexos de modo que o processo ficasse agora não-linear. Em termos de dados complexos, os

2 0 componentes reais e imaginários de um conjunto de dados SAR

são processados em um mesmo tempo. Uma vez que os componentes reais e imaginários do conjunto de dados SAR estão sendo tratados como um objeto complexo, o algoritmo de difusão anisotrópica complexo pode retirar o conteúdo da cena dos

dados extremamente ruidosos, o que melhora por sua vez a resolução da imagem SAR visualizada.

Um método implementado por computador para processar imagens SAR complexas compreende a recepção de um conjunto de dados SAR complexos para uma imagem SAR que

3 0 compreende uma pluralidade de pixels, e a aplicação do

algoritmo de difusão anisotrópica complexo à série de dados SAR complexos. 0 conjunto de dados SAR complexos compreende uma parte real e imaginária para cada pixel. Se o conjunto de dados SAR complexos for recebida no espaço da freqüência, o espaço da freqüência é convertido em espaço da imagem. 0 espaço da freqüência corresponde à fase e à intensidade para cada pixel, e o espaço da imagem corresponde à fase e à 5 amplitude para cada pixel.

0 algoritmo de difusão anisotrópica complexo também pode ser utilizado no processamento interferométrico de dados SAR, particularmente para medições do abrandamento em áreas urbanas, por exemplo. 0 abrandamento é um deslocamento do 10 terreno em que a elevação da superfície da terra está diminuindo em relação ao nível do mar. Um método implementado por computador para processamento interferométrico de imagens SAR compreende a recepção de primeiro e segundo conjuntos de dados SAR complexos de uma mesma cena, com o segundo conjunto 15 de dados SAR complexos sendo deslocado em fase com respeito ao primeiro conjunto de dados SAR complexos. Cada conjunto de dados SAR complexos pode compreender uma pluralidade de pixels. Um interferograma é formado com base no primeiro e segundo conjuntos de dados SAR complexos para prover uma 20 diferença de fase entre os mesmos. O algoritmo de difusão anisotrópica complexo é aplicado ao interferograma, em que o interferograma compreende uma parte real e imaginária para cada pixel. Um filtro de choque é aplicado ao interferograma.

0 algoritmo de difusão anisotrópica complexo mitiga 25 localmente o ruído enquanto que ao mesmo tempo preserva as descontinuidades da cena no interferograma. O filtro de choque é utilizado para remover borrões da imagem. A suavização não-linear através do algoritmo anisotrópica complexo e do realce do limite através do filtro de choque

3 0 aumenta a exatidão e a qualidade da medição da diferença de fase. Uma vez que o abrandamento é medido ao utilizar objetos persistentes na cena (isto é, edifícios) como pontos de referência, a melhora da qualidade do limite dos objetos persistentes melhora a medição do abrandamento.

0 método também pode compreender a execução de um desenrolamento de fase variacional bidimensional no interferograma após a aplicação do filtro de choque. A 5 diferença de fase entre as duas imagens SAR registradas é relacionada a uma quantidade física desejada de interesse, tal como a topografia da superfície. A diferença de fase pode ser registrada apenas no modulo 2B, e as técnicas interferométricas atuais recuperam principalmente a fase 10 absoluta (a fase desmembrada) daquela registrada (fase enrolada) ao utilizar valores distintos que têm uma tendência de suavizar os dados.

O algoritmo de desenrolamento de fase variacional de acordo com a presente invenção pode ser baseado em uma 15 função de custo para controlar a suavização. Ao invés de prover uma suavização global baseada nas propriedades dos dados, o algoritmo de desenrolamento de fase variacional deixa bordas intactas e suaviza seletivamente a área adjacente às bordas. Em conseqüência disto, o processamento 20 interferométrico dos dados SAR com base no algoritmo de difusão anisotrópica complexo, no filtro de choque e no desenrolamento de fase variacional melhoram coletivamente a qualidade do limite, o que melhora por sua vez a medição do abrandamento.

0 algoritmo de difusão anisotrópica também pode ser

utilizado na compressão e descompressão de imagens SAR. Uma vantagem da aplicação do algoritmo de difusão anisotrópica é que o tamanho da imagem SAR após a compressão resulta em um arquivo com um tamanho menor, independentemente do esquema da

3 0 compressão utilizado. Em termos de descompressão de uma imagem SAR, ao comprimir dinamicamente os dados SAR, quantificar esses dados, e então descomprimir com o algoritmo de difusão anisotrópica um arquivo de tamanho menor também é obtido. Quanto maior a escala dinâmica melhor a relação de compressão. Em conseqüência disto, o armazenamento e a transmissão das imagens SAR comprimidas e descomprimidas com base no algoritmo de difusão anisotrópica ocupam menos espaço 5 e largura de faixa.

0 algoritmo de difusão anisotrópica também pode ser utilizado na extração de elevação e no registro para imagens SAR. Um método implementado por computador para registrar imagens SAR compreende a seleção de primeira e segunda 10 imagens SAR a serem registradas, o processamento individual da primeira e segunda imagens SAR selecionadas com um algoritmo de difusão anisotrópica, e o registro da primeira e segunda imagens SAR após o processamento. O método também pode compreender a aplicação de um filtro de choque às 15 respectivas primeira e segunda imagens SAR processadas antes do registro. Este esquema propicia uma exatidão maior para o registro da imagem SAR, o que permite por sua vez que os dados da elevação sejam melhor extraídos com base nas imagens SAR registradas.

Ainda uma outra aplicação do algoritmo de difusão

anisotrópica e de um filtro de choque é com respeito à extração de vetor e estrada para a classificação material. Um método implementado por computador para a extração do vetor em imagens SAR compreende a seleção de uma imagem SAR para a 25 extração do vetor, o processamento da imagem SAR selecionada com um algoritmo de difusão anisotrópica, e a extração os dados do vetor com base na imagem SAR processada. O filtro de choque pode ser aplicado à imagem SAR processada antes da extração. Os dados da imagem da estrada podem então ser

3 0 extraídos com base nos dados de vetor extraídos. Os dados são delineados desse modo ao utilizar um esquema coerente do algoritmo de difusão anisotrópica. Este esquema propicia uma exatidão maior para a extração da estrada. A FIGURA 1 é um diagrama de blocos esquemático da coleta e processamento de imagens SAR de acordo com a presente invenção.

As FIGURAS 2 e 3 são os respectivos gráficos da distribuição de intensidades de pixel e de valores de gradiente para um conjunto de dados SAR de acordo com a presente invenção.

A FIGURA 4 é uma imagem do conjunto de dados SAR que corresponde aos gráficos mostrados nas FIGURA 2 e 3.

A FIGURA 5 é uma imagem dos gradientes que

correspondem aos gráficos mostrados nas FIGURAS 2 e 3.

A FIGURA 6 é uma imagem dos gradientes realçados com um esquema de ponderação de acordo com a técnica anterior.

A FIGURA 7 é a mesma imagem dos gradientes

mostrados na FIGURA 6 realçados com um algoritmo de difusão anisotrópica de acordo com a presente invenção.

A FIGURA 8 é uma imagem original antes da filtragem de acordo com a presente invenção.

2 0 As FIGURAS 9, 10 e 11 são imagens que correspondem

à imagem original mostrada na FIGURA 8 após a filtragem de acordo com a técnica anterior.

A FIGURA 12 é uma imagem que corresponde à imagem original mostrada na FIGURA 8 após a filtragem com um algoritmo de difusão anisotrópica de acordo com a presente invenção.

A FIGURA 13 é uma imagem original antes da filtragem de acordo com a presente invenção.

A FIGURA 14 é uma imagem que corresponde à imagem original mostrada na FIGURA 13 após a filtragem com um algoritmo de difusão anisotrópica complexo de acordo com a presente invenção.

A FIGURA 15 é um fluxograma que ilustra o processamento interferométrico não-linear de dados SAR para medições do abrandamento de acordo com a presente invenção.

A FIGURA 16 é uma imagem multi-espectral de vetor aproximado original antes da aplicação de um filtro de choque de acordo com a presente invenção.

A FIGURA 17 é uma imagem que corresponde à imagem original mostrada na FIGURA 16 após a aplicação do filtro de choque de acordo com a presente invenção.

As FIGURAS 18A e 18B são imagens bidimensionais e

tridimensionais de uma cena original antes do processamento interferométrico de acordo com a presente invenção.

As FIGURAS 19A-22B são imagens bidimensionais e tridimensionais que correspondem à cena original mostrada nas FIGURAS 18A e 18B que ilustram vários estágios da aplicação

do processamento SAR interferométrico de acordo com a presente invenção.

A FIGURA 23 é uma imagem bidimensional de cima para baixo de uma cena original antes do processamento interferométrico de acordo com a presente invenção.

As FIGURAS 24-27 são imagens de cima para baixo

bidimensionais que correspondem à cena original mostrada na FIGURA 23 ilustrando vários estágios da aplicação de processamento SAR interferométrico de acordo com a presente invenção.

2 5 A FIGURA 28 é uma imagem original antes da

filtragem de acordo com a presente invenção.

A FIGURA 2 9 é uma imagem que corresponde à imagem original mostrada na FIGURA 28 após a aplicação de um filtro de Gauss de acordo com a técnica anterior.

3 0 A FIGURA 3 0 é uma imagem que corresponde à imagem

original mostrada na FIGURA 28 após múltiplas iterações de filtragem com um algoritmo de difusão anisotrópica complexo de acordo com a presente invenção. As FIGURAS 31-33 são imagens que ilustram a compressão de dados SAR de acordo com a presente invenção.

As FIGURAS 34-36 são imagens que ilustram a descompressão de dados SAR de acordo com a presente invenção.

As FIGURAS 37-44 são imagens e gráficos que

ilustram o registro de imagens SAR de acordo com a presente invenção.

As FIGURAS 45-48 são imagens que ilustram a extração da estrada de acordo com a presente invenção.

A presente invenção será descrita agora de maneira

mais abrangente em seguida com referência aos desenhos em anexo, nos quais são mostradas as realizações preferidas da invenção. A presente invenção pode, no entanto, ser incorporada em muitas formas diferentes e não deve ser 15 interpretada como limitada às realizações aqui apresentadas. Ao invés disto, estas realizações são fornecidas de modo que esta descrição seja plena e completa, e irão indicar completamente o âmbito da invenção aos elementos versados na técnica. As referências numéricas semelhantes referem-se aos 20 elementos semelhantes por toda parte.

Com referência inicialmente à FIGURA 1, um radar de abertura sintética (SAR) 50 coleta os dados SAR e um sistema implementado por computador 60 processa os dados SAR. 0 SAR 50 é suportado por uma plataforma transportada por via aérea 25 52, tal como um avião, sobre uma área de interesse. A plataforma transportada por via aérea 52 também poderia ser baseada no espaço.

A área de interesse ilustrada é uma área urbana, tal como uma cidade, a qual inclui uma série de edifícios 54.

3 0 Na área urbana também são incluídas características relativamente pequenas tais como as árvores 56 e as estradas 58, por exemplo, em comparação aos edifícios 54. Alternativamente, a área de interesse poderia ser uma área rural, com muito poucos, quando então nenhum edifício 54.

Os elementos versados na técnica irão apreciar que uma imagem SAR é recebida primeiramente como um valor complexo antes de ser convertida em um valor real para a 5 visualização. A imagem SAR é recebida inicialmente no espaço k que inclui a fase e a intensidade. Ao tomar a transformação de Fourier inversa dos dados retornados no espaço k, um espaço da imagem é gerado. O espaço da imagem inclui a fase e a amplitude. Uma vez que o espaço da imagem inclui as 10 informações da fase e da amplitude, ele apresenta salpicos. Para ver a imagem SAR, a fase é removida do espaço da imagem ao tomar a magnitude dos dados no espaço da imagem. Isto produz uma imagem detectada ou real para a visualização.

Uma vez que as imagens SAR são coletadas, elas 15 podem ser armazenadas em uma mídia de armazenamento 70, tal como um disco magnético, por exemplo, para a transferência a um computador 62. Dentro do computador 62, as imagens SAR podem ser armazenadas como parte de um banco de dados das imagens SAR a serem processadas. Naturalmente que outros

2 0 métodos apropriados para transferir dados SAR também podem

ser utilizados, tal como apreciado prontamente pelos elementos versados na técnica. As imagens SAR coletadas podem ser avaliadas como complexas, ou reais.

Um monitor 64 é conectado ao computador 62 para 25 visualizar as imagens SAR após o processamento. Dispositivos de entrada tais como um teclado 66 e um mouse 68 também são conectados ao computador 62. De acordo com a presente invenção, o computador 62 inclui um processador 68 para processar as imagens SAR.

3 0 Um aspecto para melhorar a resolução de dados SAR

para que se pareçam mais como dados EO é baseado na modificação de uma equação de calor, que é uma equação diferencial parcial linear de segunda ordem. A equação de calor é tal como segue:

&u(x,t) .. .

-Z-- CV U (x, 0, C e

ot (I)

or ~ div(cVu(x,i)) = V · cVu(x, t)

Sob determinadas condição, uma solução fundamental da equação de calor é a função da densidade de Gauss. A equação de calor também pode ser escrita tal como segue:

= í//v(c( vi yi/)Vií(x, y, f)), ceíl

= c(a-, V, Í)V 2Ii(jf, V,f) + Vc(x, y, t) · Vu(x, y, t)

A variável c é conhecida como uma função de

convecção, e t corresponde ao tempo, e x e y formam um número complexo baseado nos dados SAR coletados. De acordo com a presente invenção, a variável c não é um valor constante. Isto permite vantajosamente que a equação de calor seja 10 adaptada aos dados SAR que estão sendo processados, uma vez que c não é um valor constante.

Em conseqüência disto, uma classe de funções pode ser criada para obter os resultados desejados, em que cada função corresponde aos dados específicos SAR que estão sendo 15 processados. Dentro de uma área urbana, os difusores em um conjunto de dados SAR podem ser diferentes dos difusores em um outro conjunto de dados SAR, de modo que os respectivos conjuntos de dados SAR sejam discrepantes.

Uma vez que os conjuntos de dados SAR não são 20 similares, a aplicação de uma mesma função (isto é, um algoritmo de difusão anisotrópica) resulta na melhora da resolução do conjunto de dados SAR mais proximamente compatíveis com a função. Para o outro conjunto de dados SAR que não é proximamente compatível com a função, a sua 25 resolução não será tão boa quanto se uma função mais proximamente compatível fosse utilizada. 0 mesmo podem ser dito sobre conjuntos de dados de áreas rurais. Ao mudar a função de convecção c para melhor compatibilizar com um conjunto de dados SAR particular, então a imagem de valor real correspondente para a visualização terá uma resolução melhor do que se c for um valor constante.

Quando a função de convecção c é um valor constante, os conjuntos de dados SAR discrepantes são tratados igualmente. Um c não-constante permite que o algoritmo de difusão anisotrópica simultaneamente borre e afine um conjunto de dados SAR. Ao ajustar matematicamente a equação de calor 10 através da função de convecção c, os conjuntos de dados SAR discrepantes podem ser melhor processados.

Para ajustar matematicamente a equação de calor, o ruído no conjunto de dados SAR precisa ser determinado. Uma abordagem para determinar o ruído é baseada na coleta de 15 estatística nos valores de gradiente em cada conjunto de dados SAR a ser processado. A estatística pode ser coletada no próprio conjunto de dados SAR real, ou pode ser predeterminada com base nos conjuntos de dados SAR similares que já foram processados. Outras abordagens para determinar o

2 0 ruído incluem um esquema de janelas de Fourier ou uma

decomposição de ondulações, tal como deve ser compreendido prontamente pelos elementos versados na técnica.

A referência é dirigida às FIGURAS 2-4 para ilustrar o cálculo do ruído em um conjunto de dados SAR com 25 base na análise estatística, que é utilizada por sua vez para ajustar matematicamente a função de convecção c dentro do algoritmo de difusão anisotrópica a ser aplicado ao conjunto de dados SAR. Uma distribuição das intensidades de pixels para o conjunto de dados SAR é fornecida na FIGURA 2, ao

3 0 passo que a distribuição dos valores de gradiente para os

pixels é fornecida na FIGURA 3. A distribuição de valores de gradiente é o número de valores de gradiente em um mesmo valor. Por exemplo, o ponto 100 que atinge seu pico em 12xl04 significa que há 120.000 gradientes que estão no mesmo valor.

A coleta de estatística para o ruído pode ser baseada em um desvio padrão dos valores de gradiente. Uma vez que o desvio padrão é determinado, uma constante 5 predeterminada pode ser adicionada a um múltiplo do desvio padrão para obter o limite de ruído k. Uma vez que o limite de ruído k tenha sido determinado para o conjunto de dados SAR, o algoritmo de difusão anisotrópica correspondente é aplicado para suavizar os valores à direita de k enquanto não 10 suaviza os valores à esquerda de k. Ao não suavizar os valores à esquerda de k, as bordas na cena são preservadas.

0 desvio padrão para a distribuição ilustrada de valores de gradiente na FIGURA 3 é indicado pela linha que corresponde à referência 102. 0 limite k é ajustado como duas 15 vezes o desvio padrão mais uma constante. Variações desta abordagem, assim como outras abordagens estatísticas, podem ser utilizadas para selecionar o limite k desejado, tal como deve ser apreciado prontamente pelos elementos versados na técnica.

2 0 O limite k é conhecido como valor de k. Os valores

elevados de gradiente correspondem aos difusores brilhantes, que devem permanecer inalterados. Portanto, o valor de k é ajustado sem borrar os difusores brilhantes. Por outro lado, os valores de gradiente parecem similares devem ser 25 suavizados. Ao ajustar o limite k, as classes diferentes de funções podem ser utilizadas para criar os resultados desejados específicos para os conjuntos de dados SAR que estão sendo processados. Uma vantagem do ajuste da função de convecção c através do limite de ruído k é que o usuário 30 final não tem que fazer uma seleção entre múltiplos algoritmos de difusão anisotrópica para escolher aquele que é mais apropriado para o processamento do conjunto de dados SAR. Ao invés disto, esta seleção é feita de maneira autônoma uma vez que o limite de ruído k tenha sido selecionado.

A distribuição de intensidades de pixels e de valores de gradiente nas FIGURAS 2 e 3 corresponde à imagem mostrada na FIGURA 4, e à imagem dos gradientes mostrada na 5 FIGURA 5. As duas imagens são muito similares. Cada imagem inclui uma escala de dB 110 que representa a quantidade de brilho para visualizar a imagem.

Para a imagem dos gradientes ilustrada na FIGURA 6, um esquema de ponderação convencional, tal como a ponderação 10 de Taylor, por exemplo, foi aplicado durante o processamento do conjunto de dados SAR. Os subtítulos ilustrados indicam que as bordas dos edifícios não estão claras, nem as árvores e as sombras podem ser facilmente determinadas.

Na FIGURA 7, o limite de ruído k para o conjunto de 15 dados SAR foi ajustado com base em uma análise estatística da distribuição dos valores de gradiente, tal como discutido acima. A filtragem adaptada, que também pode ser indicada como um filtro inteligente, combina melhor com os dados SAR que estão sendo processados ao ajustar onde a suavização deve 20 ser executada. Em conseqüência disto, os subtítulos ilustrados indicam que as bordas dos edifícios estão mais claras, e as árvores e as sombras estão mais claras, o que é vantajoso para uma melhor mensuração.

As comparações da mesma cena ao utilizar esquemas 25 filtrando diferentes serão discutidas agora com referência às FIGURAS 8-12. Os vários difusores de interesse são circundados em cada FIGURA. A cena original antes da filtragem é mostrada na FIGURA 8. Na FIGURA 9, a filtragem da cena original é baseada em uma janela de grade de 9. Cada 3 0 grupo de 3x3 pixels é calculado em média, e isto é repetido para todos os pixels no conjunto de dados SAR. Um filtro de Gauss padrão foi aplicado à cena original tal como mostrado na FIGURA 10, e um algoritmo de difusão anisotrópica padrão foi aplicado à cena original tal como mostrado na FIGURA 11.

Para melhor remover o ruído em torno dos difusores na cena original, enquanto os componentes dos difusores ficam intactos, o ruído no conjunto de dados SAR é determinado 5 primeiramente. Com base no limite de ruído k selecionado, o algoritmo de difusão anisotrópica é ajustado de maneira correspondente para prover uma imagem de maior resolução, tal como mostrado na FIGURA 12.

Mesmo que um algoritmo de difusão anisotrópica seja 10 aplicado à cena na FIGURA 11, os difusores ainda têm o que é conhecido como multi-rebote em torno dos mesmos. Com o multirebote, as ondas atingem os difusores e a interface com a terra, e em conseqüência disto rebatem todos em torno dos difusores. 0 multi-rebote se parece com o ruído, mas em 15 algumas situações pode divulgar informações úteis sobre os difusores. Na FIGURA 12, a função de convecção c foi selecionada de modo que o multi-rebote fosse removido.

A vantagem de controlar seletivamente a função de convecção c com base em uma análise estatística da

2 0 distribuição de valores de gradiente para o conjunto de dados

SAR que estão sendo processados provê uma maior resolução. A suavização entre regiões e a preservação das bordas são providas para as imagens corrompidas pelo ruído aditivo. Nos casos onde os conjuntos de dados SAR contêm salpicos, o 25 algoritmo de difusão anisotrópica com a função ajustável de convecção c produz a redução de salpico sensível a borda.

A função de convecção c seletivamente controlada pode ser vantajosamente aplicada em dados complexos crus (isto é, componentes reais e imaginários) e nas imagens

3 0 detectadas (isto é, apenas os componentes reais) ao utilizar

hardware e/ou software para melhorar a fidelidade total do conjunto de dados SAR. Isto também pode ser feito de maneira autônoma com base na seleção do limite de ruído para os conjuntos de dados SAR que estão sendo processados. Cenas do tipo EO de alta resolução podem desse modo ser criadas a partir dos conjuntos de dados SAR. Ao remover simultaneamente o ruído e ao suavizar as áreas de dados similares no conjunto 5 de dados SAR, os dados de alta freqüência são preservados. Conseqüentemente, a textura da informação e as estruturas lineares são preservadas, o que propicia uma avaliação mais exata entre os dados EO e SAR.

Um outro aspecto da invenção refere-se a um 10 algoritmo de difusão anisotrópica complexo. As equações para o algoritmo de difusão anisotrópica como tal discutido acima foram re-derivadas para dados complexos de modo que o processo ficasse agora não-linear. Em termos de dados complexos, os componentes reais e imaginários de um conjunto 15 de dados SAR são processados em um mesmo tempo.

Por outro lado, mesmo que os algoritmos de difusão anisotrópica tivessem sido aplicados aos dados complexos SAR, isto foi feito de uma forma linear. Isto significa que os componentes reais e imaginários dos dados SAR complexos são

2 0 processados separadamente, e os resultados são então combinados conjuntamente.

As equações re-derivadas para o algoritmo de difusão anisotrópica complexo são tal como segue: div(gi\'í)VI} =

ο ' <ix

õx

fçto' Sr I

<3)

Kz

' Dh V í õ<i

1 4»

v t /x J (m

(lí

( Λ -- / 'jdoϊ3νι 2ί^ΚΪ+{^Yüf.1+/*T^V| K1 m -K2 v \ õx Λ àx~ J v Sv &v } \py õxdy ) V <?*'’ Â J ) hítHfHfHÜ, V&-J lôv.l {Sy J ; Z ffh f ί õb I K: * õx‘ Õx HfHSHiH?)) . Uma vez que os componentes reais e imaginários do conjunto de dados SAR estão sendo tratados como um objeto complexo, o algoritmo de difusão anisotrópica complexo pode retirar o conteúdo da cena dos dados extremamente ruidosos, melhorando por sua vez a resolução da imagem visualizada.

Como uma comparação, a referência é dirigida à imagem original mostrada na FIGURA 13. A difusão anisotrópica complexa é aplicada ao conjunto de dados SAR que corresponde à imagem SAR original na FIGURA 13 para formar a imagem 10 difusa mostrada na FIGURA 14. Os limites e as características são visivelmente mais agudos na imagem difusa.

Além do algoritmo de difusão anisotrópica complexo que está sendo aplicado aos conjuntos de dados SAR de imagem simples, ele também pode ser aplicado a processamento interferométrico. Particularmente, o processamento

interferométrico para a medição do abrandamento para a cena urbana é particularmente benéfico ao utilizar o algoritmo de difusão anisotrópica complexo.

O processamento interferométrico de dados SAR será discutido agora em mais detalhes. Com referência ao fluxograma na FIGURA 15, o processamento interferométrico de dados SAR para medições do abrandamento em áreas urbanas será discutido como um exemplo ilustrado. 0 abrandamento é um deslocamento do terreno em que a elevação da superfície da terra está diminuindo em relação ao nível do mar.

As imagens SAR são recebidas em um Bloco do histórico da fase de aparar 120. Para finalidades de discussão, duas imagens SAR estão sendo comparadas. As duas imagens SAR são da mesma cena, mas as imagens estão 10 ligeiramente deslocadas uma da outra, tal como deve ser apreciado prontamente pelos elementos versados na técnica. Se as imagens SAR forem recebidas como dados crus, elas são convertidas de espaço da freqüência em espaço da imagem. 0 espaço da freqüência corresponde à fase e à intensidade, ao 15 passo que o espaço da imagem corresponde à fase e à amplitude. A fase e a amplitude para cada pixel na imagem SAR fornecem os componentes reais e imaginários para os dados complexos SAR a serem processados.

0 Bloco de histórico da fase de aparar 120 assegura 20 em um nível muito elevado que as duas imagens SAR são apropriadas para o processamento interferométrico. A interseção dos respectivos históricos da fase no espaço da freqüência é selecionada entre as duas imagens SAR, e tudo o mais é rejeitado. As duas imagens SAR são registradas no 25 bloco 122. O registro assegura que as características entre as duas imagens SAR fiquem alinhadas. Por exemplo, um canto de um edifício em uma determinada latitude/longitude/altura na primeira imagem SAR é registrado para corresponder a uma mesma latitude/longitude/altura na segunda imagem SAR. Em

3 0 conseqüência disto, os pixels são alinhados acima entre as duas imagens SAR.

0 interferograma é formado no bloco 124. A primeira imagem SAR é multiplicada pelo conjugado complexo da segunda imagem SAR. O resultado é uma diferença de fase entre as duas imagens SAR. O interferograma resultante é relacionado diretamente à altura. No interferograma, a fase para cada pixel é obtida ao tomar o arco tangente de sua parte 5 imaginária dividido por sua parte real. Para a primeira imagem SAR, os dados da fase para cada pixel são determinados. Do mesmo modo, os dados da fase para cada pixel são determinados para a segunda imagem SAR. Conforme será discutido em detalhes abaixo, os dados da fase para cada

pixel variam do interferograma entre menos pi e mais pi. Conseqüentemente, as fases envolvem ao redor.

Em seguida, um filtro de baixa passagem deve ser aplicado normalmente para suavizar o interferograma. No entanto, isto tem uma tendência de borrar as bordas na cena.

Para uma cena rural, o borrão é aceitável, mas para uma cena urbana em que o abrandamento está sendo medido em indicações específicas o borrão não é desejável, uma vez que isto afeta a exatidão do abrandamento medido.

No lugar de um filtro de baixa passagem, um

2 0 algoritmo de difusão anisotrópica complexo tal como discutido

acima é aplicado no bloco 126 e um filtro de choque é aplicado no bloco 128. Com o algoritmo de difusão anisotrópica complexo, as partes reais e imaginárias de cada pixel são processadas como um objeto complexo, isto é, um

processamento não-linear. Por outro lado, o processamento linear envolve o processamento separado das partes reais e imaginárias e então a combinação dos resultados conjuntamente. 0 algoritmo de difusão anisotrópica complexo mitiga localmente o ruído enquanto ao mesmo tempo preserva as

3 0 descontinuidades de cena no interferograma.

0 filtro de choque é utilizado para eliminar o borrão da imagem tal como deve ser compreendido prontamente pelos elementos versados na técnica. Em outras palavras, os limites na cena são realçados ao utilizar morfologia matemática. A equação que corresponde ao filtro de choque é tal como segue:

hiperbólica não-linear, primeira parte da equação corresponde à erosão/dilatação que é determinada por Laplace. A segunda parte da equação é uma magnitude do gradiente. Para ilustrar a aplicação do filtro de choque, uma imagem multi-espectral de vetor próximo do original é mostrada na FIGURA 16, e a 10 aplicação do filtro de choque à imagem é mostrada na FIGURA 17. Os limites são observáveis mais agudos após a aplicação do filtro de choque.

anisotrópica complexo e do realce do limite através do filtro 15 de choque aumenta a exatidão e a qualidade da medição da diferença de fase. Uma vez que o abrandamento é medido ao utilizar objetos persistentes na cena (isto é, edifícios) como pontos de referência, a melhora da qualidade do limite dos objetos persistentes melhora a medição do abrandamento.

2 0 Uma vez que a fase só pode variar entre mais pi e

menos pi, ela é chamada uma fase enrolada. Se não houver nenhum envolvimento de ambigüidade nas fases entre as duas imagens SAR, o abrandamento pode então ser medido no bloco 13 0. No entanto, se existir um envolvimento de ambigüidade, 25 como é tipicamente o caso, então um desenrolamento de fase variacional é executado no bloco 132.

interferograma, que é a diferença de fase entre a primeira e a segunda imagens SAR registradas. A diferença de fase entre

3 0 as duas imagens SAR registradas é relacionada a uma

(4)

A equação de choque é uma equação diferencial

A suavização não-linear através do algoritmo

0 desenrolamento de fase variacional é aplicado ao quantidade física desejada de interesse, tal como a topografia da superfície. A diferença de fase pode ser registrada somente no modulo 2B, e as técnicas interferométricas atuais recuperam principalmente a fase 5 absoluta (a fase desmembrada) da fase registrada (fase enrolada) utilizando valores distintos. O desenrolamento de fase atual pode ser executado por algoritmos de árvore de corte de resíduo e algoritmos dos menores quadrados, por exemplo.

Para executar o desenrolamento de fase, a fase é

determinada a partir do interferograma, que é um objeto complexo com partes reais e imaginária. O arco tangente da parte imaginária em relação à parte real fornece as respectivas fases. A amplitude é rejeitada e a fase permanece.

Uma vez que a fase só pode variar entre mais pi e menos pi, ela é chamada fase enrolada. Na realidade, no entanto, a fase vai de mais infinito a menos infinito. Isto é onde as dificuldades residem no processo interferométrico.

O objetivo é determinar o mapeamento apropriado

para ir do espaço mais/menos pi ao espaço mais/menos infinito. No entanto, as imagens finitas são limitadas pela altura do objeto mais alto na cena. Se o edifício mais alto for de 800 pés, a diferença é baseada então no nível na terra 25 e 800 pés. Em teoria, mais/menos infinito é matematicamente correto, mas na realidade a variação fica entre zero e a altura do objeto mais alto na cena. O desenrolamento de fase variacional lida com não-linearidades e descontinuidades nos dados. Existe uma ambigüidade nos envolvimentos de fase nos 30 limites de mais/menos pi antes de tirar a diferença de fase entre as duas imagens SAR. O ponto em que a fase é enrolada é conhecido como linhas limítrofes.

O algoritmo de desenrolamento de fase variacional é bidimensional. As técnicas de desenrolamento de fases unidimensionais podem ser re-derivadas para duas dimensões ao utilizar os requisitos que se aplicam especificamente ao problema de abrandamento. Outras técnicas de desenrolamento 5 de fase bidimensional que são disponíveis também podem ser adaptadas.

Na análise de imagem, a segmentação é a divisão de uma imagem digital em regiões múltiplas (conjuntos de pixels) de acordo com algum critério. 0 objetivo da segmentação é 10 tipicamente encontrar objetos de interesse. Algumas técnicas comuns para a segmentação incluem a fixação de limites, o crescimento de região e a etiquetação de componentes de conexão. Os contornos ativos também constituem um método comum.

O algoritmo de desenrolamento de fase variacional é

baseado na função de Mumford-Shah ou na função de custo, tal como indicado abaixo:

E{f.C) = fi\(f-gfdA+a { \VjfdA±Y§(k (5)

U OvC C

Ω

mc

Y^ds

ú

A equação determina que f e C irão fornecer a fase desenrolada. 0 primeiro termo é a aproximação da suavização

2 0 de pedaço f a g (a imagem) com descontinuidades ao longo de C. Esta parte da equação pode ser imaginada como um termo de fidelidade dos dados que mede a qualidade de f. 0 segundo termo da equação é o termo da lisura. Isto pode ser visto como o modelo prévio para f resultar em C. 0 terceiro corresponde ao comprimento de C. Normalmente, há uma penalidade para o comprimento de arco excessivo. A função 5 original de Mumford-Shah ou a função de custo utilizou a medição de Hausdorff para conjuntos de descontinuidades mais gerais. De acordo com a presente invenção, C é restringido para ser uma curva suave a fim de ser substituída pelo comprimento do arco.

Para minimizar a função de Mumford-Shah ou a função

de custo, uma nova função de custo é desenvolvida para melhor solucionar a descontinuidade dos dados entre as linhas limítrofes. A nova função de custo é tal como segue:

E(f<C) = pl(fx-gJ-dÁ*fij(fy-gydA*m I \Vf\*d + yÇds (6)

0 Ü OVC C

O primeiro termo expressa os gradientes entre a 15 fase enrolada e desenrolada. O segundo termo expressa o conhecimento prévio da cena a ser processada. O terceiro termo impõe limites no comprimento limítrofe máximo do interferograma não processado. A maioria dos algoritmos de desenrolamento de fase trabalha com dados suaves. Após ter 20 determinado a solução mínima para a função de custo E(f, C), uma conversão é feita em uma equação diferencial parcial (PDE). A PDE é então solucionada.

0 algoritmo de desenrolamento de fase variacional tira vantagem do fato que os dados são pré-processados com o 25 algoritmo de difusão anisotrópica complexo. 0 algoritmo de difusão anisotrópica complexo é projetado para não suavizar as descontinuidades. Conseqüentemente, a abordagem variacional para o desenrolamento de fase tira vantagem do fato que os dados serão ainda descontínuos.

3 0 0 desenrolamento de fase variacional será discutido agora com referência aos gráficos mostrados nas FIGURAS 18- 27. Uma cena original de dois edifícios 150, 152 e a terra correspondente 160 adjacente aos edifícios é fornecida nas FIGURAS 18A e 18B. A FIGURA 18B é um gráfico tridimensional 5 da imagem SAR, e a FIGURA 18A é uma vista de cima para baixo da mesma imagem SAR. Na cena original, a terra 160 é um morro que é quase tão alto quanto um dos edifícios 152.

Um interferograma enrolado da cena original é fornecido nas FIGURAS 19A e 19B. Uma vez que o interferograma 10 é um objeto complexo, a sua fase é determinada de modo que o interferograma possa ser visualizado. Conseqüentemente, a linha central do gráfico na FIGURA 19B é em radianos. Para cada pixel de x e de y há um valor da fase, que é enrolada. Isto significa que a faixa fica sempre entre mais/menos pi. 15 As bordas da faixa de mais/menos pi são as linhas limítrofes 170 e 172.

Na FIGURA 19B, o centro do morro 160 caiu. A diferença de fase não inclui nenhum ruído. 0 ruído é adicionado artificialmente à cena para simular uma coleta

2 0 real, tal como mostrado nas FIGURAS 2OA e 20B. Nas FIGURAS 21A e 21B, o ruído é mitigado com a aplicação de um filtro suavizante. Uma vez que os algoritmos de suavização normal têm uma tendência de suavizar as descontinuidades, as linhas limítrofes 170, 172 foram suavizadas. Em conseqüência disto, 25 a resolução foi reduzida. Um interferograma de difusor anisotrópico complexo é provido nas FIGURAS 22A e 22B. As linhas limítrofes 17 0, 172 são mais agudas, e as bordas do edifício são preservadas.

Um outro conjunto de exemplos será discutido agora com referência às FIGURAS 23-27. Uma vista de cima para baixo de quatro edifícios 180, 182, 184 e 186 é mostrada na FIGURA 23. Cada edifício fica em uma altura diferente tal como indicado por uma sombra diferente. Um interferograma sem ruido é mostrado na FIGURA 24. Por causa do enrolamento de fase entre mais/menos pi, os quatro edifícios 180, 182, 184 e 186 parecem estar na mesma altura que aquela indicada pela mesma sombra. Quando o ruído é adicionado ao interferograma, 5 os edifícios 180, 182, 184 e 186 ficam escondidos pelo ruído, tal como mostrado na FIGURA 25. Um interferograma de ruído anisotropicamente difuso complexo de acordo com a invenção é mostrado na FIGURA 26. Os edifícios 180, 182, 184 e 186 são extraídos do ruído após a aplicação do algoritmo de difusão 10 anisotrópica complexo. Uma vista aproximada do edifício 18 0 é fornecida na FIGURA 27 para ilustrar como as bordas limite são mantidas.

Com referência novamente ao fluxograma na FIGURA 15, a geometria de cada SAR que fornece uma respectiva imagem 15 SAR é estimada no bloco 134. É feita uma determinação a respeito de onde cada SAR foi localizado no momento em que a imagem correspondente foi tomada. Se o primeiro SAR estiver apontando a uma determinada latitude/longitude/altura, então haverá uma alta confiança nos valores de pixel contanto que a 20 latitude/longitude/altura seja para esse pixel.

A fase desenrolada em radianos é convertida em altura no bloco 136. Quando uma medição de fase desenrolada é obtida ela é em radianos. Uma conversão é feita então dos radianos na altura. Uma equação de forma fechada toma o valor 25 do radiano em altura tal como compreendido prontamente pelos elementos versados na técnica. A altura fornece a medição necessária para determinar o abrandamento entre as duas imagens SAR. Para finalidades de ilustração, uma imagem original é mostrada na FIGURA 28, a imagem original alisada 30 com um filtro de Gauss é mostrada na FIGURA 29, e a imagem original filtrada com um algoritmo de difusão anisotrópica complexo após vinte iterações é mostrada na FIGURA 30. Entre as duas imagens, as bordas limite são visivelmente mais encrespadas na FIGURA 30.

A descontinuidade é mantida desse modo pelo algoritmo de difusão anisotrópica por causa das propriedades do algoritmo. 0 algoritmo opera basicamente nos gradientes de 5 modo que fica sabendo, com base estritamente no gradiente da imagem, se deve ou não suavizar. Se o gradiente estiver abaixo do limite do ruído que é ajustado acima da parte dianteira, o algoritmo deve seguir adiante e suavizar o gradiente. Se o gradiente estiver acima do limite, o 10 algoritmo não irá suavizar o gradiente a fim de manter ou preservar uma borda de uma estrutura correspondente. Conforme observado acima, o limite é ajustado preferivelmente com base no conhecimento da cena.

O desenrolamento de fase variacional é derivado da 15 matemática de erro denominada cálculo variacional. Os princípios do cálculo variacional são utilizados para compatibilizar com um algoritmo de desenrolamento de fase que lida com descontinuidades. Os algoritmos de desenrolamento atuais têm uma tendência de suavizar os dados. Por outro

2 0 lado, o algoritmo de desenrolamento de fase variacional é

baseado em uma função de custo para controlar a suavização. Em vez de prover suavização global com base nas propriedades dos dados, o algoritmo de desenrolamento de fase variacional deixa bordas intactas e suaviza seletivamente a área adjacente às bordas.

Um outro aspecto da invenção é a aplicação do algoritmo de difusão anisotrópica ao comprimir e ao descomprimir as imagens SAR. Os conjuntos de dados SAR podem ser relativamente grandes, e quando um conjunto de dados SAR

3 0 é comprimido a resolução é geralmente reduzida durante o

processo. Este é um resultado de esquemas de compressão de pré-processamento de perda. Os algoritmos do préprocessamento de perda geralmente degradam os difusores em uma cena. Além disso, o volume de dados pode sufocar as capacidades de processamento atuais.

Os algoritmos de pré-processamento mais comuns agem como filtros de baixa passagem. Os seguintes esquemas de 5 compressão tentam agrupar os dados de uma maneira que encontre similaridades entre todos os dados: a análise de componente independente (ICA), a transformação de ondulações (filtros Gabor) e esquemas de exploração do paralelismo. Devido à faixa dinâmica dos dados SAR, é difícil limitar os 10 dados de uma maneira tal que os dados possam ser bem agrupados.

Um método implementado por computador para comprimir imagens SAR compreende a recepção de uma imagem SAR a ser comprimida, a aplicação de um algoritmo de difusão 15 anisotrópica à imagem SAR, e a compressão da imagem SAR após ter aplicado o algoritmo de difusão anisotrópica à mesma. Uma vantagem da aplicação do algoritmo de difusão anisotrópica é que o tamanho da imagem SAR após a compressão resulta em um arquivo de menor tamanho, independentemente do esquema de

2 0 compressão utilizado. Conseqüentemente, o armazenamento e a transmissão da imagem SAR comprimida ocupam menos espaço e largura de faixa.

Para finalidades da comparação, a métrica para a compressão é baseada na cena original mostrada na FIGURA 31. 25 A cena original tem um tamanho de arquivo não-Tiff comprimido de 691 kB. A compressão JPEG da cena original reduz o tamanho do arquivo para 62 kB, ao passo que a compressão Winzip da cena original reduz o tamanho do arquivo para 41 kB.

A filtragem da cena original com um filtro de Gauss é mostrada na FIGURA 32. A compressão JPEG da imagem original com filtragem de Gauss reduz o tamanho do arquivo para 55 kB, ao passo que a compressão Winzip da imagem original com filtragem de Gauss reduz o tamanho do arquivo para 33 kB. De acordo com a presente invenção, a filtragem da cena original com um algoritmo de difusão anisotrópica é mostrada na FIGURA 33. 0 tamanho da cena filtrada difusa anisotrópica é ainda do mesmo tamanho que a cena original sem 5 filtragem e com filtragem de Gauss. A compressão JPEG da cena anisotropicamente difusa reduz o tamanho do arquivo para 44 kB, ao passo que a compressão Winzip da cena anisotropicamente difusa reduz o tamanho do arquivo para 23 kB.

A TABELA 1 fornece uma comparação de lado a lado

entre as imagens diferentes. Quando a difusão anisotrópica foi aplicada a algumas das imagens SAR, uma compressão maior pode ser conseguida do que quando o algoritmo de difusão anisotrópica não foi aplicado. A imagem filtrada com difusão 15 anisotrópica tem uma compressão JPEG de 16:1 e uma relação de compressão Winzip de 30:1.

TABELA 1

Tabela de Relação de Compressão Tipo de imagem Tiff não Compressão Winzip comprimida JPEG Imagem original 1 11:1 16 : 1 Com filtragem Gauss 1 13 : 1 21: 1 Filtrada 1 16 : 1 30 : 1 anisotrópica Em termos de descompressão, a descompressão é

executada com base na difusão anisotrópica. Mais 20 particularmente, um método implementado por computador para descomprimir imagens SAR compreende a recepção de uma imagem SAR a ser descomprimida, a execução de uma compressão de faixa dinâmica na imagem SAR, a quantificação da imagem comprimida SAR, e a descompressão da imagem SAR comprimida e 25 quantificada ao aplicar um algoritmo de difusão anisotrópica à mesma. A quantificação pode ser em unidade8, por exemplo. A compressão da faixa dinâmica é um processo não-linear.

Agora é feita referência às FIGURAS 34-36 para ilustrar a descompressão. A imagem original a ser descomprimida é mostrada na FIGURA 34, e tem um tamanho de armazenamento de imagem de 65,68 MB. A compressão não-linear da faixa dinâmica é aplicada e os resultados são 5 quantificados em unidade8, tal como mostrado na FIGURA 35. 0 tamanho de armazenamento da imagem é agora de 4,11 MB. A árvore e as sombras não são bem definidas na imagem quant i f i c ada.

A descompressão da imagem quantificada com um algoritmo de difusão anisotrópica é mostrada na FIGURA 36. As árvores e as sombras são agora melhor definidas. Um filtro de choque pode até mesmo ser aplicado para realçar ainda mais a imagem SAR visualizada. Ao comprimir dinamicamente os dados SAR e então quantificar esses dados, a quantidade de dados que são requeridos durante a transmissão é reduzida significativamente. Em média, esses conjuntos de dados devem requerer 4,11/65,7 = 6,25% dos dados da cena para a transmissão. Quanto maior a faixa dinâmica, maior a relação de compressão. Para difusores muito brilhantes em uma cena, fica comprimida ainda mais. Mesmo se um usuário for provido com uma imagem comprimida/descomprimida de perda, a aplicação das abordagens de compressão/descompressão acima irá melhorar realmente a qualidade da imagem original para a visualização.

A extração/registro da elevação ao utilizar a 25 difusão anisotrópica tal como discutido acima para a formação de imagens com ruído e imagens SAR será discutida agora com referência às FIGURAS 37-44. Os dados de ruído efetuam a exatidão da correlação, do registro (o mesmo ou sensor transversal) e da extração da elevação. Atualmente, os 30 filtros de baixa passagem são utilizados para dados de ruído. A apodização é utilizada para dados SAR.

Um método implementado por computador para registrar imagens SAR compreende a seleção da primeira e segunda imagens SAR a serem registradas, o processamento individual da primeira e segunda imagens SAR selecionadas com um algoritmo de difusão anisotrópica, e o registro da primeira e segunda imagens SAR após o processamento. Um 5 filtro de choque é aplicado preferivelmente à primeira e segunda imagens SAR processadas respectivas antes do registro. Os dados da elevação podem então ser extraídos com base nas imagens SAR registradas.

Para finalidades de ilustração, duas imagens SAR não registradas são mostradas nas FIGURAS 37 e 38. Para obter a métrica em relação às vantagens de utilização de um algoritmo de difusão anisotrópica durante o registro, a correlação é entre a imagem 1 na FIGURA 37, que é a referência, e a imagem 2 na FIGURA 38, que é a imagem secundária. Um mapa de coeficiente de correlação para as duas imagens é determinado. O valor máximo do mapa de correlação é obtido, isto é, o pico. A posição do pico no mapa da correlação indica o deslocamento entre os dados, isto é, o deslocamento registrado. Toda a imagem 1 é correlacionada com toda a imagem 2. O valor máximo da correlação é obtido após o registro. A imagem é então filtrada ao utilizar o algoritmo de difusão anisotrópica.

0 pico da correlação é de 0,9081 na superfície da correlação não filtrada tal como mostrado na FIGURA 39. Após 25 a filtragem, o pico da correlação é de 0,9674 tal como mostrado na FIGURA 40. Depois que o registro é aplicado, as imagens das FIGURAS 3 7 e 38 não se movem, tal como mostrado nas FIGURAS 41 e 42. Estas imagens não são filtradas. O pico da correlação para o registro é de 0,9081 e a pós-correlação 30 de toda a imagem depois do registro é de 0,3634. Após a filtragem, as imagens correspondentes são mostradas nas FIGURAS 43-44. O pico da correlação para o registro é de

0,9674 e a pós-correlação de toda a imagem depois do registro é de 0,8094. O registro é incrementado em 9,4%, e a correlação é incrementada após o registro em 45%.

A extração de vetor/estrada ao utilizar processamento de difusão anisotrópica não-linear e filtros de 5 choque para a classificação material serão discutidos agora com referência às FIGURAS 45-48. Existe uma demanda para a extração da estrada de cenas refletivas simples. Atualmente, kernels de suavização são aplicados aos dados. Filtros morfológicos também são aplicados (dilatação/erosão) aos 10 dados. A extração de vetor/estrada pode então ser provida ao utilizar um filtro de Gauss/de cruzamento zero.

De acordo com a presente invenção, um método implementado por computador para a extração do vetor em imagens SAR compreende a seleção de uma imagem SAR para a 15 extração do vetor, o processamento da imagem SAR selecionada com um algoritmo de difusão anisotrópica, e a extração dos dados do vetor com base na imagem SAR processada. Um filtro de choque pode ser aplicado à imagem SAR processada antes da extração. Os dados da imagem da estrada podem então ser 20 extraídos com base nos dados de vetor extraídos. Os dados são delineados desse modo ao utilizar um esquema coerente de algoritmo de difusão anisotrópica. Este esquema fornece uma exatidão maior para a extração da estrada. Uma cena original é mostrada na FIGURA 45. Depois que a difusão anisotrópica é 25 aplicada, o alvo é melhor definido, tal como mostrado na FIGURA 46. Em geral, os alvos são melhor delineados para superfícies refletivas simples. Isto se presta bem para a segmentação.

O algoritmo de difusão anisotrópica também ser

3 0 aplicado a outras imagens depois que elas já foram filtradas por outros tipos de filtro. Estes tipos de filtro incluem um filtro Prewitt e um filtro Roberts, por exemplo. Uma imagem original que foi filtrada por um filtro Prewitt é mostrada na FIGURA 47. A FIGURA 48 mostra a mesma imagem após a aplicação do algoritmo de difusão anisotrópica. 0 alvo ilustrado, bem como a cerca da ligação de cadeia são melhor definidos.

Claims (8)

1. MÉTODO IMPLEMENTADO EM COMPUTADOR PARA REGISTRAR IMAGENS DE RADAR DE ABERTURA SINTÉTICA (SAR), caracterizado pelo fato de compreender: a seleção de primeira e segunda imagens de SAR a serem registradas; o processamento individual da primeira e segunda imagens de SAR selecionadas com um algoritmo de difusão anisotrópico, em que o processamento para cada imagem de SAR selecionada compreende: a determinação de ruído na imagem de SAR, a seleção de um limite de ruído para a imagem de SAR com base no ruído determinado, e o ajuste matemático do algoritmo de difusão anisotrópico com base no limite de ruído selecionado; e o registro da primeira e segunda imagens de SAR após o processamento.

2. MÉTODO IMPLEMENTADO EM COMPUTADOR, de acordo com a reivindicação 1, caracterizado pelo fato de compreender adicionalmente a aplicação de um filtro de choque às respectivas primeira e segunda imagens de SAR processadas antes do registro.

3. MÉTODO IMPLEMENTADO EM COMPUTADOR, de acordo com a reivindicação 1, caracterizado pelo fato de compreender adicionalmente a extração dos dados da elevação com base nas imagens de SAR registradas.

4. MÉTODO IMPLEMENTADO EM COMPUTADOR, de acordo com a reivindicação 1, caracterizado pelo fato de que o algoritmo de difusão anisotrópico é baseado em uma equação de calor que compreende um termo não constante; e em que o ajuste matemático do algoritmo de difusão anisotrópico compreende o ajuste do termo não constante.

5. SISTEMA COMPUTADORIZADO PARA REGISTRAR IMAGENS DE RADAR DE ABERTURA SINTÉTICA (SAR), caracterizado pelo fato de compreender: um banco de dados para armazenar as imagens de SAR a serem registradas; e um processador para registrar as imagens de SAR do dito base de dados, em que o registro compreende: a seleção da primeira e segunda imagens de SAR a serem registradas, o processamento individual da primeira e segunda imagens de SAR selecionadas com um algoritmo de difusão anisotrópico, em que o processamento para cada imagem de SAR selecionada compreende: a determinação de ruído na imagem de SAR, a seleção de um limite de ruído para a imagem de SAR com base no ruído determinado, e o ajuste matemático do algoritmo de difusão anisotrópica com base no limite de ruído selecionado, e o registro da primeira e segunda imagens de SAR após o processamento.

6. SISTEMA COMPUTADORIZADO, de acordo com a reivindicação 5, caracterizado pelo fato de que o dito processador é configurado para aplicar um filtro de choque às respectivas primeira e segunda imagens de SAR processadas antes do registro.

7. SISTEMA COMPUTADORIZADO, de acordo com a reivindicação 5, caracterizado pelo fato de que o dito processador ê configurado para extrair os dados da elevação com base nas imagens de SAR registradas.

8. SISTEMA COMPUTADORIZADO, de acordo com a reivindicação 5, caracterizado pelo fato de que o algoritmo de difusão anisotrópico é baseado em uma equação de calor que compreende um termo não constante; e em que o ajuste matemático do algoritmo de difusão anisotrópico compreende o ajuste do termo não constante.