Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [numpy]

NumPy is a Python module, adding support for large, multi-dimensional arrays and matrices, along with a large library of high-level mathematical functions to operate on these arrays.

1
vote
0answers
11 views

GDAL rastercalc giving unexpected results

I have a raster file containing the following values: data values from 0 to 8488 nodata values values equalling -3.40282347e+38 The values equalling -3.40282347e+38 are presumably supposed to be ...
0
votes
1answer
27 views

Open 4d NetCDF subdatasets with GDAL in Python

I have NetCDF files containing subdatasets with 4 dimensions: time, height, latitude, longitude. Here is an example of the output of GDAL GetSubDatasets(): ('NETCDF:"Martinique_ALT22.nc":height', '[...
-1
votes
1answer
38 views

Writing multiband raster from 3d numpy array gives an error [closed]

Can you help to fix this code? I tried to filter multiband (240 bands) raster as numpy array: import rasterio import numpy as np from scipy.signal import savgol_filter with rasterio.open('D:/...
0
votes
1answer
57 views

“No module named numpy” Python Error in ArcGIS Label Expression?

I have developed code to label points that pull formatting from a dictionary. I developed this code and it was working fine two weeks ago, today I try to open my map and I am now met with an error ...
0
votes
0answers
34 views

Broadcast Shapely Point on numpy array

I have some code that looks something like this: def get_objs(all_objs, line, my_distance): """ Gets all objs within my_distance of line. all_objs is a numpy array of shape (x, 4) ...
0
votes
0answers
18 views

Calculating Raster Percentile Values in ModelBuilder

I want to calculate the 75th percentile value of a raster dataset (floating, 32Bit) and use this value as an input to a set null expression. I have multiple rasters to perform this on so I need to ...
2
votes
1answer
40 views

NumPyArrayToFeatureClass gives RuntimeError: create table

I've tried to change polygon feautre to point by using tools: FeatureClassToNumPyArray NumPyArrayToFeatureClass. My code as below: import arcpy import numpy as np in_file='C:/ios/...
1
vote
1answer
28 views

Extracting values from array created with FeatureClassToNumPyArray and converted to DataFrame?

I cannot understand how to use the values in the arrays created with arcpy.FeaturesToNumPyArray. I want to create an array and then convert them to a pandas DataFrame. So I created an array from ...
1
vote
1answer
52 views

Array in python just full of zeros when opened using GDAL

I'm working on the below script. What I'm trying to do is create a stats.txt file output which has the Land Cover Class ID, Fractional Coverage (of that class) and the number of pixels (of that class) ...
1
vote
1answer
34 views

Classifying quantiles in raster dataset with ArcPy?

I am trying to write a Python script which uses the reclassify function in ArcPy to create 9 classes in a raster file. I want to use the quantile method to create the breaks, but so far, I can't seem ...
0
votes
0answers
70 views

Binary Raster (numpy array) to shapefile

i did a classification with a min/max threshold for stormdamaged areas. For that i calculated two indices (NSDVI) and (SAVI). I intersected two numpy arrays from those results. I recieve a binary ...
2
votes
0answers
35 views

Fill masked array values with reflect padding

I have a square matrix (image) with masked values in an area. The mask corresponds to missing values in areas close to the border (always). I want to create a new image where the masked values are ...
0
votes
1answer
97 views

How to create a GeoDataFrame points grid from a numpy.mgrid?

I'm trying to create a GeoDataFrame points grid from a numpy.mgrid. I already have created the numpy.mgrid with coordinates of future points. >>>arr = np.mgrid[-1827342.1753285176:-1827523....
0
votes
2answers
38 views

If field 1 > field 2 update field 3 in loop

I am trying to update the table above. What I want to do is if the values in the PMP fields are less than the Max fields update the check fields with PASS, if not FAIL. I started doing this trying ...
1
vote
1answer
72 views

NumPy array to Raster file (GeoTIFF)

I converted a raster file (one band grey scale) to an numpy array and set a range to it. I Did like that: import gdal import numpy ds = gdal.Open('H:/GDAL_docs/a.tif') myarray = numpy.array(ds....
2
votes
1answer
37 views

Get data from a raster using two coordinate arrays with Python

I'm trying to find a function in Python, that given two arrays with geographic coordinates, returns a list or array of raster values in those positions. that is to say: import numpy as np """geog ...
0
votes
1answer
68 views

Loading all rasters in a directory and naming the variables iteratively

Hi I'm working with prism data and I already have clipped (masked) rasters (*.bil files) of a specific area for each day of a year. I'd like to open all these rasters with rasterio, convert them into ...
1
vote
1answer
90 views

Comparing pixel values at a given point in two rasters

I have two raster images, one for land use and one for NDVI, i would like to calculate LAI based on land use. So i need to compare the pixel value in the NDVI raster and Land use raster then select ...
1
vote
0answers
47 views

Numpy import error in OSGeo4W Shell

I and my colleagues are running into a very frustrating Python import error while using the OSGeo4W Shell. Any time we attempt to use gdal_calc in a Python script, we get the following error message: ...
1
vote
1answer
108 views

Using rasterio and matplotlib how can I create a window for an NDVI image to fit into?

I am using Rasterio, Numpy, and Matplotlib. I created a window with the dimensions of 800x600 and plotted the band 5 for a Landsat 8 image, that worked just fine. I also have the nir image so then I ...
0
votes
1answer
162 views

Using Rasterio, Matplotlib, and Numpy trying to calculate and plot NDVI of a raster image getting error: IndexError: list index out of range?

First I create two variables and store my raster images in them and I use MatPlotLib to plot a subset of the original "nirband" raster image. That works fine. I am getting an IndexError: list index ...
1
vote
0answers
35 views

Using Null Values in Raster To NumPy?

I am trying a convert a raster dataset [a simple 4 x 4] into a numpy array using RasterToNumPyArray. Here is my sample raster dataset: RasterToNumPyArray however, returns zero for Null values: ...
3
votes
0answers
68 views

QGIS 3 - Save processing output to memory and turn into numpy array?

I'm currently developing a python plugin for QGIS 3. In short, the code conducts among others a viewshed algorithm, turns the output in a numpy array and calculates some spatial metrics the user can ...
0
votes
1answer
91 views

How to convert masks generated by rasterio to channel last mask

The problem is about using rasterio to generate mask from GeoJSON file which contains geographic coordinates. The function that I used is from: from rasterio.tools.mask import mask with rasterio....
1
vote
0answers
304 views

Writing numpy array to rasterio raster file generates empty file

I'm using Rasterio to read GeoTIFF from Sentinel 2 and calculating NDVI, but when I try to save as a new TIF file, the file turns to transparent and empty file. What am I missing? That's the code: ...
0
votes
0answers
98 views

numpy array coordinates to polyline in shapefile Python

I am calculating the shortest path in a numpy array which was converted from a raster (.tif). I have a list representing a path like [(1,3),(2,3)(3,3)]. How can I transform this into a polyline in ...
1
vote
1answer
111 views

Write array to CSV inside arc script tool (works in pycharm, doesn't work in arc)

I have an arcpy script in which I am writing an array to a CSV (for later use in the script to make a routeing layer). I've tried a couple iterations, and all work like a charm in PyCharm (can't ...
2
votes
0answers
179 views

Networkx shortest paths for two point layers

I have a task to know the best way from houses to subway station. There are three types of data on a picture below: orange dots are the entrances to subway stations violet dots are centerspoints of ...
1
vote
0answers
74 views

Point displacement between map (Cesium 3D) and raster (numpy array)

It is better to show than tell: On this image, the layer below the colorful shape is a trivial heat map in raster format. The colorful shape is the product of the skimage.graph.MCP_Geometric method. ...
1
vote
1answer
197 views

Writing a raster file from multi-dimension array using Python

I've been attempting to create a raster file that has multiple bands (10 or more) from a numpy array with multiple dimension(i.e. 10 or more). The script I'm currently using exports the TIFF file I ...
3
votes
1answer
168 views

Determine size of GeoTIFF header

Is there a way to find the size (in bytes) of the header information in a GeoTIFF file using Python? I am not trying to read the header (I know I can use gdalinfo for this) but rather to figure out ...
2
votes
1answer
455 views

Efficiently read large tif raster to a numpy array with GDAL?

I have been using the Python GDAL API to read tif raster files as NumPy arrays. Previously, I simply read the raster into an array directly with GDAL: ds = gdal.Open('example.tif') fullarray = np....
1
vote
1answer
200 views

Numpy array to raster with ArcPy [closed]

I have a Python script to download data in netCFD format. I have managed to extract the information I need in a masked array. This looks like: print(my_array.shape()) >> (700, 299, 217) As a ...
0
votes
0answers
283 views

Opening numpy array to rasterio?

I have numpy array that I know bbox for and can transform like this: import rasterio.plot from rasterio.transform import from_bounds t = from_bounds(minx, miny, maxx, maxy, sorted_np[3,:,:,7].shape[...
0
votes
0answers
27 views

Google Earth Engine Output - Returns Unexpected Values with np.sum()

I am exporting a CHIRPS precipitation raster from Google Earth Engine. This is part of a larger collection of data that I'll be exporting, but I'm just inspecting one image at the moment. I've ...
1
vote
2answers
471 views

NumPy memory error on large rasters

I have multiple large rasters. Their 5 percentile and 95 percentile are needed to be calculated. After loading the rasters to the ArcMap, I am using the following codes - import numpy import arcpy ...
3
votes
1answer
54 views

Raster representation in arcpy/numpy

I have a .tif raster file which contains integer values at 1, 2, 3, 4 and 5. Importing this file into ArcMap also shows a number of cells with Nodata value where the above classes are not present. I ...
2
votes
2answers
397 views

How to create a TIFF file using GDAL from a numpy array and specifying NoData value

I have some tiff files I obtained from splitting a bigger tiff file using gdal translate. The border tiff files have some nodata values, which I can see in numpy as number 15. I want to manipulate ...
1
vote
0answers
1k views

GDAL 2.3.1 is installed but the llinux terminal is using GDAL 2.2.2

I asked this question before realising this site existed: https://stackoverflow.com/questions/51476826/gdal-2-3-1-is-installed-but-the-llinux-terminal-is-using-gdal-2-2-2 I've installed GDAL 2.3.1 ...
0
votes
1answer
112 views

NumPyArrayToTable() RuntimeError: The value type is incompatible with the field type. [id]

# Create a dataframe from sorted JSON dictionary containing desired data df = pd.DataFrame(data, columns=['lat', 'lon', 'id']) # coerce dtype 'O' object values to np array native float64 type df['id']...
4
votes
1answer
533 views

Use Google Earth Engine to read MODIS NDVI into python

I'm trying to adapt this example for my needs: https://mygeoblog.com/2017/10/06/from-gee-to-numpy-to-geotiff/ The difference is that I'm using the MODIS NDVI image collection. # define the image img ...
2
votes
1answer
738 views

How to split multiband image into image tiles using Rasterio?

I have been using this answer on GIS SE to split 3-band (RGB) imagery into 256x256 3-band image tiles: import os, gdal in_path = '/path/to/indata/' input_filename = 'dtm_5.tif' out_path = '/path/to/...
2
votes
1answer
186 views

ArcGIS region group open source equivalent in numpy/python?

I am trying to get a unique ID for contiguous raster objects with the same value, like what the Region Group tool from ArcGIS does. However I need an open source solution using python/numpy/gdal/...
2
votes
0answers
321 views

LiDAR point cloud to numpy array to raster?

I'm looking for a way to create a height raster (raster of z values - i.e. not just ground poinits). I can accomplish this with gdal by creating a point vector layer from the numpy array then using ...
2
votes
1answer
180 views

NumPyArrayToTable() (SystemError: error return without exception set)

My goal for this script is to extract some information from an Overpass API query and ultimately join it with the table of another feature class I have created. So far I have completed the following ...
2
votes
1answer
155 views

Arcpy “MemoryError”; Get Classification Breaks on a TIF Raster; Jenks Breaks; Python

I have a large .tif raster with integer values. I wish to generate Jenks classification breaks using code and then use those breaks to reclassify my .tif file. I came across a block of code that I ...
1
vote
0answers
147 views

Hotspot Analysis QGIS Error: An error has occurred while executing Python code: ValueError: Illegal field requested in GetField()

I am doing a hotspot analysis for chicken spatial distribution in QGIS using Ubuntu. I have already installed the repositories used for hotspot analysis. However, an error occurred after I chose a ...
0
votes
0answers
35 views

Create polygon mask out of SAR Geom range azimuth coordinates

I have a SAR image, say s_img1, and a subset of min/max list of vertices as [(min_azimuth, min_range), (min_azimuth, max_range), (max_azimuth, min_range), (max_azimuth, max_range)] i.e they are a part ...
6
votes
1answer
1k views

numpy array to GTiff using rasterio without source raster

I have a numpy array that represents rasterized data from a LiDAR point cloud. For example: import numpy as np import rasterio arr = np.random.randint(5, size = (100, 100)) I know information about ...
0
votes
0answers
163 views

numpy array to gdal.Band

I have a numpy array that I want to convert to polygons using the gdal.Polygonize function, as outlined using this Polygonize a Raster Band approach. Unlike the example listed, however, my data is ...