Laser Sensing Display for UI interfaces in the real world
Fork of skinGames_forktest by
hardwareIO/hardwareIO.cpp@24:4e52031a495b, 2012-06-18 (annotated)
- Committer:
- mbedalvaro
- Date:
- Mon Jun 18 08:00:57 2012 +0000
- Revision:
- 24:4e52031a495b
- Parent:
- 23:bf666fcc61bc
Various modifications for show at TV Asahi
Who changed what in which revision?
User | Revision | Line number | New contents of line |
---|---|---|---|
mbedalvaro | 0:345b3bc7a0ea | 1 | #include "hardwareIO.h" |
mbedalvaro | 0:345b3bc7a0ea | 2 | |
mbedalvaro | 0:345b3bc7a0ea | 3 | HardwareIO IO; // preintantiation of cross-file global object IO |
mbedalvaro | 0:345b3bc7a0ea | 4 | |
mbedalvaro | 0:345b3bc7a0ea | 5 | // -------------------------------------- (0) SETUP ALL IO (call this in the setup() function in main program) |
mbedalvaro | 0:345b3bc7a0ea | 6 | |
mbedalvaro | 0:345b3bc7a0ea | 7 | Serial pc(USBTX, USBRX); // tx, rx |
mbedalvaro | 0:345b3bc7a0ea | 8 | LocalFileSystem local("local"); // Create the local filesystem under the name "local" |
mbedalvaro | 0:345b3bc7a0ea | 9 | |
mbedalvaro | 0:345b3bc7a0ea | 10 | SPI spiDAC(MOSI_PIN, MISO_PIN, SCK_PIN); // mosi, miso, sclk |
mbedalvaro | 0:345b3bc7a0ea | 11 | DigitalOut csDAC(CS_DAC_MIRRORS); |
mbedalvaro | 0:345b3bc7a0ea | 12 | |
mbedalvaro | 23:bf666fcc61bc | 13 | AnalogIn ain(POT_ANALOG_INPUT); |
mbedalvaro | 23:bf666fcc61bc | 14 | |
mbedalvaro | 23:bf666fcc61bc | 15 | DigitalOut Laser_Red(LASER_RED_PIN); // NOTE: this is NOT the lock in sensing laser (actually, not used yet) |
mbedalvaro | 0:345b3bc7a0ea | 16 | DigitalOut Laser_Green(LASER_GREEN_PIN); |
mbedalvaro | 0:345b3bc7a0ea | 17 | DigitalOut Laser_Blue(LASER_BLUE_PIN); |
mbedalvaro | 0:345b3bc7a0ea | 18 | |
mbedalvaro | 0:345b3bc7a0ea | 19 | void HardwareIO::init(void) { |
mbedalvaro | 23:bf666fcc61bc | 20 | Laser_Red = 0; // note: this is not the lockin-laser! |
mbedalvaro | 0:345b3bc7a0ea | 21 | Laser_Green = 0; |
mbedalvaro | 0:345b3bc7a0ea | 22 | Laser_Blue = 0; |
mbedalvaro | 0:345b3bc7a0ea | 23 | |
mbedalvaro | 0:345b3bc7a0ea | 24 | //Serial Communication setup: |
mbedalvaro | 14:0fc33a3a7b4b | 25 | pc.baud(115200);// |
mbedalvaro | 14:0fc33a3a7b4b | 26 | // pc.baud(921600);//115200);// |
mbedalvaro | 0:345b3bc7a0ea | 27 | |
mbedalvaro | 0:345b3bc7a0ea | 28 | // Setup for lock-in amplifier and pwm references: |
mbedalvaro | 0:345b3bc7a0ea | 29 | lockin.init(); |
mbedalvaro | 0:345b3bc7a0ea | 30 | |
mbedalvaro | 0:345b3bc7a0ea | 31 | // Setup the spi for 8 bit data, high steady state clock, |
mbedalvaro | 0:345b3bc7a0ea | 32 | // second edge capture, with a 10MHz clock rate |
mbedalvaro | 0:345b3bc7a0ea | 33 | csDAC = 1; |
mbedalvaro | 0:345b3bc7a0ea | 34 | spiDAC.format(16,0); |
mbedalvaro | 0:345b3bc7a0ea | 35 | spiDAC.frequency(16000000); |
mbedalvaro | 0:345b3bc7a0ea | 36 | |
mbedalvaro | 0:345b3bc7a0ea | 37 | // default initial mirror position: |
mbedalvaro | 0:345b3bc7a0ea | 38 | writeOutX(CENTER_AD_MIRROR_X); |
mbedalvaro | 0:345b3bc7a0ea | 39 | writeOutY(CENTER_AD_MIRROR_Y); |
mbedalvaro | 0:345b3bc7a0ea | 40 | |
mbedalvaro | 0:345b3bc7a0ea | 41 | // Load LUT table: |
mbedalvaro | 0:345b3bc7a0ea | 42 | setLUT(); |
mbedalvaro | 0:345b3bc7a0ea | 43 | } |
mbedalvaro | 0:345b3bc7a0ea | 44 | |
mbedalvaro | 0:345b3bc7a0ea | 45 | //write on the first DAC, output A (mirror X) |
mbedalvaro | 10:6f8e48dca1bd | 46 | void HardwareIO::writeOutX(unsigned short value){ |
mbedalvaro | 0:345b3bc7a0ea | 47 | if(value > MAX_AD_MIRRORS) value = MAX_AD_MIRRORS; |
mbedalvaro | 0:345b3bc7a0ea | 48 | if(value < MIN_AD_MIRRORS) value = MIN_AD_MIRRORS; |
mbedalvaro | 0:345b3bc7a0ea | 49 | |
mbedalvaro | 0:345b3bc7a0ea | 50 | value |= 0x7000; |
mbedalvaro | 0:345b3bc7a0ea | 51 | value &= 0x7FFF; |
mbedalvaro | 0:345b3bc7a0ea | 52 | |
mbedalvaro | 0:345b3bc7a0ea | 53 | csDAC = 0; |
mbedalvaro | 0:345b3bc7a0ea | 54 | spiDAC.write(value); |
mbedalvaro | 0:345b3bc7a0ea | 55 | csDAC = 1; |
mbedalvaro | 0:345b3bc7a0ea | 56 | } |
mbedalvaro | 0:345b3bc7a0ea | 57 | |
mbedalvaro | 0:345b3bc7a0ea | 58 | //write on the first DAC, output B (mirror Y) |
mbedalvaro | 10:6f8e48dca1bd | 59 | void HardwareIO::writeOutY(unsigned short value){ |
mbedalvaro | 0:345b3bc7a0ea | 60 | if(value > MAX_AD_MIRRORS) value = MAX_AD_MIRRORS; |
mbedalvaro | 0:345b3bc7a0ea | 61 | if(value < MIN_AD_MIRRORS) value = MIN_AD_MIRRORS; |
mbedalvaro | 0:345b3bc7a0ea | 62 | |
mbedalvaro | 0:345b3bc7a0ea | 63 | value |= 0xF000; |
mbedalvaro | 0:345b3bc7a0ea | 64 | value &= 0xFFFF; |
mbedalvaro | 0:345b3bc7a0ea | 65 | |
mbedalvaro | 0:345b3bc7a0ea | 66 | csDAC = 0; |
mbedalvaro | 0:345b3bc7a0ea | 67 | spiDAC.write(value); |
mbedalvaro | 0:345b3bc7a0ea | 68 | csDAC = 1; |
mbedalvaro | 0:345b3bc7a0ea | 69 | } |
mbedalvaro | 0:345b3bc7a0ea | 70 | |
mbedalvaro | 0:345b3bc7a0ea | 71 | void HardwareIO::setRedPower(int powerValue){ |
mbedalvaro | 0:345b3bc7a0ea | 72 | if(powerValue > 0){ |
mbedalvaro | 0:345b3bc7a0ea | 73 | lockin.setLaserPower(true); |
mbedalvaro | 0:345b3bc7a0ea | 74 | } |
mbedalvaro | 0:345b3bc7a0ea | 75 | else{ |
mbedalvaro | 0:345b3bc7a0ea | 76 | lockin.setLaserPower(false); |
mbedalvaro | 0:345b3bc7a0ea | 77 | } |
mbedalvaro | 0:345b3bc7a0ea | 78 | } |
mbedalvaro | 0:345b3bc7a0ea | 79 | void HardwareIO::setGreenPower(int powerValue){ |
mbedalvaro | 0:345b3bc7a0ea | 80 | if(powerValue > 0){ |
mbedalvaro | 0:345b3bc7a0ea | 81 | Laser_Green = 1; |
mbedalvaro | 0:345b3bc7a0ea | 82 | } |
mbedalvaro | 0:345b3bc7a0ea | 83 | else{ |
mbedalvaro | 0:345b3bc7a0ea | 84 | Laser_Green = 0; |
mbedalvaro | 0:345b3bc7a0ea | 85 | } |
mbedalvaro | 0:345b3bc7a0ea | 86 | } |
mbedalvaro | 0:345b3bc7a0ea | 87 | void HardwareIO::setBluePower(int powerValue){ |
mbedalvaro | 0:345b3bc7a0ea | 88 | if(powerValue > 0){ |
mbedalvaro | 0:345b3bc7a0ea | 89 | Laser_Blue = 1; |
mbedalvaro | 0:345b3bc7a0ea | 90 | } |
mbedalvaro | 0:345b3bc7a0ea | 91 | else{ |
mbedalvaro | 0:345b3bc7a0ea | 92 | Laser_Blue = 0; |
mbedalvaro | 0:345b3bc7a0ea | 93 | } |
mbedalvaro | 0:345b3bc7a0ea | 94 | } |
mbedalvaro | 0:345b3bc7a0ea | 95 | |
mbedalvaro | 23:bf666fcc61bc | 96 | void HardwareIO::setRGBPower(unsigned char color) { |
mbedalvaro | 23:bf666fcc61bc | 97 | lockin.setLaserPower(color&0x04>0? false : true); |
mbedalvaro | 0:345b3bc7a0ea | 98 | Laser_Green=(color&0x02)>>1; |
mbedalvaro | 23:bf666fcc61bc | 99 | Laser_Blue =color&0x01; |
mbedalvaro | 0:345b3bc7a0ea | 100 | } |
mbedalvaro | 0:345b3bc7a0ea | 101 | |
mbedalvaro | 22:d87317d7ca91 | 102 | void HardwareIO::showLimitsMirrors(int times) { |
mbedalvaro | 24:4e52031a495b | 103 | unsigned short pointsPerLine=150; |
mbedalvaro | 22:d87317d7ca91 | 104 | int shiftX = (MAX_AD_MIRRORS - MIN_AD_MIRRORS) / pointsPerLine; |
mbedalvaro | 22:d87317d7ca91 | 105 | int shiftY = (MAX_AD_MIRRORS - MIN_AD_MIRRORS) / pointsPerLine; |
mbedalvaro | 24:4e52031a495b | 106 | |
mbedalvaro | 24:4e52031a495b | 107 | Laser_Green=1; |
mbedalvaro | 24:4e52031a495b | 108 | |
mbedalvaro | 24:4e52031a495b | 109 | //for (int repeat=0; repeat<times; repeat++) { |
mbedalvaro | 24:4e52031a495b | 110 | |
mbedalvaro | 24:4e52031a495b | 111 | Timer t; |
mbedalvaro | 24:4e52031a495b | 112 | t.start(); |
mbedalvaro | 24:4e52031a495b | 113 | while(t.read_ms()<times*1000) { |
mbedalvaro | 22:d87317d7ca91 | 114 | |
mbedalvaro | 22:d87317d7ca91 | 115 | writeOutX(MIN_AD_MIRRORS);writeOutY(MIN_AD_MIRRORS); |
mbedalvaro | 24:4e52031a495b | 116 | |
mbedalvaro | 22:d87317d7ca91 | 117 | for(int j=0; j<pointsPerLine; j++){ |
mbedalvaro | 24:4e52031a495b | 118 | wait_us(200);//delay between each points |
mbedalvaro | 22:d87317d7ca91 | 119 | writeOutY(j*shiftY + MIN_AD_MIRRORS); |
mbedalvaro | 22:d87317d7ca91 | 120 | } |
mbedalvaro | 22:d87317d7ca91 | 121 | |
mbedalvaro | 22:d87317d7ca91 | 122 | writeOutX(MIN_AD_MIRRORS);writeOutY(MAX_AD_MIRRORS); |
mbedalvaro | 22:d87317d7ca91 | 123 | for(int j=0; j<pointsPerLine; j++) { |
mbedalvaro | 22:d87317d7ca91 | 124 | wait_us(200);//delay between each points |
mbedalvaro | 24:4e52031a495b | 125 | writeOutX(j*shiftX + MIN_AD_MIRRORS); |
mbedalvaro | 22:d87317d7ca91 | 126 | } |
mbedalvaro | 22:d87317d7ca91 | 127 | |
mbedalvaro | 22:d87317d7ca91 | 128 | writeOutX(MAX_AD_MIRRORS);writeOutY(MAX_AD_MIRRORS); |
mbedalvaro | 22:d87317d7ca91 | 129 | for(int j=0; j<pointsPerLine; j++) { |
mbedalvaro | 22:d87317d7ca91 | 130 | wait_us(200);//delay between each points |
mbedalvaro | 22:d87317d7ca91 | 131 | writeOutY(-j*shiftX + MAX_AD_MIRRORS); |
mbedalvaro | 22:d87317d7ca91 | 132 | } |
mbedalvaro | 22:d87317d7ca91 | 133 | |
mbedalvaro | 22:d87317d7ca91 | 134 | writeOutX(MAX_AD_MIRRORS);writeOutY(MIN_AD_MIRRORS); |
mbedalvaro | 22:d87317d7ca91 | 135 | for(int j=0; j<pointsPerLine; j++) { |
mbedalvaro | 22:d87317d7ca91 | 136 | wait_us(200);//delay between each points |
mbedalvaro | 22:d87317d7ca91 | 137 | writeOutX(-j*shiftX + MAX_AD_MIRRORS); |
mbedalvaro | 22:d87317d7ca91 | 138 | } |
mbedalvaro | 22:d87317d7ca91 | 139 | |
mbedalvaro | 22:d87317d7ca91 | 140 | } |
mbedalvaro | 24:4e52031a495b | 141 | t.stop(); |
mbedalvaro | 24:4e52031a495b | 142 | Laser_Green=0; |
mbedalvaro | 22:d87317d7ca91 | 143 | } |
mbedalvaro | 22:d87317d7ca91 | 144 | |
mbedalvaro | 10:6f8e48dca1bd | 145 | void HardwareIO::scan_serial(unsigned short pointsPerLine){ |
mbedalvaro | 0:345b3bc7a0ea | 146 | //scan the total surface with a custom resolution |
mbedalvaro | 0:345b3bc7a0ea | 147 | //send the lockin value for each point as a byte on the serial port to the PC |
mbedalvaro | 0:345b3bc7a0ea | 148 | //use "scanSLP_save" to see the data on processing |
mbedalvaro | 0:345b3bc7a0ea | 149 | |
mbedalvaro | 0:345b3bc7a0ea | 150 | |
mbedalvaro | 0:345b3bc7a0ea | 151 | int shiftX = (MAX_AD_MIRRORS - MIN_AD_MIRRORS) / pointsPerLine; |
mbedalvaro | 0:345b3bc7a0ea | 152 | int shiftY = (MAX_AD_MIRRORS - MIN_AD_MIRRORS) / pointsPerLine; |
mbedalvaro | 0:345b3bc7a0ea | 153 | |
mbedalvaro | 0:345b3bc7a0ea | 154 | for(int j=0; j<pointsPerLine; j++){ |
mbedalvaro | 0:345b3bc7a0ea | 155 | writeOutX(MIN_AD_MIRRORS); |
mbedalvaro | 0:345b3bc7a0ea | 156 | writeOutY(j*shiftY + MIN_AD_MIRRORS); |
mbedalvaro | 0:345b3bc7a0ea | 157 | |
mbedalvaro | 0:345b3bc7a0ea | 158 | wait_us(300);//begining of line delay |
mbedalvaro | 0:345b3bc7a0ea | 159 | for(int i=0; i<pointsPerLine; i++){ |
mbedalvaro | 0:345b3bc7a0ea | 160 | writeOutX(i*shiftX + MIN_AD_MIRRORS); |
mbedalvaro | 0:345b3bc7a0ea | 161 | |
mbedalvaro | 14:0fc33a3a7b4b | 162 | wait_us(200);//delay between each points |
mbedalvaro | 10:6f8e48dca1bd | 163 | |
mbedalvaro | 10:6f8e48dca1bd | 164 | // SEND A VALUE BETWEEN 0 and 255: |
mbedalvaro | 0:345b3bc7a0ea | 165 | pc.putc(int(255.0*lockin.getMedianValue()/4095));//printf("%dL",int(valueLockin*255));//pc.putc(int(lockin*255));// |
mbedalvaro | 0:345b3bc7a0ea | 166 | } |
mbedalvaro | 0:345b3bc7a0ea | 167 | } |
mbedalvaro | 0:345b3bc7a0ea | 168 | } |
mbedalvaro | 0:345b3bc7a0ea | 169 | |
mbedalvaro | 0:345b3bc7a0ea | 170 | //load Look-up Table from LUT.TXT file |
mbedalvaro | 0:345b3bc7a0ea | 171 | //or create the file with scanLUT() if not existing. |
mbedalvaro | 0:345b3bc7a0ea | 172 | void HardwareIO::setLUT(){ |
mbedalvaro | 0:345b3bc7a0ea | 173 | |
mbedalvaro | 0:345b3bc7a0ea | 174 | FILE *fp = fopen(LUT_FILENAME, "r"); // Open file on the local file system for writing |
mbedalvaro | 0:345b3bc7a0ea | 175 | if(fp){ |
mbedalvaro | 0:345b3bc7a0ea | 176 | //load the file into the lut table; keep the SAME resolution! |
mbedalvaro | 0:345b3bc7a0ea | 177 | fread(lut,sizeof(uint16),LUT_RESOLUTION*LUT_RESOLUTION,fp); |
mbedalvaro | 0:345b3bc7a0ea | 178 | fclose(fp); |
mbedalvaro | 0:345b3bc7a0ea | 179 | } |
mbedalvaro | 0:345b3bc7a0ea | 180 | else{ |
mbedalvaro | 0:345b3bc7a0ea | 181 | //fclose(fp); |
mbedalvaro | 0:345b3bc7a0ea | 182 | //if the file "LUT.TXT" doesn't exist, create one with scanLUT() |
mbedalvaro | 0:345b3bc7a0ea | 183 | lockin.setLaserPower(true); |
mbedalvaro | 0:345b3bc7a0ea | 184 | scanLUT(); |
mbedalvaro | 0:345b3bc7a0ea | 185 | } |
mbedalvaro | 0:345b3bc7a0ea | 186 | |
mbedalvaro | 0:345b3bc7a0ea | 187 | } |
mbedalvaro | 0:345b3bc7a0ea | 188 | |
mbedalvaro | 0:345b3bc7a0ea | 189 | //scan the total surface with a fixed 2^x resolution |
mbedalvaro | 0:345b3bc7a0ea | 190 | //create the Look-Up Table used to "flatten" the scan according to the position |
mbedalvaro | 0:345b3bc7a0ea | 191 | // |
mbedalvaro | 0:345b3bc7a0ea | 192 | //To Do: maybe detect high frequency to be sure the area is clean and empty? |
mbedalvaro | 0:345b3bc7a0ea | 193 | void HardwareIO::scanLUT(){ |
mbedalvaro | 0:345b3bc7a0ea | 194 | |
mbedalvaro | 0:345b3bc7a0ea | 195 | //reset lut table |
mbedalvaro | 0:345b3bc7a0ea | 196 | for(int j=0; j<LUT_RESOLUTION; j++){ |
mbedalvaro | 0:345b3bc7a0ea | 197 | for(int i=0; i<LUT_RESOLUTION; i++){ |
mbedalvaro | 0:345b3bc7a0ea | 198 | lut[i][j] =0; |
mbedalvaro | 0:345b3bc7a0ea | 199 | } |
mbedalvaro | 0:345b3bc7a0ea | 200 | } |
mbedalvaro | 0:345b3bc7a0ea | 201 | |
mbedalvaro | 0:345b3bc7a0ea | 202 | int delayScanning = 300; //in us |
mbedalvaro | 0:345b3bc7a0ea | 203 | |
mbedalvaro | 0:345b3bc7a0ea | 204 | //define the distance between each points (from 0 to 4096) and the offset (here 0) |
mbedalvaro | 0:345b3bc7a0ea | 205 | float shiftX = 1.0*(MAX_AD_MIRRORS - MIN_AD_MIRRORS) / (LUT_RESOLUTION-1); |
mbedalvaro | 0:345b3bc7a0ea | 206 | float shiftY = 1.0*(MAX_AD_MIRRORS - MIN_AD_MIRRORS) / (LUT_RESOLUTION-1); |
mbedalvaro | 0:345b3bc7a0ea | 207 | float offsetX = MIN_AD_MIRRORS; |
mbedalvaro | 0:345b3bc7a0ea | 208 | float offsetY = MIN_AD_MIRRORS; |
mbedalvaro | 0:345b3bc7a0ea | 209 | |
mbedalvaro | 0:345b3bc7a0ea | 210 | //move the mirrors to the first position |
mbedalvaro | 0:345b3bc7a0ea | 211 | writeOutX(MAX_AD_MIRRORS);writeOutY(MIN_AD_MIRRORS); |
mbedalvaro | 0:345b3bc7a0ea | 212 | wait_us(500); |
mbedalvaro | 0:345b3bc7a0ea | 213 | |
mbedalvaro | 0:345b3bc7a0ea | 214 | float x, y; |
mbedalvaro | 0:345b3bc7a0ea | 215 | |
mbedalvaro | 10:6f8e48dca1bd | 216 | //scan the surface NB_SCANS times |
mbedalvaro | 0:345b3bc7a0ea | 217 | //the total value in lut[i][j] shouldn't exceed uint16 !!! |
mbedalvaro | 10:6f8e48dca1bd | 218 | for(int loop=0; loop<NB_SCANS; loop++){ |
mbedalvaro | 0:345b3bc7a0ea | 219 | for(int j=0; j<LUT_RESOLUTION; j++){ |
mbedalvaro | 0:345b3bc7a0ea | 220 | y = shiftY*j + offsetY ; |
mbedalvaro | 0:345b3bc7a0ea | 221 | writeOutY(int(y)); |
mbedalvaro | 0:345b3bc7a0ea | 222 | //scan from right to left |
mbedalvaro | 0:345b3bc7a0ea | 223 | for(int i=LUT_RESOLUTION-1; i>=0; i--){ |
mbedalvaro | 0:345b3bc7a0ea | 224 | x = shiftX*i + offsetX; |
mbedalvaro | 0:345b3bc7a0ea | 225 | writeOutX(int(x)); |
mbedalvaro | 0:345b3bc7a0ea | 226 | wait_us(delayScanning); |
mbedalvaro | 0:345b3bc7a0ea | 227 | lut[i][j] += lockin_read(); |
mbedalvaro | 0:345b3bc7a0ea | 228 | } |
mbedalvaro | 0:345b3bc7a0ea | 229 | //re-scan from left to right |
mbedalvaro | 0:345b3bc7a0ea | 230 | for(int i=0; i<LUT_RESOLUTION; i++){ |
mbedalvaro | 0:345b3bc7a0ea | 231 | x = shiftX*i + offsetX; |
mbedalvaro | 0:345b3bc7a0ea | 232 | writeOutX(int(x)); |
mbedalvaro | 0:345b3bc7a0ea | 233 | wait_us(delayScanning); |
mbedalvaro | 0:345b3bc7a0ea | 234 | lut[i][j] += lockin_read(); |
mbedalvaro | 0:345b3bc7a0ea | 235 | } |
mbedalvaro | 0:345b3bc7a0ea | 236 | } |
mbedalvaro | 0:345b3bc7a0ea | 237 | } |
mbedalvaro | 0:345b3bc7a0ea | 238 | |
mbedalvaro | 0:345b3bc7a0ea | 239 | |
mbedalvaro | 0:345b3bc7a0ea | 240 | //save tab in file |
mbedalvaro | 0:345b3bc7a0ea | 241 | FILE *fp; |
mbedalvaro | 0:345b3bc7a0ea | 242 | #ifdef LUT_FILENAME |
mbedalvaro | 0:345b3bc7a0ea | 243 | fp = fopen(LUT_FILENAME, "w"); // Open file on the local file system for writing |
mbedalvaro | 0:345b3bc7a0ea | 244 | fwrite(lut,sizeof(uint16),LUT_RESOLUTION*LUT_RESOLUTION,fp); |
mbedalvaro | 0:345b3bc7a0ea | 245 | fclose(fp); //close the file (the mBed will appear connected again) |
mbedalvaro | 0:345b3bc7a0ea | 246 | #endif |
mbedalvaro | 0:345b3bc7a0ea | 247 | |
mbedalvaro | 0:345b3bc7a0ea | 248 | #ifdef LUT_H_FILENAME |
mbedalvaro | 22:d87317d7ca91 | 249 | //save tab in Human readable file (not used by the program, this is just for checking) |
mbedalvaro | 22:d87317d7ca91 | 250 | // NOTE: we divide the content of the lut table by NB_SCANS, for easy reading (values should be between 0-4095) |
mbedalvaro | 0:345b3bc7a0ea | 251 | fp = fopen(LUT_H_FILENAME, "w"); // Open file on the local file system for writing |
mbedalvaro | 0:345b3bc7a0ea | 252 | fprintf(fp, "scan resolution: %d x %d\r\n",LUT_RESOLUTION, LUT_RESOLUTION); |
mbedalvaro | 0:345b3bc7a0ea | 253 | for(int j=0; j<LUT_RESOLUTION; j++){ |
mbedalvaro | 0:345b3bc7a0ea | 254 | for(int i=0; i<LUT_RESOLUTION; i++){ |
mbedalvaro | 23:bf666fcc61bc | 255 | fprintf(fp, "X=%d,\tY=%d,\tI=%d\t \r\n", int(shiftX*i + offsetX), int(shiftY*j + offsetY), int(1.0*lut[i][j]/NB_SCANS) ); |
mbedalvaro | 0:345b3bc7a0ea | 256 | } |
mbedalvaro | 0:345b3bc7a0ea | 257 | } |
mbedalvaro | 0:345b3bc7a0ea | 258 | fclose(fp); //close the file (the mBed will appear connected again) |
mbedalvaro | 0:345b3bc7a0ea | 259 | #endif |
mbedalvaro | 0:345b3bc7a0ea | 260 | |
mbedalvaro | 0:345b3bc7a0ea | 261 | } |
mbedalvaro | 0:345b3bc7a0ea | 262 | |
mbedalvaro | 0:345b3bc7a0ea | 263 | |
mbedalvaro | 10:6f8e48dca1bd | 264 | //Return the lockin value "corrected with the Look-UpTable" - this means a RATIO between two reflectivities (and normally, this is <1). |
mbedalvaro | 10:6f8e48dca1bd | 265 | float HardwareIO::lockInCorrectedValue(unsigned short x, unsigned short y){ |
mbedalvaro | 0:345b3bc7a0ea | 266 | //*******Correction using DIRECT approximation |
mbedalvaro | 0:345b3bc7a0ea | 267 | #ifdef LUT_DIRECT |
mbedalvaro | 10:6f8e48dca1bd | 268 | return 2.0* NB_SCANS * lockin_read() / (lut[x >> LUT_BITS_SHIFT][y >> LUT_BITS_SHIFT]); // 2 * NB_SCANS is the number of recorded sample added to one position of the LUT (scan is performed twice: left-right and right-left) |
mbedalvaro | 22:d87317d7ca91 | 269 | #endif |
mbedalvaro | 22:d87317d7ca91 | 270 | |
mbedalvaro | 0:345b3bc7a0ea | 271 | //*******Correction using BILINEAR approximation |
mbedalvaro | 0:345b3bc7a0ea | 272 | #ifdef LUT_BILINEAR |
mbedalvaro | 14:0fc33a3a7b4b | 273 | unsigned short X = x >> LUT_BITS_SHIFT; //mirror "x" is 12bits, LUT "X" needs 4bits when lut is 17x17 |
mbedalvaro | 14:0fc33a3a7b4b | 274 | unsigned short Y = y >> LUT_BITS_SHIFT; //mirror "y" is 12bits, LUT "Y" needs 4bits when lut is 17x17 |
mbedalvaro | 0:345b3bc7a0ea | 275 | float dx = 1.0*(x & LUT_BITS_MASK)/(LUT_BITS_MASK+1); //weight to apply on X (mask with 255 and norm) |
mbedalvaro | 0:345b3bc7a0ea | 276 | float dy = 1.0*(y & LUT_BITS_MASK)/(LUT_BITS_MASK+1); //weight to apply on Y (mask with 255 and norm) |
mbedalvaro | 0:345b3bc7a0ea | 277 | |
mbedalvaro | 0:345b3bc7a0ea | 278 | //Wheighted mean approximation of the Look-Up Table at the position (x,y): |
mbedalvaro | 0:345b3bc7a0ea | 279 | float wmLUT = (1-dy)*( (1-dx)*lut[X][Y] + dx*lut[X+1][Y] ) + dy*( (1-dx)*lut[X][Y+1] + dx*lut[X+1][Y+1] ); |
mbedalvaro | 0:345b3bc7a0ea | 280 | |
mbedalvaro | 22:d87317d7ca91 | 281 | return 2.0* NB_SCANS * lockin_read() / wmLUT;// 2 * NB_SCANS is the number of recorded sample added to one position of the LUT (scan is performed twice: left-right and right-left) |
mbedalvaro | 22:d87317d7ca91 | 282 | #endif |
mbedalvaro | 22:d87317d7ca91 | 283 | |
mbedalvaro | 0:345b3bc7a0ea | 284 | //*******Correction using LINEAR approximation |
mbedalvaro | 0:345b3bc7a0ea | 285 | #ifdef LUT_LINEAR |
mbedalvaro | 14:0fc33a3a7b4b | 286 | unsigned short X = x >> LUT_BITS_SHIFT; //mirror "x" is 12bits, LUT "X" needs 4bits when lut is 17x17 |
mbedalvaro | 14:0fc33a3a7b4b | 287 | unsigned short Y = y >> LUT_BITS_SHIFT; //mirror "y" is 12bits, LUT "Y" needs 4bits when lut is 17x17 |
mbedalvaro | 0:345b3bc7a0ea | 288 | float dx = 1.0*(x & LUT_BITS_MASK)/(LUT_BITS_MASK+1); //weight to apply on X (mask with 255 and norm) |
mbedalvaro | 0:345b3bc7a0ea | 289 | float dy = 1.0*(y & LUT_BITS_MASK)/(LUT_BITS_MASK+1); //weight to apply on Y (mask with 255 and norm) |
mbedalvaro | 0:345b3bc7a0ea | 290 | float linearLUT, dzx, dzy; |
mbedalvaro | 0:345b3bc7a0ea | 291 | |
mbedalvaro | 0:345b3bc7a0ea | 292 | if(dx>dy){ //if the position is on the "top-right" triangle |
mbedalvaro | 0:345b3bc7a0ea | 293 | dzx = (lut[X+1][Y] - lut[X][Y]) * dx; |
mbedalvaro | 0:345b3bc7a0ea | 294 | dzy = (lut[X+1][Y+1] - lut[X+1][Y]) * dy; |
mbedalvaro | 0:345b3bc7a0ea | 295 | } |
mbedalvaro | 0:345b3bc7a0ea | 296 | else{ //if the position is on the "bottom-left" triangle |
mbedalvaro | 0:345b3bc7a0ea | 297 | dzy = (lut[X][Y+1] - lut[X][Y]) * dy; |
mbedalvaro | 0:345b3bc7a0ea | 298 | dzx = (lut[X+1][Y+1] - lut[X][Y+1]) * dx; |
mbedalvaro | 0:345b3bc7a0ea | 299 | } |
mbedalvaro | 0:345b3bc7a0ea | 300 | |
mbedalvaro | 0:345b3bc7a0ea | 301 | //linear approximation of the Look-Up Table at the position (x,y): |
mbedalvaro | 0:345b3bc7a0ea | 302 | linearLUT = lut[X][Y] + dzx + dzy; |
mbedalvaro | 10:6f8e48dca1bd | 303 | return 2.0* NB_SCANS * lockin_read() / linearLUT; // 2 * NB_SCANS is the number of recorded sample added to one position of the LUT (scan is performed twice: left-right and right-left) |
mbedalvaro | 0:345b3bc7a0ea | 304 | |
mbedalvaro | 22:d87317d7ca91 | 305 | #endif |
mbedalvaro | 22:d87317d7ca91 | 306 | |
mbedalvaro | 22:d87317d7ca91 | 307 | //*******No corrections, just return the value divided by 4096 (this means we are assuming that the surface is everywhere perfectly reflective - we supposedly get the max value always) |
mbedalvaro | 22:d87317d7ca91 | 308 | #ifdef NO_LUT |
mbedalvaro | 10:6f8e48dca1bd | 309 | return 1.0* lockin_read()/4096; |
mbedalvaro | 22:d87317d7ca91 | 310 | #endif |
mbedalvaro | 22:d87317d7ca91 | 311 | |
mbedalvaro | 0:345b3bc7a0ea | 312 | } |