Deze studie werd goedgekeurd door de Commissie voor Klinisch Onderzoek en Ethische Zaken van het Shenzhen People's Hospital (Goedkeuringsnummer LL-KY-2025300-01). De verplichting tot het verkrijgen van geïnformeerde toestemming werd opgeheven door de ethische commissie, aangezien dit retrospectieve onderzoek gebruikmaakte van bestaande histopathologische afbeeldingen en klinische materialen zonder directe patiëntcontacten of interventies. Alle klinische gegevens en histopathologische afbeeldingen werden geanonimiseerd vóór analyse, en er werden geen persoonsgegevens gebruikt in dit onderzoek. Alle voor dit onderzoek gebruikte gegevens werden verwerkt overeenkomstig de ethische normen van de instellingsspecifieke en nationale onderzoekscommissies.
Kaderoverzicht
De algehele architectuur van het voorgestelde UAD-kader, RareCode, is geïllustreerd in Figuur 1Het kader maakt gebruik van een parallelle dual-branch-architectuur. Voor elke 512 x 512 ingangsafbeelding worden niet-overlappende 32 x 32 en 64 x 64 patches geëxtraheerd als respectievelijk fijnschalige en grofschalige ingangen. De 64 x 64 patches worden eerst herschaald naar 32 x 32 voordat ze aan het netwerk worden doorgegeven, zodat beide branches dezelfde encoder-decoder ingangsgrootte delen terwijl ze verschillende receptieve velden behouden. Encoders in elke branch projecteren de geëxtraheerde kenmerken naar een latente ruimte, waar de MHC-module discretisatiebeperkingen oplegt. Vervolgens reconstrueren decoders de afbeelding uit gekwantiseerde kenmerken om reconstructiefouten in het ruimtelijke domein te berekenen. Het CAR-mechanisme meet vervolgens de mate van semantische afwijking door codebookactiveringsfrequenties in kaart te brengen. Tot slot combineert het model reconstructiescores en CAR-scores via een gewogen combinatie om afbeeldingsniveau afwijkingscores te verkrijgen. Het RareCode-kader wordt end-to-end getraind met uitsluitend normale samples en vereist geen annotaties van afwijkingen. De dual-branch-architectuur werd gebruikt om zowel lokale cellulaire kenmerken als bredere klierstructuren te kunnen vastleggen. De MHC-module werd gebruikt om deze kenmerken te modelleren met codebooks van verschillende capaciteit.
Encoder-decoder-architectuur
RareCode construeert structureel onafhankelijke encoders en decoders voor de fijnschalige en grofschalige takken. Elke tak gebruikt hetzelfde encoder-decoderontwerp, maar deelt geen parameters. De encoder bestaat uit vier convolutieblokken, elk bestaande uit een 3 × 3-convolutie, batchnormalisatie, ReLU-activatie en dropout. De kanaalbreedte neemt toe van 64 naar 128, vervolgens naar 256, en de uiteindelijke feature map wordt afgevlakt en geprojecteerd naar een 64-dimensionale latente representatie. De decoder weerspiegelt de encoder met een volledig verbonden projectielaag, gevolgd door vier getransponeerde convolutielagen, waarbij elke patch wordt gereconstrueerd naar de oorspronkelijke invoergrootte van het netwerk. De reconstructiefout wordt berekend als de gemiddelde kwadratische fout tussen de invoer en de gereconstrueerde patches. Door deze encoder-decoderstructuur te combineren met discrete codebookbeperkingen, leert RareCode compacte representaties van normaal weefsel en meet het afwijkingen tijdens inferentie.
Multischaal hiërarchische codeboek
Om ziekteverschijnselen vast te leggen op meerdere niveaus van abstractie—van brede weefselstructuren tot gelokaliseerde cellulaire variaties—de MHC-module maakt gebruik van K discrete codeboeken
met verschillende capaciteiten n1, n2, ..., nK. In de uiteindelijke FourScales-configuratie bevat elke tak vier codeboeken met respectievelijk 64, 128, 256 en 512 invoerposities, waarbij elke codeboekinvoerpositie een 64-dimensionale embedding heeft. Voor elke gecodeerde patch-feature wordt de dichtstbijzijnde codeboekinvoerpositie geselecteerd op basis van Euclidische afstand. De gekwantiseerde uitvoer van de verschillende codeboeken wordt vervolgens samengevoegd met behulp van leerbare gewichten die genormaliseerd zijn over de codeboeken heen. Kleinere codeboeken, vanwege de sterkere compressiebeperkingen, worden verwacht grofkorrelige, prototypische patronen te coderen die lokale variaties abstract weergeven, terwijl grotere codeboeken fijnkorreligere kenmerken behouden die specifiekere structurele eigenschappen vastleggen. Voor elk codeboek
, waarbij ek(i)
Rd geeft de i-de insluitingsvector in het k-de codeboek aan, continue kenmerken worden via opzoeking van de dichtstbijzijnde buur (Eq. 1 en 2) afgebeeld op de discrete codeboekruimte:


Om adaptieve fusie van multischaalkenmerken te realiseren, introduceren we leerbare gewichten. Na softmax-normalisatie wordt een gewogen sommatie uitgevoerd over de gekwantiseerde uitvoer van elk codeboek (verg. 3):

waar σ(·) staat voor de softmax-functie, die ervoor zorgt dat de gewichten voldoen aan Σk σ(wk) = 1. Dit ontwerp stelt het model in staat om automatisch de bijdrageverhoudingen van codeboeken op verschillende granulariteitsniveaus aan te passen op basis van de semantische inhoud van de invoerkenmerken.
Om het leerproces van de codeboekoptimalisatie te verbeteren, gebruiken we de standaard VQ-verliesfunctie (Vgl. 4):

waar sg[·] duidt de stop-gradiëntoperatie aan, en β = 0,25 is de gewichtscoëfficiënt van het commitment-verlies. De eerste term stimuleert codeboekvectoren om zich te bewegen richting de uitvoer van de encoder, terwijl de tweede term de uitvoer van de encoder stimuleert om consistent te blijven met hun overeenkomstige codeboekvectoren. Het volledige mechanisme van multi-schaal VQ is geïllustreerd in Figuur 2.
Beoordeling van de zeldzaamheid van codeboekactivatie
CAR meet semantische afwijkingen op basis van codeboek-activatiestatistieken, geschat uit normale trainingsvoorbeelden. Na het trainen van het model werden alle normale trainingsafbeeldingen zonder gegevensverrijking door de encoder en de MHC-module geleid. Voor elke tak en elk codeboek werd de toewijzingsindex van elke patch-feature geregistreerd, en het aantal toewijzingen aan elk codeboek-element werd geteld. De tellingen werden vervolgens genormaliseerd om de activatiefrequentieverdeling voor dat codeboek te verkrijgen (Vgl. 5):

Tijdens de inferentie werd elke testafbeelding verwerkt via dezelfde procedure van patch-extractie, codering en toewijzing aan het dichtstbijzijnde codewoordenboek. Voor elke toegewezen codewoordenboekvermelding wordt de zeldzaamheidsscore berekend als de negatieve logaritme van de activeringsfrequentie in de normale trainingsset (Vgl. 6):

waar ε is een kleine constante die wordt gebruikt voor numerieke stabiliteit. Een codeboekitem met lage frequentie krijgt derhalve een hogere zeldzaamheidsscore, wat suggereert dat de overeenkomstige patch-feature minder consistent is met de geleerde verdeling van normaal weefsel. Zeldzaamheidsscores op patch-niveau werden gemiddeld binnen elk beeld en over de fijne en grove schaalstructuren heen om een beeldniveau-zeldzaamheidsrespons te verkrijgen.
Om de activatiefrequentierariteit aan te vullen, berekenen we ook een percentielgebaseerde kwantisatie-afstandsscore (Eq. 7):

Deze score is afgeleid van de Euclidische afstand tussen elke gecodeerde kenmerkvector en de dichtstbijzijnde codeboekvermelding. De afstandsverdeling wordt geschat op basis van normale trainingsmonsters, en de afstanden van testmonsters worden omgezet in percentielscores. Hogere percentielscores geven aan dat het gecodeerde kenmerk moeilijker is weer te geven met behulp van de geleerde prototypen van het normale codeboek.
De uiteindelijke CAR-score combineert activatiezeldzaamheid en kwantisatie-afstand over alle codeboeken heen, gebruikmakend van de leerbare codeboekgewichten (Eq. 8):

Dezelfde genormaliseerde codebookgewichten die worden gebruikt voor multischaalfusie van kenmerken, worden toegepast op fusie op scoresniveau, waardoor consistentie wordt behouden tussen representatieleren en afwijkingsscorering. De resulterende CAR-score wordt op afbeeldingsniveau berekend en vervolgens gecombineerd met de reconstitutiescore om de uiteindelijke afwijkingsscore te verkrijgen. Het algehele CAR-scoringsmechanisme is geïllustreerd in Figuur 3.
Definitieve berekening van de anomalie-score
We combineren fouten in de ruimtelijke reconstitutie met semantische CAR-scores. Voordien worden beide scores afzonderlijk genormaliseerd met behulp van percentielgebaseerde min-max-normalisatie om de invloed van extreme uitbijters te verkleinen. De uiteindelijke anomalie-score op beeldniveau wordt vervolgens gedefinieerd als (Vgl. 9):

waar Sreconstructie geeft de genormaliseerde gemiddelde kwadratische reconstructiefout aan, die de reconstructiekwaliteit van de sample in pixelruimte weerspiegelt; SCAR is de CAR-score zoals hierboven beschreven, die graden van semantische afwijking weergeeft. De hyperparameter α
[0,1] regelt de relatieve weging tussen de twee componenten, waarbij optimale waarden worden bepaald op validatiesets.
Trainingsdoelstelling
Het RareCode-framework wordt end-to-end getraind met alleen normale samples. Het totale trainingsverlies bestaat uit de volgende componenten (Vgl. 10):

Elke term is als volgt gedefinieerd: Reconstructieverlies:
, het meten van de gemiddelde kwadratische fout tussen invoerpatches p en gereconstrueerde patches D(zq) na kwantisering. Dit verlies wordt afzonderlijk toegepast op de 32 × 32-patch tak (
) en 64 × 64 patch-vertakking (
) om effectief leren van kenmerken op beide ruimtelijke schalen te waarborgen.
Beschrijving van de dataset
Om RareCode te evalueren, zijn experimenten uitgevoerd op een klinisch geannoteerde histopathologische beeldverzameling van colorectale kanker (CRC) die beelden bevat van colorectale weefselsecties gekleurd met hematoxyline en eosine (H&E)&E), verkregen van Shenzhen People's ziekenhuis. Alle afbeeldingen in deze dataset zijn afkomstig van echte klinische gevallen. De afbeeldingen zijn gedigitaliseerd bij 40x vergroting met een oorspronkelijke resolutie van 1024 × 1024 pixels. Om de betrouwbaarheid en kwaliteit van de annotaties te waarborgen, zijn alle labels van de steekproefcategorieën gezamenlijk beoordeeld en bevestigd door twee of meer senior professionele pathologen volgens een dubbelblinde protocol.
Op basis van histopathologische kenmerken werd de dataset onderverdeeld in drie klassen: normaal weefsel (1.226 afbeeldingen), kanker (1.215 afbeeldingen) en ontsteking (1.214 afbeeldingen). In de binaire UAD-instelling werden kanker- en ontstekingsmonsters beschouwd als afwijkend, terwijl normale monsters dienden als referentieverdeling. Deze formulering weerspiegelt de beoogde triagetask om gevallen die verdere pathologische beoordeling vereisen, te scheiden van gevallen die morfologisch normaal zijn.
Alle H&E-gekleurde afbeeldingen werden herschaald naar 512 × 512 pixels voordat ze als netwerkinput worden gebruikt. Voor de evaluatie werd gestratificeerde 5-voudige kruisvalidatie toegepast. In elke vouw werden de normale monsters eerst verdeeld in trainings-, validatie- en aparte testdeelverzamelingen. RareCode werd uitsluitend getraind op de normale trainingsmonsters. De validatiedeelverzameling, die eveneens uitsluitend normale monsters bevatte, werd gebruikt voor modelselectie, afstemming van hyperparameters en selectie van het fusiegewicht αDe apart gehouden testset, bestaande uit niet eerder gezien normale monsters en afwijkende monsters waaronder kanker en ontsteking, werd uitsluitend gebruikt voor de definitieve beoordeling van de prestaties en niet gebruikt voor het trainen van het model, de selectie van hyperparameters of α selectie. De representatieve categorieën, de opbouw van de dataset en het validatieprotocol zijn weergegeven in Figuur 4.
Implementatiedetails
Het voorgestelde RareCode-framework en de vergelijkingsmodellen zijn geïmplementeerd met behulp van PyTorch. Alle experimenten zijn uitgevoerd op een werkstation uitgerust met een enkele NVIDIA GeForce RTX 4060 Ti GPU (16 GB). Voor de configuratie van de netwerkarchitectuur, om kenmerken op verschillende ruimtelijke schalen te kunnen vastleggen, verwerkt de dual-branch architectuur niet-overlappende patches op twee schalen: 32 × 32 pixels (kleiner receptief veld, lokale textuurstrepen vastleggend) en 64 × 64 pixels (groter receptief veld, waardoor een breder ruimtelijk verband wordt vastgelegd), respectievelijk. De continue dimensie van de functie-embedding die door beide takencoders wordt geproduceerd, is uniform ingesteld op 64. In de MHC-module configureren we vier codeboeken met verschillende capaciteiten voor elke tak, met discrete prototype-aantallen van respectievelijk 64, 128, 256 en 512.
Tijdens de optimalisatie werd het netwerk end-to-end getraind gedurende 50 epochs met een batchgrootte van 8. We gebruiken de AdamW-optimizer voor gewichtsupdates, met een initiële leersnelheid van 5 × en gewichtsverval om overfitting te voorkomen. Daarnaast werd, om een vlotte convergentie tijdens de training te garanderen, gebruikgemaakt van een leerfrequentie-schema met cosinus-annealing, waardoor een fijnere optimalisatie in latere trainingsfases mogelijk werd. De patch-groottes en codeboek-groottes werden geselecteerd om een balans te bieden tussen multi-schaalrepresentatie en computationele kosten. De uiteindelijke FourScales-configuratie werd ondersteund door ablatie-analyse, en α werd geselecteerd op basis van de validatieset alvorens de testset te evalueren.
Voor de afbraak van de decodercomplexiteit is een variant met een complexe decoder geïmplementeerd op basis van de FourScales-configuratie. Deze variant gebruikte dezelfde datasplitsing, codeboekgroottes, trainingsepochs, optimizer, leerfrequentie, op validatie gebaseerde α selectie en evaluatiemetingen als het FourScales-model. De basise decoder werd vervangen door een EMCAD-stijl decoder met multi-schaal depthwise convolutieblokken en Convolutional Block Attention Modules (CBAM). Na elk van de eerste drie getransponeerde-convolutie-omhoogsamplingstappen, parallel 3 × 3, 5 × 5 en 7 × 7 diepteconvoluties werden toegepast, gevolgd door 1 × 1 puntsgewijze fusie, batchnormalisatie, ReLU-activatie, een residuverbinding en CBAM-aandacht. CBAM omvatte kanaalaandacht op basis van gemiddelde- en max-poolingdescriptoren en ruimtelijke aandacht met behulp van een 7 × 7 convoluties. De laatste fase maakte gebruik van getransponeerde convolutie en sigmoid-activatie om 32 te reconstrueren × 32 patches.