Tutaj prezentujemy protokół oparty na komputerowym systemie wizyjnym (CVS) do określenia zachowania podczas topnienia wielofazowych systemów żywnościowych.
Artykuł metodologiczny
Tutaj prezentujemy protokół oparty na komputerowym systemie wizyjnym (CVS) do określenia zachowania podczas topnienia wielofazowych systemów żywnościowych.
Lody są systemem wielofazowym, w którym fazy ciekła, stała i gazowa są ściśle powiązane. Ciągła faza ciekła otacza pęcherzyki powietrza oraz kryształki lodu i zawiera częściowo skrystalizowane tłuszcze, białka koloidalne, sole, cukry (ewentualnie skrystalizowane) oraz stabilizatory. Skład lodów różni się w zależności od wymagań lokalnego rynku i ewentualnych regulacji. Choć technologia przetwarzania wpływa na charakterystykę gotowego produktu, każdy składnik odgrywa istotną rolę w określaniu jakości lodów1. Zachowanie podczas topnienia jest jednym z najważniejszych wskaźników jakości lodów, biorąc pod uwagę zjawiska zachodzące zarówno podczas konsumpcji, jak i w jamie ustnej. Wraz z penetracją ciepła w głąb lodów kryształki lodu topnieją, a woda dyfunduje i miesza się z fazą surowiczą, która może wyciekać przez pozostałą strukturę2. Produkt szybko topniejący jest niepożądany zarówno ze względu na komfort spożycia, jak i gwarancję wyższej odporności na szok termiczny. Jednak produkty topniejące powoli również mogą wskazywać na pewne wady w recepturze1. Wiadomo, że za właściwości topnienia odpowiada mikrostruktura lodów3, lecz do tej pory publikowano sprzeczne wyniki, co pokazuje, że wiedza na temat wpływu czynników mikrostrukturalnych na topnienie jest wciąż ograniczona4. W związku z tym niezbędne są dalsze badania w celu wyjaśnienia mechanizmu topnienia, co jest kluczowe również przy projektowaniu nowych receptur3.
Właściwości topnienia są zazwyczaj oceniane metodami grawimetrycznymi i wyrażane w postaci czasu rozpoczęcia oraz szybkości topnienia5. Określoną porcję lodów umieszcza się na drucianej siatce w komorze z kontrolowaną temperaturą, a następnie rejestruje się wagę stopionego produktu. Na podstawie krzywej zależność wagi od czasu można wyróżnić trzy fazy: fazę opóźnienia, podczas której następuje przenikanie ciepła; fazę szybkiego topnienia, w której rozcieńczona frakcja surowicy przepływa przez strukturę lodów z maksymalną prędkością, oraz fazę stacjonarną, w której większość produktu spłynęła2.
Metodą grawimetryczną można rozpoznać produkty wolno i szybko topniejące; jednak istotny jest również wygląd lodów podczas topnienia, ponieważ zachowanie kształtu wiąże się z dobrym składem i strukturą produktu6. Zatem procedura oparta na systemie wizji komputerowej (CVS) może wspierać istniejącą już metodologię grawimetryczną, umożliwiając badanie wyglądu produktu podczas topnienia. Systemy CVS mogą pozyskiwać liczne atrybuty żywności3 (np. rozmiar, wagę, kształt, teksturę i kolor) z dokładnością, której nie jest w stanie dostrzec ludzkie oko. Takie systemy składają się zazwyczaj z kamer cyfrowych i oprogramowania do przetwarzania obrazu7. W rzeczywistości protokół oparty na CVS obejmuje dwa główne kroki: 1) akwizycję obrazu i 2) przetwarzanie obrazu. Można zastosować różne poziomy przetwarzania obrazu7, od najprostszych po bardziej złożone, takie jak metody głębokiego uczenia w rozwoju sztucznej inteligencji8,9. W ostatnim czasie poświęcono wiele uwagi systemom CVS w sektorze spożywczym i opracowano dużą liczbę zastosowań w zakresie kontroli bezpieczeństwa żywności, monitorowania przetwarzania żywności, wykrywania ciał obcych oraz w innych dziedzinach. Są one szybkie, wydajne i nieniszczące, stanowiąc tym samym wartościowe narzędzia do dostarczania konsumentom bezpiecznej żywności wysokiej jakości10.
W dziedzinie badań nad lodami zaproponowano metodę analizy obrazu do badania rekrystalizacji lodu za pomocą mikroskopii optycznej11. W ostatnim czasie obrazy z tomografii komputerowej rentgenowskiej były przetwarzane w celu analizy trójwymiarowej mikrostruktury materiałów miękkich i porowatych, w tym lodów3. Jednak opracowanie prostych cyfrowych obrazów z urządzeń CCD (Charge Coupled Device) może oferować szereg zalet pod kątem łatwości pozyskiwania danych oraz oddania wyglądu lodów w taki sposób, w jaki są one postrzegane przez konsumentów. Niektórzy autorzy prezentują obrazy lodów podczas topnienia12, lecz według naszej najlepszej wiedzy, wyodrębnianie wskaźników liczbowych z obrazów zostało po raz pierwszy opisane przez Moriano i Alamprese13.
W związku z powyższym, proponowany tutaj protokół, oparty na pracy Moriano i Alamprese13, przedstawia prosty system CVS, który może zostać zastosowany w celu wsparcia istniejącej już metodologii grawitacyjnej do badania zachowania lodów podczas topnienia. Schemat blokowy proponowanej metody przedstawiono na Rysunku 1. Zastosowanie takiego systemu umożliwia obliczenie dwóch wskaźników topnienia związanych z utrzymaniem kształtu oraz szybkością topnienia. W szczególności w pracy po raz pierwszy szczegółowo opisano konfigurację eksperymentalną i procedurę cyfrowej akwizycji obrazów podczas topnienia lodów oraz etapy przetwarzania obrazów. Ponadto przedstawiono wyniki uzyskane dla lodów produkowanych z różnymi słodzikami (tj. sacharozą, sukromaltem i erytrytolem), aby wykazać potencjał tej metody.

Rycina 1: Schemat blokowy proponowanych metodologii. Podsumowanie ogólnych etapów dla proponowanego systemu wizji komputerowej oraz metody grawimetrycznej do badania zachowania lodów podczas topnienia. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
1. Układ doświadczalny i procedury dla prób topnienia

Rycina 2: Układ do testu topnienia. Na rysunku przedstawiono sposób przygotowania próby topnienia w szafie termostatycznej: na wadze cyfrowej ustawiono miarkowany cylinder szklany w celu zbierania i ważenia roztopionych lodów. Próbkę lodów umieszczono na metalowym sitcie z drucianej siatki w lejku zawieszonym nad cylindrem. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
2. Przetwarzanie obrazów w celu obliczenia indeksów topnienia
(1)
(2)
(3)Rycina uzupełniająca 1: Przestrzenna kalibracja obrazu. (A) Przejdź do okna Measure > Calibration > Spatial w oprogramowaniu do analizy obrazu. Wybierz New, a następnie zaznacz Image, aby otworzyć okno Scaling. Wskaż długość referencyjną w jednostce, na którą mają zostać przeliczone piksele (np. milimetry). (B>) Ostrożnie nałóż zielony pasek na fragment referencyjny odpowiadający wskazanej długości i kliknij OK. Kliknij tutaj, aby pobrać ten plik.
Rysunek uzupełniający 2: Wycinanie AOI i konwersja na skalę szarości. (A) Konwersja obszaru zainteresowania (AOI) na skalę szarości oraz (B) obraz wynikowy. Kliknij tutaj, aby pobrać ten plik.
Rycina uzupełniająca 3: Wybór parametrów do pomiaru. W oknie Select Measurement można wybrać parametry do pomiaru; w przypadku oceny topnienia lodów należy wybrać pole (area), szerokość ramki (box width) oraz wysokość ramki (box height). Proszę kliknąć tutaj, aby pobrać ten plik.
Rycina uzupełniająca 4: Segmentacja próbki lodów.W oknie „Segmentation” możliwe jest wybranie zakresów histogramu tak, aby objął on dokładnie obszar kształtu lodów. Kliknij tutaj, aby pobrać ten plik.
Rysunek uzupełniający 5: Filtrowanie obiektów i funkcja Count.Czerwone linie zaznaczają rozpoznane jasne obiekty. Po zastosowaniu funkcji count i otwarciu okna „View, Measurement data” zostaną wyświetlone wyniki dla wybranych parametrów (A). Aby odfiltrować tylko kształt loda, w oknie „Select Measurement” można wybrać zakres minimalnej i maksymalnej powierzchni, dzięki czemu policzone zostaną parametry tylko jednego obiektu (B). Kliknij tutaj, aby pobrać ten plik.
Rycina uzupełniająca 6: Indeks zachowania kształtu (R). Przedstawiono wysokość pudełka (Ht, czerwona przerywana linia) oraz szerokość pudełka (Wt, czarna ciągła linia) wykorzystywane do obliczenia indeksu zachowania kształtu (R). Proszę kliknąć tutaj, aby pobrać ten plik.

Rysunek 3: Krzywe retencji kształtu i pola. Przykład krzywych retencji kształtu (A) i pola (B) dla lodów, na których naniesiono średnie wartości Rt/R0 oraz At/A0 w funkcji czasu; słupki błędów odpowiadają wartościom odchylenia standardowego uzyskanych w powtórzeniach analizy. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.
3. Opracowanie danych grawimetrycznych
(4)
(5)
Rysunek 4: Krzywa grawimetryczna. Przykład krzywej topnienia lodów uzyskanej metodą grawimetryczną. Oryginalna krzywa jest przedstawiona na czerwono; wybrany szereg danych w części liniowej zaznaczono na zielono; wyliczona linia regresji jest przedstawiona czarnymi kropkami. Przedstawiono również równanie oraz współczynnik determinacji (R2) linii regresji. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
UWAGA: Aby uzyskać wiarygodne wyniki do analizy statystycznej, należy powtórzyć całą procedurę próby topnienia i przetwarzania obrazów co najmniej trzy razy dla każdej próbki.
Jako przykład proponowanych wyników CVS przedstawiono analizy roztopu dla trzech różnych formulacji lodów, porównane z danymi uzyskanymi metodą wagową. W szczególności zbadano zachowanie podczas topnienia lodów wykonanych z wykorzystaniem różnych słodzików (czyli sacharozy, sukromaltu i erytrytolu).
Tabela 1 oraz Rycina 5A przedstawiają wyniki indeksu zachowania kształtu (Rt/R0) dla trzech próbek lodów podczas topnienia. Do 45 minuty zaobserwowano podobne zachowanie kształtu we wszystkich próbkach; następnie próbka zawierająca sacharozę traciła kształt szybciej niż pozostałe produkty, co skutkowało istotnie niższym wskaźnikiem Rt/R0 zgodnie z analizą wariancji jednoczynnikowej (ANOVA) oraz testem najmniejszej istotnej różnicy (LSD) (p < 0,05). Tylko próbka zawierająca erythrytol zachowała swój kształt do końca badania, dając mierzalne wartości Rt/R0.
| SUKROZA | SUCROMALT | ERYTRYTOL | ||||
| Czas (min) | At/A0 | s.d. | At/A0 | s.d. | At/A0 | s.d. |
| 0 | 1.00 | - | 1.00 | - | 1.00 | - |
| 15 | 0.91a | 0.01 | 0.93a | 0.02 | 0.92a | 0.04 |
| 30 | 0.76a | 0.04 | 0.82a | 0.02 | 0.84a | 0.04 |
| 45 | 0.56a | 0.04 | 0.66ab | 0.02 | 0.74b | 0.03 |
| 60 | 0.32a | 0.01 | 0.46b | 0.02 | 0.62c | 0.01 |
| 75 | 0.13a | 0.01 | 0.24b | 0.02 | 0.52c | 0.02 |
| 90 | n.d. | - | n.d. | - | 0.40 | 0.01 |
| n.d., nie wykrywalne | ||||||
Tabela 1: Wartości średnie i odchylenia standardowe wskaźników zachowania kształtu dla trzech próbek lodów wytwarzanych z wykorzystaniem różnych słodzików. Różne litery małe w tym samym wierszu wskazują na istotnie różne wartości średnie (p < 0,05).

Rycina 5: Krzywe utrzymywania kształtu i powierzchni dla trzech próbek lodów przygotowanych z wykorzystaniem różnych słodzików. (A) Krzywe utrzymywania kształtu oraz (B) krzywe utrzymywania powierzchni. Różowa linia reprezentuje próbkę z sacharozą, traktowaną jako odniesienie; żółta linia – próbkę z sukromaltą, a niebieska linia – próbkę z erytrytolem. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.
Próbka z erytrytolem wykazała również niższą szybkość topnienia, o czym świadczą wskaźniki topnienia At/A0 (Tabela 2 oraz Rysunek 5B). Od 45 minuty próbka z erytrytolem wykazywała istotnie (p < 0,05) największe zachowanie powierzchni, a po niej lody z zawartością sukromaltu i sacharozy.
| SACHAROZA | SUCROMALT | ERYTRYTOL | ||||
| Czas (min) | Rt/R0 | s.d. | Rt/R0 | s.d. | Rt/R0 | s.d. |
| 0 | 1.00 | - | 1.00 | - | 1.00 | - |
| 15 | 0.98a | 0.04 | 0.97a | 0.01 | 0.93a | 0.02 |
| 30 | 0.94a | 0.05 | 0.94a | 0.02 | 0.89a | 0.02 |
| 45 | 0.80a | 0.03 | 0.89a | 0.02 | 0.84a | 0.02 |
| 60 | 0.64a | 0.04 | 0.78b | 0.01 | 0.76b | 0.03 |
| 75 | 0.54a | 0.03 | 0.65a | 0.01 | 0.60a | 0.07 |
| 90 | n.d. | - | n.d. | - | 0.45 | 0.06 |
| n.d., nie wykrywalne | ||||||
Tabela 2: Średnie wartości i odchylenia standardowe indeksów retencji powierzchni dla trzech próbek lodów wytwarzanych z wykorzystaniem różnych słodzików. Różne małe litery w tym samym wierszu wskazują średnie wartości istotnie różne (p < 0,05).
Wyniki potwierdzono również analizą wagową, której dane przedstawiono w Tabeli 3 jako czas rozpoczęcia topnienia (tS) oraz szybkość topnienia (gramy/minuta). Jednak w tym przypadku nie stwierdzono istotnych różnic stosując analizę wariancji ANOVA (p > 0,05).
| SUKROZA | SUKROMALT | ERYTRYTOL | ||||
| średnia | odchylenie standardowe | średnia | odchylenie standardowe | średnia | odchylenie standardowe | |
| ts (min) | 15.3a | 1.4 | 16.4a | 1.4 | 17.2a | 2.6 |
| prędkość topnienia (g/min) | 2.67a | 0.26 | 2.29a | 0.09 | 1.96a | 0.41 |
Tabela 3: Średnie wartości i odchylenia standardowe indeksów topnienia obliczone metodą wagową dla trzech próbek lodów wytwarzanych z wykorzystaniem różnych słodzików. ts, czas rozpoczęcia topnienia. Różne małe litery w tym samym wierszu wskazują istotnie różne średnie wartości (p < 0,05).
Ryc. 6 przedstawia niektóre przetworzone (krok 2.1) obrazy uzyskane podczas prób topnienia (krok 1.3) trzech próbek lodów; można zauważyć dobrą retencję kształtu wszystkich próbek oraz najniższą szybkość topnienia lodów przygotowanych z erytrytolem, które po 90 min nie uległy całkowitemu stopieniu.

Rycina 6: Rozpływ trzech próbek lodów. Obrazy pozyskane co 15 minut podczas rozpływania się trzech próbek lodów przygotowanych z wykorzystaniem różnych słodzików. Przetworzone obrazy trzech próbek lodów podczas topnienia. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.
Zaproponowany system CVS pozwala na obliczenie indeksów utrzymania kształtu i powierzchni próbek lodów podczas topnienia, a ponadto umożliwia wizualizację procesu topnienia. Może on być łączony z tradycyjną metodą grawimetryczną stosowaną do oceny zachowania lodów podczas topnienia5, aby uzyskać wyniki dotyczące aspektu wizualnego produktów. Jest to bardzo istotne, ponieważ konsumenci oceniają ich jakość również na podstawie wyglądu zewnętrznego, a zdolność do zachowania kształtu podczas topnienia jest powiązana ze składem i strukturą lodów1,3. Zastosowanie tego systemu w przemyśle spożywczym może usprawnić kontrolę jakości, a także badania i rozwój nowych receptur.
Przykład wyników uzyskanych dla topnienia trzech próbek lodów z różnymi stosowanymi słodzikami wykazał, że proponowana metoda CVS jest czuła na różnice w składzie (Tabela 1 i Tabela 2), nawet bardziej niż metoda grawimetryczna, która nie była w stanie wykazać istotnych różnic między próbkami (Tabela 3). Czułość indeksu utrzymania kształtu na recepturę lodów została wykazana również w naszych poprzednich pracach13,14,15,16. Ze względu na wysoką czułość proponowanej metody, może ona nawet zastąpić procedurę grawimetryczną.
Poza przygotowaniem próbki lodów, najbardziej krytyczne etapy w opracowanym systemie CVS wiążą się z prawidłowym przechwytywaniem i segmentacją obrazu7. Odniesienie danych Rt i At do wartości początkowych obliczonych w czasie 0 zapobiega możliwym błędnym wynikom wynikającym z różnic w ustawieniu eksperymentów między powtórzeniami. Jednakże błędy paralaksy podczas akwizycji obrazu mogą prowadzić do przeszacowania powierzchni próbki podczas topnienia, szczególnie jeśli obrazy zebrane w jednej próbie nie są spójne z powodu zmiany położenia kamery między kolejnymi ujęciami. Dlatego należy zwrócić szczególną uwagę na ustawienie procedury akwizycji obrazu, upewniając się, że metalowy ekran z siatki z próbką na górze jest idealnie wyrównany z obiektywem kamery, oraz starannie definiując i zabezpieczając pozycję kamery przed szafą termostatyczną. Podobnie, błędna segmentacja obrazu może prowadzić do nadmiernego lub niedoszacowania powierzchni lodów, a także H i W, co skutkuje bezsensownymi wskaźnikami zachowania kształtu. W tym przypadku głównym źródłem błędu są możliwe wypustki na krawędzi porcji lodów umieszczonych na siatce lub pęcherzyki powietrza pojawiające się na powierzchni lodów podczas topnienia (Rycina 7), które muszą zostać usunięte z jasnego kształtu przed pomiarem powierzchni pola oraz wartości W i H pola. Można to rozwiązać, korzystając z funkcji Split Object oprogramowania do analizy obrazu. Wysokokontrastowe tło jest również przydatne dla lepszej segmentacji obrazu7. Kolor tła należy dobrać w zależności od koloru próbki lodów, przy czym pożądane jest wysokie nasycenie barw dla uzyskania lepszego kontrastu.

Rysunek 7: Możliwe źródła błędów w segmentacji kształtu lodów. Czerwone strzałki wskazują możliwe wypustki i pęcherzyki powietrza na kształcie lodów, które mogą prowadzić do błędnego obliczenia pola powierzchni, wysokości i szerokości. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.
Jeśli obrazy są prawidłowo gromadzone i przetwarzane, metoda CVS jest jeszcze bardziej czuła i precyzyjna niż metodologia grawimetryczna. Rzeczywiście, w przedstawionych przykładach współczynniki zmienności dla analiz grawimetrycznych mieściły się w zakresie od 4% do 20%, co sprawiło, że wszystkie średnie dane nie wykazywały istotnych różnic (ANOVA, p > 0,05). Przeciwnie, dla wskaźników CVS wartości współczynnika zmienności mieściły się w zakresie od 0,1% do 13%, wykazując istotne różnice pomiędzy trzema próbkami lodów podczas topnienia (Tabela 1 i Tabela 2). Zmniejszenie szybkości topnienia spowodowane przez erytrytol jest prawdopodobnie związane z tworzeniem się amorficznych kryształów cukru podczas mrożenia lodów oraz niższą rozpuszczalnością tego słodzika w porównaniu z sacharozą i sukromaltem, co przekłada się na większą ilość zamrożonej wody w produkcie końcowym13. Inne udane zastosowania zaproponowanej tutaj metody CVS można znaleźć w naszych poprzednich badaniach13,14,15,16.
Po osiągnięciu odpowiedniej standaryzacji procedury przetwarzania obrazów, automatyczną analizę obrazu można skonfigurować za pomocą funkcji Macro-Recording (nagrywanie makr) w oprogramowaniu do analizy obrazu, co znacznie przyspiesza proces opracowania danych. Oczywiście zawsze wymagana jest dokładna weryfikacja danych przez operatora, aby uniknąć artefaktów i błędnych wyników. Ponadto można przeprowadzić analizę krzywych retencji kształtu i powierzchni w celu obliczenia parametrów kinetycznych oraz zastąpienia metody grawimetrycznej.
Protokół ten można zastosować również do innych produktów spożywczych, dla których istotna jest szybkość topnienia i wskaźnik utrzymania kształtu, takich jak musy i pianki (np. ubita śmietana mleczna i białko jaja kurzego). Najważniejszym aspektem jest zdolność produktów do utrzymania formy, natomiast temperaturę podczas prób topnienia można dostosować do konkretnego produktu.
Autorzy nie mają nic do ujawnienia.
| Nazwa | Firma | Numer katalogowy | Komentarze |
|---|---|---|---|
| Gabinet | Cavallo s.r.l. | FTX700 | Miejsce próby topnienia |
| Aparat cyfrowy | Sony Group Corp | DSC-S650 | |
| Waga cyfrowa | Gibertini Elettronica | EU-C 4002 LCD | |
| ImagePro Plus 7.0 | Media Cybernetics, Inc | N/A | Oprogramowanie do analizy obrazu |
| Microsoft Excel | Microsoft | N/A | Opracowanie danych i grafika |
| Scalecom | Gibertini Elettronica | N/A | Statyw do akwizycji oprogramowania do wagi cyfrowej |
| Manfratto | #055 |