Я использую ArcGIS 10.1 и хочу создать новый растр на основе двух ранее существующих растров. RasterToNumPyArray имеет хороший пример , который я хочу , чтобы адаптироваться.
import arcpy
import numpy
myArray = arcpy.RasterToNumPyArray('C:/data/inRaster')
myArraySum = myArray.sum(1)
myArraySum.shape = (myArray.shape[0],1)
myArrayPerc = (myArray * 1.0)/ myArraySum
newRaster = arcpy.NumPyArrayToRaster(myArrayPerc)
newRaster.save("C:/output/fgdb.gdb/PercentRaster")
Проблема в том, что он лишает пространственной привязки, а также размера ячейки. Я подумал, что это должно делать arcpy.env, но как мне установить их на основе входного растра? Я не могу разобраться.
Принимая ответ Люка, это мое предварительное решение.
Оба решения Люка правильно установили пространственную привязку, экстент и размер ячейки. Но первый метод не правильно переносил данные в массив, и выходной растр всюду заполняется ноданными. Его второй метод работает в основном, но там, где у меня большая область nodata, он заполняется блочными нулями и 255. Это может иметь отношение к тому, как я обрабатывал ячейки нодаты, и я не совсем уверен, как я это делал (хотя должен быть еще один вопрос). Я включил изображения того, о чем я говорю.
#Setting the raster properties directly
import arcpy
import numpy
inRaster0='C:/workspace/test0.tif'
inRaster1='C:/workspace/test1.tif'
outRaster='C:/workspace/test2.tif'
dsc=arcpy.Describe(inRaster0)
sr=dsc.SpatialReference
ext=dsc.Extent
ll=arcpy.Point(ext.XMin,ext.YMin)
# sorry that i modify calculation from my original Q.
# This is what I really wanted to do, taking two uint8 rasters, calculate
# the ratio, express the results as percentage and then save it as uint8 raster.
tmp = [ np.ma.masked_greater(arcpy.RasterToNumPyArray(_), 100) for _ in inRaster0, inRaster1]
tmp = [ np.ma.masked_array(_, dtype=np.float32) for _ in tmp]
tmp = ((tmp[1] ) / tmp[0] ) * 100
tmp = np.ma.array(tmp, dtype=np.uint8)
# i actually am not sure how to properly carry the nodata back to raster...
# but that's another Q
tmp = np.ma.filled(tmp, 255)
# without this, nodata cell may be filled with zero or 255?
arcpy.env.outCoordinateSystem = sr
newRaster = arcpy.NumPyArrayToRaster(myArrayPerc,ll,dsc.meanCellWidth,dsc.meanCellHeight)
newRaster.save(outRaster)
Изображение показывает результаты. В обоих случаях ноданные клетки показаны желтым цветом.
Второй метод Люка
Мой предварительный метод
источник
У меня были некоторые проблемы с получением ArcGIS для правильной обработки значений NoData с примерами, показанными здесь. Я расширил пример из блога reomtesensing.io (который более или менее похож на решения, показанные здесь), чтобы лучше обрабатывать NoData.
Очевидно, ArcGIS (10.1) нравится значение -3.40282347e + 38 как NoData. Так что я конвертирую туда и обратно между NumPy NaN и -3.40282347e + 38. Код здесь:
источник