NumPy перерасчет масива

@DarkPower · Регистрация: 20.03.2013

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

У меня тут пару вопросов:
1) Можно ли средствами NumPy перерасчетать массив?
массив raster размером 10000х10000
формула расчета:

Python
1
rater[i,j] = raster[i,j]*var_1+var2

Сам делаю обычным способом через 2 цикла:

Python
1
2
3
4
5
6
while i < 10000:
    j = 0
    while j < 10000:
        raster[j, i] = raster[i,j]*var_1+var2
        j = j +1
    i = i + 1

Но выполняется долго.
2) Как сделать перерасчет быстрее? Проц. Intel Core i5-6600K 3.5
3) Может создать 4 потока по 2500 перасчитовать? Где про это можно почитать?

4) NumPy быстрее обрабатывает массив чем способ указанный выше?
5) Можно ли как-то его тоже в 4 потока кинуть?

Jabbson · 23.11.2016, 16:28

да, можно

Python
1
raster*var_1-var2

у меня выполняется около 3 секунд при размере 10к на 10к

@DarkPower · 23.11.2016, 19:36 **[ТС]**

Сообщение от Jabbson

да, можно
PythonВыделить код
1
raster*var_1-var2
у меня выполняется около 3 секунд при размере 10к на 10к

А можете полностью написать код

Python
1
2
3
raster = numpy.array(gdalBand.ReadAsArray()) 
start_time = time.clock()
raster = K2 / (Ml*raster+Al + 1) #тут ошибку выдает TypeError: only length-1 arrays can be converted to python scalars

Jabbson · 23.11.2016, 19:47

если покажите что в переменных во всех, а то я понятия не имею, что за gdal (и кого он собственно "ждал"

)

@DarkPower · 23.11.2016, 20:03 **[ТС]**

Сообщение от Jabbson

если покажите что в переменных во всех, а то я понятия не имею, что за gdal (и кого он собственно "ждал" )

Окей.

Python
1
2
3
4
5
6
7
8
9
#gdalBand.ReadAsArray() - массив чисел типа float, [10k,10k]
raster = numpy.array(gdalBand.ReadAsArray()) #Собственно (как я понял) array to numpy.array
start_time = time.clock() #тут замеряем время
#raster = K2 / (math.log(K1/(Ml*raster+Al) + 1))
'''
K2 - тип float 1321.07 или 1201.14
K1 - тип float 774.88 или 480.88
Ml и Al тоже тип float
'''

Этот код выполняется 140сек:

Python
1
2
3
4
5
6
7
8
9
raster = numpy.array(gdalBand.ReadAsArray())
start_time = time.clock()
#raster = K2 / (math.log(K1/(Ml*raster+Al) + 1))
for col in range(xsize):
    for row in range(ysize):
        L = Ml*raster[row,col]+Al
        raster[row,col] = K2 / (math.log(K1/L + 1))
end_time = time.clock() 
print "\ttime(sec) = " + str(end_time-start_time)

Jabbson · 23.11.2016, 21:16

Сообщение от DarkPower

math.log

--> numpy.log

Добавлено через 3 минуты
у меня с массивом в 10к х 10к

Python
1
raster = K2 / (np.log(K1/(Ml*raster+Al) + 1))

выполняется от 4 до 6 секунд

@DarkPower · 23.11.2016, 21:23 **[ТС]**

Сообщение от Jabbson

math.log
--> numpy.log

Спасибо

с 130сек на 1,5сек упало =))

Jabbson · 23.11.2016, 21:56

всегда пожалуйста

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Обработчик клика мыши в браузере ПК и касания экрана в браузере на мобильном устройстве 8Observer8 02.02.2026 Содержание блога Для начала пошагово создадим рабочий пример для подготовки к экспериментам в браузере ПК и в браузере мобильного устройства. Потом напишем обработчик клика мыши и обработчик. . .	Философия технологии iceja 01.02.2026 На мой взгляд у человека в технических проектах остается роль генерального директора. Все остальное нейронки делают уже лучше человека. Они не могут нести предпринимательские риски, не могут. . .	SDL3 для Web (WebAssembly): Вывод текста со шрифтом TTF с помощью SDL3_ttf 8Observer8 01.02.2026 Содержание блога В этой пошаговой инструкции создадим с нуля веб-приложение, которое выводит текст в окне браузера. Запустим на Android на локальном сервере. Загрузим Release на бесплатный. . .	SDL3 для Web (WebAssembly): Сборка C/C++ проекта из консоли 8Observer8 30.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .
SDL3 для Web (WebAssembly): Установка Emscripten SDK (emsdk) и CMake для сборки C и C++ приложений в Wasm 8Observer8 30.01.2026 Содержание блога Для того чтобы скачать Emscripten SDK (emsdk) необходимо сначало скачать и уставить Git: Install for Windows. Следуйте стандартной процедуре установки Git через установщик. . . .	SDL3 для Android: Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 29.01.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами. Версия v3 была полностью переписана на Си, в. . .	Инструменты COM: Сохранение данный из VARIANT в файл и загрузка из файла в VARIANT bedvit 28.01.2026 Сохранение базовых типов COM и массивов (одномерных или двухмерных) любой вложенности (деревья) в файл, с возможностью выбора алгоритмов сжатия и шифрования. Часть библиотеки BedvitCOM Использованы. . .	SDL3 для Android: Загрузка PNG с альфа-каналом с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 28.01.2026 Содержание блога SDL3 имеет собственные средства для загрузки и отображения PNG-файлов с альфа-каналом и базовой работы с ними. В этой инструкции используется функция SDL_LoadPNG(), которая. . .

Jabbson 5907 / 3359 / 1036 Регистрация: 03.11.2009 Сообщений: 10,008
	23.11.2016, 19:47
	если покажите что в переменных во всех, а то я понятия не имею, что за gdal (и кого он собственно "ждал" ) 0

Jabbson 5907 / 3359 / 1036 Регистрация: 03.11.2009 Сообщений: 10,008
	23.11.2016, 21:56
	всегда пожалуйста 0

NumPy перерасчет масива

Решение