Форум программистов, компьютерный форум, киберфорум
Python: Научные вычисления
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 5.00/5: Рейтинг темы: голосов - 5, средняя оценка - 5.00
2 / 2 / 1
Регистрация: 16.04.2022
Сообщений: 1,032

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

19.07.2023, 08:35. Показов 1140. Ответов 5
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Здравствуйте!

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

Написал такой код.

Кликните здесь для просмотра всего текста
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
import numpy as np
 
height = 300
width = 600
 
img_origin=np.random.randint(0,255, (height, width, 3), np.uint8)            
extendedImage = np.expand_dims(img_origin, axis = 0)  
 
heightFragment = 50
widthFragment = 100
 
def createNewLayer(startY, startX):
    assert heightFragment + startY <= height
    assert widthFragment + startX <= width
    
    fragment = np.random.randint(0,255, (heightFragment, widthFragment, 3), np.uint8)
    
    newLayer = np.full([height, width, 3], np.nan)
 
    for i in range (heightFragment):
        for j in range (widthFragment):
            for k in range(3):
                newLayer[startY + i, startX + j, k] = fragment[i,j,k]
            
    return np.expand_dims(newLayer, axis = 0)
 
#-----На этом месте зависает!!!
for i in range(250):
    for j in range(500):
        extendedImage = np.concatenate((extendedImage, createNewLayer(startY = i, startX = j)), axis=0)
#-----
        
result = np.nanmean(extendedImage,axis=0)

В выделенном блоке кода зависает и долго висит.

Скажите, как можно оптимальнее написать средствами numpy?
0
Programming
Эксперт
39485 / 9562 / 3019
Регистрация: 12.04.2006
Сообщений: 41,671
Блог
19.07.2023, 08:35
Ответы с готовыми решениями:

Несколько изображений в одно изображение
Есть допустим несколько изображений, по высоте обычно одинаковые, по ширине разные. Хочу сделать так что бы было одно узкое изображение,...

Склейка выделенных диапазонов массива изображений в одно изображение
Добрый день,Уважаемые программисты.Есть такая проблема:bз видеофайла получаю массив bmp-кадров (сохранились в папке Windows). Далее по...

Цикл вывода изображений выводит только одно изображение
Почему цикл выводит только, первый рисунок из папки? $dir=opendir('papkasfotkami/'); while(($e=readdir($dir))!==false){ ...

5
2649 / 1660 / 267
Регистрация: 19.02.2010
Сообщений: 4,392
19.07.2023, 08:55
Цитата Сообщение от MConst Посмотреть сообщение
В выделенном блоке кода зависает и долго висит.
Там выполняются вызовы двух функций (np.concatenate() и твоя функция createNewLayer()) - разберись, какая из них висит.
Если вторая - то ты неправильно указал кривой блок кода
0
2 / 2 / 1
Регистрация: 16.04.2022
Сообщений: 1,032
19.07.2023, 17:16  [ТС]
Сделал измерения. Функция np.concatenate после каждого слоя увеличивает время работы.
0
2 / 2 / 1
Регистрация: 16.04.2022
Сообщений: 1,032
20.07.2023, 13:11  [ТС]
Удалось пока так оптимизировать
Кликните здесь для просмотра всего текста
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
import numpy as np
 
height, width, heightFragment, widthFragment, startY, startX, endY, endX   = 300, 600, 50, 100, 0, 0, 250, 1#500
 
assert heightFragment + endY <= height
assert widthFragment + endX <= width
 
img_origin=np.random.randint(0,255, (height, width, 3), np.uint8)            
 
listImages = []
listImages.append(img_origin)
 
for i in range(startY, endY, 1):
    for j in range(startX, endX, 1):
        fragment = np.random.randint(0,255, (heightFragment, widthFragment, 3), np.uint8)        
        newLayer = np.full([height, width, 3], np.nan)       
        newLayer[i:i+heightFragment, j:j+widthFragment] = fragment        
        listImages.append(newLayer)
           
result = np.nanmean(np.array(listImages),axis=0)

Если увеличивать цикл, то долго считает.

Как можно еще оптимизировать тут? Как циклы с данным расчетом можно поместить в numpy?
0
5516 / 2869 / 571
Регистрация: 07.11.2019
Сообщений: 4,760
22.07.2023, 06:52
MConst, np.nan это число с плавающей точкой, а fragment uint8. При сложении в цикле постоянно происходит конвертация.
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
import numpy as np
 
height = 300
width = 600
 
img_origin=np.random.randint(0,255, (height, width, 3), np.uint8)            
img=np.float64(img_origin)  
 
heightFragment = 50
widthFragment = 100
divisor_matrix=np.ones((height, width), dtype=np.int64)
 
for i in range(250):
    for j in range(500):
        fragment = np.random.randint(0,255, (heightFragment, widthFragment, 3))
        img[i:i+heightFragment, j:j+widthFragment]+=fragment
        divisor_matrix[i:i+heightFragment, j:j+widthFragment]+=1
     
result = np.uint8(np.divide(img, np.expand_dims(divisor_matrix, axis=-1)))
Можно и далее оптимизировать...
1
2 / 2 / 1
Регистрация: 16.04.2022
Сообщений: 1,032
22.07.2023, 13:59  [ТС]
Спасибо! Интересное и хитрое решение обойти np.nanmean, а именно конвертацию и при этом не создавать большой список с последующей обработкой, который может нагружать память компьютера.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
inter-admin
Эксперт
29715 / 6470 / 2152
Регистрация: 06.03.2009
Сообщений: 28,500
Блог
22.07.2023, 13:59
Помогаю со студенческими работами здесь

Получить 3d изображение из 3000 изображений
Возникла такая задача,что необходимо получить 3d изображение из 3000 файлов в matlab. Сами изображения постаралась загрузить в Matlab: ...

TCPDF ERROR: [Изображение] Не удалось получить изображение
Я использую TCPDF для создания PDF и заканчивая следующим сообщением об ошибке, когда я нажимаю ссылку для создания PDF: TCPDF ERROR: ...

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

Исчезло изображение слоёв в панели Слои
Хелп! Исчезло изображение в панеле слои и вместо него кисточки.. Как вернуть?!

Работа с dds файлами: получить анимированное изображение, достав из этого файлы все его части
Есть dds файл, в неё 160 маленьких изображений. К нему идет описание типа: 0 4 &quot;01微笑&quot; 66 72 78 84 Где 0 - номер...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
6
Ответ Создать тему
Новые блоги и статьи
BOINC: 22 года — и всё ещё работает
Programma_Boinc 12.03.2026
BOINC: 22 года — и всё ещё работает Дэвид Андерсон написал ретроспективу. Кратко: в 2001 году он ушёл из United Devices, где был CTO, и за несколько месяцев написал ядро BOINC — клиент, сервер,. . .
SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++
8Observer8 10.03.2026
Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:
Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE
8Observer8 10.03.2026
Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .
Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки
Neotwalker 08.03.2026
Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .
Модульная разработка через nuget packages
DevAlt 07.03.2026
Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .
Модульный подход на примере F#
DevAlt 06.03.2026
В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .
Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование
8Observer8 05.03.2026
Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .
SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D
8Observer8 04.03.2026
Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru