Borrar filtros
Borrar filtros

plot earthquake locations data on a geotiff file

5 visualizaciones (últimos 30 días)
Anitha Limann
Anitha Limann el 17 de Oct. de 2023
Comentada: Dyuman Joshi el 20 de Oct. de 2023
Hello,
I want to plot locations of earthquakes on a base map. I do have a geotiff file that i can use to create a base map but using following code i was unable to map the points on the base map. Please let me know the error in this code.
Data = load('Area_1.txt','-ascii'); % earthquake locations
lon = Data(:,1); lon = 360+lon; % longitude
lat = Data(:,2); % latitude
nv1 = importdata("volcanic.xyz"); % spreading ridge
lon1 = nv1(:,1); % longitude
lat1 = nv1(:,2); % latitude
figure(1); clf; axis equal; hold on;
basemap = readgeoraster("exportImage.tiff");
imshow(basemap);
plot(lon1, lat1,'.'); % spreading ridge
scatter(lon, lat, 15, c,'filled','o','LineWidth', 1, 'MarkerFaceAlpha', 0.9); % earthquake locations
ylabel('latitude (deg)'); xlabel('longitude (deg)');
title(['CMT Depth distribution ', Name])
xlim([min(lon) max(lon)]); ylim([min(lat) max(lat)])

Respuesta aceptada

Dyuman Joshi
Dyuman Joshi el 17 de Oct. de 2023
You should be using geoplot, geoscatter and mapshow/geoshow instead of plot(), scatter() and imshow().
If, even after making the above mentioned changes, the code is not working or the output is not as desired, please attach the data files you are working with.
  8 comentarios
Anitha Limann
Anitha Limann el 20 de Oct. de 2023
Thank you so much for your kind support. You have been extremely helpful to me. Thanks again.
Dyuman Joshi
Dyuman Joshi el 20 de Oct. de 2023
You are welcome!

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Geographic Plots en Help Center y File Exchange.

Productos

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by