Gamma ionizing radiation detection with Raspberry

[root]

Vendor page:
https://www.cooking-hacks.com/documentation/tutorials/geiger-counter-radiation-sensor-board-arduino-raspberry-pi-tutorial

Geiger board schematic
https://www.cooking-hacks.com/media/cooking/images/documentation/geiger_counter_arduino_radiation_sensor_board/radiation_board_v5.pdf

Raspberry software installing instructions:
https://circulosmeos.wordpress.com/2014/04/21/radiation-detection-with-raspberry-pi/

Demonstrative video
https://www.youtube.com/watch?v=vm97_wn2Pac

Arduino version software:
https://circulosmeos.wordpress.com/2014/10/17/radiation-detection-with-arduino-cumulative-counter-2/

Important Notes:

the geiger tube is UV light and alfa ray sensitive, to avoid false readings you need to completely cover it with some black electrical tape.
The buzzer is quite noisy (one buzz per count), if you want to reduce noise, put a piece of black electrical tape to cover the central hole.

My device images:
asterisk/software/raspberry/geiger_counter1.jpg
asterisk/software/raspberry/geiger_counter2.jpg

[root]

I modified the source code to have one reading every minute, because the original value = "1 reading every 10 seconds" does not give consistent values

su - root
cd /root/geiger_tube
nano /geiger_tube.c

EDITED on 28/01/2019

Codice: Seleziona tutto

/*
 *  ------Geiger Tube board (Raspberry Pi + arduPi library Code) Example------
 *
 *  Explanation: This example shows how to get the signal from the Geiger Tube
 *  in Arduino, we use one of the Arduino interrupt pins (PIN2).
 *  We count the time (ms) between two pulses of the Geiger tube.
 *
 *  Copyright (C) 2011 Libelium Comunicaciones Distribuidas S.L.
 *  http://www.libelium.com
 *
 *  This program is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation, either version 2 of the License, or
 *  (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see .
 *
 *  Version:      0.1
 *  Implementation:   Anartz Nuin Jiménez
 *  Version:      0.2
 *                  20140420 by circulosmeos (circulosmeos.wordpress.com)
 *                  arduPi v1.5 adapted and safe loop
 */

// include the arduPi library code:
#include "arduPi.h"

// include the LCD library code:
#include "LiquidCrystal.h"

// Threshold values for the led bar
#define TH1 75
#define TH2 150
#define TH3 250
#define TH4 400
#define TH5 600

// define geiger tube Conversion factor (cpm * conversion factor  = μSv/h)

// for J305By geiger tube
//#define CONV_FACTOR 0.00812

// for SBM-20 geiger tube (1/361)
#define CONV_FACTOR 0.00277

// Variables
int ledArray [] = {10,11,12,13,9};
int geiger_input = 2;
long count = 0;
long countPerMinute = 0;
long timePrevious = 0;
long timePreviousMeassure = 0;
long countPrevious = 0;
float radiationValue = 0.0;

void ledVar(int value);
void countPulse();

void setup(LiquidCrystal *lcd){
  pinMode(geiger_input, INPUT);
  digitalWrite(geiger_input,HIGH);
  for (int i=0;i<5;i++){
    pinMode(ledArray[i],OUTPUT);
  }

  // test led array ON
  pinMode(10, OUTPUT);
  digitalWrite(10, HIGH);
  pinMode(11, OUTPUT);
  digitalWrite(11, HIGH);
  pinMode(12, OUTPUT);
  digitalWrite(12, HIGH);
  pinMode(13, OUTPUT);
  digitalWrite(13, HIGH);
  pinMode(14, OUTPUT);
  digitalWrite(14, HIGH);
  pinMode(15, OUTPUT);
  digitalWrite(15, HIGH);
  pinMode(9, OUTPUT);
  digitalWrite(9, HIGH);

  //set up the LCD\'s number of columns and rows:
  lcd->begin(16, 2);
  lcd->clear();
  lcd->setCursor(0, 0);
  lcd->print("Radiation Sensor");
  lcd->setCursor(0,1);
  lcd->print("Board - 2019 -");
  delay(3000);

  lcd->clear();
  lcd->setCursor(0, 0);
  lcd->print("Counts*minute");
  lcd->setCursor(0,1);
  lcd->print("microSv*hour");
  delay(3000);
  //for (int i=0;i<5;i++){
  //  delay(200);
  //  lcd->scrollDisplayLeft();
  //}
  //delay(500);

  lcd->clear();
  lcd->setCursor(0, 0);
  lcd->print("- www.49v.com -");
  lcd->setCursor(0,1);
  lcd->print(" P.M. Lovisolo ");
  delay(3000);

  //led array off
  pinMode(10, OUTPUT);
  digitalWrite(10, LOW);
  pinMode(11, OUTPUT);
  digitalWrite(11, LOW);
  pinMode(12, OUTPUT);
  digitalWrite(12, LOW);
  pinMode(13, OUTPUT);
  digitalWrite(13, LOW);
  pinMode(14, OUTPUT);
  digitalWrite(14, LOW);
  pinMode(15, OUTPUT);
  digitalWrite(15, LOW);
  pinMode(9, OUTPUT);
  digitalWrite(9, LOW);

  lcd->clear();
  lcd->setCursor(0, 0);
  lcd->print("CPM=");
  lcd->setCursor(4,0);
  lcd->print(6*count);
  lcd->setCursor(8,0);
  lcd->print("counts");
  lcd->setCursor(0,1);
  lcd->print(radiationValue);
  lcd->setCursor(7,1);
  lcd->print(" uSv/h");

  attachInterrupt(2,countPulse,FALLING);

}

void loop(LiquidCrystal *lcd){
  if (millis()-timePreviousMeassure > 60000){
    detachInterrupt(2);
    countPerMinute = 1*count;
    //countPerMinute = count/10;
    radiationValue = countPerMinute * CONV_FACTOR;
    timePreviousMeassure = millis();
    //Serial.print("cpm = ");
    //Serial.print(countPerMinute,DEC);
    //Serial.print(" - ");
    //Serial.print("uSv/h = ");
    //Serial.println(radiationValue,4);
    lcd->clear();
    lcd->setCursor(0, 0);
    lcd->print("CPM=");
    lcd->setCursor(4,0);
    lcd->print(countPerMinute);
    lcd->setCursor(8,0);
    lcd->print("counts");
    lcd->setCursor(0,1);
    lcd->print(radiationValue,4);
    lcd->print((int)radiationValue);
    lcd->setCursor(7,1);
    lcd->print(" uSv/h");

    //fprintf(stderr,"CPM ");
    //fprintf(stderr,"%d\n",countPerMinute);
    fprintf(stderr,"%f",radiationValue);
    fprintf(stderr," uSv/h\n");

    //led var setting
    if(countPerMinute <= TH1) ledVar(0);
    if((countPerMinute <= TH2)&&(countPerMinute>TH1)) ledVar(1);
    if((countPerMinute <= TH3)&&(countPerMinute>TH2)) ledVar(2);
    if((countPerMinute <= TH4)&&(countPerMinute>TH3)) ledVar(3);
    if((countPerMinute <= TH5)&&(countPerMinute>TH4)) ledVar(4);
    if(countPerMinute>TH5) ledVar(5);

    count = 0;

    attachInterrupt(2,countPulse,FALLING);

    delay (60000);

  }

}

void countPulse(){
  //detachInterrupt(2);
  count++;
  //while(digitalRead(2)==0){
  //}
  //attachInterrupt(2,countPulse,FALLING);
}

void ledVar(int value){
  if (value > 0){
    for(int i=0;i<=value;i++){
      digitalWrite(ledArray[i],HIGH);
    }
    for(int i=5;i>value;i--){
      digitalWrite(ledArray[i],LOW);
    }
  }
  else {
    for(int i=5;i>=0;i--){
      digitalWrite(ledArray[i],LOW);
    }
  }
}


int main (){

    // initialize the library with the numbers of the interface pins
    // AFTER WirePi has been initialized (!)
    LiquidCrystal lcd(3,4,5,6,7,8);

    setup(&lcd);

    while(1){
        loop(&lcd);
    }
    return (0);
}



Note:
modified for the SBM-20 tube, if you want to use the J305B tube, you need to comment-out the SBM-20 CONV_FACTOR section and remove comment in the J305 section this way:

Codice: Seleziona tutto

// for J305By geiger tube
#define CONV_FACTOR 0.00812
// for SBM-20 geiger tube (1/361)
//#define CONV_FACTOR 0.00277

COMPILATION:
===========

Codice: Seleziona tutto

g++ -lrt -lpthread geiger_tube.c arduPi.cpp -o geiger_tube

setting permissions:

Codice: Seleziona tutto

chmod 755 geiger_tube

[root]

Script to put the readings into a file to be shown in a web page:

print_rad2file.sh

Codice: Seleziona tutto

/root/geiger_tube/geiger_tube 2>&1 | awk '{ print strftime("%c: "), $0; fflush(); }' | tee /var/www/rad/$( date '+%d-%m-%Y' ) &

note: the name of the generated file is the current date, ex: /var/www/rad/12-01-2019

the script starts and goes to background, to change the date in the file name with the current date (otherwise it keeps the date of the first run), you should stop and restart the script at midnight every day, by doing the following:

crontab -e

Codice: Seleziona tutto

0 0 * * * killall geiger_tube
0 1 * * * /bin/sh /root/print_rad2file.sh

then configure Apache to point to "/var/www/rad".

example of what is shown in a browser by pointing it to http://localhost/rad:

Index of /rad
Name Last modified Size Description
[<] Parent Directory -
[?] 11-01-2019 2019-01-11 20:36 166
[?] 12-01-2019 2019-01-12 23:36 466
[?] 13-01-2019 2019-01-13 23:36 366
-------------------------------------------------------------------
Apache/2.4.10 (Raspbian) Server at 192.168.1.83 Port 80

click on file name (date) to open it:

ven 11 gen 2019 20:45:23 CET: CPM=9
ven 11 gen 2019 20:45:23 CET: 0.073080 uSv/h
ven 11 gen 2019 20:46:23 CET: CPM=5
ven 11 gen 2019 20:46:24 CET: 0.040600 uSv/h
ven 11 gen 2019 20:47:24 CET: CPM=5
ven 11 gen 2019 20:47:24 CET: 0.040600 uSv/h
ven 11 gen 2019 20:48:24 CET: CPM=6
ven 11 gen 2019 20:48:24 CET: 0.048720 uSv/h
ven 11 gen 2019 20:49:24 CET: CPM=11
ven 11 gen 2019 20:49:24 CET: 0.089320 uSv/h
ven 11 gen 2019 20:50:24 CET: CPM=6
ven 11 gen 2019 20:50:24 CET: 0.048720 uSv/h
ven 11 gen 2019 20:51:24 CET: CPM=7
ven 11 gen 2019 20:51:24 CET: 0.056840 uSv/h
ven 11 gen 2019 20:52:24 CET: CPM=9
ven 11 gen 2019 20:52:24 CET: 0.073080 uSv/h
ven 11 gen 2019 20:53:24 CET: CPM=11
ven 11 gen 2019 20:53:24 CET: 0.089320 uSv/h
ven 11 gen 2019 20:54:24 CET: CPM=2
ven 11 gen 2019 20:54:24 CET: 0.016240 uSv/h
ven 11 gen 2019 20:55:24 CET: CPM=7
ven 11 gen 2019 20:55:24 CET: 0.056840 uSv/h
ven 11 gen 2019 20:56:24 CET: CPM=5
ven 11 gen 2019 20:56:24 CET: 0.040600 uSv/h
ven 11 gen 2019 20:57:24 CET: CPM=9
ven 11 gen 2019 20:57:24 CET: 0.073080 uSv/h
ven 11 gen 2019 20:58:24 CET: CPM=10
ven 11 gen 2019 20:58:24 CET: 0.081200 uSv/h
ven 11 gen 2019 20:59:24 CET: CPM=9
ven 11 gen 2019 20:59:24 CET: 0.073080 uSv/h
ven 11 gen 2019 21:00:24 CET: CPM=8
ven 11 gen 2019 21:00:24 CET: 0.064960 uSv/h

Note:
CPM = counts per minute
uSv/h = microSievert per hour

[root]

To redirect output to csv file in "comma-delimited" format [unix-time,date,time,value,unit of measure]

example:

Codice: Seleziona tutto

1547278704,2019-01-12,08:38:24,0.024360,uSv/h

note: unix-time (unique value) is needed to avoid duplicates when importing data to mysql (see next post "example3")


First of all I modified the source code to put a "comma" between data and uSv/h, also I commented out the CPM section, since I don't need the CPM value,
than I modfied the script to put "commas" instead of spaces in the output.

modified source code:

su - root
cd /root/geiger_tube
nano measure.c

Codice: Seleziona tutto

/*
 *  ------Geiger Tube board (Arduino Code) Example--------
 *
 *  Explanation: This example shows how to get the signal from the Geiger Tube
 *  in Arduino, we use one of the Arduino interrupt pins (PIN2).
 *  We count the time (ms) between two pulses of the Geiger tube.
 *
 *  Copyright (C) 2011 Libelium Comunicaciones Distribuidas S.L.
 *  http://www.libelium.com
 *
 *  This program is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation, either version 2 of the License, or
 *  (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see <http://www.gnu.org/licenses/>.
 *
 *  Version:      0.3
 *  Design:      Marcos Yarza, David Gascon
 *  Implementation:   Marcos Yarza
 *
 *  Version     0.4
 *              20140420 by circulosmeos for RPi compatibility
 *              circulosmeos.wordpress.com
 *
 *  modified by P.M. Lovisolo Italy - http://forum.49v.com
 */

// include the arduPi library code:
#include "arduPi.h"

// include the LCD library code:
#include "LiquidCrystal.h"

// Threshold values for the led bar
#define TH1 75
#define TH2 150
#define TH3 250
#define TH4 400
#define TH5 600

// define geiger tube Conversion factor (cpm * conversion factor  = μSv/h)

// for J305By geiger tube
//#define CONV_FACTOR 0.00812

// for SBM-20 geiger tube (1/361)
#define CONV_FACTOR 0.00277

// Variables
int ledArray [] = {10,11,12,13,9};
int geiger_input = 2;
long count = 0;
long countPerMinute = 0;
long timePrevious = 0;
long timePreviousMeassure = 0;
long countPrevious = 0;
float radiationValue = 0.0;

void ledVar(int value);
void countPulse();

void setup(LiquidCrystal *lcd){
  pinMode(geiger_input, INPUT);
  digitalWrite(geiger_input,HIGH);
  for (int i=0;i<5;i++){
    pinMode(ledArray[i],OUTPUT);
  }

  // test led array ON
  pinMode(10, OUTPUT);
  digitalWrite(10, HIGH);
  pinMode(11, OUTPUT);
  digitalWrite(11, HIGH);
  pinMode(12, OUTPUT);
  digitalWrite(12, HIGH);
  pinMode(13, OUTPUT);
  digitalWrite(13, HIGH);
  pinMode(14, OUTPUT);
  digitalWrite(14, HIGH);
  pinMode(15, OUTPUT);
  digitalWrite(15, HIGH);
  pinMode(9, OUTPUT);
  digitalWrite(9, HIGH);

  //set up the LCD\'s number of columns and rows:
  lcd->begin(16, 2);
  lcd->clear();
  lcd->setCursor(0, 0);
  lcd->print("Radiation Sensor");
  lcd->setCursor(0,1);
  lcd->print("Board - 2019 -");
  delay(3000);

  lcd->clear();
  lcd->setCursor(0, 0);
  lcd->print("Counts*minute");
  lcd->setCursor(0,1);
  lcd->print("microSv*hour");
  delay(3000);
  //for (int i=0;i<5;i++){
  //  delay(200);
  //  lcd->scrollDisplayLeft();
  //}
  //delay(500);

  lcd->clear();
  lcd->setCursor(0, 0);
  lcd->print("- www.49v.com -");
  lcd->setCursor(0,1);
  lcd->print(" P.M. Lovisolo ");
  delay(3000);

  //led array off
  pinMode(10, OUTPUT);
  digitalWrite(10, LOW);
  pinMode(11, OUTPUT);
  digitalWrite(11, LOW);
  pinMode(12, OUTPUT);
  digitalWrite(12, LOW);
  pinMode(13, OUTPUT);
  digitalWrite(13, LOW);
  pinMode(14, OUTPUT);
  digitalWrite(14, LOW);
  pinMode(15, OUTPUT);
  digitalWrite(15, LOW);
  pinMode(9, OUTPUT);
  digitalWrite(9, LOW);

  lcd->clear();
  lcd->setCursor(0, 0);
  lcd->print("CPM=");
  lcd->setCursor(4,0);
  lcd->print(6*count);
  lcd->setCursor(8,0);
  lcd->print("counts");
  lcd->setCursor(0,1);
  lcd->print(radiationValue);
  lcd->setCursor(7,1);
  lcd->print(" uSv/h");

  attachInterrupt(2,countPulse,FALLING);

}

void loop(LiquidCrystal *lcd){
  if (millis()-timePreviousMeassure > 60000){
    detachInterrupt(2);
    countPerMinute = 1*count;
    //countPerMinute = count/10;
    radiationValue = countPerMinute * CONV_FACTOR;
    timePreviousMeassure = millis();
    //Serial.print("cpm = ");
    //Serial.print(countPerMinute,DEC);
    //Serial.print(" - ");
    //Serial.print("uSv/h = ");
    //Serial.println(radiationValue,4);
    lcd->clear();
    lcd->setCursor(0, 0);
    lcd->print("CPM=");
    lcd->setCursor(4,0);
    lcd->print(countPerMinute);
    lcd->setCursor(8,0);
    lcd->print("counts");
    lcd->setCursor(0,1);
    lcd->print(radiationValue,4);
    lcd->print((int)radiationValue);
    lcd->setCursor(7,1);
    lcd->print(" uSv/h");

    //fprintf(stderr,"CPM ");
    //fprintf(stderr,"%d\n",countPerMinute);
    fprintf(stderr,"%f",radiationValue);
    fprintf(stderr," uSv/h\n");

    //led var setting
    if(countPerMinute <= TH1) ledVar(0);
    if((countPerMinute <= TH2)&&(countPerMinute>TH1)) ledVar(1);
    if((countPerMinute <= TH3)&&(countPerMinute>TH2)) ledVar(2);
    if((countPerMinute <= TH4)&&(countPerMinute>TH3)) ledVar(3);
    if((countPerMinute <= TH5)&&(countPerMinute>TH4)) ledVar(4);
    if(countPerMinute>TH5) ledVar(5);

    count = 0;

    attachInterrupt(2,countPulse,FALLING);

    delay (60000);

  }

}

void countPulse(){
  //detachInterrupt(2);
  count++;
  //while(digitalRead(2)==0){
  //}
  //attachInterrupt(2,countPulse,FALLING);
}

void ledVar(int value){
  if (value > 0){
    for(int i=0;i<=value;i++){
      digitalWrite(ledArray[i],HIGH);
    }
    for(int i=5;i>value;i--){
      digitalWrite(ledArray[i],LOW);
    }
  }
  else {
    for(int i=5;i>=0;i--){
      digitalWrite(ledArray[i],LOW);
    }
  }
}


int main (){

    // initialize the library with the numbers of the interface pins
    // AFTER WirePi has been initialized (!)
    LiquidCrystal lcd(3,4,5,6,7,8);

    setup(&lcd);

    while(1){
        loop(&lcd);
    }
    return (0);
}


compile:

Codice: Seleziona tutto

g++ -lrt -lpthread measure.c arduPi.cpp -o measure

setting permissions:

Codice: Seleziona tutto

chmod 755 measure

new script to be executed:

Codice: Seleziona tutto

measure 2>&1 | awk '{ print strftime("%s %Y-%m-%d %H:%M:%S"), $0; fflush(); }'  | sed -e 's/\s\+/,/g;w /var/www/rad/readings.csv' &

[root]

Create database name: "data" > table name: "letture_radiazioni" with the following coloumns:

unix_time (int11)
date (date)
time (time)
value (decimal(10,6)
unit (text)

UNIQUE KEY = 'unix_time` to avoid importing duplicates

or you can just import my msql_dump :

data.sql

Codice: Seleziona tutto

-- MySQL dump 10.13  Distrib 5.5.62, for debian-linux-gnu (armv8l)
--
-- Host: localhost    Database: data
-- ------------------------------------------------------
-- Server version   5.5.62-0+deb8u1

/*!40101 SET @OLD_CHARACTER_SET_CLIENT=@@CHARACTER_SET_CLIENT */;
/*!40101 SET @OLD_CHARACTER_SET_RESULTS=@@CHARACTER_SET_RESULTS */;
/*!40101 SET @OLD_COLLATION_CONNECTION=@@COLLATION_CONNECTION */;
/*!40101 SET NAMES utf8 */;
/*!40103 SET @OLD_TIME_ZONE=@@TIME_ZONE */;
/*!40103 SET TIME_ZONE='+00:00' */;
/*!40014 SET @OLD_UNIQUE_CHECKS=@@UNIQUE_CHECKS, UNIQUE_CHECKS=0 */;
/*!40014 SET @OLD_FOREIGN_KEY_CHECKS=@@FOREIGN_KEY_CHECKS, FOREIGN_KEY_CHECKS=0 */;
/*!40101 SET @OLD_SQL_MODE=@@SQL_MODE, SQL_MODE='NO_AUTO_VALUE_ON_ZERO' */;
/*!40111 SET @OLD_SQL_NOTES=@@SQL_NOTES, SQL_NOTES=0 */;

--
-- Table structure for table `letture_radiazioni`
--

DROP TABLE IF EXISTS `letture_radiazioni`;
/*!40101 SET @saved_cs_client     = @@character_set_client */;
/*!40101 SET character_set_client = utf8 */;
CREATE TABLE `letture_radiazioni` (
  `unix_time` int(11) DEFAULT NULL,
  `date` date DEFAULT NULL,
  `time` time DEFAULT NULL,
  `value` decimal(10,6) DEFAULT NULL,
  `unit` text,
  UNIQUE KEY `unix_time` (`unix_time`)
) ENGINE=InnoDB DEFAULT CHARSET=latin1;
/*!40101 SET character_set_client = @saved_cs_client */;

--
-- Dumping data for table `letture_radiazioni`
--

LOCK TABLES `letture_radiazioni` WRITE;
/*!40000 ALTER TABLE `letture_radiazioni` DISABLE KEYS */;
INSERT INTO `letture_radiazioni` VALUES (1547276291,'2019-01-12','07:58:11',0.040600,'uSv/h'),(1547276351,'2019-01-12','07:59:11',0.081200,'uSv/h'),(1547276411,'2019-01-12','08:00:11',0.089320,'uSv/h'),(1547276471,'2019-01-12','08:01:11',0.056840,'uSv/h'),(1547276531,'2019-01-12','08:02:11',0.064960,'uSv/h'),(1547276591,'2019-01-12','08:03:11',0.073080,'uSv/h'),(1547276651,'2019-01-12','08:04:11',0.056840,'uSv/h'),(1547276711,'2019-01-12','08:05:11',0.073080,'uSv/h'),(1547276771,'2019-01-12','08:06:11',0.032480,'uSv/h'),(1547276831,'2019-01-12','08:07:11',0.089320,'uSv/h'),(1547276891,'2019-01-12','08:08:11',0.097440,'uSv/h'),(1547276951,'2019-01-12','08:09:11',0.064960,'uSv/h'),(1547277011,'2019-01-12','08:10:11',0.064960,'uSv/h'),(1547277071,'2019-01-12','08:11:11',0.064960,'uSv/h'),(1547277131,'2019-01-12','08:12:11',0.073080,'uSv/h'),(1547277191,'2019-01-12','08:13:11',0.064960,'uSv/h'),(1547277251,'2019-01-12','08:14:11',0.073080,'uSv/h'),(1547277311,'2019-01-12','08:15:11',0.040600,'uSv/h'),(1547277371,'2019-01-12','08:16:11',0.064960,'uSv/h'),(1547277431,'2019-01-12','08:17:11',0.064960,'uSv/h'),(1547277491,'2019-01-12','08:18:11',0.064960,'uSv/h'),(1547277551,'2019-01-12','08:19:11',0.056840,'uSv/h'),(1547277611,'2019-01-12','08:20:11',0.056840,'uSv/h'),(1547277671,'2019-01-12','08:21:11',0.089320,'uSv/h'),(1547277731,'2019-01-12','08:22:11',0.048720,'uSv/h'),(1547277791,'2019-01-12','08:23:11',0.064960,'uSv/h'),(1547277851,'2019-01-12','08:24:11',0.056840,'uSv/h'),(1547277911,'2019-01-12','08:25:11',0.016240,'uSv/h'),(1547277971,'2019-01-12','08:26:11',0.048720,'uSv/h'),(1547278031,'2019-01-12','08:27:11',0.064960,'uSv/h'),(1547278091,'2019-01-12','08:28:11',0.089320,'uSv/h'),(1547278151,'2019-01-12','08:29:11',0.040600,'uSv/h'),(1547278211,'2019-01-12','08:30:11',0.073080,'uSv/h'),(1547278271,'2019-01-12','08:31:11',0.105560,'uSv/h'),(1547278704,'2019-01-12','08:38:24',0.040600,'uSv/h'),(1547278764,'2019-01-12','08:39:24',0.048720,'uSv/h'),(1547278824,'2019-01-12','08:40:24',0.064960,'uSv/h');
/*!40000 ALTER TABLE `letture_radiazioni` ENABLE KEYS */;
UNLOCK TABLES;

--
-- Dumping events for database 'data'
--

--
-- Dumping routines for database 'data'
--
/*!40103 SET TIME_ZONE=@OLD_TIME_ZONE */;

/*!40101 SET SQL_MODE=@OLD_SQL_MODE */;
/*!40014 SET FOREIGN_KEY_CHECKS=@OLD_FOREIGN_KEY_CHECKS */;
/*!40014 SET UNIQUE_CHECKS=@OLD_UNIQUE_CHECKS */;
/*!40101 SET CHARACTER_SET_CLIENT=@OLD_CHARACTER_SET_CLIENT */;
/*!40101 SET CHARACTER_SET_RESULTS=@OLD_CHARACTER_SET_RESULTS */;
/*!40101 SET COLLATION_CONNECTION=@OLD_COLLATION_CONNECTION */;
/*!40111 SET SQL_NOTES=@OLD_SQL_NOTES */;

-- Dump completed on 2019-01-12  8:41:11

script to read the "readings.csv" file and insert data into mysql:

import_csv2sql.sh

Codice: Seleziona tutto

mysql -uUSER -pPASSWOD --local_infile=1 data -e "LOAD DATA LOCAL INFILE '/var/www/rad/readings.csv' INTO TABLE letture_radiazioni FIELDS TERMINATED BY ',' enclosed by '\"'"

note: change USER and PASSWORD with your db credentials


launch the above script every minute:
crontab -e

Codice: Seleziona tutto

*/1 * * * * /bin/sh /root/import_csv2sql.sh

[root]

Work in progress ...

It will be ready very soon

[root]

First test to plot data in a web page, I used python, pandas and ploty.

https://plot.ly/~pilovis/6/utext-uusvh- ... ix-time/#/

X axis = Unix time
Y axis = Radiation value (uSv/h)

That's not what I want, but it is a good starting point

the code I used:

plot_graph.py

Codice: Seleziona tutto

import MySQLdb
conn = MySQLdb.connect(host="localhost", user="pilovis", passwd="SECRET", db="data")
cursor = conn.cursor()
cursor.execute('select unix_time, date, time, value, unit from letture_radiazioni');
rows = cursor.fetchall()
import pandas as pd
df = pd.DataFrame( [[ij for ij in i] for i in rows] )
df.rename(columns={0: 'unix_time', 1: 'date', 2: 'time', 3: 'value', 4: 'unit'}, inplace=True);
df = df.sort(['value'], ascending=[1]);
import plotly.plotly as py
from plotly.graph_objs import *
py.sign_in('pilovis', 'SECRET')
unit = df['unit']

trace1 = Scatter(
     x=df['unix_time'],
     y=df['value'],
     text=unit,
     mode='markers'
)
layout = Layout(
     xaxis=XAxis( type= 'linear', title='unix_time' ),
     yaxis=YAxis( type= 'linear', title='uSv/h' )
)
data = Data([trace1])

fig = Figure(data=data, layout=layout)
py.iplot(fig, filename='radiation')

[root]

After two full days of testing I decided to increase the measurement time from 1 minute to 10, this gives more consistent data results.

Codice: Seleziona tutto

/*
 *  ------Geiger Tube board (Raspberry Pi + arduPi library Code) Example------
 *
 *  Explanation: This example shows how to get the signal from the Geiger Tube
 *  in Arduino, we use one of the Arduino interrupt pins (PIN2).
 *  We count the time (ms) between two pulses of the Geiger tube.
 *
 *  Copyright (C) 2011 Libelium Comunicaciones Distribuidas S.L.
 *  http://www.libelium.com
 *
 *  This program is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation, either version 2 of the License, or
 *  (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see .
 *
 *  Version:      0.1
 *  Implementation:   Anartz Nuin Jiménez
 *  Version:      0.2
 *                  20140420 by circulosmeos (circulosmeos.wordpress.com)
 *                  arduPi v1.5 adapted and safe loop
 */

#include "arduPi.h"

// Not needed since arduPi v1.5
//SerialPi Serial;
//WirePi Wire;
//SPIPi SPI;

// Conversion factor - CPM to uSV/h
#define CONV_FACTOR 0.00812

int geiger_input = 2;
long count = 0;
long timePreviousMeassure = 0;
long countPerMinute = 0;
float radiationValue = 0.0;

void countPulse();

void setup(){
   pinMode(geiger_input, INPUT);
   digitalWrite(geiger_input,HIGH);
   attachInterrupt(2,countPulse, FALLING);
}

void loop(){

  if (millis()-timePreviousMeassure > 600000){

    timePreviousMeassure = millis();

   //detachInterrupt for accesing securely to count variable
   detachInterrupt(2);

   countPerMinute = count/10;
   radiationValue = countPerMinute * CONV_FACTOR;

//   fprintf(stderr,"CPM ");
//   fprintf(stderr,"%d\n",countPerMinute);
   fprintf(stderr,"%f",radiationValue);
   fprintf(stderr," uSv/h\n");

   count = 0;

   attachInterrupt(2,countPulse,FALLING);
   
   delay(600000);
  }
}

void countPulse(){
  count++;
}

int main (){
   setup();
   while(1){
      loop();
   }
   return (0);
}

[root]

To have just the "date time readingvalue unit" printed out, substituting commas with spaces and readingvalue with just three decimals,

example: 2019-01-17 21:18:11 0.040 uSv/h

this the modified file:

nano printout-html.sh

Codice: Seleziona tutto

#!/bin/bash
while :
do
/usr/bin/tail -n 1 /var/www/rad/readings.csv | cut --complement -c  1-11 | cut --complement -c  37-39 | sed -e 's/,/ /g' > /var/www/misura.txt
/bin/sleep 10
done

[root]

My Arduino code to measure uSv/h, send it to stdout and to activate LCD display and bargraph

Prerequisites: arduPi library

measure.c

Codice: Seleziona tutto

/*
 *  ------Geiger Tube board (Arduino Code) Example--------
 *
 *  Explanation: This example shows how to get the signal from the Geiger Tube
 *  in Arduino, we use one of the Arduino interrupt pins (PIN2).
 *  We count the time (ms) between two pulses of the Geiger tube.
 *
 *  Copyright (C) 2011 Libelium Comunicaciones Distribuidas S.L.
 *  http://www.libelium.com
 *
 *  This program is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation, either version 2 of the License, or
 *  (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see <http://www.gnu.org/licenses/>.
 *
 *  Version:      0.3
 *  Design:      Marcos Yarza, David Gascon
 *  Implementation:   Marcos Yarza
 *
 *  Version     0.4
 *              20140420 by circulosmeos for RPi compatibility
 *              circulosmeos.wordpress.com
 *
 *  modified by P.M. Lovisolo Italy - http://forum.49v.com
 */

// include the arduPi library code:
#include "arduPi.h"

// include the LCD library code:
#include "LiquidCrystal.h"

// Threshold values for the led bar
#define TH1 35
#define TH2 75
#define TH3 180
#define TH4 350
#define TH5 500

// Conversion factor - CPM to uSV/h
#define CONV_FACTOR 0.00812

// Variables
int ledArray [] = {10,11,12,13,9};
int geiger_input = 2;
long count = 0;
long countPerMinute = 0;
long timePrevious = 0;
long timePreviousMeassure = 0;
long countPrevious = 0;
float radiationValue = 0.0;

void ledVar(int value);
void countPulse();

void setup(LiquidCrystal *lcd){
  pinMode(geiger_input, INPUT);
  digitalWrite(geiger_input,HIGH);
  for (int i=0;i<5;i++){
    pinMode(ledArray[i],OUTPUT);
  }

  // test led array ON
  pinMode(10, OUTPUT);
  digitalWrite(10, HIGH);
  pinMode(11, OUTPUT);
  digitalWrite(11, HIGH);
  pinMode(12, OUTPUT);
  digitalWrite(12, HIGH);
  pinMode(13, OUTPUT);
  digitalWrite(13, HIGH);
  pinMode(14, OUTPUT);
  digitalWrite(14, HIGH);
  pinMode(15, OUTPUT);
  digitalWrite(15, HIGH);
  pinMode(9, OUTPUT);
  digitalWrite(9, HIGH);

  //set up the LCD\'s number of columns and rows:
  lcd->begin(16, 2);
  lcd->clear();
  lcd->setCursor(0, 0);
  lcd->print("Radiation Sensor");
  lcd->setCursor(0,1);
  lcd->print("Board - 2019 -"); 
  delay(2000);

  lcd->clear();
  lcd->setCursor(0,1); 
  lcd->print("Geiger Muller tube mod J305BY");
  delay(500);
  for (int i=0;i<5;i++){
    delay(200); 
    lcd->scrollDisplayLeft();
  }
  delay(500);

  lcd->clear();
  lcd->setCursor(0, 0);
  lcd->print("- www.49v.com -");
  lcd->setCursor(0,1);
  lcd->print(" P.M. Lovisolo ");   
  delay(3500);

  //led array off
  pinMode(10, OUTPUT);
  digitalWrite(10, LOW);
  pinMode(11, OUTPUT);
  digitalWrite(11, LOW);
  pinMode(12, OUTPUT);
  digitalWrite(12, LOW);
  pinMode(13, OUTPUT);
  digitalWrite(13, LOW);
  pinMode(14, OUTPUT);
  digitalWrite(14, LOW);
  pinMode(15, OUTPUT);
  digitalWrite(15, LOW);
  pinMode(9, OUTPUT);
  digitalWrite(9, LOW);

  lcd->clear();
  lcd->setCursor(0, 0);
  lcd->print("CPM=");
  lcd->setCursor(4,0);
  lcd->print(6*count);
  lcd->setCursor(8,0);
  lcd->print("counts");
  lcd->setCursor(0,1);
  lcd->print(radiationValue);
  lcd->setCursor(7,1);
  lcd->print(" uSv/h");

  attachInterrupt(2,countPulse,FALLING);

}

void loop(LiquidCrystal *lcd){
  if (millis()-timePreviousMeassure > 60000){
    detachInterrupt(2);
    countPerMinute = 1*count;
    radiationValue = countPerMinute * CONV_FACTOR;
    timePreviousMeassure = millis();
    //Serial.print("cpm = ");
    //Serial.print(countPerMinute,DEC);
    //Serial.print(" - ");
    //Serial.print("uSv/h = ");
    //Serial.println(radiationValue,4);
    lcd->clear();
    lcd->setCursor(0, 0);
    lcd->print("CPM=");
    lcd->setCursor(4,0);
    lcd->print(countPerMinute);
    lcd->setCursor(8,0);
    lcd->print("counts");
    lcd->setCursor(0,1);
    lcd->print(radiationValue,4);
    lcd->print((int)radiationValue);
    lcd->setCursor(7,1);
    lcd->print(" uSv/h");

    //fprintf(stderr,"CPM ");
    //fprintf(stderr,"%d\n",countPerMinute);
    fprintf(stderr,"%f",radiationValue);
    fprintf(stderr," uSv/h\n");

    //led var setting
    if(countPerMinute <= TH1) ledVar(0);
    if((countPerMinute <= TH2)&&(countPerMinute>TH1)) ledVar(1);
    if((countPerMinute <= TH3)&&(countPerMinute>TH2)) ledVar(2);
    if((countPerMinute <= TH4)&&(countPerMinute>TH3)) ledVar(3);
    if((countPerMinute <= TH5)&&(countPerMinute>TH4)) ledVar(4);
    if(countPerMinute>TH5) ledVar(5);

    count = 0;

    attachInterrupt(2,countPulse,FALLING);

    delay (60000);

  }

}

void countPulse(){
  //detachInterrupt(2);
  count++;
  //while(digitalRead(2)==0){
  //}
  //attachInterrupt(2,countPulse,FALLING);
}

void ledVar(int value){
  if (value > 0){
    for(int i=0;i<=value;i++){
      digitalWrite(ledArray[i],HIGH);
    }
    for(int i=5;i>value;i--){
      digitalWrite(ledArray[i],LOW);
    }
  }
  else {
    for(int i=5;i>=0;i--){
      digitalWrite(ledArray[i],LOW);
    }
  }
}


int main (){

    // initialize the library with the numbers of the interface pins
    // AFTER WirePi has been initialized (!)
    LiquidCrystal lcd(3,4,5,6,7,8);

    setup(&lcd);

    while(1){
        loop(&lcd);
    }
    return (0);
}

Note: my data readings range from 3 to 9 CPM, my tube (North Optic J305β) conversion factor is 0.00812, I set one reading for 6000 ms (measure time = 1 minute).

compile the code for Raspberry:

Codice: Seleziona tutto

g++ -lrt -lpthread measure.c LiquidCrystal.cpp Print.cpp arduPi.cpp -o geiger_measure

setting permission to file:

Codice: Seleziona tutto

chmod 755 measure