projects / BK-2020-03.git / blame_incremental
| Commit | Line | Data |
|---|---|---|
| 1 | #!/usr/bin/env python3 | |
| 2 | ||
| 3 | # Desc: Expand or reduce list of int/floats to specified length. | |
| 4 | ||
| 5 | import numpy as np | |
| 6 | ||
| 7 | def medianSubset(listIn: list = [], listOutLen: int = 0) -> list: | |
| 8 | # Input: int: listOutLen: quantity of elements in output list | |
| 9 | # list: listIn: input list consisting of integers or floats | |
| 10 | # Output: list: ints/floats of specified size | |
| 11 | # Ref/Attrib: PEP 3107 typing https://stackoverflow.com/a/21384492 | |
| 12 | # Version: 0.1.0 | |
| 13 | #print('DEBUG:listOutLen:' + str(listOutLen)) | |
| 14 | #print('DEBUG:listIn:' + str(listIn)) | |
| 15 | ||
| 16 | # Exit for invalid input | |
| 17 | if not isinstance(listOutLen, int): | |
| 18 | raise ValueError('ERROR:Not a valid int:' + str(listOutLen)) | |
| 19 | else: | |
| 20 | if not listOutLen > 0: | |
| 21 | raise ValueError('ERROR:Invalid value:' + str(listOutLen)) | |
| 22 | if not isinstance(listIn, list): | |
| 23 | raise ValueError('ERROR:Not a valid list:' + str(listOutLen)) | |
| 24 | if not all([( (isinstance(x,int)) or (isinstance(x,float)) ) for x in listIn]): | |
| 25 | raise ValueError('ERROR:Input list contains something besides integers or floating point numbers.') | |
| 26 | ||
| 27 | # Initialize listOut | |
| 28 | listOut = [None] * listOutLen | |
| 29 | #print('DEBUG:listOut:' + str(listOut)) | |
| 30 | ||
| 31 | # Calc listIn length | |
| 32 | listInLen = len(listIn) | |
| 33 | #print('DEBUG:listInLen:' + str(listInLen)) | |
| 34 | ||
| 35 | # Calc subset length float | |
| 36 | subsetLenFloat = ( (max([listInLen,listOutLen]) - 1) /min([listInLen,listOutLen])) | |
| 37 | subsetIndRatio = ( (listInLen)/(listOutLen) ) | |
| 38 | #print('DEBUG:subsetLenFloat: %.5f' % subsetLenFloat) | |
| 39 | #print('DEBUG:subsetLenFloat2: %.5f' % subsetIndRatio) | |
| 40 | ||
| 41 | # Iterate for each element in listOut | |
| 42 | for i_out in range(listOutLen): | |
| 43 | #print('DEBUG:i_out:' + str(i_out)) | |
| 44 | ## Decide to expand or reduce listIn to produce listOut | |
| 45 | if listInLen > listOutLen: | |
| 46 | ### reduce listIn to listOut | |
| 47 | #print('DEBUG:listOutLen:' + str(listOutLen)) | |
| 48 | #print('DEBUG:listInLen:' + str(listInLen)) | |
| 49 | if i_out == 0: | |
| 50 | #### Initialize subsetIndLo in first loop | |
| 51 | subsetIndLo = int(0) | |
| 52 | #print('DEBUG:subsetIndLo:' + str(subsetIndLo)) | |
| 53 | #print('DEBUG:i_out:' + str(i_out)) | |
| 54 | #### Calc indices of i_out'th subset of listIn | |
| 55 | subsetIndHi = (listInLen - 1) * (i_out + 1) // listOutLen | |
| 56 | subsetLen = subsetIndHi - subsetIndLo + 1 | |
| 57 | #print('DEBUG:subsetIndLo:' + str(subsetIndLo)) | |
| 58 | #print('DEBUG:subsetIndHi:' + str(subsetIndHi)) | |
| 59 | #print('DEBUG:subsetLen:' + str(subsetLen)) | |
| 60 | #### Extract subset from listIn using indices inclusively | |
| 61 | subset = listIn[ int(subsetIndLo) : int(subsetIndHi)+1 ] | |
| 62 | #print('DEBUG:subset:' + str(subset)) | |
| 63 | #### Calculate median for subset | |
| 64 | subsetMedian = np.median(subset) | |
| 65 | #print('DEBUG:subset median:' + str(subsetMedian)) | |
| 66 | #### Set listOut element | |
| 67 | listOut[i_out] = subsetMedian | |
| 68 | #### Housekeeping | |
| 69 | ##### Update subsetIndLo for next loop | |
| 70 | subsetIndLo = subsetIndHi + 1 | |
| 71 | #print('DEBUG:Updated subsetIndLo:' + str(subsetIndLo)) | |
| 72 | elif listOutLen > listInLen: | |
| 73 | ### Expand listIn to listOut | |
| 74 | #print('DEBUG:listOutLen:' + str(listOutLen)) | |
| 75 | #print('DEBUG:listInLen:' + str(listInLen)) | |
| 76 | #### Identify index list of lists mapping listIn to ListOut | |
| 77 | expandIndex = int(i_out / subsetLenFloat) | |
| 78 | expandIndex = min([expandIndex,(listInLen - 1)]) | |
| 79 | #print('DEBUG:expandIndex:' + str(expandIndex)) | |
| 80 | listOut[i_out] = listIn[expandIndex] | |
| 81 | #print('DEBUG:listOut[i_out]:' + str(listOut[i_out])) | |
| 82 | elif listOutLen == listInLen: | |
| 83 | listOut = listIn | |
| 84 | #print('DEBUG:end for loop===========') | |
| 85 | return listOut | |
| 86 | ||
| 87 | example_list1 = [738,253,104,543,501,352,770,945,992,920,977,435,783,592,261,176,888] | |
| 88 | example_list2 = [641,679,371] | |
| 89 | example_list3 = [738,253,104,543,501,352,770,945,992,920,977,435,783,592,261,176,888,767,169,495,382,291,116,480,8,52,761,455,743,592,54,490,596,8,816,134,641,705,2,675,414,511,220,60,116,563,19,316,605,735,196,210,107,848,630,392,700,94,824,329,307,693,374,511,17,488,626,730,606,273,875,466,790,40,681,954,947,981,775,498,49,327,268,492,799,273,403,839,895,200,329,373,819,302,456,23,944,543,997,605,338,919,385,783,339,825,457,761,994,530,862,504,690,720,10,102,662,579,835,747,337,471,409,904,321,928,908,573,984,853,244,499,940,468,480,668,128,570,584,325,940,230,738,624,672,601,384,242,800,172,125,150,124,958,242,560,829,22,614,707,677,413,800,536,269,748,209,821,212,939,196,2,830,815,235,847,456,181,920,413,153,935,626,215,821,813,847,875,550,656,237,760,718,179,549,63,839,204,631,917,95,733,341,881,539,561,227,492,549,423,565,856,384,376,809,920,116,609,205,644,791,187,15,55,776,376,882,905,937,304,626,312,263,16,163,651,422,656,116,495,158,156,981,167,786,536,987,697,88,354,335,281,718,562,169,783,483,16,314,344,760,189,965,494,68,314,806,26,866,741,283,302,761,518,713,264,570,239,410,629,47,547,270,94,228,303,811,795,658,102,96,438,157,62,826,27,382,847,341,953,671,147,746,408,201,536,809,235,778,41,266,682,686,620,279,167,183,533,471,173,145,676,774,191,885,583,694,1000,484,115,159,31,704,248,852,441,901,705,829,496,784,163,88,1,503,986,929,631,839,124,55,136,907,342,863,328,417,674,881,567,876,30,812,515,606,469,725,478,100,149,881,12,221,839,138,354,712,550,698,931,914,610,704,513,477,527,388,341,373,77,673,970,612,384,613,321,411,382,273,100,55,760,246,161,268,348,22,533,851,515,807,76,396,254,554,503,322,292,821,345,147,303,225,389,356,119,530,502,407,474,824,900,561,624,292,893,597,60,393,856,675,27,257,535,5,879,809,944,904,327,209,626,440,102,528,970,839,77,0,987,422,414,337,499,266,640,194,454,414,517,221,83,4,473,84,39,831,640,557,685,45,91,410,327,734,470,928,22,837,583,570,671,773,94,317,935,365,319,965,465,723,471,77,420,683,581,604,17,373,767,722,15,742,533,497,847,8,84,842,28,526,591,466,396,807,676,564,695,967,47,886,266,638,767,58,763,775,778,25,699,895,64,616,11,369,521,856,277,726,165,281,692,736,256,545,912,954,325,620,220,407,304,191,363,884,818,138,783,559,164,516,880,791,884,246,963,522,547,356,475,302,582,687,399,862,109,524,570,867,961,780,611,491,677,977,879,803,651,41,400,634,415,750,129,89,861,265,669,833,833,333,221,102,505,353,264,31,782,595,838,365,141,475,697,46,503,71,205,164,811,46,57,278,278,12,297,33,326,106,916,352,528,616,21,129,534,485,753,410,952,758,26,152,123,578,42,33,881,145,71,662,902,782,634,256,598,265,953,135,778,407,990,643,309,448,381,620,559,208,648,830,332,59,420,456,558,439,503,182,728,486,568,986,172,652,176,380,542,92,551,335,602,854,485,543,752,860,163,120,607,207,710,298,711,32,908,228,41,481,313,581,874,443,273,204,265,567,792,326,173,29,762,398,923,18,312,144,965,764,342,884,395,624,853,709,7,424,11,836,677,210,717,950,525,782,39,770,972,993,517,495,323,772,133,941,300,809,517,574,243,664,772,721,232,304,429,864,135,667,871,371,894,485,86,42,445,426,556,588,205,533,392,900,962,112,837,961,617,302,779,793,969,199,199,591,747,321,453,861,145,451,149,887,726,414,696,421,260,92,159,56,528,879,844,408,568,603,182,492,796,656,876,566,468,237,400,629,366,868,145,150,151,561,738,241,891,813,235,400,476,224,241,169,157,969,540,901,925,723,536,421,369,24,560,882,817,677,234,16,241,464,923,849,300,144,261,916,529,586,905,268,950,673,697,664,840,725,926,603,96,845,410,321,9,996,844,952,380,272,125,140,719,402,383,499,9,48,955,429,458,965,584,381,739,984,10,282,165,486,33,353,616,304,215,293,362,148,542,815,659,874,890,502,369,499,533,144,876,169,375,951,784,320,782,763,473,38,871,790,933,522,88,702,387,705,544,575,719,757,867,470,807,628,914,719,764,536,748,766,419,42,604,678,742,321,280,147,975,918,982,440,842,723,910,893,111,560,922,837,116,858,374,262,878,354,397,638,338,853] | |
| 90 | example_list4 = [641,679,371,3.141592] | |
| 91 | example_list5 = [641,679,371,3.141592, 'jan'] | |
| 92 | #print(medianSubset(example_list,3)) | |
| 93 | #print(medianSubset(example_list2,17)) | |
| 94 | #print(medianSubset(example_list3,3)) | |
| 95 | #print(medianSubset(example_list3,1)) | |
| 96 | #print(medianSubset(example_list4,1)) | |
| 97 | print(medianSubset(example_list3,47)) |