Привет, Хабр!
Идея этого поста родилась при работе с геммологами, которые занимаются оценкой качества драгоценных изделий и драгоценных камней, оценкой стоимости и выявлением подделок. Исследования были направлены на то, чтобы выяснить, что может рассказать о внутренней структуре драгоценностей наша программа Smart Tomo Engine (ознакомиться с которой можно у нас на сайте Smart Engines).
Существует как минимум три типа жемчуга: искусственный, культивированный и натуральный. Искусственный жемчуг - жемчуг действительно искусственно созданный человеком. В его появлении моллюск не участвует. В производстве искусственного жемчуга принимает участие только человек от начала и до конца. Шарики из пластмассы или стекла покрывают многократно различными перламутровыми красками. Искусственный жемчуг не может содержать такое разнообразие цветовых оттенков, как натуральный жемчуг, такой жемчуг гораздо дешевле. Его геммологи могут без труда определить. Процесс образования природного и культивированного жемчуга практически одинаков. Отличие в случае образования природного жемчуга инородным телом являются мелкие паразиты, а для образования культивированного жемчуга человек сам начинает процесс созревания жемчужины, вставляя инородное тело в раковину моллюска. Стоимость культивированного жемчуга намного ниже, но вот определить подделку не заглядывая внутрь изделия уже трудно. Поэтому первой задачей, которая стояла перед нами, было определение размера ядра культивированного жемчуга с помощью метода томографии.
Выполненное на нашем микротомографе исследование позволило построить трёхмерную модель жемчужины.
![Трёхмерная модель культивированного жемчуга, выращенного ядерным способом. Трёхмерная модель культивированного жемчуга, выращенного ядерным способом.](https://webcf.waybackmachine.org/web/20240322003949/https://habrastorage.org/getpro/habr/upload_files/63b/fd4/d82/63bfd4d824dfd8891d2c8203ebae43d1.gif)
Но уже просто увидев проекции жемчужины мы подумали, а нельзя ли измерить размер ядра просто по проекциям, ведь визуально мы видим окружность. И действительно радиус удалось измерить “на глазок”, и он оказался равным 6.08 мм.
![Определение радиуса ядра искусственного жемчуга по проекциям (Радиус сферы, мм: 6.08) Определение радиуса ядра искусственного жемчуга по проекциям (Радиус сферы, мм: 6.08)](https://webcf.waybackmachine.org/web/20240322003949/https://habrastorage.org/getpro/habr/upload_files/1ae/423/733/1ae423733b88cfc56eb449d2b32575fa.png)
Но что нам скажут трёхмерные алгоритмы сегментации? Для этого нами был реализован следующий пайплайн сегментации и измерения линейных размеров жемчужины:
Начальное приближение для первого активного контура
Итеративный поиск контуров с инициализацией по соседнему слою
Выбор диапазона слоев по визуальной оценке качества
Подбор 4-х параметров сферы (радиус, координаты центра) методом градиентного спуска по точкам, полученным в п.3
![Пайплайн метода определения границ ядра жемчужины Пайплайн метода определения границ ядра жемчужины](https://webcf.waybackmachine.org/web/20240322003949/https://habrastorage.org/getpro/habr/upload_files/d88/dd0/2f7/d88dd02f74d2bad86afcbda7b2305012.png)
Разработанный нами метод позволил вычислить параметры жемчужины:
Радиус ядра жемчужины, мм: 6.12
Суммарный объем жемчужины, мм3: 1777
Суммарный объем ядра, мм3: 962
Процентное соотношение ядра к объему: 54.11%
![Трёхмерная модель жемчужины Трёхмерная модель жемчужины](https://webcf.waybackmachine.org/web/20240322003949/https://habrastorage.org/getpro/habr/upload_files/d9a/10d/ec7/d9a10dec7fbfadbc2630e2732aad3928.png)
![Результаты определения границ ядра жемчужины по томографической реконструкции. Результаты определения границ ядра жемчужины по томографической реконструкции.](https://webcf.waybackmachine.org/web/20240322003949/https://habrastorage.org/getpro/habr/upload_files/51d/908/e93/51d908e93f8e1d2af7056afb455e7286.gif)
Результаты сегментации можно посмотреть на гифке. Так же на правой панели гифки для сравнения приведен результат работы классического детектора границ Канни.
Полученные результаты показали, что в случае простых объектов можно вычислять их размеры просто по рентгеновским изображениям, тем более что, например, жемчуг может изменить окрас и потемнеть при долгом воздействии рентгена на него.
А как быть с более сложными несимметричными камнями? Вот например у нас есть сапфир, трещины в котором были реставрированы методом термообработки с заполнением трещин свинцовым стеклом. Как посчитать сколько там стекла? Вот тут без томографии не обойтись.
![Рентгеновская проекция корунда (сапфира) Рентгеновская проекция корунда (сапфира)](https://webcf.waybackmachine.org/web/20240322003949/https://habrastorage.org/getpro/habr/upload_files/40c/7ad/224/40c7ad22436dc9ee58f8b1b0421a65f8.png)
![Корунд (сапфир), облагороженный методом термообработки с заполнением трещин свинцовым стеклом. Корунд (сапфир), облагороженный методом термообработки с заполнением трещин свинцовым стеклом.](https://webcf.waybackmachine.org/web/20240322003949/https://habrastorage.org/getpro/habr/upload_files/a40/5fa/b8e/a405fab8e61ef092653ceebd8a663a0d.gif)
Нами с помощью инструментов Smart Tomo Engine была выполнена томографическая реконструкция, проведена чистка от шумов и сегментация. Проведенное исследование показало:
Суммарный объем камня, мм3: 454
Суммарный объем трещин, мм3: 69
Процентное соотношение трещин к объему: 15.18%
![Пайплайн метода сегментации Пайплайн метода сегментации](https://webcf.waybackmachine.org/web/20240322003949/https://habrastorage.org/getpro/habr/upload_files/5ff/8db/a56/5ff8dba56f45e467b0890843c78cb17c.png)
![Результаты сегментации Результаты сегментации](https://webcf.waybackmachine.org/web/20240322003949/https://habrastorage.org/getpro/habr/upload_files/e43/f01/fa7/e43f01fa7ad0d52617d5c906da782d34.gif)
А какой вывод можно сделать из нашей статьи? Он прост. Микротомографию можно использовать для создания “цифрового паспорта” ювелирных изделий. Причем вопросов к алгоритмам реконструкции уже нет (как мы показали, алгоритмически можно точно определить все особенности драгоценных камней). Проблема, как это часто бывает, остается только лишь в юридическом (законодательном) русле.