Questions tagged [python]
Python is an open source interpreted programming language used in many GIS programs.
12,256 questions
1
vote
1
answer
216
views
TROPOMI NO2 and tropospheric NO2 fields bugs
I downloaded the monthly NO2 data for the city of Bogotá using the GEE platform. I am analyzing the monthly average for April 2020 of both total NO2 and tropospheric NO2.
A large amount of missing ...
CommunityBot
- 1
0
votes
1
answer
282
views
Using cloud confidence to create cloud mask from Landsat 8 QA_PIXEL
I want to use the Landsat QA_PIXEl from this dataset to create a cloud mask in GEE. I only want to mask out the areas where there is high confidence that cloud cover is present. Looking at the ...
CommunityBot
- 1
0
votes
1
answer
735
views
Renaming attachments from a geodatabase
I am completing a task that requires me to collect info and pictures of gravestones in cemeteries using the Esri Collector app. There are thousands of attachments within the GDB that I would like to ...
CommunityBot
- 1
2
votes
1
answer
192
views
Grab/Index the raster bands before and after a unique band for each pixel (depth-wise) in a raster stack
I'm working with a time series stack of 28 raster bands. Each band is an image for a specific date, spanning all of the year 2020, and they are stacked in order. The stack (change_point_raster) is in ...
CommunityBot
- 1
0
votes
2
answers
2k
views
Problem with gdal.WarpOptions concerning reference
I have a tiff file that I want to crop down to a shp cutline and reproject.
If the code is:
clip_spec=gdal.WarpOptions(format='GTiff',cutlineDSName=os.getcwd()+'/moksha/DEM/cutline2.shp',cropToCutline=...
CommunityBot
- 1
-2
votes
0
answers
43
views
Which library to use for calculating flow accumulation from DEM raster?
I have acquired the DEM raster for a watershed area. I want to first calculate the flow accumulation using D8(Deterministic 8-Node) algorithm. Earlier in our lab people have used GRASS in Python to ...
- 65.5k
1
vote
2
answers
3k
views
Converting Raster to Shapefile using filtered values of raster
I have multiple raster files (RGB, 255) and I want to convert those raster files to shapefiles but I want only that area as a polygon where there is data, which means I want to exclude the zero values
...
CommunityBot
- 1
3
votes
1
answer
312
views
What is the syntax for assigning a feature attribute value to a variable in the QGIS Advanced Python Field Calculator?
Using the QGIS Advanced Python Field Calculator, I am able to assign the new field a value from a feature an existing field using value = <Attribute_1>. I run into a problem when I am using ...
CommunityBot
- 1
3
votes
1
answer
909
views
Image collection monthly averages using geemap package
I am trying to convert some code snippets from GEE JavaScript to Python. This code computes the monthly averages of NDVI images given a range of years and a range of months.
This is my code in ...
CommunityBot
- 1
1
vote
1
answer
4k
views
Reprojecting unequal number of latitude and longitude: pyproj.exceptions.ProjError: x, y, z, and time must be same size
I have a rectangular matrix of geographic climate data with a 1d row of latitude coordinates (size=834) and a 1d column with longitude coordinates (size=1115). They are both in unit of m (EPSG:32630) ...
CommunityBot
- 1
1
vote
1
answer
28
views
PDAL: open stream for ' ' with error 'No error'
I'm trying to run a pdal pipeline - which has worked before (on a different computer, different data around 3 years ago) - on a singular las dataset - however it keeps throwing up an error: PDAL: ...
- 75.7k
0
votes
1
answer
750
views
Performance difference between xarray open_mfdataset and open_dataset
I'm dealing with CMIP6 Scenariomip data, I have eight ESMs precipitation grid that each one is composed by multiple netcdf files. There is one ESMs's netcdf files are seperately reserved every decade, ...
CommunityBot
- 1
1
vote
2
answers
2k
views
Python imshow will not overlay -- only display last line's map
I have the following code that when run separately, displays two maps:
Map1:
f = h5py.File(filename[0], 'r')
group_id='Soil_Moisture_Retrieval_Data_AM'
var_id = 'soil_moisture'
a = f[group_id][var_id][...
- 83.1k
0
votes
1
answer
700
views
getInfo use in earthengine-api (Python)
When I write the an algorithm in Google Earth Engine using JavaScript by using the Code Editor, it works fast. After thatif I use the Map.addLayer()， the result of algorithm can show in the map ...
CommunityBot
- 1
2
votes
2
answers
522
views
Filtering points appearing in list of convex hull polygons
I have two DataFrames df & df1:
df: containing latitude, longitude, and binary cluster (0 or 1)
latitude Longitude cluster
30.400091 -100.722830 0
30.500091 -100.762830 0
30.600091 ...
CommunityBot
- 1