Разработка специализированной компьютерной системы для определения кислотно-основного состояния жидкой среды

Токарев Д. В. (Руководитель – проф. Герасимов И. Г.)

ДонНТУ, кафедра АСУ

Abstract

Тоkаrеv D.V. Development of specialized computer system for determination environment acid-base balance. The algorithm of definition environment acid-base balance under the color image is submitted in this job, the methods of a filtration of the digital image are described, taking into account their specificity, the optimum color model is chosen, the equations of linear regress between pH and parameters of the image are received.

Введение

Водородный показатель имеет важное общебиологическое значение, в связи с чем в процессе эволюции у большинства живых организмов выработался ряд механизмов, обеспечивающих относительное постоянство этого показателя в клетке. Роль этого фактора определяется в первую очередь его влиянием на активность ферментов и состояние других белковых молекул. Кроме того, поскольку большинство реакций в клетках протекает в водной среде, избыток или недостаток ионов Н⁺ может существенно влиять на протекание также различных неферментативных реакций. Внутрижелудочная рН метрия нужна именно для того, чтобы ответить на вопрос: когда, где и сколько у Вас вырабатывается кислоты? Какой она концентрации? Как она распределена по стенкам пищевода, желудка и двенадцатиперстной кишки? Как эффективно работают клетки, вырабатывающие щелочной секрет, необходимый для нейтрализации кислоты в нижней части желудка? Еще измерение кислотности нужно для того, чтобы после соответствующего исследования подобрать лекарства, которые лучше всего подойдут для лечения именно Вас, с учетом особенностей вашего организма. Поэтому, человек нуждается в быстрых методах определении pH жидких сред организма.

Современная медицина и техника предоставляют несколько способов определения кислотности, исследуя пробы жидких сред организма человека, однако моя магистерская посвящена разработке компьютерной системы для определения pH по цветному изображению соответствующих сред (без непосредственного анализа образцов), поэтому наши с научным руководителем исследования можно назвать исследованиями “переднего фронта”.

Понятие кислотности среды

Показатель водородных ионов (H⁺) — это физическая величина, считающаяся количественной мерой активности водородных ионов в растворе (кислотности или щелочности раствора). Для удобства пользуются величиной, равной отрицательному десятичному логарифму от молярной концентрации соответствующих ионов (р-функция, рН). Тогда величина рН = -Lg[НзО⁺]. (1)

Величина рН (-Lg[НзО⁺]) в этом случае равна 7. В кислой среде [НзО⁺]>[ОН^-], а рН<7. Для щелочных сред, наоборот, [НзО⁺]<[ОН^-], а рН>7. Состоянию нормы соответствует узкий интервал значений рН — от 7,37 до 7,44.

Обзор выполненных работ

Методы определения pH среды и существующие системы

Для определения рН используют главным образом индикаторный и ионометрический методы. Индикаторный метод используется в основном в тех случаях, когда необходимо быстро приблизительно определить величину рН исследуемого раствора. Этот метод, как правило, затруднительно применять при исследовании мутных и окрашенных растворов. Ионометрический метод дает возможность определения рН любых водных растворов с гораздо большей точностью (±0,03 ед. рН).

Индикаторный метод

Индикаторный метод основан на применении кислотно-основных индикаторов. Такие индикаторы обычно являются органическими соединениями, проявляющими свойства слабых кислот или оснований. Реакция диссоциации или ассоциации индикаторов приводит к изменению окраски. Специфическое свойство индикаторов состоит в том, что их исходная и прореагировавшая формы значительно отличаются по окраске.

Типичную реакцию кислотно-основного индикатора (условно обозначаемого как Hind) можно представить в следующем виде:
(кислая форма) или (2)
(щелочная форма). (3)
При этом кислота Hind и сопряженное с ней основание Ind⁺ различаются окраской.

Равновесие диссоциации индикаторов характеризуется соответствующими константами Kind — эта константа связана с кислотностью среды следующим уравнением:

рН=pKind+lg[lnd⁺]/[Hind], где (4)

рН — р-функция концентрации ионов водорода;
pKind — показатель индикатора (-LgKind);
[Ind⁺] и [Hind] — концентрации щелочной и кислотной форм индикатора соответственно.

Раствор, содержащий индикатор, непрерывно изменяет свою окраску при изменении рН, т.к. согласно вышеприведенному уравнению меняется соотношение между концентрациями щелочной и кислотной форм индикатора, имеющими, как уже было сказано, различную окраску. Человеческий глаз, однако, не очень чувствителен к таким изменениям. Обычно требуется 5-10-кратный избыток одной из форм, чтобы только ее окраска воспринималась наблюдателем; дальнейшее увеличение избытка этой формы трудно заметить. Поэтому интервал перехода окраски индикатора укладывается, как правило, в интервал, определяемый следующим соотношением: рН=pKind ± l. Наиболее отчетливо изменение цвета будет наблюдаться примерно в середине интервала при значении рН, близком к показателю индикатора.

Методы ионометрии

Методы прямой потенциометрии

1. Метод концентрационного элемента основан на измерении разности потенциалов между двумя идентичными ИСЭ, находящимися в растворе известного состава и в анализируемом растворе с неизвестным содержанием определяемого иона. Вычисления активной концентрации ведутся с использованием вышеприведенных уравнений. Метод эффективен в автоматических методах анализа.

2. Метод градуированного графика основан на измерении потенциала индикаторного электрода в растворе с неизвестной концентрацией определяемого иона и расчете этой концентрации по уравнению регрессии, найденному по серии градуировочных растворов с известной концентрацией этого же иона. Для снижения погрешности анализа градуировочный график строят по серии растворов, состав которых (концентрация инертного электролита и рН) максимально приближен к составу анализируемого раствора. Как правило, для этого во все измеряемые растворы вводят специальные буферные смеси, обеспечивающие постоянство ионной силы, рН и устраняющие мешающее влияние ионов, сопутствующих определяемому. Метод используют также в автоматизированных системах анализа. В настоящее время выпускают иономеры, в которых отсчет активности (обычно p-функция) любого иона можно производить непосредственно по шкале измерительного прибора после его настройки по стандартным растворам. При работе с такими приборами нет необходимости в построении калибровочного графика.

3. Метод стандартных добавок основан на введении в анализируемый раствор объемом V_X, содержащий неизвестную концентрацию потенциалопределяющего иона С_X, точного объема V_доб. раствора этого же иона с известной концентрацией С_доб. и вычислении С_X по измеренному значению изменения потенциала индикаторного электрода E по формуле: (5)

Основное достоинство этого метода заключается в том, что все измерения проводятся в присутствии всех компонентов пробы. Поэтому хорошей точности определения можно достичь, даже если существенная часть определяемого иона находится в комплексной форме.

Косвенные методы анализа

Применение этих методов наиболее эффективно в анализе полизарядных ионов. Наибольшее распространение получили методы косвенного определения ионов, образующих с потенциалопределяющим ионом комплексные или труднорастворимые соединения, и методы, где потенциалопределяющий ион используется для индикации конечной точки титрования.

Существующие системы

На современном этапе научного прогресса известно много специализированных приборов для определения кислотности среды. Ведущим разработчиком по странам СНГ в этой области является государственное научно-производственное предприятие “Исток-Система” (г.Фрязино Московской обл.).

Среди их разработок можно выделить pH-метры: микропроцессорный ацидогастрометр “Гастротест” и компьютерная система внутрижелудочной рН-метрии “Гастроскан”, а также ацидогастрометром АГМ-01, суточный ацидогастромонитор “Гастроскан-24”, и появилась модификация “Гастроскана” – “Гастроскан-5” [5].

Основные достоинства pH-метров следующие:
          1. Быстрота получения результата.
          2. Точность (средняя точность — 0.05 — 0.1 единиц значения рН).
          3. Возможность постоянного контроля за изменением рН и получения управляющего сигнала. Некоторые могут быть подключены к компьютеру.
         4. Измерения можно проводить в непрозрачных, мутных и окрашенных средах и вязких пастах.
          5. Характерен широкий диапазон измерений.

К основным недостаткам pH-метров следует отнести следующие:
          1. Высокая стоимость.
          2. Селективность основной части электродов не так велика, чтобы производить непосредственное измерение активной концентрации анализируемого иона в любой среде.
         3. Срок службы электрода — главной части такого измерителя примерно около года для непрерывной работы (конечно зависит от типа электрода).

Преимущества индикаторного метода:
1. Низкая стоимость.
2. Быстрота получения результата.

Недостатки индикаторного метода:
          1. Большая погрешность полученного результата (этот метод применяют для приблизительного определения рН).
          2. Метод затруднительно применять при исследовании мутных и окрашенных растворов.
         3. Не универсальность.
         4. Метод не автоматизирован.

Цель работы, практическая ценность, научная новизна

Целью магистерской работы является автоматизация индикаторного метода, что позволит существенно снизить стоимость проведения эксперимента (отпадет необходимость использовать дорогое оборудование (pH-зонды и pH-метры)). В современных компьютерных системах для анализа необходимо предварительно получение экспериментального образца (пробы крови, желудочного сока, желчи и т.д.) или наличие доступной выращенной культуры (на практике возможны случаи, когда исследуемая культура содержится в специальных условиях и попытки измерения в ней pH с помощью pH-метров приводят к нарушению структуры культуры и, соответственно к ее гибели). В разрабатываемой системе в качестве исходных данных для определения кислотности среды используются цветные снимки жидких сред с индикатором в формате bmp, полученные с помощью микроскопа и видеосистемы, что исключает возможность гибели исследуемой культуры. Т.е. научная новизна заключается в том, что разрабатываемая система не будет нуждаться в непосредственном анализе образца, а будет производить исследования с изображением культуры. Это является основным отличием разрабатываемой системы от существующих, что делает ее уникальной, а направление исследования проб по изображению – принципиально новым на сегодняшнем этапе научных разработок в этой области.

Перечень решаемых в работе задач

Основными задачами в процессе программной реализации для решения поставленной задачи являются:
       1. Вычесть “фон” из изображения (среда без индикатора).
       2. Изучить существующие методы повышения качества изображения для решения поставленной задачи, учитывая специфику исходных изображений.
       3. Произвести фильтрацию изображения или удалить артефакты, используя оптимальный метод.

       4. Получить свертки снимка (параметры цветного снимка).
       5. Выбрать и реализовать оптимальную цветовую модель для определения как можно большего числа параметров изображения.
       6. Найти зависимость между pH среды и параметрами изображения, получить уравнение регрессии.

В результате анализа теоретических данных (изучение методов фильтрации изображений, изучение цветовых моделей, регрессионный анализ) было создано программное обеспечение, реализующее следующие функции:
       1. Удаление артефактов (выполняется на основе математического аппарата теории вероятностей).
       2. Вычитание “фона” из изображения.
       3. Повышение качества изображения (фильтрация изображения, манипуляция с яркостью). Реализованы следующие фильтры: усредняющий, пороговый, медианный, Гауссово размывание, линейные фильтры (с использованием курсовых масок наложения, оператора Лапласа, других масок для сглаживания и размытия шумов на изображении).
       4. Реализация цветовых моделей HSB, HLS, HSV и CMYK (исходные изображения представлены в RGB-модели).
       5. Получение параметров изображения для каждой цветовой модели (например, красный, синий и зеленый составляющие цвета в модели RGB, яркость, оттенок цвета, освещенность, насыщенность в HLS-модели).
       6. Построение характеристических гистограмм для всех параметров изображения во всех цветовых моделях.
       7. Получение параметров гистограмм (координаты пика, площадь изображения без артефактов, спектральная интенсивность).
       8. Анализ полученных результатов: нахождение зависимостей между pH среды и результатами, отбор параметров изображения, удовлетворяющих некоторой зависимости.
       9. Нахождение математических зависимостей между кислотностью и отобранными параметрами посредством нахождения уравнений линейной и квадратичной регрессии, определение коэффициентов корреляции.
       10. Построение графиков полученных зависимостей для удобства восприятия.

Алгоритмы и методы

Удаление фона

Фон – среда без индикатора. Удаление фона из изображения сводится к удалению тех пикселей, значения цвета которых в исследуемом изображении и изображении фона одинаковы.

Манипуляции с контрастностью

Разработанное программное обеспечение выполняет манипуляции с контрастностью путем “равнопроцентного” изменения красной, зеленой и синей составляющих RGB-цвета.

Удаление артефактов

          Удаление артефактов реализуется с помощью математического аппарата теории вероятностей:
          На первом этапе обработки изображения производится его разложение по шкале Grayscale, т.е. происходит преобразование цветного изображения в серое, в котором градации серого соответствуют уровню ярости.
          Исходное изображение bmp представлено в модели RGB. Известно, что цветовая модель RGB — аддитивная модель. В ней цвета воспроизводятся путём сложения основных цветов красного (R), зелёного (G) и синего (B). Тогда перевод в серое изображение происходит по формуле:
       Яркость = A*Red+B*Green+C*Blue   (6), где Red, Green и Blue – красная, зеленая и синяя составляющие цвета в диапазоне [0..255], коэффициенты A,B,C определяют долю содержания каждой составляющей цвета и рассчитываются по формулам:
       A = Red/( Red + Green + Blue)
       B = Green/( Red + Green + Blue)
       C = Blue/( Red + Green + Blue)              (7)
       для каждой точки изображения.
       На втором этапе обработки производится расчет средневзвешенного значения яркости и СКО по формуле:

, (8)

      где — средневзвешенное значение яркости;
       b – яркость текущей точки; 0<=b<=L-1;
       P(b) = P{Y[i,j]=b} – вероятность того, что [i,j] –я точка имеет яркость b.
        ,где N(b) – число пикселей, имеющих яркость b.                       (9)
       M – общее число пикселей.
       СКО рассчитывается по формуле: = .                             (10)
       На третьем этапе рассчитываются нижний и верхний предел значения яркости:
   PredeL1 = -                                                                             (11)

PredeL2 = + (12)

и оставляются пикселей, имеющих яркость: PredeL1 < b < PredeL2 (13)
Остальные пиксели удаляются.

Вторичное удаление артефактов происходит по такому же принципу, но без учета уже удаленных точек и призвано для повышения качества удаления.

Фильтрация изображения

Усредняющий фильтр

Т.к. исходные изображения содержат импульсный шум, т.е. значения цвета некоторых пикселей искажаются настолько, что теряют всякую информацию о их неискаженных значениях, актуально произвести фильтрацию изображения.

В программе реализованы следующие методы фильтрации изображения: линейные фильтры, основанные на использовании масок, служащих для сглаживания и размывания изображения, усредняющий фильтр, пороговый, гауссово размывание.

Большое количество операций в области обработки изображений можно представить в виде линейной системы:

(14)

где f (x,y) и h(x,y) — соответственно, начальное и конечное изображение, а g(x,y) — импульсная реакция системы. Другими словами, g(x,y) — это операция над f(x,y), которая создает h(x,y), Для такой системы результат h(x,y) есть свертка (конволюция) f(x,y) с импульсной реакцией g(x,y), определенной в дискретных терминах для изображения NxM:

(15)

Если f и h — изображения, то свертка становится вычислением взвешенных сумм пикселей изображения. Это вычисление делается с учетом сколь угодно большой квадратной таблицы значений под названием "маска свертки". Это вычисление и есть то, что делает фильтрация.

Один из самых простых фильтров реализуется операцией локального усреднения, где значение каждого пикселя заменяется на среднее всех значений в некоторой окрестности, определенной при помощи маски свертки. Например, возьмем окрестность 3x3 для пикселя (i,j), получим:

(16)

Если g[i,j] = 1/9 для каждого [i,j] в маске свертки, то операция свертывания упрощается до локального усреднения сетки из 3х3 пикселей с центром в [i,j]. Заметим, что для 9 пикселей, представленных в этой операции, сумма весов равна 1 (9*1/9=1). Когда для фильтра усреднения используется маска свертки NxN, значение N контролирует "количество" фильтрации. С его увеличением уменьшается количество шума на изображении, но вместе с тем теряются и детали. Так что нужно искать компромисс между выбором конкретного значения N, а выбор маски свертки будет зависеть от нужного уровня фильтрации и конечной детализации изображения.

Пороговый фильтр

Пороговый фильтр является модификацией усредняющего, и отличие заключается в том, что замена значения яркости на среднее производится только в том случае, если разность между значением яркости и полученным средним превышает установленный порог.

Линейный фильтр с использованием масок наложения

Наиболее простые алгоритмы подавления шума на изображении – линейные. Это значит, что каждый пиксел обработанного изображения получается линейной комбинацией нескольких пикселей исходного изображения. Линейная комбинация любых величин представляет собой сумму этих величин, каждая из которых умножена на свой постоянный коэффициент. Таким образом, новое значение пиксела J определяется через старые значения нескольких соседних пикселей I:

(17)

Коэффициенты называются весами и не зависят от положения на изображении точки (i,j) и значений пикселей . Деление на сумму всех весов необходимо для получения единичного коэффициента передачи, чтобы процедура подавления шума не вызывала смещения средней яркости обработанного изображения. Область, в которой определены коэффициенты и в которой производится суммирование, называется маской или апертурой. Обычно маска имеет форму квадрата, и весовые коэффициенты записываются в виде матрицы. Например, для часто встречающейся квадратной маски 3 на 3 пикселя матрица весовых коэффициентов имеет вид:

(18)

В моей работе используются фильтры, сглаживающие шум с маской размером 3х3:

(19)

Гауссово размывание

Эмпирическим путем установлено, что наиболее эффективным методом фильтрации, учитывая специфику исходных изображений, является гауссово размывание.

Фильтр Гаусса похож на усредняющий фильтр, но значения в маске выбираются согласно форме функции Гаусса. Например, функция Гаусса с нулевым средним в одном измерении выглядит так:

, (20)

где дисперсия (разброс функции) определяет ее ширину.

В обработке изображений используется двумерная дискретная функция Гаусса с нулевым средним:

(21)

Она служит для сглаживания. Фильтрация при помощи функции Гаусса обладает определенными преимуществами, которые делают ее особенно удобной именно для сглаживания. При вычислении функции Гаусса выполняются манипуляции с параметром (установленное значение дисперсии обычно намного больше вычисленного, что увеличивает область и “силу” размывания изображения).

Во-первых, функция Гаусса симметрична. Другими словами, она не имеет никаких "приоритетов по движению" во время сглаживания, что особенно удобно, когда обрабатываемые зоны ориентированы произвольным (или заранее неизвестным) образом, и нет причин сглаживать в каком-то определенном направлении.

Во-вторых, функция Гаусса имеет только одну выпуклость, поэтому фильтр на ее основе заменяет каждый пиксель взвешенным средним значением соседних пикселей (как и в случае усредняющего фильтра), поэтому вес пикселя монотонно уменьшается с удалением от центрального пикселя.

Фильтр центрируется на одном пикселе (i,j). Последний модифицируется при помощи: 1) Умножения каждого окружающего пикселя, включая и центральный, на его соответствующий вес из фильтра и суммирования всех результатов 2) Деления суммы, полученной в шаге 1, на сумму весов из фильтра. В результате получим новое значение для пикселя (i,j). Это позволяет локальным деталям остаться в фильтрованном изображении. Наконец, ширина (и таким образом, уровень сглаживания) напрямую связывается с s, так что если s увеличивается, увеличивается и степень сглаживания. Можно контролировать этот параметр для балансирования между количественной составляющей сглаживания и размыванием результирующего изображения.

Оптимальная цветовая модель HSB

Т.к. магистерская работа основана на непосредственной работе с цветом, то для получения параметров изображения целесообразно использовать цветовую модель HSB, поскольку она имеет цветовой обхват, значительно больший, чем модель RGB. (Диапазон цветов, который может быть воспроизведен, зафиксирован или описан каким-либо способом, называется цветовым охватом).

Трехканальная модель HSB получила название по первым буквам английских слов Hue (Оттенок), Saturation (Насыщенность) и Brightness (Яркость). Это очень понятная и удобная в работе модель. Она, теоретически, описывает как аддитивные, так и субтрактивные цвета. Модель HSB основана на модели RGB, но имеет другую систему координат. Любой цвет в модели HSB определяется своим цветовым тоном (оттенком), насыщенностью (то есть процентом добавленной к цвету белой краски) и яркостью (процентом добавленной черной краски).

Цветовой тон характеризует положение энергетического максимума излучения в спектре. Все цветовые тона часто графически изображают в виде круга. Красный цвет принят за нулевое значение, положение остальных цветов характеризуется величиной угла между данным тоном и красным и может изменяться в пределах от 0 до 360. Для чистых спектральных цветов цветовой тон (Hue) является исчерпывающей характеристикой. Однако цвет может быть осветлен (уменьшена его насыщенность) или затемнен (уменьшена яркость).

Насыщенность является мерой чистоты цвета. Чем ближе излучение к монохроматическому, тем выше его насыщенность. Для чистых спектральных цветов параметр Saturation имеет максимальное значение. Насыщенность может изменяться от максимума (100 или 255 градаций) до нуля. Чем меньше насыщенность, тем светлее цвет. При нулевой насыщенности чистый спектральный цвет становится белым.

Яркость (Brightness, Intensity, Luminance) — это объективный (измеряемый) параметр излучаемого цвета, определяющий освещенность или затемненность цвета. Яркость характеризует энергию излучения. Его субъективный аналог это светлота. Все цвета рассмотренного выше цветового круга имеют максимальную яркость (100%), и ярче уже быть не могут. Как и в случае с насыщенностью, остается только уменьшать яркость до минимума (0%), чтобы получить черный цвет. Уменьшение яркости цвета означает зачернение этого цвета.

Поскольку в основе модели находится цветовой круг; в котором начало и конец совпадают, цветовое пространство HSB удобнее всего представить в виде цилиндра (Рисунок 1). Спектральные цвета расположены по верхнему радиусу. К центру круга убывает насыщенность цвета. По высоте цилиндра убывает яркость цветов. Нижняя плоскость цилиндра — черная. Каждый срез цилиндра — это спектральный круг с уменьшающейся яркостью. Серая шкала — линия, соединяющая середины верхнего и нижнего оснований цилиндра.

Рисунок 1. Цветовая модель HSB

Параметры цвета

Основными “цветовыми” параметрами изображения являются: RGB-составляющие цвета, оттенок, насыщенность и яркость в модели HSB, их средневзвешенные значения.

Уравнение линейной регрессии

Регрессионный анализ позволяет приближенно определить форму связи между результативным и факторными признаками, а также решить вопрос о том, значима ли эта связь. Вид функции, с помощью которой приближенно выражается форма связи, выбирают заранее, исходя из содержательных соображений или визуального анализа данных. Математическое решение задачи основано на методе наименьших квадратов. Для нахождения уравнения линейной регрессии использовалась следующая формула:

, (23)

где

(24) (25)

В этом уравнении параметр — свободный член; графически он представляет отрезок ординаты (у) в системе прямоугольных координат. Этот коэффициент, называемый коэффициентом регрессии, показывает, на какую величину в среднем изменяется результативный признак у при изменении факторного признака х на единицу. Параметр называется коэффициентом регрессии. С точки зрения аналитической геометрии — угловой коэффициент, определяющий наклон линии регрессии по отношению к осям координат,
— дисперсия случайной величины x, определяется:

, (26)
где n — число наблюдений.

В качестве независимой переменной x выступает значение pH, функцией y(x) является один из параметров изображения.

Заключение

В процессе выполнения магистерской работы получены уравнения регрессии между pH среды и параметрами изображения, что соответствует поставленным задачам. Например, уравнение регрессии между модой синей составляющей цвета модели RGB, полученной по гистограмме цветов, и pH среды представлено на рисунке 2.

Рисунок 2. Уравнение регрессии между модой синей составляющей цвета модели RGB и pH среды (получено с использованием Delphi 6)

Следующей задачей является – выбор из набора полученных зависимости с максимальным уровнем корреляции и минимальной дисперсией (минимальным разбросом точек из статистической выборки параметров (набор значений параметров определенного типа)) от графика соответствующего уравнения регрессии. Для получения достоверных данных также необходимо выполнить ряд экспериментальных исследований:

       1. Провести ряд экспериментов с целью увеличения объема статистической выборки для подтверждения и уточнения полученных уравнений регрессии, построенных на основе минимального вариационного ряда.
       2. С целью уточнения регрессионных прямых этого необходимо проанализировать фрагменты исходных изображений, не содержащие артефакты, получить для них уравнения регрессии, ввести поправочные коэффициенты.
       3. Провести ряд экспериментов с целью выявления параметров и условий проведения видеосъемки среды, оказывающих влияние на параметры изображения (толщина пробы, предметного стекла, освещенность и т.д.), выполнить анализ этого влияния.

Литература

Конспект лекций “Теория вероятностей и математическая статистика”, Махно С.Я., 2002.
Конспект лекций “Обработка сигналов и изображений”, Скобцов Ю.А., 2003.
Статья. Химические основы лабораторных технологий. Кислотность среды, измерение кислотности среды, индикаторы. http://www.labinfo.ru/metod/basa/002.htm

4. Статья. Химические основы лабораторных технологий. Ионоселективные электроды. Ионометрия. http://www.labinfo.ru/metod/basa/003.htm

5. Основной сайт предприятия "Исток-Система" – основного российского разработчика и производителя приборов для исследования кислотности. http://www.i-s.ru/

    6. Статья. Базовая обработка изображения и сглаживание высот. http://www.gamasutra.com/resource_guide/20010716/chen_01.htm
    7. Ковальченко И.Д., Бородкин Л.И. и др. Линейная корреляция.  http://www.yartel.ru/stat/vzaim2.html
    8. Стефан Стефанов. Цвет и цветовоспроизведение в полиграфии. http://www.aqualon.ru/
    9. Шавелева Л.И. Работа в PageMaker. Ч. V. Цвет. Гл. 19. Определение цветов.  http://graph.power.nstu.ru