User Tools

Site Tools


neuroimagen:pipe05

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
neuroimagen:pipe05 [2022/01/10 10:45]
osotolongo [Bajar proyecto de XNAT]
neuroimagen:pipe05 [2022/05/02 08:34] (current)
osotolongo [Dependencias]
Line 1: Line 1:
 ====== Pipeline de Procesamiento de FACE v0.5 ====== ====== Pipeline de Procesamiento de FACE v0.5 ======
  
-Nueva version en etapa de planificacion. 
  
 ===== Porque una nueva version? ===== ===== Porque una nueva version? =====
  
-  * // Introducir containers. Preparar el pipeline para utilizarlo en Microsoft Azure, AWS o similar. Todas las llamadas a FSL, Freesurfer o ANTs deben sustituirse por una llamada a containers. // **Nota: De momento no parece ser necesario y los containers de FSL o Fs pueden ser mas complicados de lo que parecen** DELETEME 
   * <del>Modularizacion de las llamadas a SLURM. Los script sbatch han de crearse y lanzarse a traves del modulo SLURM.pm. Deben reescribirse los scripts.</del> DONE ([[neuroimagen:neuro4.pm|Modulos del pipeline]])   * <del>Modularizacion de las llamadas a SLURM. Los script sbatch han de crearse y lanzarse a traves del modulo SLURM.pm. Deben reescribirse los scripts.</del> DONE ([[neuroimagen:neuro4.pm|Modulos del pipeline]])
   * Reparalelizar los procesos de registro, y extraccion de metricas PET en modulos mas pequeños :?:   * Reparalelizar los procesos de registro, y extraccion de metricas PET en modulos mas pequeños :?:
Line 13: Line 11:
   * <del>Usar Freesurfer 7.2 [[https://surfer.nmr.mgh.harvard.edu/fswiki/ReleaseNotes|Release Notes]]</del> DONE ([[neuroimagen:freesurfer7|Update Freesurfer 7.2]])   * <del>Usar Freesurfer 7.2 [[https://surfer.nmr.mgh.harvard.edu/fswiki/ReleaseNotes|Release Notes]]</del> DONE ([[neuroimagen:freesurfer7|Update Freesurfer 7.2]])
  
-===== TODO =====  
  
-  - acenip Containers? DELETEME 
- 
-{{ :neuroimagen:map_container_build.png?600 |}} 
  
 ===== Dependencias ===== ===== Dependencias =====
Line 28: Line 22:
   * fmriprep: https://fmriprep.org/en/stable/   * fmriprep: https://fmriprep.org/en/stable/
   * dicom3tools: https://www.dclunie.com/dicom3tools.html (//no estrictamente necesario pero conveniente//)   * dicom3tools: https://www.dclunie.com/dicom3tools.html (//no estrictamente necesario pero conveniente//)
-  * xnatapic: https://github.com/asqwerty666/xnatapic (para las operaciones con XNAT)+  * xnatapic: https://github.com/asqwerty666/xnatapic (para las operaciones con XNAT). Intentando eliminar estas dependencias
  
 ===== DB ===== ===== DB =====
Line 502: Line 496:
 </code> </code>
  
-===== Bajar proyecto de XNAT =====+==== Freesurfer QC ==== 
 + 
 +[[neuroimagen:visualqc_freesurfer_xnat|Realizar el QC a la segmentacion de Freesurfer de los sujetos analizados en XNAT]] 
 + 
 +==== Bajar proyecto de XNAT ====
  
 Una a vez procesada la segmentación de Freesurfer, si queremos hacer algo mas complejo, hemos de bajar ordenadamente el directorio de procesamiento de cada sujeto. El primer paso sería crear el proyecto local, pero no es necesario (al menos en principio) identificar la posicion de las imagenes RAW. Es decir que practicamente se puede hacer a mano a partir del listado de XNAT. Una a vez procesada la segmentación de Freesurfer, si queremos hacer algo mas complejo, hemos de bajar ordenadamente el directorio de procesamiento de cada sujeto. El primer paso sería crear el proyecto local, pero no es necesario (al menos en principio) identificar la posicion de las imagenes RAW. Es decir que practicamente se puede hacer a mano a partir del listado de XNAT.
Line 527: Line 525:
 [osotolongo@brick03 mri_face]$ mkdir working [osotolongo@brick03 mri_face]$ mkdir working
 [osotolongo@brick03 mri_face]$ cat ~/.config/neuro/mriface.cfg [osotolongo@brick03 mri_face]$ cat ~/.config/neuro/mriface.cfg
-DATA = /old_nas/mriface+DATA = /old_nas/mri_face
 SRC = /old_nas/mri_face/raw SRC = /old_nas/mri_face/raw
-WORKING = /old_nas/mriface/working +WORKING = /old_nas/mri_face/working 
-BIDS = /old_nas/mriface/bids+BIDS = /old_nas/mri_face/bids
 </code> </code>
  
Line 566: Line 564:
 </code> </code>
  
 +y ahora,
 +
 +<code bash>
 +[osotolongo@brick03 mri_face]$ xnat_getfs.pl -p mriface -x unidad
 +Getting XNAT subject list
 +xnatapic list_subjects --project_id unidad --label > /old_nas/mri_face/xnat_subjects.list
 +Getting experiments
 +Sorting, joining, keeping shit together
 +OK, now I'm going on
 +Downloading and extracting
 +xnatapic get_fsresults --experiment_id XNAT05_E00046 --all-tgz /tmp/fstar.ikl0N3/mriface_0001.tar.gz
 +...
 +...
 +...
 +</code>
 +
 +y queda todo almacenado como si se hubiera ejecutado FS localmente.
  
  
Line 572: Line 587:
 ===== Hippocampal Subfields ===== ===== Hippocampal Subfields =====
  
-A esta version se le ha añadido la capacidad de extraer la segmentacion del hipocampo y la amigadala+(detalles en [[neuroimagen:hsf|Segmentacion del hipocampo]]) 
 + 
 +A esta version se le ha añadido la capacidad de extraer la segmentacion del hipocampo y la amigdala. Primero debemos tener el proyecto creado y la segmentacion de FS en local. 
 + 
 +Ahora solo hacemos, 
 + 
 +<code bash> 
 +[osotolongo@brick03 mri_face]$ hsfrecon.pl mriface  
 +[osotolongo@brick03 acenip]$ queue 
 +   JOBID PARTITION                                 NAME     USER ST       TIME  NODES NODELIST(REASON) 
 +  218945      fast                    hsf_recon_mriface osotolon  R       2:16      1 brick02 
 +  218946      fast                    hsf_recon_mriface osotolon  R       2:06      1 brick02 
 +  218947      fast                    hsf_recon_mriface osotolon  R       1:56      1 brick02 
 +  218948      fast                    hsf_recon_mriface osotolon  R       1:46      1 brick02 
 +  218949      fast                    hsf_recon_mriface osotolon  R       1:36      1 brick02 
 +  218950      fast                    hsf_recon_mriface osotolon  R       1:26      1 brick02 
 +  218951      fast                    hsf_recon_mriface osotolon  R       1:16      1 brick02 
 +  218952      fast                    hsf_recon_mriface osotolon  R       1:06      1 brick02 
 +  218953      fast                    hsf_recon_mriface osotolon  R       0:56      1 brick02 
 +  218954      fast                    hsf_recon_mriface osotolon  R       0:46      1 brick02 
 +  218955      fast                    hsf_recon_mriface osotolon  R       0:36      1 brick02 
 +  218956      fast                    hsf_recon_mriface osotolon  R       0:26      1 brick02 
 +  218957      fast                    hsf_recon_mriface osotolon  R       0:16      1 brick02 
 +  218958      fast                    hsf_recon_mriface osotolon  R       0:06      1 brick02 
 +... 
 +... 
 +... 
 + 
 +</code> 
 + 
 +y, tras recibir el aviso de finalización, 
 + 
 +{{:neuroimagen:ab6f090c-66a1-4e44-a140-ead38526f604.jpeg?400|}} 
 + 
 +se sacan las métricas, 
 + 
 +<code bash> 
 +[osotolongo@brick03 mri_face]$ hsf_metrics.pl mriface 
 +</code> 
 + con lo que los resultados quedan en //mriface_hsf_metrics.csv//
 + 
 +<code bash> 
 +[osotolongo@brick03 mri_face]$ head mriface_hsf_metrics.csv 
 +Subject,lh.CA1-body,lh.CA1-head,lh.CA3-body,lh.CA3-head,lh.CA4-body,lh.CA4-head,lh.GC-ML-DG-body,lh.GC-ML-DG-head,lh.HATA,lh.Hippocampal_tail,lh.Whole_hippocampal_body,lh.Whole_hippocampal_head,lh.Whole_hippocampus,lh.fimbria,lh.hippocampal-fissure,lh.molecular_layer_HP-body,lh.molecular_layer_HP-head,lh.parasubiculum,lh.presubiculum-body,lh.presubiculum-head,lh.subiculum-body,lh.subiculum-head,rh.CA1-body,rh.CA1-head,rh.CA3-body,rh.CA3-head,rh.CA4-body,rh.CA4-head,rh.GC-ML-DG-body,rh.GC-ML-DG-head,rh.HATA,rh.Hippocampal_tail,rh.Whole_hippocampal_body,rh.Whole_hippocampal_head,rh.Whole_hippocampus,rh.fimbria,rh.hippocampal-fissure,rh.molecular_layer_HP-body,rh.molecular_layer_HP-head,rh.parasubiculum,rh.presubiculum-body,rh.presubiculum-head,rh.subiculum-body,rh.subiculum-head,whole.CA1-body,whole.CA1-head,whole.CA3-body,whole.CA3-head,whole.CA4-body,whole.CA4-head,whole.GC-ML-DG-body,whole.GC-ML-DG-head,whole.HATA,whole.Hippocampal_tail,whole.Whole_hippocampal_body,whole.Whole_hippocampal_head,whole.Whole_hippocampus,whole.fimbria,whole.hippocampal-fissure,whole.molecular_layer_HP-body,whole.molecular_layer_HP-head,whole.parasubiculum,whole.presubiculum-body,whole.presubiculum-head,whole.subiculum-body,whole.subiculum-head 
 +mriface_0001,160.140423,643.111512,121.733859,181.450963,136.089629,176.885647,142.768001,214.595877,75.325365,675.274931,1293.282294,2090.659074,4059.216299,74.752421,188.162283,249.214902,399.634974,65.916415,147.868612,140.671683,260.714447,193.066639,160.417880,557.584677,144.919215,150.610933,158.435914,146.932857,167.188920,174.063578,58.895612,646.290298,1294.041853,1785.289454,3725.621605,64.672166,186.310011,244.291453,344.342486,61.677945,112.838806,121.130778,241.277501,170.050586,320.558303,1200.696189,266.653074,332.061896,294.525543,323.818504,309.956921,388.659455,134.220977,1321.565229,2587.324147,3875.948528,7784.837904,139.424587,374.472294,493.506355,743.97746,127.59436,260.707418,261.802461,501.991948,363.117225 
 +mriface_0002,103.287906,482.663069,70.404188,136.628715,108.229725,130.864657,123.395396,154.621026,59.824302,483.804995,1049.338826,1653.070021,3186.213842,77.842256,135.308845,188.309856,312.397416,75.537688,143.016055,140.327036,234.853444,160.206112,94.023260,490.761234,91.329668,133.129240,126.997684,132.734252,141.050346,152.248955,60.223939,448.709529,1104.744740,1655.016815,3208.471084,61.685305,166.993028,202.158720,317.133128,73.087146,146.284817,132.921255,241.214940,162.777666,197.311166,973.424303,161.733856,269.757955,235.227409,263.598909,264.445742,306.869981,120.048241,932.514524,2154.083566,3308.086836,6394.684926,139.527561,302.301873,390.468576,629.530544,148.624834,289.300872,273.248291,476.068384,322.983778 
 +mriface_0003,93.472053,292.211426,63.399442,75.138635,83.350543,79.036618,89.423755,86.168844,28.511534,351.012237,744.342056,991.830003,2087.184295,27.602650,125.530329,142.094642,188.495860,53.322443,103.487495,80.480452,141.511475,108.464190,106.828836,342.747996,71.336504,75.118041,90.448081,86.619239,95.673876,95.907348,32.929040,404.474760,846.803326,1187.683959,2438.962045,30.090636,164.192786,160.769831,228.038129,68.750950,119.683304,102.904924,171.972258,154.668293,200.300889,634.959422,134.735946,150.256676,173.798624,165.655857,185.097631,182.076192,61.440574,755.486997,1591.145382,2179.513962,4526.14634,57.693286,289.723115,302.864473,416.533989,122.073393,223.170799,183.385376,313.483733,263.132483 
 +mriface_0004,90.025085,425.035298,61.411596,95.384809,85.516784,100.256754,93.233374,112.652146,55.431434,376.308526,872.067514,1349.532440,2597.908481,41.677259,159.021456,159.268126,242.911788,69.730972,134.091369,110.008671,206.843921,138.120570,110.674897,370.394924,62.026268,95.396101,92.756480,97.983884,102.218936,105.636656,52.204311,385.544462,849.524182,1210.263846,2445.332491,34.370706,176.590676,159.272561,218.078316,58.590188,103.594471,85.746939,184.609864,126.232527,200.699982,795.430222,123.437864,190.78091,178.273264,198.240638,195.45231,218.288802,107.635745,761.852988,1721.591696,2559.796286,5043.240972,76.047965,335.612132,318.540687,460.990104,128.32116,237.68584,195.75561,391.453785,264.353097 
 +mriface_0005,93.859758,324.806652,77.916251,83.807912,89.245013,87.995024,94.582035,96.679923,49.024978,414.336672,810.804979,1157.485099,2382.626750,37.909804,141.689681,154.416912,208.037965,68.748935,98.145165,107.828875,164.730041,130.554836,119.825141,406.823919,85.275754,95.169639,94.333964,107.044569,102.725566,121.103714,48.089261,506.876257,921.452891,1405.736164,2834.065313,37.298695,155.083729,182.466672,265.513688,81.957020,102.535701,121.564807,196.991399,158.469547,213.684899,731.630571,163.192005,178.977551,183.578977,195.039593,197.307601,217.783637,97.114239,921.212929,1732.25787,2563.221263,5216.692063,75.208499,296.77341,336.883584,473.551653,150.705955,200.680866,229.393682,361.72144,289.024383 
 +</code>
  
 ===== Preprocesamiento SBM ===== ===== Preprocesamiento SBM =====
  
-<code> +El analisis SBM debe prepararse y hacerse independientemente del flujo de proyecto ya que requiere preparar una DB independiente con datos externos a neuroimagen. No obstante, pueden preprocesarse las MRI dentro de la estructura general del proyecto.  
-pqcache.pl+ 
 +Para esto basta con ejecutar un //recon_all -qcache -subjid project_subjectid//. Asi que es tan simple como recorrer la DB y ejecutar esto en paralelo con, 
 + 
 + 
 +<code bash
 +pqcache.pl proyecto
 </code> </code>
-==== Freesurfer QC ==== 
  
-[[neuroimagen:visualqc_freesurfer_xnat|Realizar el QC a la segmentacion de Freesurfer de los sujetos analizados en XNAT]]+y ahi quedan todas las medias preparadas para distintos valores de FWHM, 
 + 
 +<code bash> 
 +[osotolongo@brick03 ~]$ ls /old_nas/subjects/bioface_0001/surf/*fwhm* 
 +/old_nas/subjects/bioface_0001/surf/lh.area.fwhm0.fsaverage.mgh             /old_nas/subjects/bioface_0001/surf/lh.volume.fwhm25.fsaverage.mgh          /old_nas/subjects/bioface_0001/surf/rh.jacobian_white.fwhm15.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.area.fwhm10.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/lh.volume.fwhm5.fsaverage.mgh           /old_nas/subjects/bioface_0001/surf/rh.jacobian_white.fwhm20.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.area.fwhm15.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/lh.w-g.pct.mgh.fwhm0.fsaverage.mgh      /old_nas/subjects/bioface_0001/surf/rh.jacobian_white.fwhm25.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.area.fwhm20.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/lh.w-g.pct.mgh.fwhm10.fsaverage.mgh     /old_nas/subjects/bioface_0001/surf/rh.jacobian_white.fwhm5.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.area.fwhm25.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/lh.w-g.pct.mgh.fwhm15.fsaverage.mgh     /old_nas/subjects/bioface_0001/surf/rh.sulc.fwhm0.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.area.fwhm5.fsaverage.mgh             /old_nas/subjects/bioface_0001/surf/lh.w-g.pct.mgh.fwhm20.fsaverage.mgh     /old_nas/subjects/bioface_0001/surf/rh.sulc.fwhm10.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.area.pial.fwhm0.fsaverage.mgh        /old_nas/subjects/bioface_0001/surf/lh.w-g.pct.mgh.fwhm25.fsaverage.mgh     /old_nas/subjects/bioface_0001/surf/rh.sulc.fwhm15.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.area.pial.fwhm10.fsaverage.mgh       /old_nas/subjects/bioface_0001/surf/lh.w-g.pct.mgh.fwhm5.fsaverage.mgh      /old_nas/subjects/bioface_0001/surf/rh.sulc.fwhm20.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.area.pial.fwhm15.fsaverage.mgh       /old_nas/subjects/bioface_0001/surf/lh.white.H.fwhm0.fsaverage.mgh          /old_nas/subjects/bioface_0001/surf/rh.sulc.fwhm25.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.area.pial.fwhm20.fsaverage.mgh       /old_nas/subjects/bioface_0001/surf/lh.white.H.fwhm10.fsaverage.mgh         /old_nas/subjects/bioface_0001/surf/rh.sulc.fwhm5.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.area.pial.fwhm25.fsaverage.mgh       /old_nas/subjects/bioface_0001/surf/lh.white.H.fwhm15.fsaverage.mgh         /old_nas/subjects/bioface_0001/surf/rh.thickness.fwhm0.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.area.pial.fwhm5.fsaverage.mgh        /old_nas/subjects/bioface_0001/surf/lh.white.H.fwhm20.fsaverage.mgh         /old_nas/subjects/bioface_0001/surf/rh.thickness.fwhm10.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.curv.fwhm0.fsaverage.mgh             /old_nas/subjects/bioface_0001/surf/lh.white.H.fwhm25.fsaverage.mgh         /old_nas/subjects/bioface_0001/surf/rh.thickness.fwhm15.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.curv.fwhm10.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/lh.white.H.fwhm5.fsaverage.mgh          /old_nas/subjects/bioface_0001/surf/rh.thickness.fwhm20.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.curv.fwhm15.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/lh.white.K.fwhm0.fsaverage.mgh          /old_nas/subjects/bioface_0001/surf/rh.thickness.fwhm25.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.curv.fwhm20.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/lh.white.K.fwhm10.fsaverage.mgh         /old_nas/subjects/bioface_0001/surf/rh.thickness.fwhm5.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.curv.fwhm25.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/lh.white.K.fwhm15.fsaverage.mgh         /old_nas/subjects/bioface_0001/surf/rh.volume.fwhm0.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.curv.fwhm5.fsaverage.mgh             /old_nas/subjects/bioface_0001/surf/lh.white.K.fwhm20.fsaverage.mgh         /old_nas/subjects/bioface_0001/surf/rh.volume.fwhm10.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.jacobian_white.fwhm0.fsaverage.mgh   /old_nas/subjects/bioface_0001/surf/lh.white.K.fwhm25.fsaverage.mgh         /old_nas/subjects/bioface_0001/surf/rh.volume.fwhm15.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.jacobian_white.fwhm10.fsaverage.mgh  /old_nas/subjects/bioface_0001/surf/lh.white.K.fwhm5.fsaverage.mgh          /old_nas/subjects/bioface_0001/surf/rh.volume.fwhm20.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.jacobian_white.fwhm15.fsaverage.mgh  /old_nas/subjects/bioface_0001/surf/rh.area.fwhm0.fsaverage.mgh             /old_nas/subjects/bioface_0001/surf/rh.volume.fwhm25.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.jacobian_white.fwhm20.fsaverage.mgh  /old_nas/subjects/bioface_0001/surf/rh.area.fwhm10.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/rh.volume.fwhm5.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.jacobian_white.fwhm25.fsaverage.mgh  /old_nas/subjects/bioface_0001/surf/rh.area.fwhm15.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/rh.w-g.pct.mgh.fwhm0.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.jacobian_white.fwhm5.fsaverage.mgh   /old_nas/subjects/bioface_0001/surf/rh.area.fwhm20.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/rh.w-g.pct.mgh.fwhm10.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.sulc.fwhm0.fsaverage.mgh             /old_nas/subjects/bioface_0001/surf/rh.area.fwhm25.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/rh.w-g.pct.mgh.fwhm15.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.sulc.fwhm10.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/rh.area.fwhm5.fsaverage.mgh             /old_nas/subjects/bioface_0001/surf/rh.w-g.pct.mgh.fwhm20.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.sulc.fwhm15.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/rh.area.pial.fwhm0.fsaverage.mgh        /old_nas/subjects/bioface_0001/surf/rh.w-g.pct.mgh.fwhm25.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.sulc.fwhm20.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/rh.area.pial.fwhm10.fsaverage.mgh       /old_nas/subjects/bioface_0001/surf/rh.w-g.pct.mgh.fwhm5.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.sulc.fwhm25.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/rh.area.pial.fwhm15.fsaverage.mgh       /old_nas/subjects/bioface_0001/surf/rh.white.H.fwhm0.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.sulc.fwhm5.fsaverage.mgh             /old_nas/subjects/bioface_0001/surf/rh.area.pial.fwhm20.fsaverage.mgh       /old_nas/subjects/bioface_0001/surf/rh.white.H.fwhm10.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.thickness.fwhm0.fsaverage.mgh        /old_nas/subjects/bioface_0001/surf/rh.area.pial.fwhm25.fsaverage.mgh       /old_nas/subjects/bioface_0001/surf/rh.white.H.fwhm15.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.thickness.fwhm10.fsaverage.mgh       /old_nas/subjects/bioface_0001/surf/rh.area.pial.fwhm5.fsaverage.mgh        /old_nas/subjects/bioface_0001/surf/rh.white.H.fwhm20.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.thickness.fwhm15.fsaverage.mgh       /old_nas/subjects/bioface_0001/surf/rh.curv.fwhm0.fsaverage.mgh             /old_nas/subjects/bioface_0001/surf/rh.white.H.fwhm25.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.thickness.fwhm20.fsaverage.mgh       /old_nas/subjects/bioface_0001/surf/rh.curv.fwhm10.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/rh.white.H.fwhm5.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.thickness.fwhm25.fsaverage.mgh       /old_nas/subjects/bioface_0001/surf/rh.curv.fwhm15.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/rh.white.K.fwhm0.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.thickness.fwhm5.fsaverage.mgh        /old_nas/subjects/bioface_0001/surf/rh.curv.fwhm20.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/rh.white.K.fwhm10.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.volume.fwhm0.fsaverage.mgh           /old_nas/subjects/bioface_0001/surf/rh.curv.fwhm25.fsaverage.mgh            /old_nas/subjects/bioface_0001/surf/rh.white.K.fwhm15.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.volume.fwhm10.fsaverage.mgh          /old_nas/subjects/bioface_0001/surf/rh.curv.fwhm5.fsaverage.mgh             /old_nas/subjects/bioface_0001/surf/rh.white.K.fwhm20.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.volume.fwhm15.fsaverage.mgh          /old_nas/subjects/bioface_0001/surf/rh.jacobian_white.fwhm0.fsaverage.mgh   /old_nas/subjects/bioface_0001/surf/rh.white.K.fwhm25.fsaverage.mgh 
 +/old_nas/subjects/bioface_0001/surf/lh.volume.fwhm20.fsaverage.mgh          /old_nas/subjects/bioface_0001/surf/rh.jacobian_white.fwhm10.fsaverage.mgh  /old_nas/subjects/bioface_0001/surf/rh.white.K.fwhm5.fsaverage.mgh 
 +</code> 
 + 
 +Ver https://github.com/asqwerty666/acenip/blob/main/bin/pqcache.pl 
  
 ===== DTI preproc ===== ===== DTI preproc =====
Line 1971: Line 2086:
 [osotolongo@brick03 facehbi]$ for x in `cat xnat_pet_experiments.list`; do xnatapic run_pipeline --project_id facehbi --pipeline RegisterPETwithMRImatch --experiment_id ${x}; done [osotolongo@brick03 facehbi]$ for x in `cat xnat_pet_experiments.list`; do xnatapic run_pipeline --project_id facehbi --pipeline RegisterPETwithMRImatch --experiment_id ${x}; done
 </code> </code>
 +
 +==== Extraer resultados ====
 +
 +La idea fundamental aqui es obtener los valores de SUVR y centiloide calculados en XNAT. Estos pueden obtenerse con,
 +
 +<code bash>
 +$ xnatapic get_registration_report --project_id f5cehbi --output xnat_pet_results.csv
 +</code>
 +
 +y esto vamos a limpiarlo un poco,
 +
 +<code bash>
 +$ awk -F"," '{print $2","$6","$7}' xnat_pet_results.csv | sed 's/"//g' | tail -n +2 | sort -t, -k 1 > pet.results
 +</code>
 +
 +La lista de sujetos del proyecto se obtiene con,
 +
 +<code bash>
 +$  xnatapic list_subjects --project_id f5cehbi --label > xnat_subjects.list
 +</code>
 +
 +Esto se puede ordenar rapidamente haciendo,
 +
 +<code bash>
 +$ sort -t, -k 1 xnat_subjects.list > xnat_subjects_sorted.list
 +</code>
 +
 +Y entonces podemos integrar los ID de proyecto con los resultaos, a traves de los ID de XNAT, y seleccionar los datos que queremos guardar,
 +
 +<code bash>
 +$ join -t, xnat_subjects_sorted.list pet.results | sort -t, -k 2 | awk -F"," '{if ($3) print $2";"$3";"$4}' | sed '1iSubject;SUVR;Centiloid'"
 +</code>
 +
 +Todo esto se resume en el script [[https://github.com/asqwerty666/acenip/blob/main/bin/xnat_pullcl.pl | xnat_pullcl.pl]]. Los archivos intermedios se borran una vez que se obtienen los resultados. Las opciones son,
 +
 +  * -p : nombre del proyecto en XNAT, se ignora si existe //-x//. Solo existe por compatibilidad.
 +  * -x : nombre del proyecto en XNAT
 +  * -o : archivo de salida (default STDOUT)
 +
  
  
Line 2167: Line 2321:
 ---- ----
  
-==== Extraer resultados ==== 
- 
-La idea fundamental aqui es obtener los valores de SUVR y centiloide calculados en XNAT. Estos pueden obtenerse con, 
- 
-<code bash> 
-$ xnatapic get_registration_report --project_id f5cehbi --output xnat_pet_results.csv 
-</code> 
- 
-y esto vamos a limpiarlo un poco, 
- 
-<code bash> 
-$ awk -F"," '{print $2","$6","$7}' xnat_pet_results.csv | sed 's/"//g' | tail -n +2 | sort -t, -k 1 > pet.results 
-</code> 
- 
-La lista de sujetos del proyecto se obtiene con, 
- 
-<code bash> 
-$  xnatapic list_subjects --project_id f5cehbi --label > xnat_subjects.list 
-</code> 
- 
-Esto se puede ordenar rapidamente haciendo, 
- 
-<code bash> 
-$ sort -t, -k 1 xnat_subjects.list > xnat_subjects_sorted.list 
-</code> 
- 
-Y entonces podemos integrar los ID de proyecto con los resultaos, a traves de los ID de XNAT, y seleccionar los datos que queremos guardar, 
- 
-<code bash> 
-$ join -t, xnat_subjects_sorted.list pet.results | sort -t, -k 2 | awk -F"," '{if ($3) print $2";"$3";"$4}' | sed '1iSubject;SUVR;Centiloid'" 
-</code> 
- 
-Todo esto se resume en el script [[https://github.com/asqwerty666/acenip/blob/main/bin/xnat_pullcl.pl | xnat_pullcl.pl]]. Los archivos intermedios se borran una vez que se obtienen los resultados. Las opciones son, 
- 
-  * -p : nombre del proyecto en XNAT, se ignora si existe //-x//. Solo existe por compatibilidad. 
-  * -x : nombre del proyecto en XNAT 
-  * -o : archivo de salida (default STDOUT) 
  
 ===== PET-tau ===== ===== PET-tau =====
Line 2264: Line 2381:
 </code> </code>
  
-y para cada experiento gurdamos los resultados del analisis de FS en un directorio,+y para cada experiento guardamos los resultados del analisis de FS en un directorio,
  
 <code perl> <code perl>
Line 2280: Line 2397:
 my $fsout = $fsoutput.'/fsrecon'; my $fsout = $fsoutput.'/fsrecon';
 make_path $fsout; make_path $fsout;
-</code perl>+</code>
  
 podemos utilizar las sintaxis de la funcion //fs_file_metrics()// almacenadas en [[neuroimagen:neuro4.pm#fsmetricspm|FSMetrics.pm]]. //L ultima version actualizada, siempre [[https://github.com/asqwerty666/acenip/blob/main/lib/FSMetrics.pm|en github]]//. podemos utilizar las sintaxis de la funcion //fs_file_metrics()// almacenadas en [[neuroimagen:neuro4.pm#fsmetricspm|FSMetrics.pm]]. //L ultima version actualizada, siempre [[https://github.com/asqwerty666/acenip/blob/main/lib/FSMetrics.pm|en github]]//.
neuroimagen/pipe05.1641811551.txt.gz · Last modified: 2022/01/10 10:45 by osotolongo