Paso 1: Cargar imágenes

Este ejemplo实用功能imágenes三维空间de la cabeza del mismo paciente。关于登记的问题，关于图像的问题，关于图像的问题，关于图像的问题，关于图像的问题，关于运动的问题。线性记录的目标，运动的意象，和幻想的意象。在这里，在我的脑海里，在我的脑海里，在我的脑海里RMN ponderada T1。这是一种运动的形象，它是一种文化的形象。关于档案利用程序和资料回顾evaluación登记imágenes (RIRE)。图象的二进制档案和连续数据利用。multibandreadUtilice la función para obtener los metadata asociados a cada imagen。helperReadHeaderRIREPuede实用工具el siguiente enlace para encontrar más información档案馆的严肃格式:数据格式

fixedHeader = helperReadHeaderRIRE(“rirePatient007MRT1.header”）;movingHeader = helperReadHeaderRIRE(“rirePatient007CT.header”）;fixedVolume = multibandread(“rirePatient007MRT1.bin”，.．.[fixedHeader。行,fixedHeader。列、fixedHeader.Slices]…'int16=>single'， 0， 'bsq'， ' ee-be');movingVolume = multibandread('rirePatient007CT.bin'，…[movingHeader。行,movingHeader。列、movingHeader.Slices]…'int16=>single'， 0， 'bsq'， ' ee-be');

La función es una función辅助比例para ayudar a橘子眼La calidad de los resulttos de registro 3D。helperVolumeRegistrationPuede girar interactivamente la vista y ambos ejes permanecerán sincronizados。

helperVolumeRegistration (fixedVolume movingVolume);

También在地球上的未来空间和个人空间volúmenes在地球上的运动方式和思想alineación在地球上的将军volúmenes。imshowpairEn la imagenen superpuesta de, las áreas灰色对应áreas que tienen intensidades similares, mientras que las áreas品红和绿色的muestran lugares donde una imagenes más灿烂的la otra。imshowpair对不正确的观测记录的利用volúmenes图像的半分割的轴向番茄través中心的卷。imshowpair

centerFixed = size(fixedVolume)/2;centerMoving = size(movingVolume)/2;数字imshowpair (movingVolume (:,:, centerMoving (3)), fixedVolume (:,:, centerFixed (3)))；标题(“未注册轴向切片”）

Paso 2: Configurar el Registro Inicial

La función facilita La elección del optimizador y La configuración de métricas correctos para usar con。imregconfigimregisterEstas dos imágenes不同的程序，la RMN y la TC, por lo que la opción“多模态”es adecuada。

[optimizer,metric] = imregconfig(“多通道”）;

El algoritmo utilizado por convergerá para obtener mejores resulttos más rápidamente cuando se especque información de referencia space清醒la resolución y/o ubicación de las imágenes de entrada。imregisterEn este caso, la resolución de los conjuntos de datos CT y MRI se define En los元数据de la imagen。运用元数据和裁判对象，参考空间和通行原则，和记事的分项论证。imref3d

Rfixed = imref3d(size(fixedVolume)，fixedHeader.PixelSize(2)，fixedHeader.PixelSize(1)，fixedHeader.SliceThickness);Rmoving = imref3d(size(movingVolume)，movingHeader.PixelSize(2)，movingHeader.PixelSize(1)，movingHeader.SliceThickness);

空间定义参考对象的如意算盘dónde空间定义参考对象的如意算盘volúmenes宇宙协调系统的如意算盘cuál空间定义参考对象的如意算盘extensión空间定义参考对象的如意算盘dimensión。La丙达XWorldLimits de Rmoving定义La posición del volumen móvil en La dimensión X. La丙达PixelExtentInWorld定义el tamaño de cada píxel en unidades del mundo en La dimensión X(一个lo largo de las columnas)。地球外运动体积0.3269 a 334.97坐标系统X地球中心píxel tiene una extensión de 0.6536mm。Las unidades están en milímetros porque la información de encabezado utilizada para conconir la referencia spatial established en milímetros。La丙烯，图像，范围，世界，决定，La extensión，完整的世界，图像，活动，世界，完整的世界。

Rmoving。XWorldLimits

ans =1×20.3268 - 334.9674

Rmoving。PixelExtentInWorldX

Ans = 0.6536

Rmoving。ImageExtentInWorldX

Ans = 334.6406

La desalineación entre los dos volúmenes incluye La conversión, el escalado y La rotación。使用una transformación de similud para registrar las imágenes。

计算计算计算，计算计算，计算计算，计算计算，计算计算。imregister

特别的configuración没有预先决定，对一个初始化半径，对一个优化半径，对一个主要的，收敛的结果，对登记。

优化器。InitialRadius = 0.004;movingRegisteredVolume = imregister(movingVolume,Rmoving, fixedVolume,Rfixed，“刚性”，优化器，度量);

使用新记忆程序进行检查alineación不分段轴向tomado través在中心位置volúmenes记录者的新记忆方法在中心位置记录者的想法。imshowpair

数字imshowpair (movingRegisteredVolume (:,:, centerFixed (3)), fixedVolume (:,:, centerFixed (3)))；标题(“登记体积轴向切片”）

前侧轴向，注册中心realizó正确。Utilícelo登记新记录，登记新记录，登记新记录，登记新记录，登记新记录，登记新记录，登记新记录，登记新记录，登记新记录。helperVolumeRegistration

helperVolumeRegistration (fixedVolume movingRegisteredVolume);

第3节:观察者transformación geométrica三维的运动与幻想。

La función se puede效用，cuando está interesado enla estimación de transformación geométrica que se效用，para formar La imagen de salida registrada。使用错误的算法和错误的论证。imregtformimregisterimregtformimregisterimregister我们在一起inspección我们在一起realizó我们在一起transformación geométrica我们在一起transformación geométrica我们在一起。imregisterimregtform

geomtform = imregtform(movingVolume,Rmoving, fixedVolume,Rfixed，“刚性”，优化器，度量)

几何形式= affine3d与属性:维度:3 T: [4x4双]

结果不符合要求transformación geométrica。Este对象包括una丙，T, que定义la matriz de transformación afín 3D。

geomtform。T

ans =4×40.9704 -0.0143 -0.2410 0 0.0228 0.9992 0.0324 0 0.2404 -0.0369 0.9700 0 -15.8648 -17.5692 29.1830 1.0000

El método transformPointsForward de la transformación geométrica se puede utizar para determinar dónde se puede determinar dónde se puede utizar un punto [u,v,w] en los mapas de imágenes en movimiento como resultado del registro。关于宇宙空间的参考文献和特定形式的数据，transformación geométrica关于宇宙空间的协调系统的数据。El método transformPointsForward se utiliza a continuación para determinar la ubicación transformada del centro de la imagen en movimiento en El sistema de coordenadas universe。

centerXWorld = mean(remove . xworldlimits);centerYWorld = mean(删除。yworldlimits);centerZWorld = mean(remove . zworldlimits);[xWorld,yWorld,zWorld] = transformPointsForward(geomtform,centerXWorld,centerYWorld,centerZWorld);

Puede效用método worldto下标de Rfixed para determinar元素del volume fijo que se alinea和centro del volume en moviento。

[r,c,p] = worldToSubscript(Rfixed,xWorld,yWorld,zWorld)

R = 116

C = 132

P = 13

第4段:欣赏estimación de transformación geométrica画面和动作的体积。

La función se puede utilzar para aplicar La estimación de transformación geométrica de imregtform un volumen 3D。imwarpEl nombre/valor 'OutputView' se效用para definir un论证de referencia决定性空间límites世界和la resolución de la imagen remuestreada de salida。在世界上产生了许多错误的结果，中间的参考对象，空间和意象。imregister为了创造一个在我们心中的灵魂límites在我们心中的世界resolución在我们心中的世界和我们的生命。在我们的世界límites世界resolución我们的行动volúmenes我们的世界，我们的通信píxeles píxeles我们的行动volúmenes我们的行动。

movingRegisteredVolume = imwarp(movingVolume,Rmoving,geomtform，“双三次的”，“OutputView”, Rfixed);

利用imshowpair新para un扇形轴向través登记生产体积中心。imwarp

数字imshowpair (movingRegisteredVolume (:,:, centerFixed (3)), fixedVolume (:,:, centerFixed (3)))；标题(“登记体积轴向切片”）