Questions tagged [code-review]

Questions requesting a review of working code should be asked at the Code Review Stack Exchange

Filter by
Sorted by
Tagged with
0 votes
0 answers
33 views

Should exporting image take days?

New to Earth Engine on this end, I am using trying to export specific bands of an image collection, filtered by date, and summarized as shown below. It takes at least two days for the image to ...
user avatar
1 vote
0 answers
40 views

A faster way to extract geometry xy of geodataframe

I need a faster way to extract geometry xy of a large dataset which is in a geodataframe format. I have developed a lambda function as below to do that but I need a faster way. The following is a ...
user avatar
2 votes
0 answers
46 views

Get file path of Geopackage with PyQGIS [duplicate]

I want to get the Geopackage file path where my layer is stored in. Layer = QgsProject.instance().mapLayersByName('DAT')[0] print(Layer.source()) But this only gives me the path with layer pointer: /...
user avatar
  • 457
1 vote
0 answers
99 views

Optimizing code to run faster in Google Earth Engine

I am trying to query properties (~10) from Earth Engine datasets (Temperature, Precipitation, GLDAS outputs, etc.) for a given coordinate point and time range in my dataset (90,000 points). I’ve done ...
user avatar
1 vote
0 answers
145 views

Finding distance from point to closest LineString given multiple LineStrings efficiently using Python

I have a dataset with many Points and I also have a GeoJSON file consisting of road networks using LineStrings. I want to estimate traffic by seeing which roads have the largest points on or beside ...
user avatar
2 votes
1 answer
65 views

Reducing number of for loops in ArcPy code

I want to clean this code in such a way that there are not so many for loops. For loops in the following code are placed to trap FCs in multiple group layers and do some work on them. My Data in TOC ...
user avatar
1 vote
0 answers
21 views

TPI calculation by PDAL [closed]

Is there a convenient way to calculate the topographic position index by PDAL? I have LAS file and what I do is iteratively for every point calculate the TPI. The coordinate information I take from ...
user avatar
0 votes
0 answers
218 views

Combining image collection from different years using Google Earth Engine

I'm trying to estimate superficial runoff for urban areas, but aparently having problems to generate for 2020. The NDVI collection goes to 2020, but the urban collection only goes until 2019. In line ...
user avatar
0 votes
1 answer
100 views

Export image collecton with name using Google Earth Engine

I am using GEE for a new collection of images and I need to export it by automating the name of each image, for example, that month one is called January, and so on, how can I do it? My code is here: /...
user avatar
  • 5
0 votes
0 answers
46 views

Loop over dataframes and extract specific values to generate multiple rasters/plots

I would seek help to revise this script. I need to generate multiple plots. For example, the following code works well ggplot and gsave but I would prefer to use raster, and I tried something such as ...
user avatar
  • 39
0 votes
0 answers
85 views

Split GEE Image into spatially discrete patches & intersect with another layer

I have an Image in GEE that corresponds to a species' spatial distribution (its 'Area of Habitat' = 'AOH'). Ultimately, what I need is to identify how many spatially discrete patches of habitat there ...
user avatar
  • 107
0 votes
0 answers
91 views

Optimizing minimum distance calculations between multiple polygons in R

I have a list of polygons and would like to create a minimum distance matrix by calculating the minimum distance between each possible pair. I am currently using the code below, but this is running ...
user avatar
0 votes
0 answers
67 views

Speeding up my Python code for reading from Excel, performing geospatial analysis and writing to Excel files

I have written a Python snippet that reads lat and long stored in an Excel file. Converts them to a point which is then used to perform multiple geospatial analysis including buffer, intersection, and ...
user avatar
1 vote
0 answers
272 views

Downloading and merging 0.5 minute WorldClim data for (long) list of coordinates using R

I want to download temperature and precipitation data from WorldClim for a huge list of coordinates. I am working in R, and have a solution using getData() from raster(), but it's going to take ...
user avatar
  • 111
-3 votes
1 answer
85 views

Reducing number of lines / more efficient coding method for field maps [closed]

I am looking for ways to reduce the number of lines when making field maps in ArcPy (note, I am using ArcPy with ArcMap 10.6). I cannot seem to figure out a way, and am stuck with a significant amount ...
user avatar
0 votes
1 answer
101 views

Calculating Rumple_Index for clipped Inventory Plot LAS files

I have successfully implemented the Rumple_Index example provided in the LidR book. I am now trying to calculate the same Rumple_Index for my forest field plots. I am trying to integrate the ...
user avatar
  • 11
1 vote
1 answer
102 views

Improving PyQGIS distance matrix efficiency

I have written a custom processing tool, which I use within a large graphical model in QGIS. I am looking for a way to speed it up. First step of this processing tool is to connect every point of one ...
user avatar
-1 votes
1 answer
1k views

Converting TIFF file and changing CRS

Is there a fastest way to: Load a .tif file; Convert it into GeoDataFrame; and Change Coordinate Reference System? I am using Python, and what I tried to do is: Use rasterio to check if the Source ...
user avatar
-1 votes
1 answer
251 views

Problem Visualizing and Downloading LANDSAT 5 Images With Google Earth Engine

I am new to GEE. I'm working with Landsat images. First, visualizing them in the map, then downloading them and working with them afterwards in QGis 3.6 to obtain certain indexes (NDVI, NDWI, Wetness) ...
user avatar
1 vote
1 answer
283 views

Easier Way to Record Yearly Mean NDVI

For my Google Earth Engine project, I am finding the mean NDVI for each year between 2000 and 2020 in protected areas a part of the World Database on Protected Areas (WDPA). As an example, I have ...
user avatar
0 votes
1 answer
60 views

How do I add the soil moisture data I am trying to get onto my Google drive?

Here is my code How do I add the soil moisture data I am trying to get onto my Google drive? var dataset = ee.ImageCollection('NASA_USDA/HSL/SMAP_soil_moisture') .filter(ee.Filter....
user avatar
1 vote
2 answers
139 views

Making ArcPy code run faster

I have an ArcPy code that I am using to find various characteristics about certain watersheds in Indiana. I start with a shapefile of all 800 watersheds and go through it row by row and perform a ...
user avatar
0 votes
0 answers
89 views

Error when displaying progress bar

My initial code read an image, apply a median filter to it, and save the smoothed image. It looked like this : #Image reading rasterpath = r"C:/Users/SIEXT5646/Desktop/Donnes/2018_gris.tif" ds = gdal....
user avatar
1 vote
1 answer
148 views

Arcpy Delete Multiple Layers

I'm writing an add-in for Arcmap and at the end of it I am left with multiple layers from buffers, joins, and selections. I am trying to delete these layers before the add-in finishes so to keep the ...
user avatar
4 votes
2 answers
2k views

Filling dictionary with list of row attributes using arcpy.SearchCursor

I have a dBase database file and I'm viewing and editing the table in ArcGIS Pro. I want to use the arcpy.da.SearchCursor function to fill a dictionary and create a key-value pair by iterating ...
user avatar
2 votes
1 answer
253 views

Creating NDVI from raster using ArcPy

I have to create NDVI's for multiple raster images, and so want to find out the most efficient way to do it. The code below works, but I wondered if there was a more efficient way to do it? # ...
user avatar
5 votes
1 answer
248 views

Optimizing Intersection Area Loops on big layers using PyQGIS

Using 3.10.1-A Coruña on Windows 10 I am looping over some (10) national shapefiles and finding intersection and areas but the process is very slow. I have two loops - one loop for each postcode and ...
user avatar
0 votes
0 answers
2k views

Tool script for GeoJSON to Shapefile

I am creating an ArcMap tool script that takes user input and then takes a GeoJSON file from a specified folder and then opens it in ArcMap as a shapefile. The issue I am running into is that I am ...
user avatar
3 votes
0 answers
379 views

Alternative to Rasterio Mask module for clipping a raster using masked geojson

I am clipping a raster (.tif file) using a mask vector (.geojson) file. with rasterio.open(url) as src: out_image, out_transform = rasterio.mask.mask(src, shapes, crop=True) But, this ...
user avatar
  • 719
1 vote
1 answer
59 views

Improving ArcPy script for ArcMap to only return rows of selected polygons instead of all polygons in layer

I am debugging a code provided to me that merges two polygons and alters the fields of "FID1' and "CALC_AREA." I have realized that the current code searches every layer in the map document for ...
user avatar
0 votes
1 answer
47 views

Creating dictionary of link and its adjacent links (having common vertex) in ArcPy

How do I create a dictionary that has the linkID of a link as key and linkID of touching links as values? I have used this code in ArcPy but it takes a very long time. arcpy.CopyFeatures_management(...
user avatar
  • 755
0 votes
1 answer
539 views

Google Earth Engine : error: Cannot use 'in' operator to search for 'type' in {"type":"LinearRing","coordinates":

I'm trying to export Sentinel-1 bands with this code in Google Earth Engine but I have recently this error: Cannot use 'in' operator to search for 'type' in {"type":"LinearRing","coordinates": I ...
user avatar
1 vote
0 answers
143 views

Optimizing Google Earth Engine cloud cover for Sentinel 5p

I am trying to generate cloud free Sentinel 5p images, I have a basic cloud detection code working but it runs very slowly (most likely due to nested map and reduce features). Is there a way to do ...
user avatar
2 votes
0 answers
166 views

Making Dijkstra faster

I would like to compute shortest path between 2,5 millions pairs. I created table with id, origin, the nearest node of origin, destination, the nearest node of destination. I also created index on ...
user avatar
  • 195
1 vote
0 answers
148 views

QGIS 3.8 Python Code for calculate area

Is this the the best solution to calculate the area of the features? parent = iface.mainWindow() mc = iface.mapCanvas() fldname = "m2" #field which has the values for squaremeters myLyr = mc....
user avatar
0 votes
0 answers
28 views

Getting Exported bands must have compatible data types; found inconsistent types: Float64 and Float32 using Google Earth Engine [duplicate]

I'm trying to export sentinel-1 bands with this code in Google Earth Engine but I have recently this error. I didn't have this problem in the past. So here is the code. The code: var start = new ...
user avatar
  • 9
2 votes
1 answer
235 views

Optimising Large Spatial Queries in SQL

I'm sure similar answers have been given elsewhere but no matter what I have tried I am unable to optimise this query. Details: TableA has 3.5mill records TableB Had 57000 records but I have ...
user avatar
0 votes
1 answer
63 views

Getting same result with less code with PythonWin ArcPy? [closed]

I have finished writing code and it does what I need without having any error. However, I am wondering if there is any option to get the same result in fewer lines, pecially the for loop inside the ...
user avatar
  • 29
0 votes
1 answer
76 views

Automate Selection/Copy Features from Selection

I am new to Python and ArcPy, but I would like to automate this process and/or make my code more efficient here. I am basically selecting a layer by attribute, then copying features from the (selected)...
user avatar
  • 11
1 vote
1 answer
48 views

Speeding up script to flag nodes within 100 feet of segments using PostGIS?

I'm working with two datasets, attempting to flag when an intersection is within 100 feet of a network subset. There are over 40,000 intersections and my network is ~ 1400 segments. The problem with ...
user avatar
  • 325
2 votes
1 answer
464 views

Populating layer field with dict.values() in PyQGIS 3?

I've coded the script below in PyQGIS 3 to get the azimuth of starting or ending points for a line layer that touches a border. It works fine, but I want to populate the azimuth field layer with the ...
user avatar
1 vote
0 answers
523 views

Calculating nearest distance between points and Spatial Polygons in R? [closed]

So right now I have pretty large datasets, about 200,000 (lat,long) locations and a Large SpatialPolygonsDataFrame with a polygons list of length 428386 spread in a single city of around 450 square ...
user avatar
0 votes
1 answer
299 views

Making looping through line segments for vertex points faster?

I am trying to speed up some current code that I have for looping through lines, getting their segments and ultimately the vertex points in Python. Can anyone recommend a faster method than what I ...
user avatar
11 votes
1 answer
6k views

Counting points in polygons with sf package of R?

I'm quite new to R and especially GIS with R. So recently I had to proceed to a quite simple analysis: counting health centers (points) per administrative unit (polygons). Quite common process in ...
user avatar
  • 353
0 votes
0 answers
120 views

Why is internal memory not overwritten?

I have the following code in Python for ArcMap to first select a property polygon and then with this selection to intersect the layer which contains lines (sewers) (which pass through or are within ...
user avatar
  • 1
1 vote
1 answer
1k views

I can not export the image to the drive in Google Earth Engine (GEE)

Whenever I try to export to the drive the image gives me the following error: Error: Invalid argument: 'image' must be of type Image. This is my code: var mod11 = ee.ImageCollection("MODIS/MOD11A2"...
user avatar
1 vote
1 answer
438 views

Using accumulate function for temperature units?

I am working on a piece of code that is designed to aggregate growing degree days for a particular region. To start, the equation for growing degree days is: Daily GDD (°F) = Daily Average ...
user avatar
4 votes
2 answers
314 views

ArcPy UpdateCursor deleting duplicate records?

I am working with .dbf files. I have two columns of data that have duplicates. A column of parcel numbers with duplicates and a column of letters that are duplicated. I want to delete all parcel ...
user avatar
  • 41
1 vote
1 answer
524 views

Reducing hourly rainfall data to monthly data for particular period of time and making a timeseries chart

I am working with GSMaP image collection (hourly rainfall data) in Google Earth Engine. I would like to make a script which is calculating and making a chart of sum rainfall in particular location ...
user avatar
2 votes
4 answers
91 views

More Efficient Conditional Statements - ArcGIS Field Calculator

I am new to Python and have successfully created the following statement that writes a change period value (01 or 02) based on the municipality and date fields. I am wondering if anyone could point ...
user avatar