$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
O conjunto de dados inclui um total de 1.190 fatias de tomografia computorizada de 110 casos, categorizadas em três classes: normais (55 casos), benignas (15 casos) e malignas (40 casos). Cada caso compreende múltiplas fatias de tomografia computorizada (variando de aproximadamente 80 a 200 fatias por caixa), oferecendo diversas vistas axiais da região torácica. As imagens de tomografia computorizada foram obtidas no DICOM usando o scanner SOMATOM com parâmetros padrão de imagem: tensão do tubo = 120 kV, espessura do corte = 1 mm, largura da janela = 350–1.200 Unidades de Hounsfield (HU) e valores centrais da janela = 50–600 HU, durante uma suspensão total de inspiração.
Todas as imagens foram totalmente desidentificadas antes da análise para remover qualquer informação pessoal identificável (PII). O conjunto de dados foi eticamente aprovado pelos conselhos institucionais de revisão dos centros médicos participantes, com consentimento por escrito renunciado pelo conselho de supervisão. Os casos incluíam indivíduos de origens variadas, como funcionários públicos, agricultores e moradores de várias províncias iraquianas (incluindo Bagdá, Wasit, Diyala, Salahuddin e Babilônia) com diversidade em gênero, idade, nível educacional e status de vida.
Como o conjunto de dados está disponível publicamente e foi desidentificado, não foi necessária aprovação ética adicional para seu uso neste estudo. O conjunto de dados segue os termos e condições especificados por seus colaboradores originais e pela plataforma de hospedagem.
A metodologia desta pesquisa utiliza um processo em múltiplas etapas projetado para criar um modelo de predição com a ajuda do conjunto de dados de câncer depulmão IQ-OTH/NCCD 21. Essa metodologia cria um terreno forte para dispositivos de saúde focados no consumidor, onde é necessária menos latência, o que pode proporcionar maior precisão. A Figura 2 mostra o mecanismo de funcionamento do modelo proposto.

Figura 2: Diagrama de fluxo de trabalho do modelo proposto representando pré-processamento, aumento, classificação do Transformador de Visão e etapas de avaliação. Por favor, clique aqui para ver uma versão maior desta figura.
1. Descrição do conjunto de dados
O conjunto de dados de câncer de pulmão IQ-OTH/NCCD foi coletado no Hospital Universitário de Oncologia do Iraque/Centro Nacional para Doenças Cancerígenas durante o outono de 2019. A parte de treinamento do conjunto de dados utilizado continha 1.097 imagens cuja descrição é dada na Tabela 2, que ilustra a contagem de exemplos de diferentes classes.
| Classe | Número de imagens |
| Normal | 416 |
| Benigno | 120 |
| Maligno | 561 |
Tabela 2: Amostras de imagens de tomografia computadorizada normais, benignas e malignas do conjunto de dados.
O conjunto de dados de câncer de pulmão IQ-OTH/NCCD foi escolhido por sua representação abrangente das tomografias normais, benignas e malignas. Embora outros conjuntos de dados, como LIDC-IDRI e NSCLC-Radiomics, estejam disponíveis, eles focam principalmente na detecção de nódulos ou carecem de amostras suficientes de casos benignos. O conjunto de dados IQ-OTH/NCCD é particularmente adequado para nosso estudo, pois permite que o Transformador Visual aprenda características discriminativas entre as três classes, possibilitando uma classificação robusta de padrões sutis em imagens de tomografia computorizada pulmonar.
2. Pré-processamento de dados
O pré-processamento é uma etapa importante para melhorar o desempenho do modelo por meio da consistência e dos ajustes apropriados nos dados que serão processados pela rede neural. Para garantir consistência no tamanho de entrada da rede neural, escalamos todas as imagens para a mesma dimensão de 256 x 256 pixels, e normalizamos os dados para um valor em pixel entre 0 e 1, na esperança de melhorar o treinamento e a convergência do modelo. A Tabela 3 contém os valores da técnica de aumento.
| Técnica | Valor |
| Normalização | - |
| Redimensionamento | image_size x image_size |
| RandomRotation | fator=0,02 |
| RandomZoom | height_factor=0,2, width_factor=0,2 |
Tabela 3: Técnicas de aumento aplicadas com intervalos de parâmetros.
Para aumentar a robustez do modelo contra o sobreajuste e melhorar sua capacidade de generalização, técnicas de aumento de dados como rotações aleatórias e zooms são aplicadas, simulando vários ângulos e tamanhos das manifestações do câncer de pulmão conforme podem aparecer em diferentes pacientes. Rotações aleatórias com ângulos amostrados uniformemente a partir de 0,02 radianos e transformações aleatórias de zoom com fatores variados de altura e largura foram aplicadas neste estudo. As imagens após o aumento são exibidas na Figura 3.

Figura 3: Imagens de tomografia após aumento demonstrando rotação, zoom e normalização para melhorar a generalização do modelo. Por favor, clique aqui para ver uma versão ampliada desta figura.
Essas técnicas de pré-processamento são necessárias, e a ampliação tem sido usada em todas as classes para garantir a robustez do modelo.
3. Arquitetura do modelo
Adaptando o inovador framework Vision Transformer (ViT) para lidar eficazmente com dados de imagens complexos, a arquitetura do modelo tem como objetivo categorizar as tomografias em três categorias: normais, benignas e malignas. Com esse método, são processadas imagens de entrada de 256 x 256 pixels. Para garantir consistência e promover um aprendizado de modelos mais eficiente, cada imagem passa por várias operações de pré-processamento, como redimensionamento e normalização. Para melhorar a resiliência do modelo a mudanças na escala e orientação da imagem, o projeto começa com uma camada de aumento de dados que utiliza métodos como normalização, redimensionamento, rotações aleatórias de 0,02 radianos e zooms aleatórios de até 20%. Após a ampliação, o modelo retira 16 x 16 patches de pixel de cada imagem, assim cada imagem tem 256 patches.
O Algoritmo 1 define o fluxo de trabalho do modelo proposto, como ele é construído e o ajuste fino realizado em relação ao modelo.
Algoritmo 1: Classificação do câncer de pulmão usando Transformadores de Visão
Entrada: Conjunto de imagens CT I do conjunto de dados IQ-OTH/NCCD
Resultado: Resultados de classificação C em categorias: Normal, Benigno, Maligno
Pré-processamento:
for each image i in I do
i <- Resize(i, 256 x 256) // Normalize and resize images
i <- Normalize(i)
i <- DataAugmentation(i) // Apply random rotations and zooms
end for
Configuração do modelo:
Initialize Vision Transformer (ViT) Model
Set number of patches P = 16 x 16
Set number of heads H = 6 in multi-head attention
Formação:
for epoch = 1 to MaxEpochs do
for each batch b in I do
Patches <- ExtractPatches(b, P)
EncodedPatches <- PatchEncoder(Patches)
AttentionWeights <- MultiHeadAttention(EncodedPatches, H)
ClassLogits <- TransformerEncoder(AttentionWeights)
Loss <- ComputeLoss(ClassLogits, TrueLabels)
UpdateModelWeights(Loss)
end for
if EarlyStoppingCriteriaMet (Val_Loss No Improvement Patience = 5 Epochs) then
break
end if
end for
Validação:
Split data into TrainingSet (80%) and ValidationSet (20%)
ComputeValidationMetrics(ValidationSet)
Avaliação:
C <- Classify(I)
ComputePerformanceMetrics(C)
Return C
Para preservar as relações espaciais entre patches, esses patches são achatados (equação 1) em vetores de 768 dimensões (porque cada patch tem 16 x 16 x 3 = 768) e então enviados por uma camada de codificação de patch, que projeta cada vetor para um espaço de 64 dimensões integrando embeddings posicionais. O núcleo da arquitetura é composto por oito camadas transformadoras, cada uma com um mecanismo de atenção multi-cabeças com seis cabeças, permitindo que o modelo foque em diferentes partes da imagem simultaneamente e capture uma ampla gama de características. Ele é representado usando as equações 2, 3 e 4.

Onde μ é a média e σ2 é a variância da entrada x, e ε é uma pequena constante para evitar divisão por zero.

Onde Q é a matriz de consulta, K é a matriz chave, V é a matriz de valores e dk é a dimensão dos vetores chave.


Onde WiQ, WiK e WiV são as matrizes de pesos aprendidas para consultas, chaves e valores, respectivamente, e WO é a matriz de pesos de saída.
O diagrama em blocos do modelo proposto foi mostrado na Figura 4.

Figura 4: Diagrama em blocos do modelo Vision Transformer com cabeça MLP, detalhando as principais camadas de processamento e arquitetura. Por favor, clique aqui para ver uma versão ampliada desta figura.
Cada camada transformadora inclui uma percepção multi-camada (MLP) com unidades ocultas inicialmente escaladas para 128 e depois para 64, expandindo e comprimindo efetivamente o fluxo de informações para capturar dependências complexas.
O dropout é aplicado estrategicamente a uma taxa de 0,1 nos mecanismos de atenção e MLPs para evitar o sobreajuste. A saída do codificador do transformador é processada por uma cabeça MLP composta por duas camadas densas com 2.048 e 1.024 unidades, respectivamente, empregando Unidades Lineares de Erro Gaussiano (GELU) para ativação, que tem demonstrado bom desempenho em aplicações de aprendizado profundo. Isso é alcançado usando a equação 5.

Onde W1 e W2 são matrizes de peso, b1 e b2 são vieses, e GELU é a função de ativação utilizada.
Uma queda de 0,1 foi usada nas camadas transformadoras para evitar sobreajuste sem perder informações importantes de características. A função de ativação GELU foi utilizada por suas características suaves e não lineares, que promovem o fluxo de gradiente e a convergência em modelos baseados em transformadores. O abandono é regularizado usando a equação 6.
Onde p é a taxa de desistência (por exemplo, 0,1 ou 0,5), e a camada de desistência define aleatoriamente uma fração p das unidades de entrada a zero durante o treinamento.
Uma taxa de dropout de 0,5 nessas camadas ajuda ainda mais a mitigar o sobreajuste. A camada final do modelo é uma camada logits (equação 7) que gera três logits de classe correspondentes às categorias normal, benigna e maligna, usando uma função de perda de entropia cruzada categórica esparsa (equação 8) apropriada para esse cenário de classificação multiclasse.


Onde zi é o vetor logits para a classe i, SoftMax(z i)) representa as probabilidades previstas, e yi é o verdadeiro rótulo da classe.
O modelo é compilado com o otimizador AdamW (equações 9, 10, 11, 12 e 13), uma extensão do otimizador Adam que lida com decaimento de peso de forma mais eficaz, definido com uma taxa de aprendizado de 0,001 e decaimento de peso de 0,0001, otimizando tanto a velocidade de aprendizado quanto a estabilidade do modelo.





Onde mt e vt são o primeiro e o segundo momento dos gradientes,
e
são os momentos corrigidos por viés, η é a taxa de aprendizado, e ε é uma pequena constante para evitar divisão por zero.
O modelo, quando treinado, utiliza um lote de 32 para aproveitar a aceleração da GPU e facilitar cálculos mais rápidos e é configurado para treinar por 100 épocas.
Mas o treinamento possui um mecanismo de parada precoce na perda de validação para limitar o treinamento quando o modelo para de melhorar, economizando recursos computacionais e evitando sobretreinamento de 5 épocas consecutivas. O desempenho do modelo é acompanhado com métricas como precisão categórica esparsa e precisão top 2, que informam sobre a capacidade de classificação do modelo nas categorias previstas mais prováveis. Callbacks para treinamento de modelos incluem checkpoints, que salvam o modelo de maior desempenho de acordo com a precisão da validação, garantindo que a versão mais bem-sucedida do modelo seja mantida após o término do processo de treinamento. Essa arquitetura robusta e bem planejada aproveita as capacidades dos transformers para processar dados de imagem médica com uma abordagem altamente avançada para utilizar métodos contemporâneos de IA em grandes aplicações em diagnósticos em saúde.
4. Treinamento e validação
O modelo foi treinado em um ambiente Kaggle usando uma GPU NVIDIA P100, e uma mini-descida de gradiente em lote com tamanho de lote de 32 foi usada durante o processo de treinamento. O otimizador usado durante o treinamento foi o AdamW, com taxa de aprendizado de 0,001 e decaimento de peso de 0,0001; Isso era um bom equilíbrio entre uma convergência rápida e alguma regularização. O modelo foi treinado por 100 épocas, porém, foi utilizado o parar precocemente na perda de validação com paciência de 5 épocas. Simplificando, se o treinamento não melhorasse a perda de validação por 5 épocas contínuas, o treinamento pararia para por sobreajuste e eficiência computacional. Na maioria dos casos, esse modelo restaurava pesos da época com menor perda de validação, indicando que estava funcionando de forma ótima e não precisava rodar as 100 épocas completas. Um total de cerca de 32 minutos para treinar o modelo.
Durante o treinamento, as métricas incluíram precisão categórica esparsa e precisão top-2, monitoradas tanto em conjuntos de treinamento quanto de validação. Essas métricas fornecem insights sobre com que frequência um modelo prevê a classe correta e se a classe correta está entre as duas principais previsões, o que em aplicações médicas é importante porque classificações erradas de alta confiança podem ter consequências significativas. As curvas de aprendizado para perda e precisão são visualizadas na Figura 5.
.
Figura 5: Curvas de precisão e perdas de treinamento e validação ao longo de 50 épocas mostrando aprendizado estável e sobreajuste mínimo. Por favor, clique aqui para ver uma versão maior desta figura.
Nenhuma validação cruzada foi realizada no presente estudo. Essa configuração permitiu experimentação eficiente com a arquitetura do modelo, incluindo ajustes nas camadas dos transformadores e nos tamanhos das cabeças MLP, garantindo que o modelo se generalizasse bem em dados de validação não vistos.
5. Análise estatística
O desempenho do modelo Vision Transformer foi analisado estatisticamente com base no conjunto de dados de câncer de pulmão IQ-OTH/NCCD. Precisão, recordação, pontuação F1 e precisão foram calculadas usando as equações 14, 15, 16 e 17, respectivamente. Essas quatro são consideradas as medidas convencionais para tarefas de classificação, pois podem revelar informações sobre o desempenho do modelo na classificação e classificação por classe para cada classe.




A recordação é uma medida da capacidade de um modelo de identificar todas as instâncias relevantes em uma classe, enquanto a precisão é a proporção de identificações positivas que estavam realmente corretas. A pontuação F1 equilibra recordação e precisão quando ambas são importantes.
As medidas de desempenho usadas para esta tarefa foram a equação 18 do Coeficiente de Correlação de Matthews (MCC) e a equação de Kappa 19 de Cohen.


Onde P0 é o acordo observado, enquanto Pe é o acordo esperado.
O MCC é uma medida abrangente do desempenho na classificação, levando em conta tanto verdadeiros positivos quanto negativos, falsos positivos e falsos negativos. Seu valor pode estar entre -1 e 1, sendo 1 a previsão perfeita, 0 a previsão aleatória e -1 a total discordância entre os rótulos previsto e verdadeiro. O MCC foi valioso para comparações entre diferentes desempenhos de modelos quando aplicado a conjuntos de dados desequilibrados e fornecendo uma medida equilibrada independentemente dos tamanhos de turma diferentes.
Como parte de uma análise estatística adicional, a pontuação F2, priorizando a recordação, foi calculada conforme mostrado pela equação 20.

O desempenho do modelo também foi quantificado usando Erro Quadrático Médio (MSE), Erro Quadrático Médio Raiz (RMSE) e Erro Absoluto Médio (MAE) - que são tipicamente medições para tarefas de regressão, mas modificadas aqui para a tarefa de classificação para quantificar o quão grande era em média o erro, sem considerar a direção.
Para determinar se seria prático integrar ViTs em dispositivos de saúde voltados para o consumidor, realizamos resultados extensos de avaliação de desempenho focando exclusivamente na eficiência temporal do modelo. Criamos funções para relatar o tempo necessário para prever uma ou simplesmente múltiplas imagens, pois isso se torna particularmente relevante para aplicações em tempo real. Para cada imagem, antes de começar a prever, os dados são primeiro pré-processados para estarem no formato da imagem de entrada desejada pelo modelo. O tempo em que a previsão começou e o tempo em que a previsão foi concluída, nós registrávamos para obter os resultados do tempo e da latência. O tempo e a latência média foram calculados usando as equações 21 e 22, respectivamente.


Onde:
- N é o número de imagens (amostras).
- TEND é o tempo registrado após prever a i-ésima imagem.
- Tstart é o tempo registrado antes de prever a i-ésima imagem.
Experimentos adicionais foram realizados para avaliar a robustez do modelo Vision Transformer que propusemos, introduzindo sistematicamente ruído e desfoque adicionais às imagens. Ruído gaussiano era adicionado a cada pixel com fator de ruído de 0,4 enquanto os valores dos pixels se cortavam na faixa de [0,1]. O fator de desfoque foi diminuído por meio do desfoque Gaussiano com tamanho de grão de 45× 45. Esses experimentos são projetados para avaliar o modelo de forma abrangente sob degradações simuladas da qualidade perceptiva da imagem.