Ufdjrw

Project Euler #1 in C++

Sum letters are not two different

How to write the following sign?

I am having problem understanding the behavior of below code in JavaScript

Do I really need to have a message in a novel to appeal to readers?

What would you call this weird metallic apparatus that allows you to lift people?

What does it mean that physics no longer uses mechanical models to describe phenomena?

How much damage would a cupful of neutron star matter do to the Earth?

What order were files/directories outputted in dir?

Why does it sometimes sound good to play a grace note as a lead in to a note in a melody?

An adverb for when you're not exaggerating

How could we fake a moon landing now?

What was the first language to use conditional keywords?

Question about debouncing - delay of state change

Crossing US/Canada Border for less than 24 hours

Converted a Scalar function to a TVF function for parallel execution-Still running in Serial mode

Why does the remaining Rebel fleet at the end of Rogue One seem dramatically larger than the one in A New Hope?

Selecting user stories during sprint planning

What is the meaning of 'breadth' in breadth first search?

Is grep documentation about ignoring case wrong, since it doesn't ignore case in filenames?

Why do we bend a book to keep it straight?

Hangman Game with C++

What is this clumpy 20-30cm high yellow-flowered plant?

Time to Settle Down!

Plot a 3-D surface from a table of coordinates in Python

.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty height:90px;width:728px;box-sizing:border-box;

0

I haven't found an answer to this yet: I have a grid defined in a text file with four columns: (lon,lat,depth,slip). Each row is a grid point.

I can generate a scatter plot of these points using the following simple code:

# Main imports:
import numpy as np
from pylab import *
from mpl_toolkits.mplot3d import Axes3D

# Read the grid:
points = np.loadtxt("grid.txt")

# Retrieve parameters from the grid:
lon = points[:,0]
lat = points[:,1]
depth = points[:,2]
slip = points[:,3]

# 3-D plot of the model:
fig = figure(1)
ax = fig.add_subplot(111, projection='3d')
p = ax.scatter(lon, lat, depth, c=slip, vmin=0, vmax=max(slip), s=30, edgecolor='none', marker='o')
fig.colorbar(p)
title("Published finite fault in 3-D")
ax.set_xlabel("Longitude [degrees]")
ax.set_ylabel("Latitude [degrees]")
ax.set_zlabel("Depth [km]")
ax.invert_zaxis()
jet()
grid()
show()

And I get the following figure:

What I want to do is to be able to interpolate those points to create a "continuous" surface grid and plot it in both 2-D and 3-D plots. Therefore, somehow I've to consider all (lon,lat,depth,slip) in the interpolation. I'd appreciate your suggestions. Thanks in advance!

python python-3.x matplotlib

share|improve this question

asked Mar 8 at 20:41

Carlos HerreraCarlos Herrera

6016

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  Try tri-surface-plots
  
  – Sheldore
  Mar 8 at 20:42
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  It could work, but how can I account for the values in the fourth column (slip)? I'm trying to do that right now, still unsuccessfully.
  
  – Carlos Herrera
  Mar 8 at 20:58
  
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  I’m sure there are many examples of trisurf out there... you just have to search for them
  
  – Sheldore
  Mar 8 at 21:05
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  Check this question. It depends a bit on how your data is ordered (which is unknown here).
  
  – ImportanceOfBeingErnest
  Mar 8 at 21:22
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  In my data file, none of the columns are sorted in any way. I've checked some examples of trisurf. For a 2-D plot, the 'slip' values of my data could be used in 'z'. But in a 3-D plot, the grid is a plane with an azimuth and dip, and each point has a value of 'slip' that I want to interpolate and color on the plane. That is the main difficulty.
  
  – Carlos Herrera
  Mar 8 at 21:35
  

  
  
  
  

 | 
show 1 more comment

0

I haven't found an answer to this yet: I have a grid defined in a text file with four columns: (lon,lat,depth,slip). Each row is a grid point.

I can generate a scatter plot of these points using the following simple code:

# Main imports:
import numpy as np
from pylab import *
from mpl_toolkits.mplot3d import Axes3D

# Read the grid:
points = np.loadtxt("grid.txt")

# Retrieve parameters from the grid:
lon = points[:,0]
lat = points[:,1]
depth = points[:,2]
slip = points[:,3]

# 3-D plot of the model:
fig = figure(1)
ax = fig.add_subplot(111, projection='3d')
p = ax.scatter(lon, lat, depth, c=slip, vmin=0, vmax=max(slip), s=30, edgecolor='none', marker='o')
fig.colorbar(p)
title("Published finite fault in 3-D")
ax.set_xlabel("Longitude [degrees]")
ax.set_ylabel("Latitude [degrees]")
ax.set_zlabel("Depth [km]")
ax.invert_zaxis()
jet()
grid()
show()

And I get the following figure:

What I want to do is to be able to interpolate those points to create a "continuous" surface grid and plot it in both 2-D and 3-D plots. Therefore, somehow I've to consider all (lon,lat,depth,slip) in the interpolation. I'd appreciate your suggestions. Thanks in advance!

python python-3.x matplotlib

share|improve this question

asked Mar 8 at 20:41

Carlos HerreraCarlos Herrera

6016

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  Try tri-surface-plots
  
  – Sheldore
  Mar 8 at 20:42
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  It could work, but how can I account for the values in the fourth column (slip)? I'm trying to do that right now, still unsuccessfully.
  
  – Carlos Herrera
  Mar 8 at 20:58
  
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  I’m sure there are many examples of trisurf out there... you just have to search for them
  
  – Sheldore
  Mar 8 at 21:05
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  Check this question. It depends a bit on how your data is ordered (which is unknown here).
  
  – ImportanceOfBeingErnest
  Mar 8 at 21:22
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  In my data file, none of the columns are sorted in any way. I've checked some examples of trisurf. For a 2-D plot, the 'slip' values of my data could be used in 'z'. But in a 3-D plot, the grid is a plane with an azimuth and dip, and each point has a value of 'slip' that I want to interpolate and color on the plane. That is the main difficulty.
  
  – Carlos Herrera
  Mar 8 at 21:35
  

  
  
  
  

 | 
show 1 more comment

I haven't found an answer to this yet: I have a grid defined in a text file with four columns: (lon,lat,depth,slip). Each row is a grid point.

I can generate a scatter plot of these points using the following simple code:

# Main imports:
import numpy as np
from pylab import *
from mpl_toolkits.mplot3d import Axes3D

# Read the grid:
points = np.loadtxt("grid.txt")

# Retrieve parameters from the grid:
lon = points[:,0]
lat = points[:,1]
depth = points[:,2]
slip = points[:,3]

# 3-D plot of the model:
fig = figure(1)
ax = fig.add_subplot(111, projection='3d')
p = ax.scatter(lon, lat, depth, c=slip, vmin=0, vmax=max(slip), s=30, edgecolor='none', marker='o')
fig.colorbar(p)
title("Published finite fault in 3-D")
ax.set_xlabel("Longitude [degrees]")
ax.set_ylabel("Latitude [degrees]")
ax.set_zlabel("Depth [km]")
ax.invert_zaxis()
jet()
grid()
show()

And I get the following figure:

What I want to do is to be able to interpolate those points to create a "continuous" surface grid and plot it in both 2-D and 3-D plots. Therefore, somehow I've to consider all (lon,lat,depth,slip) in the interpolation. I'd appreciate your suggestions. Thanks in advance!

python python-3.x matplotlib

share|improve this question

asked Mar 8 at 20:41

Carlos HerreraCarlos Herrera

6016

# Main imports:
import numpy as np
from pylab import *
from mpl_toolkits.mplot3d import Axes3D

# Read the grid:
points = np.loadtxt("grid.txt")

# Retrieve parameters from the grid:
lon = points[:,0]
lat = points[:,1]
depth = points[:,2]
slip = points[:,3]

# 3-D plot of the model:
fig = figure(1)
ax = fig.add_subplot(111, projection='3d')
p = ax.scatter(lon, lat, depth, c=slip, vmin=0, vmax=max(slip), s=30, edgecolor='none', marker='o')
fig.colorbar(p)
title("Published finite fault in 3-D")
ax.set_xlabel("Longitude [degrees]")
ax.set_ylabel("Latitude [degrees]")
ax.set_zlabel("Depth [km]")
ax.invert_zaxis()
jet()
grid()
show()

share|improve this question

asked Mar 8 at 20:41

Carlos HerreraCarlos Herrera

6016

share|improve this question

asked Mar 8 at 20:41

Carlos HerreraCarlos Herrera

6016

asked Mar 8 at 20:41

Carlos HerreraCarlos Herrera

6016

asked Mar 8 at 20:41

Carlos HerreraCarlos Herrera

6016