Questions tagged [code-review]
Questions requesting a review of working code should be asked at the Code Review Stack Exchange
58
questions
0
votes
1answer
65 views
Reducing number of lines / more efficient coding method for field maps
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 ...
-1
votes
0answers
22 views
Speed up PyQGIS proces for attribute modification (100k features) with expressions
Context: I am currently in the process of using PyQGIS to automate the data-input required for external software. This software requires certain attributes which can be created by using QGIS to modify ...
0
votes
0answers
34 views
Visualizing and downloading Landsat 5 NDVI Composite images
I'm working with GEE (no so experienced), trying to download a series of Landsat 5 NDVI composite images. The thing is the code I have isn't showing me the images that I want to downlad afterwards.
...
0
votes
1answer
46 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 ...
1
vote
1answer
78 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 ...
-1
votes
1answer
225 views
Converting TIFF file and changing CRS
Is there a fastest way to:
Load a .tif file
Convert it into GeoDataFrame
Change Coordinate Reference System
I am using Python, and what I tried to do is:
Use rasterio to check if the Source CRS is ...
-1
votes
1answer
85 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) ...
0
votes
1answer
72 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 ...
0
votes
1answer
36 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....
1
vote
2answers
89 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 ...
0
votes
0answers
60 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....
1
vote
1answer
51 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 ...
4
votes
2answers
876 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 ...
2
votes
1answer
117 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?
# ...
5
votes
1answer
130 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 ...
0
votes
0answers
1k 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 ...
2
votes
0answers
225 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 ...
1
vote
1answer
53 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 ...
0
votes
1answer
43 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(...
0
votes
1answer
241 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 ...
1
vote
0answers
109 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 ...
1
vote
0answers
76 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....
0
votes
0answers
24 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 ...
2
votes
1answer
133 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 ...
0
votes
1answer
61 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 ...
0
votes
1answer
42 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)...
1
vote
1answer
47 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 ...
2
votes
1answer
308 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 ...
1
vote
0answers
303 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 ...
-1
votes
1answer
150 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 ...
5
votes
1answer
3k views
Counting points in polygons with sf package of R? [closed]
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 ...
0
votes
0answers
114 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 ...
1
vote
1answer
972 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"...
1
vote
1answer
227 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 ...
3
votes
2answers
157 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 ...
1
vote
1answer
339 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 ...
2
votes
4answers
88 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 ...
2
votes
0answers
451 views
Moving Feature Class to Different Feature Dataset within Same File Geodatabase using ArcPy?
I want to progammatically move a feature class into a different feature dataset within the same file geodatabase. I'm wondering if there is a way of doing this without including code to rename and ...
-1
votes
1answer
63 views
More Pythonic way to loop through UpdateCursor?
I have several sets of fields to populate based on other fields. The code below works, but I am wondering if there is a more Pythonic way to march through the i,j pairs? (Actual fieldnames are not ...
0
votes
2answers
50 views
Changing values in column based on another column in arcpy
I created new field in a table named 'Cimp'. Now I want to put there values that depends on values from column "COS_Cod" for example, if COS_Cod starts with 1.1 then Cimp is equal to 1. COS_Cod is a ...
0
votes
1answer
674 views
Fastest way to know in which section of line is a point
I have a line and a point. I want to locate where is this point in the line according to imaginary sections I make (any way I want). The output should be in which number of section is the point (...
1
vote
1answer
3k views
Fastest way to extract a raster in R (improve the time of my reproducible code)
I'm wondering if I have maximized the speed at which a mean of an area buffered around a point in a raster can be extracted.
Fastest way I have ever found to extract a raster, with the pre-cropping ...
1
vote
0answers
69 views
Unable to get points in buffer on selected feature in QGIS using Python?
I have to create a 1.5 km buffer on each latitude/longitude of layer A. I would like to know how many points of line B will be inside the buffer of layer A. I used the following code but it takes too ...
0
votes
0answers
726 views
Speeding-up processing of pair-wise distance calculation in Python?
I am trying to do a cluster analysis using DBSCAN for my time series NDVI image in Python. I am using distance time warping (DTW) to measure distances between my time series. With this, I am computing ...
2
votes
1answer
64 views
Counting only data points which are within a specified radius of other data points?
I wanted to know if I am on the right track here with what I am trying to do. I have two shapefiles, the first contains point data for all the businesses in DC and the other contains point data all ...
0
votes
1answer
229 views
Overview of available images per image collection?
I´m trying to make an overview in a table or graph of the available images in Landsat, Sentinel and Modis series for a certain area of interest. Although I could solve it as follows, this is obviously ...
1
vote
0answers
408 views
Finding a more efficient way to execute a long-run python script [closed]
I'm executing a python script which does these steps:
iterates through 19 features (makes pairs "each one with each other", order doesn't matter: A-B = B-A)
Creates Closest Facility Layer
Adds ...
1
vote
1answer
400 views
Improving poor performance of ArcPy update cursor with ArcSDE and Editor? [closed]
I am using ArcGIS 10.5 in an ArcSDE environment.
I have two tables, they are the exact same table schemas, but one is being edited and those edits are being pushed to the main table. I am a little ...
1
vote
0answers
43 views
Improve my arcpy script for simple querying [closed]
I am quite new to coding python in arcmap so I would like to get my code reviewed and hopefully improved.
I have two Shapefiles, one contain polygons of cities borders and the other is points ...
-1
votes
1answer
310 views
Exporting attributes table from csv to spatialite file using PyQGIS?
I have this python script which allows me to convert delimited text file to spatialite . I used a loop to export attributes table to the output file and it takes over 2 hours when I execute the code ...
