Abbreviation for Digital Terrain Map. This time, we will use the DEM data acquired by a camera called HiRISE mounted on the Mars probe MRO.
One of the cameras equipped with MRO. Insanely high resolution. The image above seems to be one of them. (This is a pseudo color image)
Download: HiRISE Archive Site ESP and PSP are the modes when the camera acquires the image. (Any)
This time I will read it in Python. Please use gdal or something to read it. (I will skip it.)
dtm.py
xx=int(len(data_array))+1
yy=int(len(data_array[0]))+1
data_array = np.where(data_array<=-8000, data_array*(-0.0), data_array)
x = np.arange(0, xx-1, 1)
y = np.arange(0, yy-1, 1)
Z = data_array
Y, X = np.meshgrid(y, x)
fig = plt.figure()
ax = Axes3D(fig)
ax.plot_wireframe(X, Y, Z)
plt.show()
Rough direct pasting. Get the length of the array. -Elevation below 8000 is impossible ... use it. The invalid value is 0 for the time being. Specify max and min. Plot with a mesh. By the way, there were 270 lines including comment out to see the value. Even though the contents are so thin ...
For the time being, if you do it like this
You can do this.
Now that the display is complete, let's clean it up ... I don't know what this is used for w I wonder if it can be used for something ... if it's clean!
After that, I will summarize some interesting things and what I did in the extra time of my research. I have to do research and job hunting ... I want to play games ...
Recommended Posts