studen/dynamicSPECT @ 1440d682a0bea92b2824b7807a9bae8e4d4dac57

table.ipynb12 KB
					
					文件历史
					原始文件
				
import sys
import os
import SimpleITK
import numpy
import matplotlib.pyplot
sys.path.append(os.path.join(os.path.expanduser('~'),'software','src','nixSuite','wrapper'))
import nixWrapper
nixWrapper.loadLibrary('labkeyInterface')
import labkeyInterface
net=labkeyInterface.labkeyInterface()
fconfig=os.path.join(os.path.expanduser('~'),'.labkey','network.json')
net.init(fconfig)
net.getCSRF()
import labkeyFileBrowser
sys.path.append(os.getcwd())
import downloadPatient
fb=labkeyFileBrowser.labkeyFileBrowser(net)
project='dinamic_spect/Patients'
locDir=os.path.join(os.path.expanduser('~'),'temp','dynamicSPECT')
import labkeyDatabaseBrowser
db=labkeyDatabaseBrowser.labkeyDB(net)

/usr/lib/python3/dist-packages/urllib3/connection.py:391: SubjectAltNameWarning: Certificate for merlin.fmf.uni-lj.si has no `subjectAltName`, falling back to check for a `commonName` for now. This feature is being removed by major browsers and deprecated by RFC 2818. (See https://github.com/urllib3/urllib3/issues/497 for details.)
  warnings.warn(
User: andrej studen CSRF: 4d08323bbd1605f86108a60a9d40db20
def averageK1(db,regionId):
    regionSelect={'variable':'regionId','value':str(regionId),'oper':'eq'}
    classSelect={'variable':'nclass','value':str(c),'oper':'eq'}
    ds=db.selectRows(project,'study','ImagingData',[classSelect,regionSelect])
    patients=list(set([r['PatientId'] for r in ds['rows']]))

    k1MAll=[]
    k1OAll=[]
    for p in patients:
        k1M=[]
        k1O=[]

        patientFilter={'variable':'PatientId','value':p,'oper':'eq'}
        ds=db.selectRows(project,'study','ImagingData',[patientFilter,regionSelect])
        for r in ds['rows']:
            k1M.append(r['k1MIRMedian'])
            k1O.append(r['k1OBRMedian'])
        vM=numpy.median(k1M)
        #print('{} {}'.format(vM,k1M))
        #print('{} {}'.format(numpy.median(k1M),numpy.std(k1M)))
        vO=numpy.median(k1O)
        print('[{}/{}] {} {}'.format(p,regionId,vO,k1O))
        #print('{} {}'.format(numpy.median(k1O),numpy.std(k1O)))
        k1MAll.append(vM)
        k1OAll.append(vO)
    return k1MAll,k1OAll
    
f = open("/home/studen/temp/dynamicSPECT/table.txt", "w")
for r in numpy.arange(1,10):
    k1MAll,k1OAll=averageK1(db,r)
    print(k1MAll)
    print(k1OAll)
    
    f.write('{:.2f},{:.2f},{:.2f},{:.2f},{:.2f},{:.2f}\n'.format(60*numpy.mean(k1MAll),60*numpy.median(k1MAll),60*numpy.std(k1MAll), 60*numpy.mean(k1OAll), 60*numpy.median(k1OAll),60*numpy.std(k1OAll)))
f.close()
[7TM/1] 0.008045887022546614 [0.008045887022546614, 0.004483481665888348, 2.1286405607228596]
[10KF/1] 0.0039462171230090675 [0.0039462171230090675, 0.00415070155114901, 0.0031446075367390323]
[11ZM/1] 0.0025240576597139624 [0.00312436925868039, 0.0025240576597139624, 0.002142209448394668]
[2SB/1] 0.010642740867635697 [0.03745862838837244, 0.010642740867635697, 0.010028953831119597]
[5MI/1] 0.00888906341063412 [0.008304194640398289, 0.00888906341063412, 0.0091936110897083]
[3ZF/1] 0.26867647983040466 [0.38591675611767157, 0.26867647983040466, 0.17709045381913846]
[0.005894665958762989, 0.004668591303813753, 0.002949742733330378, 0.011751648013763579, 0.008162106616002189, 0.00511268572208257]
[0.008045887022546614, 0.0039462171230090675, 0.0025240576597139624, 0.010642740867635697, 0.00888906341063412, 0.26867647983040466]
[7TM/2] 0.0017372501091852122 [0.0017372501091852122, 0.001675933130699284, 0.19584171785850896]
[10KF/2] 8.259213311648866e-05 [9.282359403286088e-06, 0.002221889457124819, 8.259213311648866e-05]
[11ZM/2] 0.0021371243829022488 [0.0021470477401390694, 0.0021371243829022488, 0.0017499690427985907]
[2SB/2] 0.003306368767070012 [0.003866263540421745, 0.003306368767070012, 0.0029943594109714824]
[5MI/2] 0.005177206176626878 [0.005017774180493417, 0.005601704036412694, 0.005177206176626878]
[3ZF/2] 0.004304534377883605 [0.003885393241151492, 0.005233206242833382, 0.004304534377883605]
[0.0032749453662799595, 0.0028943546543249484, 0.001871871695075595, 0.003855121335050874, 0.004535002996260104, 0.002431672607075657]
[0.0017372501091852122, 8.259213311648866e-05, 0.0021371243829022488, 0.003306368767070012, 0.005177206176626878, 0.004304534377883605]
[7TM/3] 0.0024026106630740326 [0.0024026106630740326, 0.0020285709174733724, 0.7314645007711776]
[10KF/3] 0.0015306808132398398 [0.0018703063484711772, 0.0015306808132398398, 0.0012136425939998793]
[11ZM/3] 0.003917784910452867 [0.003917784910452867, 0.004195671684134263, 0.0030221084430817322]
[2SB/3] 0.0018769599381238148 [0.0018769599381238148, 0.0019454598760022081, 0.0017322417901848254]
[5MI/3] 0.004714936487121288 [0.005021416534152844, 0.00453359720801415, 0.004714936487121288]
[3ZF/3] 0.005698838193463921 [0.005698838193463921, 0.006372350188361305, 0.004733389017227235]
[0.003006674656825954, 0.0024901765401445894, 0.002640441247904802, 0.002556852090616107, 0.003922843192938324, 0.0012724135802516157]
[0.0024026106630740326, 0.0015306808132398398, 0.003917784910452867, 0.0018769599381238148, 0.004714936487121288, 0.005698838193463921]
[7TM/4] 0.004195682980252796 [0.003950188996175878, 0.004195682980252796, 2.0863785597545665]
[10KF/4] 0.005017210604686076 [0.0058695581350023925, 0.005017210604686076, 0.00382322050790849]
[11ZM/4] 0.004373077263385557 [0.005304248942517635, 0.004373077263385557, 0.0038680517245407303]
[2SB/4] 0.004780494766631807 [0.004452715331717112, 0.004836906908797207, 0.004780494766631807]
[5MI/4] 0.004522149143858034 [0.004378097716940858, 0.004730024546274227, 0.004522149143858034]
[3ZF/4] 0.5730251968275395 [0.5730251968275395, 0.6151303069429822, 0.17997801559145749]
[0.004801354184919535, 0.001981110883473747, 0.00311768616798197, 0.002664618039935267, 0.008329222245675397, 0.002445709148188661]
[0.004195682980252796, 0.005017210604686076, 0.004373077263385557, 0.004780494766631807, 0.004522149143858034, 0.5730251968275395]
[7TM/5] 0.0030026127913968777 [0.0030026127913968777, 0.002966314233129647, 1.0919607438562668]
[10KF/5] 0.0018728630968653172 [0.0018728630968653172, 0.00254711409860054, 0.0017867094196874452]
[11ZM/5] 0.001744985966157502 [0.0020969180529696483, 0.001744985966157502, 0.0016118588408116715]
[2SB/5] 0.004780494766631807 [0.004452715331717112, 0.0048586456286596755, 0.004780494766631807]
[5MI/5] 0.006707032545252601 [0.006170975452619039, 0.006707032545252601, 0.007460460364806493]
[3ZF/5] 0.00858308393943204 [0.00858308393943204, 0.011845597833970808, 0.007380003330274289]
[0.004055782128720719, 0.00328245378291892, 0.0014083849974710172, 0.0029858684539020392, 0.007845910485130312, 0.002445709148188661]
[0.0030026127913968777, 0.0018728630968653172, 0.001744985966157502, 0.004780494766631807, 0.006707032545252601, 0.00858308393943204]
[7TM/6] 0.002837339604720477 [0.002837339604720477, 0.002449581166322905, 0.41387917879987024]
[10KF/6] 8.259213311648866e-05 [9.282359403286088e-06, 0.002221889457124819, 8.259213311648866e-05]
[11ZM/6] 0.0012243542725250457 [0.0014820056066890714, 0.0012243542725250457, 0.0010976608964243597]
[2SB/6] 0.0023446856820961777 [0.0024867978271161787, 0.0023446856820961777, 0.0022332197164982065]
[5MI/6] 0.0031172956941653137 [0.0029513503451490193, 0.0033183656384836225, 0.0031172956941653137]
[3ZF/6] 0.003781674236994135 [0.003781674236994135, 0.004071966210877487, 0.0034494591909788983]
[0.003138866692139931, 0.0025325379278861595, 0.0009215932677006866, 0.002621400221835252, 0.001954550346671144, 0.0015189554927926444]
[0.002837339604720477, 8.259213311648866e-05, 0.0012243542725250457, 0.0023446856820961777, 0.0031172956941653137, 0.003781674236994135]
[7TM/7] 0.0023828627975731745 [0.0023828627975731745, 0.0017047654247853699, 0.18591699611163967]
[10KF/7] 0.0019421932362865028 [0.0019421932362865028, 0.002120856579389987, 0.0018114613219635]
[11ZM/7] 0.0014335428486863065 [0.0014432654864380337, 0.0014335428486863065, 0.0014000118489706172]
[2SB/7] 0.001781126021641757 [0.0024867978271161787, 0.001781126021641757, 0.0016906990021623499]
[5MI/7] 0.0037122197226202246 [0.0037122197226202246, 0.003776907003804963, 0.003495911908165303]
[3ZF/7] 0.007206213864647251 [0.007206213864647251, 0.007899440203861475, 0.006211245602005458]
[0.0036923065171510427, 0.0020527140377196566, 0.0017667001699923747, 0.0015729443133339796, 0.0019636983396832306, 0.0013583430415635878]
[0.0023828627975731745, 0.0019421932362865028, 0.0014335428486863065, 0.001781126021641757, 0.0037122197226202246, 0.007206213864647251]
[7TM/8] 0.0030026127913968777 [0.0030026127913968777, 0.002912860211894843, 0.8457710628002826]
[10KF/8] 0.003937374142034068 [0.004542457733622267, 0.003937374142034068, 0.0034600892248924293]
[11ZM/8] 0.004351273874444507 [0.0050625173773602185, 0.004351273874444507, 0.0035935856109804007]
[2SB/8] 0.004326823649970777 [0.004452715331717112, 0.004325268011938922, 0.004326823649970777]
[5MI/8] 0.004743434827195546 [0.004743434827195546, 0.004492161915006733, 0.004857684829972067]
[3ZF/8] 0.01592586039329981 [0.015304785485526388, 0.09280392746997639, 0.01592586039329981]
[0.0032756892467843577, 0.0020949474383190852, 0.0028992798255622135, 0.0027488460983543573, 0.00416371874346105, 0.0022741875289843377]
[0.0030026127913968777, 0.003937374142034068, 0.004351273874444507, 0.004326823649970777, 0.004743434827195546, 0.01592586039329981]
[7TM/9] 0.0017372501091852122 [0.0017372501091852122, 0.0016790618998124395, 0.05339344195537287]
[10KF/9] 0.0011993853157435785 [0.0011370496488211962, 0.0013632774988599744, 0.0011993853157435785]
[11ZM/9] 0.0010565458881223932 [0.0009607852127521912, 0.0011070435774301734, 0.0010565458881223932]
[2SB/9] 0.0017336108126774022 [0.0018769599381238148, 0.0017336108126774022, 0.0015180304828638965]
[5MI/9] 0.0022169750747788675 [0.00278484013499442, 0.0022169750747788675, 0.001833834797518541]
[3ZF/9] 0.005896802413281909 [0.005896802413281909, 0.007601946240230148, 0.00578983776202438]
[0.0016943661978618529, 0.0014651121228743532, 0.0011333849162417412, 0.002449166151660834, 0.0018263583349822485, 0.001067181006343166]
[0.0017372501091852122, 0.0011993853157435785, 0.0010565458881223932, 0.0017336108126774022, 0.0022169750747788675, 0.005896802413281909]
table.ipynb 12 KB 文件历史 原始文件

table.ipynb 12 KB

文件历史原始文件
