Im new to python and I'm trying to write a function that will take a numpy array from a netcdf file with dimensions [time,height,longitude,latitude] and interpolate the function to a specified lat and lon. I have looked into scipy.interpolate but am still not sure where to go from there. any help?
grid point interpolation of wrf data to lat lon in python with numpy and scipy
4.8k Views Asked by Barry Baker At
2
There are 2 best solutions below
Related Questions in PYTHON
- How to store a date/time in sqlite (or something similar to a date)
- Instagrapi recently showing HTTPError and UnknownError
- How to Retrieve Data from an MySQL Database and Display it in a GUI?
- How to create a regular expression to partition a string that terminates in either ": 45" or ",", without the ": "
- Python Geopandas unable to convert latitude longitude to points
- Influence of Unused FFN on Model Accuracy in PyTorch
- Seeking Python Libraries for Removing Extraneous Characters and Spaces in Text
- Writes to child subprocess.Popen.stdin don't work from within process group?
- Conda has two different python binarys (python and python3) with the same version for a single environment. Why?
- Problem with add new attribute in table with BOTO3 on python
- Can't install packages in python conda environment
- Setting diagonal of a matrix to zero
- List of numbers converted to list of strings to iterate over it. But receiving TypeError messages
- Basic Python Question: Shortening If Statements
- Python and regex, can't understand why some words are left out of the match
Related Questions in NUMPY
- Why numpy.vectorize calls vectorized function more times than elements in the vector?
- Producing filtered random samples which can be replicated using the same seed
- Numpy array methods are faster than numpy functions?
- When I create a series of spectrograms from a long audio file, the colour intesities vary noticably
- How do I fix a NumPy ValueError for an inhomogeneous array shape?
- How should I troubleshoot "RuntimeWarning: invalid value encountered in arccos" in NumPy?
- Unravel by multi-index/group
- Calculating IRR Using Numpy
- Integrating with an array of upper limits without sacrificing time efficiency
- Why doesn't this code work? - Backpropagation algorithm
- How to remove integers from a mixed numpy array containing sub-arrays and integers?
- How to transfer object dataframe in sklearn.ensemble methods
- Rust cannot borrow as mutable
- Why does the following code detect this matrix as a non-singular matrix?
- How to detect the exact boundary of a Sudoku using OpenCV when there are multiple external boundaries?
Related Questions in SCIPY
- Integrating with an array of upper limits without sacrificing time efficiency
- Hourly and annual optimization problem over matrix
- Checking Event in solve_ivp
- Is polar decomposition commutative for diagonal matrices?
- loading variables from a python script in matlab
- cov2corr() for scipy sparse matrices
- Solving double integration in python using scipy.integrate
- Saving a scipy.sparse matrix directly as a regular txt file
- How can I apply scipy.interpolate.RBFInterpolator on an image / ndarray?
- Truncated normal distribution doesn't agree with untruncated normal distribution?
- Incorrect world-to-camera matrix with euler angles (yaw, pitch, roll)
- Scipy solve_ivp extremely slow/freezing
- Interpolating 3D volumetric data with scipy's RegularGridInterpolator
- Hilbert envelope in scipy matches original signal
- How to fill the area for segmentation using fill_holes function?
Related Questions in SPATIAL-INTERPOLATION
- How to interpolate missing values given a SpatRaster object in R Terra package?
- Interpolate points on curve having start/end points coordinates and curvature values (Python)
- Add padding to CubeSpline interpolation curve
- IDW Interpolation of point rainfall data using Python
- Interpolate CSV data and mask with SHP file
- Rasterize overlapping sf objects in R and specify which values to use for pixels
- First call to interpolator constructed by LinearNDInterpolator (with precomputed triangulation) slow regardless of grid size or values
- Is there a way in R i can create a new polygon (representing a certain constant value) based on two polygons representing 2 different constant values
- Vectorized method to insert arrays of different length at different spots in an ndarray
- Inverse Distance Weighting with gstat
- Area Weighted Average Grid Interpolation
- Spatial Interpolation Map Temperature using ggplot in R
- 2D slice from 3D volume (interpolation)
- Interpolation with tempdisagg package in R
- st_interpolate_aw returns error - replacement has x rows, data has y
Trending Questions
- UIImageView Frame Doesn't Reflect Constraints
- Is it possible to use adb commands to click on a view by finding its ID?
- How to create a new web character symbol recognizable by html/javascript?
- Why isn't my CSS3 animation smooth in Google Chrome (but very smooth on other browsers)?
- Heap Gives Page Fault
- Connect ffmpeg to Visual Studio 2008
- Both Object- and ValueAnimator jumps when Duration is set above API LvL 24
- How to avoid default initialization of objects in std::vector?
- second argument of the command line arguments in a format other than char** argv or char* argv[]
- How to improve efficiency of algorithm which generates next lexicographic permutation?
- Navigating to the another actvity app getting crash in android
- How to read the particular message format in android and store in sqlite database?
- Resetting inventory status after order is cancelled
- Efficiently compute powers of X in SSE/AVX
- Insert into an external database using ajax and php : POST 500 (Internal Server Error)
Popular # Hahtags
Popular Questions
- How do I undo the most recent local commits in Git?
- How can I remove a specific item from an array in JavaScript?
- How do I delete a Git branch locally and remotely?
- Find all files containing a specific text (string) on Linux?
- How do I revert a Git repository to a previous commit?
- How do I create an HTML button that acts like a link?
- How do I check out a remote Git branch?
- How do I force "git pull" to overwrite local files?
- How do I list all files of a directory?
- How to check whether a string contains a substring in JavaScript?
- How do I redirect to another webpage?
- How can I iterate over rows in a Pandas DataFrame?
- How do I convert a String to an int in Java?
- Does Python have a string 'contains' substring method?
- How do I check if a string contains a specific word?
I can't tell from the question if you are looking to (a) interpolate the entire data set to a new set of latitude/longitude coordinates or (b) get the value at a single latitude/longitude location.
Either way, if you're not married to doing this in python, it would be a lot easier to use remapgrid in Climate Data Operators (make sure you install with netCDF support). This does require the netCDF files to follow CF-Conventions (see this post on making your WRF file CF-compliant).
Documentation for remapgrid is found here. See in particular section 1.3 (pp. 9-12) for info on grid options and section 2.10 (starts p.106) for info on interpolation techniques. There are several options for how to remap, and the command you use depends on your application.
(a)Here's an example remapping using bilinear interpolation to a gaussian (128x64) grid
(b) Here's an example remapping using nearest neighbor mapping to one latitude,longitude point: